| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description | 
| RA | NC_000913 | 2,028,160 | G→C | 20.0% | intergenic (‑143/+28) | yedQ ← / ← yodC | putative membrane‑anchored diguanylate cyclase/uncharacterized protein | 
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,028,160 | 0 | G | C | 20.0% | 77.0 / 13.7 | 30 | intergenic (‑143/+28) | yedQ/yodC | putative membrane‑anchored diguanylate cyclase/uncharacterized protein | 
| Reads supporting (aligned to +/- strand): ref base G (12/12); new base C (3/3); total (15/15) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.49e-01 | |||||||||||
CCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCC  >  NC_000913/2028028‑2028307                                                                                                                                    |                                                                                                                                                   cccAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATTTGATCCgg                                                                                                                                               <  1:967470/139‑1 (MQ=255)cccAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCgg                                                                                                                                               >  2:967470/1‑139 (MQ=255)  cAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGAt                                                                                                                                                   >  4:656765/1‑133 (MQ=255)  cAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGAt                                                                                                                                                   <  3:656765/133‑1 (MQ=255)                                                    aaCAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACa                                                                                  >  1:288773/1‑148 (MQ=255)                                                                tttATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGcc                                                                                                                                             <  3:570885/77‑1 (MQ=39)                                                                tttATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGcc                                                                                                                                             >  4:570885/1‑77 (MQ=38)                                                                               tGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATg                                                                                                                                  >  1:989661/1‑73 (MQ=255)                                                                               tGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATg                                                                                                                                  <  2:989661/73‑1 (MQ=255)                                                                                 aTAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTc                                                                                                                        <  4:66117/81‑1 (MQ=255)                                                                                 aTAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTc                                                                                                                        >  3:66117/1‑81 (MQ=255)                                                                                          gACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTa                                            <  2:288773/148‑1 (MQ=255)                                                                                                  ttcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGAt                                                                                                                                >  2:663062/1‑56 (MQ=255)                                                                                                  ttcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGAt                                                                                                                                >  4:433894/1‑56 (MQ=255)                                                                                                  ttcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGAt                                                                                                                                <  1:663062/56‑1 (MQ=255)                                                                                                  ttcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGAt                                                                                                                                <  3:433894/56‑1 (MQ=255)                                                                                                   tcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTCAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACa                                                                                  <  4:427432/101‑1 (MQ=255)                                                                                                   tcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTCAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACa                                                                                  >  3:427432/1‑101 (MQ=255)                                                                                                     attGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCt                                                      <  2:457827/127‑1 (MQ=255)                                                                                                     attGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCt                                                      >  1:457827/1‑127 (MQ=255)                                                                                                               ttAGCACACCAGGCCGGATGTGATCCGGCCTGATg                                                                                                                                        >  2:65486/1‑35 (MQ=255)                                                                                                               ttAGCACACCAGGCCGGATGTGATCCGGCCTGATg                                                                                                                                        <  1:65486/35‑1 (MQ=255)                                                                                                                       tcAGGCCGGATCACATCCGGCCTGGTGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTc                                                                       <  4:10033/91‑1 (MQ=38)                                                                                                                       tcAGGCCGGATCACATCCGGCCTGGTGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTc                                                                       >  3:10033/2‑92 (MQ=38)                                                                                                                       tcAGGCCGGATCACATCCGGCCTGGTGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATAc                                  >  3:463616/2‑129 (MQ=255)                                                                                                                       tcAGGCCGGATCACATCCGGCCTGGTGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATAc                                  <  4:463616/128‑1 (MQ=255)                                                                                                                       ccAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGaa                                                                                                       <  3:138048/60‑1 (MQ=255)                                                                                                                       ccAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGaa                                                                                                       >  4:138048/1‑60 (MQ=255)                                                                                                                          ggCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGaaa                                                                 >  3:254062/1‑95 (MQ=255)                                                                                                                          ggCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGaaa                                                                 <  4:254062/95‑1 (MQ=255)                                                                                                                                    gATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATcc  >  3:656386/1‑148 (MQ=255)                                                                                                                                    |                                                                                                                                                   CCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCC  >  NC_000913/2028028‑2028307 | 
| Alignment Legend | 
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG | 
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑ |