Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,352,815 G→A 100% A92V (GCA→GTA)  rpsI ← ribosomal protein S9

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,352,8150GA99.8% 74.1 / ‑7.3 27A92V (GCA→GTA) rpsIribosomal protein S9
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (16/10);  minor base T (1/0);  total (17/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CGACAATTAACGTTTTGAGAATTGAGGTGAACGACGAGCTGCTTTAAGACCTGGTTTTTTACGTTCTTTCATACGTGGGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTGCTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAACCACCACCATGAACGTTAACTAAAACATCATAGTTACCTTTAGTTTCAGTTACATCAAATGGTTGG  >  CP000730/2352685‑2352937
                                                                                                                                  |                                                                                                                          
cGACAATTAACGTTTTTAGAATTGAGGTGAACGACGATTTGCTTTAAGACCTTGTTTTTTACGTTCTTTCATACGTTGGTCACGAGTAAGTAATCCATCGCGTTTTAAAGAACCTCTTTATTCTTGATCTACTTCtaataa                                                                                                                  <  2:113133/141‑1 (MQ=37)
      ttAACGTTTTGAGAATTGAGGTGAACGACGAGCTGCTTTAAGACCTGGTTTTTTACGTTCTTTCATACGTGGGCAACACACAAAGAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGAACTACTTCTAATAATGCACg                                                                                                            >  1:388932/1‑141 (MQ=37)
       tAACGTTTTGAGAATTGAGGTGAACGACGAGCTGCTTTAAGACCTGGTTTTTTACGTTCTTTCATACGTGGGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGCTCTACTTCTAATAATGCACGa                                                                                                           >  1:193814/1‑141 (MQ=255)
             tttGAGAATTGAGGTGAACGACGAGCTGCTTTAAGACCTGGTTTTTTTTTTTTTTTTTTACGTGGGAAACGAGTAAGTGATTCACCGCGTTTTAAAGGATTTCTGTATTTGGGCTTTTCTTCTaataataccagtgtgatt                                                                                                     >  2:94555/1‑129 (MQ=16)
                     ttGAGGTGAACGACGAGCTGCTTTAAGACCTGGTTTTTTACGTTCTTTCATACGTGGGTCACGAGTAAGTAATCCAGCGCGGGCGGTAGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACg                                                                                             >  1:287326/1‑141 (MQ=255)
                                   gaGCTGCTTTAAGACCTGGTTTTTTACGTTCTTTCATACGTGGGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATagcttgagctt                                                                               <  1:57658/141‑1 (MQ=255)
                                     gctgctTTAAGACCTGGTTTTTTACGTTCTTTCATACGTGGGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGt                                                                             <  1:253302/141‑1 (MQ=255)
                                      ctgctTTAAGACCTGGTTTTTTACGTTCTTTCATACGTGGGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTc                                                                            >  1:261481/1‑141 (MQ=255)
                                                    ggTTTTTTACGTTCTTTCATACGTGGGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCCGTGAAaccccca                                                              >  2:261216/1‑141 (MQ=255)
                                                    ggTTTTTTACGTTCGTTCATACGTGGGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAaccacca                                                              >  2:205760/1‑141 (MQ=255)
                                                      gtttttCCGTGCTTTCATACGTGGGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAaccaccacc                                                            <  2:57195/140‑1 (MQ=255)
                                                       tttttACGTTCTTTCATACGTGGGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTACTGCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCTAGTGAGAccaccacca                                                           >  2:887/1‑141 (MQ=255)
                                                        ttttACGTTCTTTCATACGTGGGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAACCACGCCCAt                                                          >  2:64001/1‑141 (MQ=255)
                                                        ttttACGTTCTTTCATACGTGGGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAACCACCACCAt                                                          <  2:309209/141‑1 (MQ=255)
                                                        ttttACGTGCTTTCATACGTGGGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATGCCCGCGCGCTCCCTTGACGGATAGCTTGAGCTTGTCCAGTGAAACCACCACCAt                                                          <  2:361051/141‑1 (MQ=255)
                                                                  tttCATACGTGGGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAACCACCACCATGAACGTTAAc                                                <  2:214237/141‑1 (MQ=255)
                                                                     cATACGTGGGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAACCACCACCATGAACGTTAACTaa                                             >  1:317074/1‑141 (MQ=255)
                                                                                  cGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTAGTCCAGTGGAACCACCAGCATGAACGATACCTGAAACATCCTCGttc                                >  2:136977/1‑140 (MQ=255)
                                                                                      aaagtaaTCGAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAACCACCACCATGAACGTTAACTAAAACATCATAGTTACCtt                            <  2:59366/140‑1 (MQ=255)
                                                                                                   gcgTTTTAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAACCACCACCATGAACGTTAACTAAAACATCATAGTTACCTTTAGTTTCAGTTAc               >  1:323233/1‑141 (MQ=255)
                                                                                                        ttAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAACCACCACCATGAAAGTTAACTAAAACATCATAGTTCCCTTTGGGTTCCGTTGCATCaa          >  2:357721/1‑141 (MQ=255)
                                                                                                         tAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAACCACCACCATGAACGTTAACTAAAACATCATAGTTACCTTTAGTTTCAGTTACATCaaa         <  2:8159/141‑1 (MQ=255)
                                                                                                         tAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAACCACCACCATGAACGTTAACTAAAACATCATAGTTACCTTTAGTTTCAGTTACATCaaa         >  1:56542/1‑141 (MQ=255)
                                                                                                         tAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAACCACCACCATGAACGTTAAATAAAACATCATAGTTACCTTTAGTTTCAGTTACATCaat         >  1:328434/1‑140 (MQ=255)
                                                                                                               aCCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAACCACCACCATGAACGTTAACTAAAACATCATAGTTACCTTTAGTTTCAGTTACATCAAAtggttg   >  1:161239/1‑141 (MQ=255)
                                                                                                                ccTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAACCACCACCATGAACGTTAACTAAAACATCATAGTTACCTTTAGTTTCAGTTACATCAAttggttgg  >  2:150898/1‑141 (MQ=255)
                                                                                                                ccTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAACCACCACCATGAACGTTAACTAAAACATCATAGTTACCTTTAGTTTCAGTTACATCAAAtggttgg  <  2:287326/141‑1 (MQ=255)
                                                                                                                                  |                                                                                                                          
CGACAATTAACGTTTTGAGAATTGAGGTGAACGACGAGCTGCTTTAAGACCTGGTTTTTTACGTTCTTTCATACGTGGGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTGCTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAACCACCACCATGAACGTTAACTAAAACATCATAGTTACCTTTAGTTTCAGTTACATCAAATGGTTGG  >  CP000730/2352685‑2352937

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: