Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,691,173 A→T 8.7% intergenic (‑75/+483) glrK ← / ← purL sensor protein kinase regulating glmY sRNA in two‑component system with response regulator GlrR/phosphoribosylformyl‑glycineamide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,691,1730AT 8.7% 195.1 / 4.6 64intergenic (‑75/+483)glrK/purLsensor protein kinase regulating glmY sRNA in two‑component system with response regulator GlrR/phosphoribosylformyl‑glycineamide synthetase
Reads supporting (aligned to +/- strand):  ref base A (29/29);  new base T (4/2);  total (33/31)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

TGACGGTAACTACGCTCCATCTCCAGCGCCGCGTTGGTCATCGCTTCACTGCGCCGGGCATCAATAAGCGTAGTGCGGTTAACCAGCGCCGCCTGATCGCTAAGCGCATTCAGGCTTTGCCAGGCTTGCCATGCCAGCACCAACAGGGGCAGCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATAAGTCGGTGAATGAGCCACTGGTTACTATTATGCAGTAACTGTGCCAATAAAGAAAATAGTTTGGTAACGTACTGATTATACGTTGCTTTGAGGGGTTTATGTCTCCTCCGCTGTTTGAATGTAAATCAGCCACTGTGTCGCTAAAAAGAGAC  >  NC_000913/2690878‑2691473
                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                                            
tGACGGTAACTACGCTCCATCTCCAGCGCCGCGTTGGTCATCGCTTCACTGCGCCGGGCATCAATAAGCGTAGTGCGGTTAACCAGCGCCGCCTGATCGCTAAGCGCATTCAGGCTTTGCCAGGCTTGCCATGCCAGCACCAACAGGGGCAGCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAATTCAAATACAACAAAGCCGGGTAATTCCCgg                                                                                                                                                                                                                                                                                                         >  2:434998/1‑301 (MQ=255)
      aaccTACGCTCCCTCTCCGGGGCGGCGTTTGTCATCGCTTACCTGCGCCGGGCAGTAAGAAGCGTAGGGCGGTTAACCAGCGCCGCCTGATCGCTAAGCGCAGTCAGGCTTTGCCAGGCTTGCCATGCCAGAACCAACAGGGGCAGCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGt                                                                                                                                                                                                                                                                                                   <  1:612027/298‑1 (MQ=255)
       tttgctttactgctctctcttccgtgtggttCATCGCTTCACTGCTGCGGGCAGCAGTCGGCGTAGTCAGGTTAATCTGCGCCACCTGATCGCTAAGCCCATTCAGGCTTTGCCAGGCTTGCCATAACAGCACCAACATGCGCAGCCGAATCGGCAATAATCCCAGCATTACCGGTTGTCGTAATGAGCGGGAAAAAACGGGCCCGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCAGATGCTACCCGGGGAGAAGTTCAGATACAACAAAGCCGGGTAATTCCCGGCTTTGtt                                                                                                                                                                                                                                                                                                  <  1:303417/271‑1 (MQ=255)
       aaCTACGCTCCATCTCCAGCGCCGCGTTGGTCATCGCTTCACTGCGCCGGGCATCAATAAGCGTAGTGCGGTTAACCAGCGCCGCCTGATCGCTAAGCGCATTCAGGCTTTGCCAGGCTTGCCATGCCAGCACCAACAGGGGCAGCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGtt                                                                                                                                                                                                                                                                                                  >  2:367823/1‑301 (MQ=255)
          tgactggagttcagacgtgtgctcttccgaTCTCTTCACTGCGCCGGGCATCAATAGCGTAGTGCGGTTAACCAGCGCCGCCTGATCGCTAAGCGCATTCAGGCTTTGCCAGGCTTGCCATGCCAGCACCAACAGGGGCAGCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATgg                                                                                                                                                                                                                                                                                              <  2:413535/272‑1 (MQ=255)
                           gccgcGTTGGTCATCGCTTCACTGCGCCGGGCATCAATAAGCGTAGTGCGGTTAACCAGCGCCGCCTGATCGCTAAGCGCATTCAGGCTTTGCCAGGCTTGCCATGCCAGCACCAACAGGGGCAGCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTa                                                                                                                                                                                                                                                                              >  2:528163/1‑301 (MQ=255)
                             cgcgTTGGTCATCGCTTCACTGCGCCGGGCATCAATAAGCGTAGTGCGGTTAACCAGCGCCGCCTGATCGCTAAGCGCATTCAGGCTTTGCCAGGCTTGCCATGCCAGCACCAACAGGGGCAGCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTAACCGGCTTTGATATGGAATAAGGCGGTGCCTAAc                                                                                                                                                                                                                                                                            >  2:134714/1‑301 (MQ=255)
                                 ggggTCATCGCTTCACTGCGCCGGGCATCAATAAGCGTAGTGCGGTTAACCAGCGCCGCCTGATCGCTAAGCGCATTCAGGCTTTGCCAGGCTTGCCATGCCAGCACCAACAGGGGCAGCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGTAATTCCCGGCTTTGttgtatctgaacttcccctcggttagca                                                                                                                                                                                                                                                                        <  1:434434/299‑27 (MQ=255)
                                   ggTCATCGCTTCACTGCGCCGGGCATCAATAAGCGTAGTGCGGTTAACCAGCGCCGCCTGATCGCTAAGCGCATTCAGGCTTTGCCAGGCTTGCCATGCCAGCACCAACAGGGGCAGCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGGGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGACGGGTAATTCCCGGCTTTGTTagaccggaagaaaacaagacagaacacc                                                                                                                                                                                                                                                                      >  2:303821/1‑274 (MQ=255)
                                        tCGCTTCACTGCGCCGGGCATCAATAAGCGTAGTGCGGTTAACCAGCGCCGCCTGATCGCTAAGCGCATTCAGGCTTTGCCAGGCTTGCCATGCCAGCACCAACAGGGGCAGCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCg                                                                                                                                                                                                                                                                 <  2:527316/301‑1 (MQ=255)
                                                  gcgcGGGGAGTCAATAAGGGTAGTGGGGTTAACCGGGGCCGCCGGATCGCTAAGCGCATTCAGGCTTGGCCAGGCTGGCCAGGCCAGAACCAACAGGGGCAGCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGGTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGtt                                                                                                                                                                                                                                                       <  1:54213/301‑1 (MQ=255)
                                                               aTAAGCGTAGTGCGGTTAACCAGCGCCGCCTGATCGCTAAGCGCATTCAGGCTTTGCCAGGCTTGCCATGCCAGCACCAACATGGGCAGCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGAAGGTTGATATAGCTAccc                                                                                                                                                                                                                                          >  1:144930/1‑299 (MQ=255)
                                                               aTAAGCGTAGTGCGGTTAACCAGCGCCGCCTGATCGCTAAGCGCATTCAGGCTTTGCCAGGCTTGCCATGCCAGCACCAACAGGGGCAGCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACCAGACATTTCGCCCGACGGTTGCTATAGCTAccc                                                                                                                                                                                                                                          >  1:420509/1‑299 (MQ=255)
                                                                 aaGCGTAGGTGGGTTTACCAGGGCCGCCTGATGGCTAAGCGCATTCAGGATTTGCCAGGCTGGCCATGCCAGCACCAACAGGGGCAGCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGGTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTg                                                                                                                                                                                                                                        <  1:527490/301‑1 (MQ=255)
                                                                  aGCGTAGTGCGGTTAACCAGCGCCGCCTGATCGCTAAGCGCATTCAGGCTTTGCCAGGCTTGCCATGCCAGCACCAACAGGGGCAGCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGa                                                                                                                                                                                                                                       >  2:158924/1‑301 (MQ=255)
                                                                        gTGCGGTTAACCAGCGCCGCCTGATCGCTAAGCGCATTCAGGCTTTGCCAGGCTTGCCATGCCAGCACCAACAGGGGCAGCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAg                                                                                                                                                                                                                                 <  2:84375/301‑1 (MQ=255)
                                                                              ttAACCAGCGCCGCCGGATCGCTAAGCTCAGTCAGGCTTTGCCAGGCTTGCCATGCCAGCACCAACAGGGGCAGCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAAGGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTg                                                                                                                                                                                                                           <  1:367339/301‑1 (MQ=255)
                                                                                 aCCAGCGCCGCCTGATCGCTAAGCGCATTCAGGCTTTGCCAGGCTTGCCATGCCAGCACCAACAGGGGCAGCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGaa                                                                                                                                                                                                                        >  2:118404/1‑300 (MQ=255)
                                                                                       gccgccTGATCGCTAAGCGCATTCAGGCTTTGCCAGGCTTGCCATGCCAGCACCAACAGGGGCAGCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACggcagg                                                                                                                                                                                                                  >  2:357023/1‑301 (MQ=255)
                                                                                         ctccTGATCTCTAATCGCATTCAGGCTTTGCCAGGCTTGCCATGCCAGCACCAACAGGGGCAGCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACggcaggca                                                                                                                                                                                                                <  2:81463/299‑1 (MQ=255)
                                                                                           ccTGATCGCTAAGCGCATTCAGGCTTTGCCAGGCTTGCCATGCCAGCACCAACAGGGGCAGCAGAATCAGCAAAAATGCCAGCATTACCAGTGGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCAcc                                                                                                                                                                                                              <  1:141001/301‑1 (MQ=255)
                                                                                                     aggcgcttgctggttttGCAATGTTTGCCATGCCAGGACCTACAGGGGCAGCAGAATCAGCAAAAATGCCAGAATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGGGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCg                                                                                                                                                                                                    <  1:148539/287‑1 (MQ=255)
                                                                                                      agggcATTCAGGCTTTGTCAGGCTTGCCATGCCAGCACCAACAGGGGCAGCATAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGt                                                                                                                                                                                                   <  1:158773/298‑1 (MQ=255)
                                                                                                                ggCTTTGCCAGGCTTGCCATGCCAGCACCAACAGGGGCAGCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTTTGCGTCATTCGTAtt                                                                                                                                                                                         >  2:328696/1‑301 (MQ=255)
                                                                                                                  gtgtGCAAGGGTTGTCATGTCAGCATCAAAAGGGGCAGCAGAATCAGGAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGGCGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTAtttt                                                                                                                                                                                       <  1:445621/298‑1 (MQ=255)
                                                                                                                               gccatgccaGCACCAACAGGGGCAGCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTTCGTCATTCGTATTTTATGTAACCCGTcc                                                                                                                                                                          >  1:531113/1‑301 (MQ=255)
                                                                                                                                       aGCACCAACAGGGGCAGCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGtcagccccttcgggacgt                                                                                                                                                                  >  1:410389/1‑284 (MQ=255)
                                                                                                                                               acgcgacgctgttccgatttaaaaTGCCAGCATTACCATTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGGTTCAAGGTGTTACTCTCGTCAGACGCGAAAAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATa                                                                                                                                                          <  1:465478/281‑1 (MQ=255)
                                                                                                                                                  gggCAGCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGtcagccccttcgggacgtgctacataaaa                                                                                                                                                       <  2:410952/301‑29 (MQ=255)
                                                                                                                                                      agcagAATCAGCAAAAATGCCAGCATTCCCAGTTGTCGTAATGAGCGGGGAAAAACGCGCCAGCGTTCCAAGGCGTTACTCTCGTCAGACGCTAATAGCCTGATGCTAACCGCGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCCACGTTTCGCCCGATGGTTTCTATAGCTACGCTGATATCAGACGTTGGACGACACGCCCCTTGTCGCAGCTCATCCGTATTTTTTGTAGCCAGTCCCGCAGGCGCTggccttcatcggt                                                                                                                                                   >  1:328478/1‑289 (MQ=255)
                                                                                                                                                               aGCAAAAATGCCAGCATTACCAGTTGTCTGAATGAGCGGGGAAAAACGGGCCAGCGTGTCAAGGGGTTACTCGCGTCAGACGCGAAGAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATAAGTCGGTGAATGAGcc                                                                                                                                          <  1:463602/301‑1 (MQ=255)
                                                                                                                                                               aGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTTCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATAAGTCGGTGAATGAAcc                                                                                                                                          >  2:187506/1‑301 (MQ=255)
                                                                                                                                                                   aaaaTGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATAAgaacggaagagcacacgtc                                                                                                                                      >  2:466042/1‑283 (MQ=255)
                                                                                                                                                                         ccAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGCATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATAAGTCGGTGAATGAGCCACTGGTTAca                                                                                                                                >  2:325599/1‑300 (MQ=255)
                                                                                                                                                                          cAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATAAGTCGGTGAATGAGCCACTGGTTACta                                                                                                                               >  2:366020/1‑301 (MQ=255)
                                                                                                                                                                             gaTTACCGGTTGTCGTGAGGGGGGGGGAAAAAGGGGCCGGCGTTTGAAGGTGTTACTGTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAGGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATAAGTCGGTGAATGAGCCACTGGTTACtatta                                                                                                                            <  1:631277/300‑1 (MQ=255)
                                                                                                                                                                                   aggtggtggaaagggggggggAAAAAGGGGCCAGGGTTTCATGGGGTTACGTGCGTAAGAGGCGAATAGCCTGATGCTAACCGGGGGGAAGTTCAGATACAACAAAGCCGGGAAGTACCCGGCTTTGGTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGGGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATAAGTCGGTGAATGAGCCACTGGTTACTATTATGCAGt                                                                                                                      <  1:325195/284‑1 (MQ=255)
                                                                                                                                                                                       gggcggaaTGGGCGGGGAAAAACGGGCCAGGGTTTGAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATAAGTCGGTGAATGAGCCACTGGTTACTATTATGCAGTAACt                                                                                                                  <  1:134589/295‑1 (MQ=255)
                                                                                                                                                                                            ccgttctcGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATAAGTCGGTGAATGAGCCACTGGTTACTATTATGCAGTAACTGTGcc                                                                                                             <  2:238790/294‑1 (MQ=255)
                                                                                                                                                                                               tGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATAAGTCGGTGAATGAGCCACTGGTTACTATTATGCAGTAACTGTGCCaaa                                                                                                          >  2:456939/1‑300 (MQ=255)
                                                                                                                                                                                                 aGCGTGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATAAGTCGGTGAATGAGCCACTGGTTACTATTATGCAGTAACTGTGCCAATaa                                                                                                        <  2:421072/301‑1 (MQ=255)
                                                                                                                                                                                                   gggggAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATAAGTCGGTGAATGAGCCACTGGTTACTATTATGCAGTAACTGTGCCAATAAAg                                                                                                      <  2:531786/300‑1 (MQ=255)
                                                                                                                                                                                                   cgtgtgctcttccgatctcGTTTCAAGTTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATAAGTCGGTGAATGTGCCACTGGTTACTATTATGCAGTAACTGTGCCAATAAAg                                                                                                      <  2:234737/283‑1 (MQ=255)
                                                                                                                                                                                                   cGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATAAGTCGGTGAATGAGACACTGGTTACTATTATGCAGTAACTGTGCCagaccgg                                                                                                      >  1:238523/1‑295 (MQ=255)
                                                                                                                                                                                                     gggAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATAAGTCGGTGAATGAGCCACTGGTTACTATTATGCAGTAACTGTGCCAATAAAGaa                                                                                                    <  2:145056/301‑1 (MQ=255)
                                                                                                                                                                                                         aaaaCGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTTCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATAAGGCGGTGAAAGAGCCACTGTTTACTATTAAGCAGTAACTGTGCCAATAAACAACATa                                                                                                >  1:483264/1‑301 (MQ=255)
                                                                                                                                                                                                                     cGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATAAGTCGGTGAATGTGCCACTGGTTAATATTATGAAGTAACTGTGCCAATCAaaaacacggaagaaggtcgg                                                                                    >  1:234471/1‑285 (MQ=255)
                                                                                                                                                                                                                      gTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATAAGTCGGTGAATGAGCCACTGGTTACTATTATGCAGTAACTGTGCCAAAAAAGAAAATAGTTAGGTAACGTa                                                                                   >  2:635485/1‑301 (MQ=255)
                                                                                                                                                                                                                        ttCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATAAGTCGGTGAATGAGCCACTGGTTACTATTATGCAGTAACTGTGCCAATAAAGAAAATAGTTTGGTAACGTACt                                                                                 >  1:137975/1‑301 (MQ=255)
                                                                                                                                                                                                                        ttCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATAAGTCGGTGAATGAGCCACTGGTTACTATTATGCAGTAACTGTGCCAATAAAGAAAATAGTTTGGTAACGTACt                                                                                 <  2:188322/301‑1 (MQ=255)
                                                                                                                                                                                                                            aGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGACGGGGCTGACATAAGTCGGTGAATGAGCCACTGGTTACTATTCTGAAGTACCTGTGCCAATAAAGAAAATAGTTTGGTAACGTACTGAtt                                                                             >  1:300745/1‑301 (MQ=255)
                                                                                                                                                                                                                                         ggcatgcgtgAATACGCTGATGCTAACGGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATAAGTCGGTGAATGAGCCACTGGTTACTATTATGCAGTAACTGTGCCAATAAAGAAAATAGTTTGGTAACGTACTGATTATACGTTGCTTTg                                                                <  1:187264/292‑1 (MQ=255)
                                                                                                                                                                                                                                                  gAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATCTAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCACGAAGGGGCTGACATAAGTCGGTGAATGAGCAACTGTTTACTATTATTCAGTAACTGTGCCACAAACGAAAATAGTATGGTACAGTACTTATTATACATTGCATTGCAGGGtttt                                                       >  1:451225/1‑300 (MQ=255)
                                                                                                                                                                                                                                                    aTAGCCTGAGGGTAACCGGGGGGAAGTTCAGATACAACAAAGCCGGGAATTACGCGGCTTTGTGATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATAAGTCGGTGAATGAGCCACTGGTTACTATTATGCAGTAACTGTGCCAATAAAGAAAATAGTTTGGTAACGTACTGATTATACGTTGCTTTGAGGGGTTTATg                                                     <  1:356539/301‑1 (MQ=255)
                                                                                                                                                                                                                                                      aGCCTGTTGTTAACCGAGGGGAAGTTCAGATACAACATAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATAAGTCGGTGAATGAGCCACTGGTTACTATTATGCAGTAACTGTGCCAATAAAGAAAATAGTTTGGTAACGTACTGATTATACGTTGCTTTGAGGGGTTTATGtc                                                   <  1:118280/301‑1 (MQ=255)
                                                                                                                                                                                                                                                        ccTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATAAGTCGGTGAATGAGCCACTGGTTACTATTATGCAGTAACTGTGCCAATAAAGAAAATAGTTTGGTAACGTACTGATTATACGTTGCTTTGAGGGGTTTATGtctc                                                 >  2:566892/1‑301 (MQ=255)
                                                                                                                                                                                                                                                        ccTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATAAGTCGGTGAATGAGCCACTGGTTACTATTATGCAGTAACTGTGCCAATAAAGAAAATAGTTTGGTAACGTACTGATTATACGTTGCTTTGAGGGGTTTATGtctc                                                 >  2:619426/1‑301 (MQ=255)
                                                                                                                                                                                                                                                           gATGCTAACCGGGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTGTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATAAGTCGGTGAATGAGCCACTGGTTACTATTATGCAGTAACTGTGCCAATAAAGAAAATAGTTTGGTAACGTACTGATTATACGTTGGTTTGAGGGGTTTATGTctcctc                                              <  2:328882/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                   ccGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATAAGTCGGTGAATGAGCCACTGGTTACTATTATGCAGTAACTGTGCCAATAAAGAAAATAGTTTGGTAACGTACTGATTATACGTTGCTTTGAGGGGTTTATGTCTCCTCCGCTGttt                                      >  1:259578/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                            aGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATAAGTCGGTGAATGAGCCACTGGTTACTATTATGCAGTAACTGTGCCAATAAAGAAAATAGTTTGGTAACGTACTGATTATACGTTGCTTTGAGGGGTTTATGTCTCCTCCGCTGTTTGAATGTaaa                             >  2:59238/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                acgagacaactAAGCCGGGAATTACCCGGCTTTTTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATAAGTCGGTGAATGAGCCACTGGTTACTATTATGCAGTAACTGTGCCAATAAAGAAAATAGTTTGGTAACGTACTGATTATACGTTGCTTTGAGGGGTTTATGTCTCCTCCGCTGTTTGAATGTAAATCAg                         <  1:365536/296‑1 (MQ=255)
                                                                                                                                                                                                                                                                                          aaaGCCGGGTAATTCCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGACGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATACGTCGGTGAATGAGCCACTGGTTACTATTATGCAGTAACTGTGCCAATAAAGAAACTAGTTTGGTAACGTACTGATTATACGTTGCTTTGAAGGGTTTATGTCTCCTCCGCTGTTTGAATGTAAATCAGCCACTttgtt               >  1:329680/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                          aaaGCCGGGTAATTCCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATAAGTCGGTGAACGAGCCACTGGTTACTATTATGCAGTAACTGTGCCAATAAAGAAAATAGTTTGGTAACGTACTGATTATACGTTGCTTTGAGGGGTTTATGTCTCCTCCGCTGTTTGAATGTAAATCAGCCACTgtgtt               >  1:97153/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                            aGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATAAGTCGGTGAATGAGCCACTGGTTACTATTATGCAGTAACTGTGCCAATAAAGAAAATAGTTTGGTAACGTACTGATTATACGTTGCTTTGAGGGGTTTATGTCTCCTCCGCTGTTTGATTGTAAATCAGCCACTGTGTCGc             >  2:93740/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                       aCCCGGCTTTGTTATGGAATAAGGCTGTTCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCGGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATAAGTCGGTGAATGAGCCACTGGTTACTATTATGCAGTAACTGTGCCAATAAAGAAAATAGTTTGGTAACGTACTGATTATACGTTGCTTTGAGGGGTTTATGTCTCCTCCGCTGTTTGAATGTAAATCAGCCACTGTGTCGCTAAAAAGAGAc  <  1:456375/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                       aCCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATAAGTCGGTGAATGAGCCACTGGTTACTATTATGCAGTAACTGTGCCAATAAAGAAAATAGTTTGGTAACGTACTGATTATACGTTGCTTTGAGGGGTTTATGTCTCCTCCGCTGTTTGAATGTAAATCAGCCACTGTGTCGCTAAAAAgagaa  >  1:405402/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                                            
TGACGGTAACTACGCTCCATCTCCAGCGCCGCGTTGGTCATCGCTTCACTGCGCCGGGCATCAATAAGCGTAGTGCGGTTAACCAGCGCCGCCTGATCGCTAAGCGCATTCAGGCTTTGCCAGGCTTGCCATGCCAGCACCAACAGGGGCAGCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATAAGTCGGTGAATGAGCCACTGGTTACTATTATGCAGTAACTGTGCCAATAAAGAAAATAGTTTGGTAACGTACTGATTATACGTTGCTTTGAGGGGTTTATGTCTCCTCCGCTGTTTGAATGTAAATCAGCCACTGTGTCGCTAAAAAGAGAC  >  NC_000913/2690878‑2691473

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