Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,678,024 C→T 100% intergenic (‑121/‑403) pntA ← / → ydgH pyridine nucleotide transhydrogenase, alpha subunit/DUF1471 family periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,678,0240CT98.3% 203.0 / ‑4.4 57intergenic (‑121/‑403)pntA/ydgHpyridine nucleotide transhydrogenase, alpha subunit/DUF1471 family periplasmic protein
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (23/33);  total (23/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01
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.

CAGCGCCCGCTTGCACAAACGCTTTATCGTCAAAACTTGCCAGTTGACCCGCGCCGCTCTCTACCGCGACGGTAAAACCCAGTTTCAGCAGCTGTTCCACTGTTTTTGGCGTTGCTGCAACACGGGTTTCATTGGTTAACCGTTCTCTTGGTATGCCAATTCGCATGATATTCCCTTCCATCGGTTTTATTGATGATGGTTTGCCTGTGTCAGGAGCCACACAAGCTGCTCATGTACGAGCTAAATGTTACTCCGTTAAAATAAATTAGTAACAAACGCCCTATAACGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTAACGATTCAGCCGTGATAGTGGGATAAACACCTTAGAACGCCGGATAAAGACTGATAATTGTCTTCGACGGTCGGGTAAAACGAGACATCGCCCCGGCACGAATCACTACTTAACATTAAATTAACTTATACAATTCAGTTGCTTCAGTAGTAATGATGCTGATACGGCTGTTTTTTAAGCATAGACGGTCATTTGAGCAGGATTAAAATT  >  NC_000913/1677738‑1678302
                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                      
cAGCGCCCGCTTGCACAAACGCTTTATCGTCAAAACTTGCCAGTTGACCCGCGCCGCTCTCTACCGCGACGGTAAAACCCAGTTTCAGCAGCTGTTCCACTGTTTTTGGCGTTGCTGCAACACGGGTTTCATTGGTTAACCGTTCTCTTGGTATGCCAATTCGCATGATATTCCCTTCCATCGGTTTTATTGATGATGGTTTGCCTGTGTCAGGAGCCACACAAGCTGCTCATGTACGAGCTAAATGTTACTCCGTTAAAATAAATTAGTAACAAACGCCCTATAATGt                                                                                                                                                                                                                                                                                      >  2:522269/1‑289 (MQ=255)
       cGCTTGCACAAACGCTTTATCGTCAAAACTTGCCAGTTGACCCGCGCCGCTCTCTACCGCGACGGTAAAACCCAGTTTCAGCAGCTGTTCCACTGTTTTTGGCGTTGCTGCAACACGGGTTTCATTGGTTAACCGTTCTCTTGGTATGCCAATTCGCATGATATTCCCTTCCATCGGTTTTATTGATGATGGTTTGCCTGTGTCAGGAGCCACACAAGCTGCTCATGTACGAGCTAAATGTTACTCCGTTAAAATAAATTAGTAACAAACGCCCTATAATGTACTGaaa                                                                                                                                                                                                                                                                               >  2:397781/1‑289 (MQ=255)
          ttGCACAAACGCTTTATCGTCAAAACTTGCCAGTTGACCCGCGCCGCTCTCTACCGCGACGGTAAAACCCAGTTTCAGCAGCTGTTCCACTGTTTTTGGCGTTGCTGCAACACGGGTTTCATTGGTTAACCGTTCTCTTGGTATGACAATTCGCATGATATTCCCTTCCATCGGTTTTATTGATGATGGTTTGCCTGTGTCAGGAGCCACACAAGCTGCTCATGTACGAGCTAAATGTTACTCCGTTAAAATAAATTAGTAACAAACGCCCTATAATGTACTGAAAAtt                                                                                                                                                                                                                                                                            <  2:597158/289‑1 (MQ=255)
                              cAAAACTTGCCAGTTGACCCGCGCCGCTCTCTACCGCGACGGTAAAACCCAGTTTCAGCAGCTGTTCCACTGTTTTTGGCGTTGCTGCAACACGGGTTTCATTGGTTAACCGTTCTCTTGGTATGCCAATTCGCATGATATTCCCTTCCATCGGTTTTATTGATGATGGTTTGCCTGTGTCAGGAGCCACACAAGCTGCTCATGTACGAGCTAAATGTTACTCCGTTAAAATAAATTAGTAACAAACGCCCTATAATGTACTGAAAATTATGCCTGTGATCTAGCGCCa                                                                                                                                                                                                                                                        <  2:219467/289‑1 (MQ=255)
                                aaaCTTGCCAGTTGACCCGCGCCGCTCTCTACCGCGACGGTAAAACCCAGTTTCAGCAGCTGTTCCACTGTTTTTGGCGTTGCTGCAACACGGGTTTCATTGGTTAACCGTTCTCTTGGTATGCCAATTCGCATGATATTCCCTTCCATCGGTTTTATTGATGATGGTTTGCCTGTGTCAGGAGCCACACAAGCTGCTCATGTACGAGCTAAATGTTACTCCGTTAAAATAAATTAGTAACAAACGCCCTATAATGTACTGAAAATTATGCCt                                                                                                                                                                                                                                                                      <  2:594605/273‑1 (MQ=255)
                                aaaCTTGCCAGTTGACCCGCGCCGCTCTCTACCGCGACGGTAAAACCCAGTTTCAGCAGCTGTTCCACTGTTTTTGGCGTTGCTGCAACACGGGTTTCATTGGTTAACCGTTCTCTTGGTATGCCAATTCGCATGATATTCCCTTCCATCGGTTTTATTGATGATGGTTTGCCTGTGTCAGGAGCCACACAAGCTGCTCATGTACGAGCTAAATGTTACTCCGTTAAAATAAATTAGTAACAAACGCCCTATAATGTACTGAAAATTATGCCt                                                                                                                                                                                                                                                                      >  1:594605/1‑273 (MQ=255)
                                      gCCAGTTGACCCGCGCCGCTCTCTACCGCGACGGTAAAACCCAGTTTCAGCAGCTGTTCCACTGTTTTTGGCGTTGCTGCAACACGGGTTTCATTGGTTAACCGTTCTCTTGGTATGCCAATTCGCATGATATTCCCTTCCATCGGTTTTATTGATGATGGTTTGCCTGTGTCAGGAGCCACACAAGCTGCTCATGTACGAGCTAAATGTTACTCCGTTAAAATAAATTAGTAACAAACGCCCTATAATGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAg                                                                                                                                                                                                                                                >  2:288166/1‑289 (MQ=255)
                                                                  cgACGGTAAAACCCAGTTTCAGCAGCTGTTCCACTGTTTTTGGCGTTGCTGCAACACGGGTTTCATTGGTTAACCGTTCTCTTGGTATGCCAATTCGCATGATATTCCCTTCCATCGGTTTTATTGATGATGGTTTGCCTGTGTCAGGAGCCACACAAGCTGCTCATGTACGCGCTAAATGTTACTCCGTTAAAATAAATTAGTAACAAACGCCCTATAATGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTa                                                                                                                                                                                                                    >  1:678062/1‑289 (MQ=255)
                                                                  cgACGGTAAAACCCAGTTTCAGCAGCTGTTCCACTGTTTTTGGCGTTGCTGCAACACGGGTTTCATTGGTTAACCGTTCTCTTGGTATGCCAATTCGCATGATATTCCCTTCCATCGGTTTTATTGATGATGGTTTGCCTGTGTCAGGAGCCACACAAGCTGCTCATGTACGAGCTAAATGTTACTCCGTTAAAATAAATTAGTAACAAACGCCCTATAATGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTa                                                                                                                                                                                                                    >  2:195947/1‑289 (MQ=255)
                                                                  cgACGGTAAAACCCAGTTTCAGCAGCTGTTCCACTGTTTTTGGCGTTGCTGCAACACGGGTTTCATTGGTTAACCGTTCTCTTGGTATGCCAATTCGCATGATATTCCCTTCCATCGGTTTTATTGATGATGGTTTGCCTGTGTCAGGAGCCACACAAGCTGCTCATGTACGAGCTAAATGTTACTCCGTTAAAATAAATTAGTAACAAACGCCCTATAATGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTa                                                                                                                                                                                                                    >  1:268066/1‑289 (MQ=255)
                                                                   gACGGTAAAGCCCAGTTTCAGAAGCTGTTCCACTGTTTGTGGCGTTGCTGCAACACGGGTTTCATTGGTTAACCGTTCTCTTTGTATGCCAATTCGCATGATATTCCCTTCCATCGGTTTTATTGATGATGGTTTGCCTGTGTCAGGAGCCACACAAGCTGCTCATGTACGAGCTAAATGTTACTCCGTTAAAATAAATTAGTAACAAACGCCCTATAATGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTaa                                                                                                                                                                                                                   <  2:646842/289‑1 (MQ=255)
                                                                                  tttCAGCAGCTGTTCCACTGTTTTTGGCGTTGCTGCAACACGGGTTTCATTGGTTAACCGTTCTCTTGGTATGCCAATTCGCATGATATTCCCTTCCATCGGTTTTATTGATGATGGTTTGCCTGTGTCAGGAGCCACCCAAGCTGCTCATGTACGAGCTAAATGTTACTCCGTTAAAATAAATTAGTAACAAACGCCCTATAATGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTAACGATTCAGCCGTGAt                                                                                                                                                                                                    >  2:22166/1‑289 (MQ=255)
                                                                                   ttCAGCAGCTGTTCCACTGTTTTTGGCGTTGCTGCAACACGGGTTTCATTGGTTAACCGTTCTCTTGGTATGCCAATTCGCATGATATTCCCTTCCATCGGTTTTATTGATGATGGTTTGCCTGTGTCAGGAGCCACACAAGCTGCTCATGTACGAGCTAAATGTTACTCCGTTAAAATAAATTAGTAACAAACGCCCTATAATGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTAACGATTCAGCCGTGa                                                                                                                                                                                                     <  1:595145/287‑1 (MQ=255)
                                                                                    tCAGCAGCTGTTCCACTGTTTTTGGCGTTGCTGCAACACGGGTTTCATTGGTTAACCGTTCTCTTGGTATGCCAATTCGCATGATATTCCCTTCCATCGGTTTTATTGATGATGGTTTGCCTGTGTCAGGAGCCACACAAGCTGCTCATGTACGAGCTAAATGTTACTCCGTTAAAATAAATTAGTAACAAACGCCCTATAATGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTAACGATTCAGCCGTGATAg                                                                                                                                                                                                  <  1:767631/289‑1 (MQ=255)
                                                                                     cagcagCTGTTCCACTGTTTTTGGCGTTGCTGCAACACGGGTTTCATTGGTTAACCGTTCTCTTGGTATGCCAATTCGCATGATATTCCCTTCCATCGGTTTTATTGATGATGGTTTGCCTGTGTCAGGAGCCACACAAGCTGCTCATGTACGAGCTAAATGTTACTCCGTTAAAATAAATTAGTAACAAACGCCCTATAATGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTAACGATTCAGCCGTGATAGt                                                                                                                                                                                                 <  2:420646/289‑1 (MQ=255)
                                                                                             gTTCCACTGTTTTTGGCGTTGCTGCAACACGGGTTTCATTGGTTAACCGTTCTCTTGGTATGCCAATTCGCATGATATTCCCTTCCATCGGTTTTATTGATGATGGTTTGCCTGTGTCAGGAGCCACACAAGCTGCTCATGTACGAGCTAAATGTTACTCCGTTAAAATAAATTAGTAACAAACGCCCTATAATGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTAACGATTCAGCCGTGATAGTGGGATaaa                                                                                                                                                                                         <  1:288166/289‑1 (MQ=255)
                                                                                              gtCCACGGGTTTGTGCGTGGCTGCAACACGGGTTTCATTGGTTAACCGTTCTCTTGGTATGCCAATTCGCATGAGATTCCCTTCCGTCGGTTTTATTGATGATGGTGTGCCTGTGTCAGGGGCCACACAAGCTGCTCATGTACGGGCTAAATGTTACTCCGTTAAAATAGATTAGTAACAAACGCCCGATAATGTACTGAAAATTATGGCGGTGATCTAGCGCGAAAAATCAGTATGTGGGCGGGAACTCGCGAAATATTAACGATTCAGCGGTGATAGTGGGATaaa                                                                                                                                                                                         <  2:678062/287‑1 (MQ=255)
                                                                                                        tttGGCGTTGCTGCAACACGGGTTTCATTGGTTAACCGTTCTCTTGGTATGCCAATTCGCATGATATTCCCTTCCATCGGTTTTATTGATGATGGTTTGCCTGTGTCAGGAGCCACACAAGCTGCTCATGTACGAGCTAAATGTTACTCCGTTAAAATAAATTAGTAACAAACGCCCTATAATGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTAACGATTCAGCCGTGATAGTGGGATAAACACCTTAGAAc                                                                                                                                                                              >  2:706534/1‑289 (MQ=255)
                                                                                                          tGGCGTTGCTGCAACACGGGTTTCATTGGTTAACCGTTCTCTTGGTATGCCAATTCGCATGATATTCCCTTCCATCGGTTTTATTGATGATGGTTTGCCTGTGTCAGGAGCCACACAAGCTGCTCATGTACGAGCTAAATGTTACTCCGTTAAAATAAATTAGTAACAAACGCCCTATAATGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTAACGATTCAGCCGTGATAGTGGGATAAACACCTTAGAACGc                                                                                                                                                                            <  1:47451/289‑1 (MQ=255)
                                                                                                          tGGCGTTGCTGCAACACGGGTTTCATTGGTTAACCGTTCTCTTGGTATGCCAATTCGCATGATATTCCCTTCCATCGGTTTTATTGATGATGGTTTGCCTGTGTCAGGAGCCACACAAGCTGCTCATGTACGAGCTAAATGTTACTCCGTTAAAATAAATTAGTAACAAACGCCCTATAATGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTAACGATTCAGCCGTGATAGTGGGATAAACACCTTAGAACGc                                                                                                                                                                            <  1:397781/289‑1 (MQ=255)
                                                                                                          tGGCGTTGCTGCAACACGGGTTTCATTGGTTAACCGTTCTCTTGGTATGCCAATTCGCATGATATTCCCTTCCATCGGTTTTATTGATGATGGTTTGCCTGTGTCAGGAGCCACACAAGCTGCTCATGTACGAGCTAAATGTTACTCCGTTAAAATAAATTAGTAACAAACGCCCTATAATGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTAACGATTCAGCCGTGATAGTGGGATAAACACCTTAGAACGc                                                                                                                                                                            <  1:518798/289‑1 (MQ=255)
                                                                                                          tGGCGTTGCTGCAACACGGGTTTCATTAGTTAACCGTTCTCTTGGTATGCCAATTCGCATGATATTCCCTTCCATCGGTTTTATTGATGATGGTTTGCCTGTGTCAGGAGCCACACAAGCTGCTCATGTACGAGCTAAATGTTACTCCGTTAAAATAAATTATTAACAAACGCCCTATAATGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTAACGATTCAGCCGTGATAGTGGGATAAACACCTTAGAACGc                                                                                                                                                                            <  1:550216/289‑1 (MQ=255)
                                                                                                                tgctgcAACACGGGTTTCATTGGTTAACCGTTCTCTTGGTATGCCAATTCGCATGATATTCCCTTCCATCGGTTTTATTGATGATGGTTTGCCTGTGTCAGGAGCCACACAAGCTGCTCATGTACGAGCTAAATGTTACTCCGTTAAAATAAATTAGTAACAAACGCCCTATAATGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTAACGATTCAGCCGTGATAGTGGGATAAACACCTTAGAACGCCGGATa                                                                                                                                                                      >  2:150671/1‑289 (MQ=255)
                                                                                                                tgctgcAACACGGGTTTCATTGGTTAACCGTTCTCTTGGTATGCCAATTCGCATGATATTCCCTTCCATCGGTTTTATTGATGATGGTTTGCCTGTGTCAGGAGCCACACAAGCTGCTCATGTACGAGCTAAATGTTACTCCGTTAAAATAAATTAGTAACAAACGCCCTATAATGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTAACGATTCAGCCGTGATAGTGGGATAAACACCTTAGAACGCCGGATa                                                                                                                                                                      <  1:260835/289‑1 (MQ=255)
                                                                                                                                  aTTGGTTAACCGTTCTCTTGGTATGCCAATTCGCATGATATTCCCTTCCATCGGTTTTATTGATGATGGTTTGCCTGTGTCAGGAGCCACACAAGCTGCTCATGTACGAGCTAAATGTTACTCCGTTAAAATAAATTAGTAACAAACGCCCTATACTGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTAACGATTCAGCCGTGATAGTGGGATAAACACCTTAGAAAGCCGTATAAAGACTGATAATTGTAtt                                                                                                                                                    >  2:123197/1‑289 (MQ=255)
                                                                                                                                     ggTTAACCGTTCTCTTGGTATGCCAATTCGCATGATATTCCCTTCCATCGGTTTTATTGATGATGGTTTGCCTGTGTCAGGAGCCACACAAGCTGCTCATGTACGAGCTAAATGTTACTCCGTTAAAATAAATTAGTAACAAACGCCCTATAATGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTAACGATTCAGCCGTGATAGTGGGATAAACACCTTAGAACGCCGGATAAAGACTGATAATTGTCTTCGa                                                                                                                                                 <  1:522269/289‑1 (MQ=255)
                                                                                                                                              ttCTCTTGGTATGCCAATTCGCATGATATTCCCTTCCATCGGTTTTATTGATGATGGTTTGCCTGTGTCAGGAGCCACACAAGCTGCTCATGTACGAGCTAAATGTTACTCCGTTAAAATAAATTAGTAACAAACGCCCTATAATGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTAACGATTCAGCCGTGATAGTGGGATAAACACCTTAGAACGCCGGATAAAGACTGATAATTGTCTTCGACGGTCGGGt                                                                                                                                        >  1:100131/1‑289 (MQ=255)
                                                                                                                                                ctctTGGTATGCCAATTCGCATGATATTCCCTTCCATCGGTTTTATTGATGATGGTTTGCCTGTGTCAGGAGCCACACAAGCTGCTCATGTACGAGCTAAATGTTACTCCGTTAAAATAAATTAGTAACAAACGCCCTATAATGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTAACGATTCAGCCGTGATAGTGGGATAAACACCTTAGAACGCCGGATAAAGACTGATAATTGTCTTCGACGGTCGGGTaa                                                                                                                                      >  1:62557/1‑289 (MQ=255)
                                                                                                                                                     ggTATGCCAATTCGCATGATATTCCCTTCCATCGGTTTTATTGATGATGGTTTGCCTGTGTCAGGAGCCACACAAGCTGCTCATGTACGAGCTAAATGTTACTCCGTTAAAATAAATTAGTAACAAACGCCCTATAATGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTAACGATTCAGCCGTGATAGTGGGATAAACACCTTAGAACGCCGGATAAAGACTGATAATTGTCTTCGACGGTCGGGTAAAACga                                                                                                                                 >  1:200457/1‑289 (MQ=255)
                                                                                                                                                              aTTCGAATGTTATTCCCTTCCATCGGTTTTATGGAGGATGGTTTGCCTGTGTCAGGAGCCACACAAGCTGCTCATGTACGAGCTAAATGTTACTCCGTTAAAATAAATTAGTAACAAACGCCCTATAATGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTAACGATTCAGCCGTGATAGTGGGATAAACACCTTAGAACGCCGGATAAAGACTGATAATTGTCTTCGACGGTCGGGTAAAACGAGACATCGcc                                                                                                                        <  2:5368/289‑1 (MQ=255)
                                                                                                                                                                    aTGATATTCCCTTCCATCGGTTTTATTGATGATGGTTTGCCTGTGTCAGGAGCCACACAAGCTGCTCATGTACGAGCTAAATGTTACTCCGTTAAAATAAATTAGTAACAAACGCCCTATAATGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTAACGATTCAGCCGTGATAGTGGGATAAACACCTTAGAACGCCGGATAAAGACTGATAATTGTCTTCGACGGTCGGGTAAAACGAGACATCGCCCCGGCa                                                                                                                  <  2:100131/289‑1 (MQ=255)
                                                                                                                                                                    aTGATATTCCCTTCCATCGGTTTTATTGATGATGGTTTGCCTGTGTCAGGAGCCACACAAGCTGCTCATGTACGAGCTAAATGTTACTCCGTTAAAATAAATTAGTAACAAACGCCCTATAATGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTAACGATTCAGCCGTGATAGTGGGATAAACACCTTAGAACGCCGGATAAAGACTGATAATTGTCTTCGACGGTCGGGTAAAACGAGACATCGCCCCGGCa                                                                                                                  <  2:62557/289‑1 (MQ=255)
                                                                                                                                                                         atTCCCTTCCATCGGTTTTATTGATGATGGTTTGCCTGTGTCAGGAGCCACACAAGCTGCTCATGTACGAGCTAAATGTTACTCCGTTAAAATAAATTAGTAACAAACGCCCTATAATGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTAACGATTCAGCCGTGATAGTGGGATAAACACCTTAGAACGCCGGATAAAGACTGATAATTGTCTTCGACGGTCGGGTAAAACGAGACATCGCCCCGGCACGAAt                                                                                                             <  2:395815/289‑1 (MQ=255)
                                                                                                                                                                                          ttATTGATGATGGTTTGCCTGTGTCAGGAGCCACACAAGCTGCTCATGTACGAGCTAAATGTTACTCCGTTAAAATAAATTAGTAACAAACGCCCTATAATGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTAACGATTCAGCCGTGATAGTGGGATAAACACCTTAGAACGCCGGATAAAGACTGATAATTGTCTTCGACGGTCGGGTAAAACGAGACATCGCCCCGGCACGAATCACTACTTAACattaaa                                                                                            <  1:706534/289‑1 (MQ=255)
                                                                                                                                                                                          ttATTGATGATGGTTTGCCTGTGTCAGGAGCCACACAAGCTGCTCATGTACGAGCTAAATGTTACTCCGTTAAAATAAATTAGTAACAAACGCCCTATAATGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTAACGATTCAGCCGTGATAGTGGGATAAACACCTTAGAACGCCGGATAAAGACTGATAATTGTCTTCGACGGTCGGGTAAAACGAGACATCGCCCCGGCACGAATCACTACTTAACattaaa                                                                                            <  1:123197/289‑1 (MQ=255)
                                                                                                                                                                                          ttATTGATGATGGTTTGCCTGTGTCAGGAGCCACACAAGATGCTCATGTACGAGCTAAATGTTACTCCGTTAAAATAAATTAGTAACAAACGCCCTATAATGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTAACGATTCAGCCGTGATAGTGGGATAAACACCTTAGAACGCCGGATAAAGACTGATAATTGTCTTCGACGGTCGGGTAAAACGAGACATCGCCCCGGCACGAATCACTACTTAACattaaa                                                                                            <  2:333764/289‑1 (MQ=255)
                                                                                                                                                                                                  gatgGTTTGCCTGTGTCAGGAGCCACACAAGCTGCTCATGTACGAGCTAAATGTTACTCCGTTAAAATAAATTAGTAACAAACGCCCTATAATGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTAACGATTCAGCCGTGATAGTGGGATAAACACCTTAGAACGCCGGATAAAGACTGATAATTGTCTTCGACGGTCGGGTAAAACGAGACATCGCCCCGGCACGAAt                                                                                                             >  2:549466/1‑264 (MQ=255)
                                                                                                                                                                                                  gatgGTTTGCCTGTGTCAGGAGCCACACAAGCTGCTCATGTACGAGCTAAATGTTACTCCGTTAAAATAAATTAGTAACAAACGCCCTATAATGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTAACGATTCAGCCGTGATAGTGGGATAAACACCTTAGAACGCCGGATAAAGACTGATAATTGTCTTCGACGGTCGGGTAAAACGAGACATCGCCCCGGCACGAAt                                                                                                             <  1:549466/264‑1 (MQ=255)
                                                                                                                                                                                                  gatgGTTTGCCTGTGTCAGGAGCCACACAAGCTGCTCATGTACGAGCTAAATGTTACTCCGTTAAAATAAATTAGTAACAAACGCCCTATAATGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTAACGATTCAGCCGTGATAGTGGGATAAACACCTTAGAACGCCGGATAAAGACTGATAATTGTCTTCGACGGTCGGGTAAAACGAGACATCGCCCCGGCACGAATCACTACTTAACATTAAATTAACTta                                                                                    >  2:86481/1‑289 (MQ=255)
                                                                                                                                                                                                           ccTGTGTCAGGAGCCACACAAGCTGCTCATGTACGAGCTAAATGTTACTCCGTTAAAATAAATTAGTAACAAACGCCCTATAATGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTAACGATTCAGCCGTGATAGTGGGATAAACACCTTAGAACGc                                                                                                                                                                            <  2:595671/192‑1 (MQ=255)
                                                                                                                                                                                                           ccTGTGTCAGGAGCCACACAAGCTGCTCATGTACGAGCTAAATGTTACTCCGTTAAAATAAATTAGTAACAAACGCCCTATAATGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTAACGATTCAGCCGTGATAGTGGGATAAACACCTTAGAACGc                                                                                                                                                                            >  1:595671/1‑192 (MQ=255)
                                                                                                                                                                                                                           cacaAGCTGCTCATGTACGAGCTAAATGTTACTCCGTTAAAATAAATTAGTAACAAACGCCCTATAATGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTAACGATTCAGCCGTGATAGTGGGATAAACACCTTAGAACGCCGGATAAAGACTGATAATTGTCTTCGACGGTCGGGTAAAACGAGACATCGCCCCGGCACGAATCACTACTTAACATTAAATTAACTTATACAATTCAGt                                                                         >  1:632950/1‑275 (MQ=255)
                                                                                                                                                                                                                           cacaAGCTGCTCATGTACGAGCTAAATGTTACTCCGTTAAAATAAATTAGTAACAAACGCCCTATAATGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTAACGATTCAGCCGTGATAGTGGGATAAACACCTTAGAACGCCGGATAAAGACTGATAATTGTCTTCGACGGTCGGGTAAAACGAGACATCGCCCCGGCACGAATCACTACTTAACATTAAATTAACTTATACAATTCAGt                                                                         <  2:632950/275‑1 (MQ=255)
                                                                                                                                                                                                                                                         aCTCCGTTAAAATAAATTAGTAACAAACGCCCTATAATGTACTGAAAATTATGCCTGTGATCTAgcg                                                                                                                                                                                                                                                           >  2:104658/1‑67 (MQ=255)
                                                                                                                                                                                                                                                         aCTCCGTTAAAATAAATTAGTAACAAACGCCCTATAATGTACTGAAAATTATGCCTGTGATCTAgcg                                                                                                                                                                                                                                                           <  1:104658/67‑1 (MQ=255)
                                                                                                                                                                                                                                                             cGTTAAAATAAATTAGTAACAAACGCCCTATAATGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTAACGATTCAGCCGTGATAGTGGGATAAACACCTTAGAACGCCGGATa                                                                                                                                                                      <  2:330837/148‑1 (MQ=255)
                                                                                                                                                                                                                                                             cGTTAAAATAAATTAGTAACAAACGCCCTATAATGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTAACGATTCAGCCGTGATAGTGGGATAAACACCTTAGAACGCCGGATa                                                                                                                                                                      >  1:330837/1‑148 (MQ=255)
                                                                                                                                                                                                                                                             cGTTAAAATAAATTAGTAACAAACGCCCTATAATGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTAACGATTCAGCCGTGATAGTGGGATAAACACCTTAGAACGCCGGATAAAGACTGATAATTGTCTTCGACGGTCGGGTAAAACGAGACATCGCCCCGGCACGAATCACTACTTAACATTAAATTAACTTATACAATTCAGTTGCTTCAGTAGTAATGATGCTGATACGGCTGTTTTTTAAGCATAGACg                         <  1:330340/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                       aatTAGTAACAAACGCCCTATAATGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGc                                                                                                                                                                                                                             >  2:131915/1‑83 (MQ=255)
                                                                                                                                                                                                                                                                       aatTAGTAACAAACGCCCTATAATGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGc                                                                                                                                                                                                                             <  1:131915/83‑1 (MQ=255)
                                                                                                                                                                                                                                                                       aatTAGTAACAAACGCCCTATAATGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTAACGATTCAGCCGTGATAGTGGGATAAACACCTTAGAACGCCGGATAAAGACTGATAATTGTCTTCGACGGTCGGGTAAAACGAGACATCGCCCCGGCACGAATCACTACTTAACATTAAATTAACTTATACAATTCAGTTGCTTCAGTAGTAATGATGCTGATACGGCTGTTTTTTAAGCATAGACGGTCATTTGAg               <  1:150671/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                        atTAGTAACAAACGCCCTATAACATACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTAACGATTCAGCCGTGATAGTGGGATAAACACCTTAGAACGCCGGATAAAGACTGATAATTGTCTTCGACGGTCGGGTAAAACGAGACATCGCCCCGGCACGAATCACTACTTAACATTAAATTAACTTATACAATTCAGTTGCTTCAGTAGTAATGATGCTGATACGGCTGTTTTTTAAGCATAGACGGTCATTTGAGc              <  2:106130/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                              accAAAAGCCCTATAATTTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACACGCAAAATATTAACGATTCAGCCGTGATAGTGGGATAAACACCTTAGAACGCCGGATAAAGACTGATAATTGTCTTCGACGGTCGGGTAAAACGAGACATCGCCCCGGCACGAATCACTACTTAACATTAAATTAACTTATACAATTCAGTTGCTTCAGTAGTAATGATGCTGATACGGCTGTTTTTTAAGCATAGACGGTCATTTGAGCAGGa          <  2:268066/285‑1 (MQ=255)
                                                                                                                                                                                                                                                                                   aCGCCCTATAATGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTAACGATTCAGCCGTGATAGTGGGATAAACACCTTAGAACGCCGGATAAAGACTGATAATTGTCTTCGACGGTCGGGTAAAACGAGACATCGCCCCGGCACGAATCACTACTTAACATTAAATTAACTTATACAATTCAGTTGCTTCAGTAGTAATGATGCTGATACGGCTGTTTTTTAAGCATAGACGGTCATTTGAGCAGGATTAAAAt   >  1:107332/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                                    cGCCCTATAATGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTAACGATTCAGCCGTGATAGTGGGATAAACACCTTAGAACGCCGGATAAAGACTGATAATTGTCTTCGACGGTCGGGTAAAACGAGACATCGCCCCGGCACGAATCACTACTTAACATTAAATTAACTTATACAATTCAGTTGCTTCAGTAGTAATGATGCTGATACGGCTGTTTTTTAAGCATAGACGGTCATTTGAGCAGGATTAAAAtt  <  2:200457/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                                         tataATGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTAACGATTCAGCCGTGATAGTGGGATAAACACCTTAGAACGCCGGATAAAGACTGATAATTGTCTTCGACGGTCGGGTAAAACGAGACATCGCCCCGGCACGAATCa                                                                                                           >  1:647103/1‑179 (MQ=255)
                                                                                                                                                                                                                                                                                         tataATGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTAACGATTCAGCCGTGATAGTGGGATAAACACCTTAGAACGCCGGATAAAGACTGATAATTGTCTTCGACGGTCGGGTAAAACGAGACATCGCCCCGGCACGAATCa                                                                                                           <  2:647103/179‑1 (MQ=255)
                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                      
CAGCGCCCGCTTGCACAAACGCTTTATCGTCAAAACTTGCCAGTTGACCCGCGCCGCTCTCTACCGCGACGGTAAAACCCAGTTTCAGCAGCTGTTCCACTGTTTTTGGCGTTGCTGCAACACGGGTTTCATTGGTTAACCGTTCTCTTGGTATGCCAATTCGCATGATATTCCCTTCCATCGGTTTTATTGATGATGGTTTGCCTGTGTCAGGAGCCACACAAGCTGCTCATGTACGAGCTAAATGTTACTCCGTTAAAATAAATTAGTAACAAACGCCCTATAACGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTAACGATTCAGCCGTGATAGTGGGATAAACACCTTAGAACGCCGGATAAAGACTGATAATTGTCTTCGACGGTCGGGTAAAACGAGACATCGCCCCGGCACGAATCACTACTTAACATTAAATTAACTTATACAATTCAGTTGCTTCAGTAGTAATGATGCTGATACGGCTGTTTTTTAAGCATAGACGGTCATTTGAGCAGGATTAAAATT  >  NC_000913/1677738‑1678302

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: