Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,160,854 C→T intergenic (‑64/‑270) sthA ← / → fabR pyridine nucleotide transhydrogenase, soluble/transcriptional repressor of fabA and fabB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,160,8540CT99.2% 434.3 / ‑6.3 122intergenic (‑64/‑270)sthA/fabRpyridine nucleotide transhydrogenase, soluble/transcriptional repressor of fabA and fabB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (58/62);  minor base G (0/1);  total (58/64)
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.04e-01

Only 100 of 122 total aligned reads displayed.

GTCGGGAATGGTCGCTGTAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGCCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAG  >  NC_000913/4160559‑4161136
                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                          
gTCGGGAATGGTCGCTGTAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAAtt                                                                                                                                                                                                                                                                                       >  1:823635/1‑301 (MQ=255)
 tCGGGAATGGTGGCTGTAAAGTGGGTTTTGATTGAATTCTAGAATGCGGCTGACGGCGTGACGGAGAGCTTTCGGCGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGACAATTg                                                                                                                                                                                                                                                                                      <  1:766289/301‑1 (MQ=255)
  cGGGAATGGTCGCTGTAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGa                                                                                                                                                                                                                                                                                     <  1:799814/301‑1 (MQ=255)
     gAATGGTCGGTGTAAAGTTGGGTTTGAGTGAATTCTGTAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCGGCAAACATTTTGATAACGCTCGATAACTGCAACGCGCGCACCTTGCTTAACCAGGCCTATTGCAGCGTCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTAAATTGATc                                                                                                                                                                                                                                                                                   <  1:452207/300‑1 (MQ=255)
         ggTCGCTGTAAAGTGGGGTTTGAGTTGATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCg                                                                                                                                                                                                                                                                               <  1:541830/300‑1 (MQ=255)
              cTGTAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCAttt                                                                                                                                                                                                                                                                         >  1:568607/1‑301 (MQ=255)
              cTGTAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCAttt                                                                                                                                                                                                                                                                         >  1:309208/1‑301 (MQ=255)
                     gTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGTTAACATTTTTTCGCGGGTAAGTCAATTGATCCCCGTCATTTTTCTgg                                                                                                                                                                                                                                                                   >  1:228080/1‑300 (MQ=255)
                            ttgattgcATGCTATGAGGCGGCGGACGGCGTGACGGAGAGCTTGCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACGGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATcc                                                                                                                                                                                                                                                           <  1:592253/301‑1 (MQ=255)
                            ttgattgaATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATcc                                                                                                                                                                                                                                                           <  1:101079/301‑1 (MQ=255)
                            ttgattgaATTCTATAAGGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATcc                                                                                                                                                                                                                                                           <  1:660695/301‑1 (MQ=255)
                            ttgattgaATTCTATAAGGCGGCTGACGGCGTGACGGAGAGCTGTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATcc                                                                                                                                                                                                                                                           <  1:99825/301‑1 (MQ=255)
                                    aTTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCACCGTCATTTATCTGGCTATATCCTGAGCgg                                                                                                                                                                                                                                                    >  1:942741/1‑300 (MQ=255)
                                       cTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCctt                                                                                                                                                                                                                                                >  1:719960/1‑301 (MQ=255)
                                       cTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCctt                                                                                                                                                                                                                                                >  1:543504/1‑301 (MQ=255)
                                       cTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCctt                                                                                                                                                                                                                                                >  1:216974/1‑301 (MQ=255)
                                       cTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCATTCGCCGACGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTCTACATAATAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCATGAGCGGCctt                                                                                                                                                                                                                                                >  1:429236/1‑301 (MQ=255)
                                                             ccGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATACTGAGCGGCCTTTGCTTTGTCTGTTTCTTACttt                                                                                                                                                                                                                          >  1:140676/2‑301 (MQ=255)
                                                                                   gTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCCTGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGTCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGCtt                                                                                                                                                                                                    >  1:5144/1‑301 (MQ=255)
                                                                                     nccccAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGa                                                                                                                                                                                                                                                                                     <  2:823597/217‑1 (MQ=255)
                                                                                     gCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTAc                                                                                                                                                                                                                             >  1:659271/1‑274 (MQ=255)
                                                                                     gCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGAtttt                                                                                                                                                                                                  >  1:746041/1‑301 (MQ=255)
                                                                                     aaagcAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCCCGCACCTTGCTTAACCAGGCCCATGGCAGCGCCTTCGCAGCCGGGGCCGGAACCGATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTgt                                                                                                                                                                                                                                                                                     <  2:655277/214‑2 (MQ=255)
                                                                                      ccccAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTg                                                                                                                                                                                                                                                                                      <  2:866821/216‑1 (MQ=255)
                                                                                         cAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTCTTCTTATACATAAAAGCAACAGTATGGTAACATTTTATCGCGGGTAAGTCAATTgt                                                                                                                                                                                                                                                                                     >  1:655300/1‑213 (MQ=255)
                                                                                                    aCCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATc                                                                                                                                                                                                                                                                      <  2:568589/218‑1 (MQ=255)
                                                                                                               cATTTTGATAACGCTTGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTatattc                                                                                                                                                                                                                                                           <  2:25673/218‑3 (MQ=255)
                                                                                                                          cGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGc                                                                                                                                                                                                                                                   >  2:848900/1‑215 (MQ=255)
                                                                                                                            cTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCg                                                                                                                                                             >  1:740883/1‑299 (MQ=255)
                                                                                                                                aTAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCctttgct                                                                                                                                                                                                                                            <  2:719932/216‑1 (MQ=255)
                                                                                                                                aTAACGGCGACGCGCGCACCTTGCTTAGACATGCACATTGCAGCGCCTTCGCCGCCGGGGCAGCAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGTTAGGGCTTAGCTGTTCTTATACATAAAAGCAACAGAATGGTAAAATTTGATCGCGGGTAAGTCAATTGTTCCCCGTCATTTGTCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACGCTCCGTTTAAGCGGCAGGGTTCCGCTGTAc                                                                                                                                                       <  1:942507/301‑1 (MQ=255)
                                                                                                                                 taccTGCGACGCGCGAACCTTGCTTAACCAGGCCCATTGCAGCGACTGCGCCGCCGGGGACGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCctttgctttg                                                                                                                                                                                                                                         <  2:228075/215‑1 (MQ=255)
                                                                                                                                          cgcgcgTACCATGATGAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCTGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAAc                                                                                                                                             <  1:785544/301‑1 (MQ=255)
                                                                                                                                              cgcACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTAc                                                                                                                                                                                                                             <  2:659248/217‑1 (MQ=255)
                                                                                                                                                  ccTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTg                                                                                                                                                                                                                        >  2:533221/1‑218 (MQ=255)
                                                                                                                                                     tGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGc                                                                                                                                                                                                                       >  2:771366/1‑216 (MQ=255)
                                                                                                                                                          ntccAGGCCCATTGCAGGTCCGGCGCCGCCGGGGGTGGAACCTATTACTATGGCATCGTAATCGGAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATGTTATCGCGGGTAAGTCAATGGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACg                                                                                                                                                                                                                <  2:723348/216‑1 (MQ=255)
                                                                                                                                                               ggCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTaaa                                                                                                                        >  1:978654/1‑301 (MQ=255)
                                                                                                                                                                gCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTan                                                                                                                                                                                                          >  2:679693/1‑217 (MQ=255)
                                                                                                                                                                gCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTa                                                                                                                                                                                                           <  2:216969/217‑1 (MQ=255)
                                                                                                                                                                gCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATAGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTa                                                                                                                                                                                                           <  2:424592/217‑1 (MQ=255)
                                                                                                                                                                   cATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTGCTTATACATAAAAGCAACAGAATGGTAACATTGTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATg                                                                                                                    <  1:394931/301‑1 (MQ=255)
                                                                                                                                                                           cgcCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTAt                                                                                                              <  1:679718/299‑1 (MQ=255)
                                                                                                                                                                           cgcCTTCGCCACCGGTGCCGGAACCTATTACTATGACAGAGTAATCGTAGGAATGTGGCATGGTAGTGTTTACCTGTTCTTATACATAAAAGCAACAAAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTg                                                                                                            <  1:621484/301‑1 (MQ=255)
                                                                                                                                                                                            ccGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGtt                                                                                                                                                                              <  2:543486/218‑1 (MQ=255)
                                                                                                                                                                                            ccGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGttcattca                                                                                            <  1:178285/300‑1 (MQ=255)
                                                                                                                                                                                                 aCCTATTACTTTGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGc                                                                                                                                                            <  1:971743/231‑1 (MQ=255)
                                                                                                                                                                                                 aCCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAgn                                                                                                                                                                         >  2:971695/1‑217 (MQ=255)
                                                                                                                                                                                                 aCCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAgn                                                                                                                                                                         >  2:952541/1‑217 (MQ=255)
                                                                                                                                                                                                 aCCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTaaa                                                                                                                                                  >  1:474135/1‑241 (MQ=255)
                                                                                                                                                                                                  ccTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACt                                                                                     <  1:842261/301‑1 (MQ=255)
                                                                                                                                                                                                      ttACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAg                                                                                 >  1:790825/1‑301 (MQ=255)
                                                                                                                                                                                                       tACTATGGAATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGc                                                                                <  1:533239/301‑1 (MQ=255)
                                                                                                                                                                                                         cTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTCCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGtt                                                                                                                                                                 >  2:301463/1‑218 (MQ=255)
                                                                                                                                                                                                         cTAGGGCATCGTAATCGTAGGAATGGTGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTGATCGCGGGTAAGTCAAGTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTGCTTACTTTTGCCATGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAg                                                                                 <  1:301471/298‑1 (MQ=255)
                                                                                                                                                                                                                       tcgtaGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTaa                                                                                                                                                   >  2:270891/1‑218 (MQ=255)
                                                                                                                                                                                                                       tcgtaGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTaa                                                                                                                                                   >  2:556767/1‑218 (MQ=255)
                                                                                                                                                                                                                        cgtaGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTaaa                                                                                                                                                  <  2:474119/218‑1 (MQ=255)
                                                                                                                                                                                                                             naaTGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAAc                                                                                                                                             <  2:746012/217‑1 (MQ=255)
                                                                                                                                                                                                                                  tgtcATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATTGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGtat                                                      <  1:398943/297‑1 (MQ=255)
                                                                                                                                                                                                                                  tgGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGtat                                                      <  1:952588/300‑1 (MQ=255)
                                                                                                                                                                                                                                    gCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTgg                                                                                   >  1:570934/1‑269 (MQ=255)
                                                                                                                                                                                                                                       tGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTCAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGa                                                                                                                                   >  2:636640/1‑218 (MQ=255)
                                                                                                                                                                                                                                       tGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGc                                                 >  1:250945/1‑300 (MQ=255)
                                                                                                                                                                                                                                            gtgCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAg                                                                                                                                                <  2:32322/198‑1 (MQ=255)
                                                                                                                                                                                                                                            gggCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTaaa                                                                                                                                                  >  1:939161/1‑198 (MQ=255)
                                                                                                                                                                                                                                            gggCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTaaa                                                                                                                                                  <  2:939114/198‑1 (MQ=255)
                                                                                                                                                                                                                                             ggCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTc                                                                                                                                                                                                                                 <  1:323513/118‑1 (MQ=255)
                                                                                                                                                                                                                                             ggCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTc                                                                                                                                                                                                                                 >  2:323504/1‑118 (MQ=255)
                                                                                                                                                                                                                                             ggCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCcagca                                            <  1:67336/299‑1 (MQ=255)
                                                                                                                                                                                                                                                    ntGTTCTTATAGATAAAAGCAACAGAAGGGTAACATTTTATCGCGGGTAAGTTAATTTATCCCCGTCATTTATCTGGCTATATCCTGAGCGTCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACTTAAAAGAACCGGCCAAAGAATTGCAGTAAata                                                                                                                      <  2:888020/217‑1 (MQ=255)
                                                                                                                                                                                                                                                    ntGTTCTTATACATAAAAGCAACAGATTGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAata                                                                                                                      <  2:797391/217‑1 (MQ=255)
                                                                                                                                                                                                                                                    cTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAata                                                                                                                      <  2:110121/218‑1 (MQ=255)
                                                                                                                                                                                                                                                      gTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGtt                                                                                                                  >  1:912528/1‑220 (MQ=255)
                                                                                                                                                                                                                                                      attCTTATACATAAAAGCAACAGAATGGTAACATTTCATCGCGGGTAAGTCAATTAATCCCCGTCCTTTATCTGGCTCCATCCTGTGCGGCCTTCGCTTTGTGTGTTTCTTACTTTTTCCCTGACGATTTATTGGATTTTTATCGACGATATTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCCGTAAataan                                                                                                                    >  2:810424/2‑216 (MQ=255)
                                                                                                                                                                                                                                                        tCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGtt                                                                                                                  <  2:309200/218‑1 (MQ=255)
                                                                                                                                                                                                                                                        ncTTATACATAAAAGCACCAGAATGGTAACAGTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACGTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGtt                                                                                                                  <  2:912482/217‑1 (MQ=255)
                                                                                                                                                                                                                                                         cTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGttn                                                                                                                 >  2:951841/1‑217 (MQ=255)
                                                                                                                                                                                                                                                         cTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAg                              >  1:568935/1‑301 (MQ=255)
                                                                                                                                                                                                                                                              aCATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTg                                                                                                            <  2:417739/218‑1 (MQ=255)
                                                                                                                                                                                                                                                                  aaaaGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTTGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGAGTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGAc                     <  1:402487/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                       cAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGCCGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCATTAAATATGTTTTATTGCGTTACCgn                                                                                                   >  2:840000/1‑217 (MQ=255)
                                                                                                                                                                                                                                                                       cAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGt                                                                                                   >  2:563367/1‑218 (MQ=255)
                                                                                                                                                                                                                                                                       aaaCAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGt                <  1:771395/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                        aaCAGAAAGGTAACAAATGATGTAGGGTAAGTCAATAGATACCAGTCATTTATATGGCTATATCATGAGCGGCCTTTGCTTGGTATGGTTCTTACTGTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTNCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTg               <  1:1020165/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                            gAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCGATCCAGTAt                                                                      <  1:698249/242‑1 (MQ=255)
                                                                                                                                                                                                                                                                               ttgTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACATTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTAt                                                                                                              <  2:449692/197‑1 (MQ=255)
                                                                                                                                                                                                                                                                               tGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACATTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTAt                                                                                                              >  1:449708/1‑199 (MQ=255)
                                                                                                                                                                                                                                                                                 gTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTCCGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACTGTTCATTCACaa                                                                                         >  2:156658/1‑218 (MQ=255)
                                                                                                                                                                                                                                                                                 gTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTa      >  1:153359/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                  tAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAAt                                                                                        <  2:484911/218‑1 (MQ=255)
                                                                                                                                                                                                                                                                                     cATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTTCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCGGGTTTCCGCTGTACGTAAACGAACCGGCAAAAGAATTGCAGTAAATATGTTTTCTTTCGTTTCCGTTCCTTCACCATa                                                                                       >  2:485468/1‑216 (MQ=255)
                                                                                                                                                                                                                                                                                     cATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTAACGTTCATTCACCATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGTAGTACTTTTGGCAAGGATTCAGACATCGCGATGGGCGTAagag  >  1:935879/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                      aTTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGttcattc                                                                                             <  1:439753/209‑1 (MQ=255)
                                                                                                                                                                                                                                                                                      aTTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGttcattc                                                                                             >  2:439738/1‑209 (MQ=255)
                                                                                                                                                                                                                                                                                       ttttATCTGGGGTAAGTCAATTGCTCTCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTgg                                                                                   <  2:570915/218‑1 (MQ=255)
                                                                                                                                                                                                                                                                                       ttttATCACGGTTAAGTCACTTTATCCCCGTCATTTATCTGGCTATTTCCTGGCCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTGTATTTGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTGTTATTGCGTTACCGTTCATTCACAATACTgg                                                                                   <  2:222182/218‑1 (MQ=255)
                                                                                                                                                                                                                                                                                        ttgacgtcgGGTAAGGTGATTGATCGCCGTCTTGTAGCGGGCTATAGCCTGAGGGGCATTTGCGTTGTCTGTTTATTACTTGTGCCCTGACGGGTTATTGGATTTTTATCGACGATACTCTCCGTATAAGCGTCAGGTTTCCGCTGTCCGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATAc                                                                                      <  2:744699/207‑1 (MQ=255)
                                                                                                                                                                                                                                                                                         ntATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAg                                                                                 <  2:739470/217‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                gggTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTCTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTACTTCACAATACTGGAGCAATCCa                                                                          >  2:394407/1‑218 (MQ=255)
                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                          
GTCGGGAATGGTCGCTGTAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGCCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAG  >  NC_000913/4160559‑4161136

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