Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,987,442 Δ1 bp 5.9% intergenic (+14/+65) ftnB → / ← yecJ ferritin B, putative ferrous iron reservoir/DUF2766 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,987,4420T. 5.9% 297.8 / 12.1 85intergenic (+14/+65)ftnB/yecJferritin B, putative ferrous iron reservoir/DUF2766 family protein
Reads supporting (aligned to +/- strand):  ref base T (37/42);  new base . (2/3);  total (40/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.75e-02

CTGTTCCAAAAAACGATGGAAGAATACGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAATGATGATTCAACCGTCAATTTCCTGCGCGATCTGGAAAAAGAACAGCAGCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTGATAACCAACTGGCAAGCAA  >  NC_000913/1987180‑1987712
                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                               
cTGTTCCAAAAAACGATGGAAGAATACGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAATGATGATTCAACCGTCAATTTCCTGCGCGATCTGGAACAAGAACAGCAGCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCCGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCCGCTGCCTTCATCACCATCGGCGCTAATGCATTgcga                                                                                                                                                                                                                                                                    >  4:386770/1‑274 (MQ=255)
      cAAAAAACGATGGAAGAATACGAGCAACGTTCTAGTACGTAAGCACAGTTTGCCGATGAAGCGAAAGACCTGAATGATGATTCAACCGTCAATTTCCTGCGCGATCTGGAAAAAGAACAGCAGCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTCCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGa                                                                                                                                                                                                                                                              <  2:450429/275‑1 (MQ=255)
          aaaCGATGGAAGAATACGAGCAAAGTTCTAGTACGTTGGAACAGTTAGCCGATGAAGCGAAAGAACTGAATGATGAGTCAACCATCAATTTCCTGCGCGATCTGGAAAAAGAACAGCAGCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGatgat                                                                                                                                                                                                                                                          <  2:504564/275‑5 (MQ=255)
                             gCAACGTTCTAGTACGTTGTCACAGTTAGCCGATGAAGCGAAAGAACTGAATGATGATTCAACCGTCAATTTCCTGCGCGTTCTGGAAAAACAACAGCAGCATGATGGTCTGTTGCTGCAAACCTTTCTCGATAAAGTGCGCAGTGCTAAACTTGCGGGTCTGTTCCCTGTGCAGACCGACCAACTTGTGCTGAACGTCCTGTCACCCCAGCTTCATTGACCATCAACGGCGCCAATTCATTGCTCCGATGAAGGTTTTGAGAAACCGATGCCTc                                                                                                                                                                                                                                       >  4:316873/1‑275 (MQ=255)
                                aCGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAATGATGATTCAACCGTCAATTTCCTGCGCGATCTGGAAAAAGAACAGCAGCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTATCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCAt                                                                                                                                                                                                                                     <  4:125967/274‑1 (MQ=255)
                                             ttGGCACAGTTAGCCGATGAAGCGAAAGAACTGAATGATGATTCAACCGTCAATTTCCTGCGCGATCTGGAAAAAGAACAGCAGCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCTCAGTGCGACACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGCTGAAGGTTCTGAGAAACCGCTGCCTCATCTGTTTGAagcagc                                                                                                                                                                                                                       >  2:330304/1‑275 (MQ=255)
                                              tGGCACAGTTAGCCGATGAAGCGAAAGAACTGAATGATGATTCAACCGTCAATTTCCTGCGCGATCTGGAAAAAGAACAGCAGCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCg                                                                                                                                                                                                                      >  2:384570/1‑275 (MQ=255)
                                                          ccGATGAAGCGAAAGAACTGAATGATGATTCAACCGTCAATTTCCTGCGCGATCTGGAAAAAGAACAGCAGCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATg                                                                                                                                                                                                          >  4:515828/1‑275 (MQ=255)
                                                              tGAAGCGAAAGAACTGAATGATGATTCAACCGTCAATTTCCTGCGCGATCTGGAAAAAGAACAGCAGCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGTCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGAt                                                                                                                                                                                                      >  2:396273/1‑275 (MQ=255)
                                                              tGAAGCGAAAGAACTGAATGATGATTCAACCGTCAATTTCCTGCGCGATCTGGAAAAAGAACAGCAGCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGCATGTCGTGTCACACCAGCTGCATTGATCCTCATCGGCGATAATGCATTGCGCCGATGAAGTTTTTGCGAAACCGCTGCCCAATCTGTTTGCAGCAGCGGTTTTTTTAATGGGAt                                                                                                                                                                                                      >  4:343039/1‑275 (MQ=255)
                                                                aaGCGAAAGAACTGAATGATGATTCAACCGTCAATTTCCTGCGCGATCTGGAAAAAGAACAGCAGCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTc                                                                                                                                                                                                    <  2:544493/275‑1 (MQ=255)
                                                                    gAAAGAACTGAATGATGATTTAACCGGCAAATTCCTGCGCGATCTGGAAAAAGAACAGCAGCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCAcc                                                                                                                                                                                                 <  4:203158/274‑1 (MQ=255)
                                                                     aaagaaagtaaTGATGATTCAACCGTCAATTTCCTGCTCGATCTGGAAAAAGAACAGCAGCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCATACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCa                                                                                                                                                                                                   <  1:93395/262‑1 (MQ=255)
                                                                     aaaGAACTGAATGATGATTCAACCGTCAATTTCCTGCGCGATCTGGAAAAAGAACAGCAGCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCAcc                                                                                                                                                                                                 <  2:66110/273‑1 (MQ=255)
                                                                              aaTGATGATTCAACCGTCAATTTCCTGCGCGATCTGGAAAAAGAACAGCAGCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTa                                                                                                                                                                                       <  2:257078/274‑1 (MQ=255)
                                                                                     atTCAACCGTCAATTTCCTGCGCGATCTGGAAAAAGAACAGCAGCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAg                                                                                                                                                                               <  4:384569/275‑1 (MQ=255)
                                                                                     atTCAACCGTCAATTTCCTGCGCGATCTGGAAAAAGAACAGCAGCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGACATGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTGTGAGAAACCGCTGCCTAAACTGTTTGAAGCAGCGGTTGTTTTAATGGCATTCACCATGTGGGGTAAACtgtag                                                                                                                                                                               >  2:537272/1‑271 (MQ=255)
                                                                                           ccGTCAATTTCCTGCGCGATCTGGAAAAAGAACAGCAGCAAGATGGTCTGTTGCTGCAAACCATGCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATATCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAGTGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAAt                                                                                                                                                                         <  2:316874/275‑1 (MQ=255)
                                                                                                   ttCCTGCGCGATCTGGAAAAAGAACAGCAGCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGTCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGa                                                                                                                                                                 <  4:396272/275‑1 (MQ=255)
                                                                                                    tCCTGCGCGATCTGGAAAAAGAACAGCAGCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGa                                                                                                                                                                 <  2:293232/274‑1 (MQ=255)
                                                                                                    tCCTGCGCGATCTGGAAAAAGAACAGCAGCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGa                                                                                                                                                                 <  2:195183/274‑1 (MQ=255)
                                                                                                     ccTGCGCGACCTGGAAAAAGAACAGCATCATGATGGTCTGTTGCTGCAAATCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATg                                                                                                                                                               <  2:520597/275‑1 (MQ=255)
                                                                                                                tGGAAAAAGAACAGCAGCATGATTGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGAt                                                                                                                                                    <  2:160908/275‑1 (MQ=255)
                                                                                                                tGGAAAAAGAACAGCAGCATGATGGTCTGTTGCTGCAAACCATGCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGAt                                                                                                                                                    <  2:470134/275‑1 (MQ=255)
                                                                                                                   aaaaaGAACAGCAGCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCt                                                                                                                                                                                                                                                                                 >  7:122209/1‑148 (MQ=255)
                                                                                                                     aaaGAACAGCAGCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTaa                                                                                                                                                                                                                                                                               <  8:122209/148‑1 (MQ=255)
                                                                                                                            agcagcATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTg                                                                                                                                                                                                                                                                       <  5:603948/148‑1 (MQ=255)
                                                                                                                             gcagcaTGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAAACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTTCGCCGAAGAAGGTTTTGAGTAACCGCTGCCTCATCTGCTTGAATCAGCGGTTTTTTTATTGGGATTCACCCTGTGGGGTAAACTGGAGGTCAATAAGGCGGAGGGTTTTTAGGATTCACCGTTTGGTg                                                                                                                                       >  4:411809/1‑275 (MQ=255)
                                                                                                                                gcaTGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGAc                                                                                                                                     <  2:521736/274‑1 (MQ=255)
                                                                                                                                gaatgatgTGCTGTTGTGGTAAACCCTGTGTGCTGAAGGGCGCAGTCGGAAACTTGCGGGTAGGTGCCCTGTGCAGACCGACCAACATGTTCTGAAGGTCGTGGCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTTCCCGCATGGTGAc                                                                                                                                     <  4:119047/272‑1 (MQ=255)
                                                                                                                                 caTGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCAATGCATTAGCGCCg                                                                                                                                                                                                                                                                  <  7:402624/148‑1 (MQ=255)
                                                                                                                                       ggTCTGTTGCTGCAAACCATTCTTGTTGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTaa                                                                                                                                                                                                                                                                               >  7:438947/1‑130 (MQ=255)
                                                                                                                                       ggTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTaa                                                                                                                                                                                                                                                                               <  8:438947/130‑1 (MQ=255)
                                                                                                                                       ggTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGTGTAGACTTGAGTTCAATACGCTCGCTGGCTGTTTGGATTGCCCGCATGGTGACCGGGTCTg                                                                                                                             >  1:44983/1‑275 (MQ=255)
                                                                                                                                       ggTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTg                                                                                                                             <  2:380264/275‑1 (MQ=255)
                                                                                                                                           tGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTTCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTgg                                                                                                                                                                                           >  2:234883/1‑209 (MQ=255)
                                                                                                                                           tGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCGGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTgcggc                                                                                                                         >  4:503407/1‑275 (MQ=255)
                                                                                                                                             ttGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTg                                                                                                                                                                                                                                                      <  5:103502/148‑1 (MQ=255)
                                                                                                                                             ttGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCAATGCATTAGCGCCg                                                                                                                                                                                                                                                                  >  8:330021/1‑136 (MQ=255)
                                                                                                                                             ttGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCAATGCATTAGCGCCg                                                                                                                                                                                                                                                                  <  7:330021/136‑1 (MQ=255)
                                                                                                                                                   cAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTCATGGGATGCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGAt                                                                                                                 >  2:188054/1‑275 (MQ=255)
                                                                                                                                                   cAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCGCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATCCACCCTGCGGGGTAACCTGGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCGGCAGGGTGGCCGGGTCTGGGGTGGCGGGAt                                                                                                                 >  4:258239/1‑275 (MQ=255)
                                                                                                                                                      aCCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTaa                                                                                                                                                                                                                                                                               >  6:629272/1‑115 (MQ=255)
                                                                                                                                                      aCCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTaa                                                                                                                                                                                                                                                                               <  5:629272/115‑1 (MQ=255)
                                                                                                                                                      aCCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCAATGCAttagcgc                                                                                                                                                                                                                                                                    <  8:1074673/125‑6 (MQ=255)
                                                                                                                                                      aCCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCAATGCAttagcgc                                                                                                                                                                                                                                                                    >  7:1074673/1‑120 (MQ=255)
                                                                                                                                                       ccATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGaa                                                                                                                                                                                                                                                             <  6:113415/131‑1 (MQ=255)
                                                                                                                                                       ccATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGaa                                                                                                                                                                                                                                                             >  5:113415/1‑131 (MQ=255)
                                                                                                                                                        cATTATTGAGGTAGTGCGCGTTCCTAANNTTGCGGTTCTGTGCCCCTTGCAAACCTCCCNACATTTTCTGAATGTCTTTTCACACCAGCTGCATTGATCTTCATCTGCGCTAATGCATTGCGCCGATGAAGGTTNTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGt                                                                                                              <  3:57933/273‑1 (MQ=255)
                                                                                                                                                          ttCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGttt                                                                                                                                                                                                                                                        >  7:318080/1‑133 (MQ=255)
                                                                                                                                                          ttCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGttt                                                                                                                                                                                                                                                        <  8:318080/133‑1 (MQ=255)
                                                                                                                                                            cTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTc                                                                                                                                                                                                                                       >  6:467285/1‑148 (MQ=255)
                                                                                                                                                                atgaAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCg                                                                                                                                                                                                                                                                  <  7:405695/117‑1 (MQ=255)
                                                                                                                                                                atgaAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCg                                                                                                                                                                                                                                                                  <  5:805382/117‑1 (MQ=255)
                                                                                                                                                                atgaAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCg                                                                                                                                                                                                                                                                  >  6:805382/1‑117 (MQ=255)
                                                                                                                                                                atgaAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCg                                                                                                                                                                                                                                                                  >  8:405695/1‑117 (MQ=255)
                                                                                                                                                                  gaAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGAGGCGGGATGGTTAGTAAAGTCg                                                                                                   >  2:518780/1‑274 (MQ=255)
                                                                                                                                                                      tGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGa                                                                                                                                                                                                                             <  5:467285/148‑1 (MQ=255)
                                                                                                                                                                          cAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCa                                                                                                                                                                                                                                                                           <  7:539661/99‑1 (MQ=255)
                                                                                                                                                                          cAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCa                                                                                                                                                                                                                                                                           >  8:539661/1‑99 (MQ=255)
                                                                                                                                                                                      tGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGtt                                                                                                                                                                                                                                >  6:918851/1‑129 (MQ=255)
                                                                                                                                                                                      tGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGtt                                                                                                                                                                                                                                <  5:918851/129‑1 (MQ=255)
                                                                                                                                                                                       gCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTa                                                                                                         >  4:220405/1‑247 (MQ=255)
                                                                                                                                                                                           gTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGt                                                                                                     <  2:220406/247‑1 (MQ=255)
                                                                                                                                                                                              tgtgCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGcc                                                                                                                                                >  4:183080/1‑201 (MQ=255)
                                                                                                                                                                                                 gCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGttt                                                                                                                                                                                                                                                        <  5:888118/94‑1 (MQ=255)
                                                                                                                                                                                                 gCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGttt                                                                                                                                                                                                                                                        >  6:888118/1‑94 (MQ=255)
                                                                                                                                                                                                  cccTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCa                                                                                                                                            <  2:183081/201‑1 (MQ=255)
                                                                                                                                                                                                    cTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACTCCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGTGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAAc                                                                <  2:258240/275‑1 (MQ=255)
                                                                                                                                                                                                         cAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTaa                                                                                                                                                                                                                                                                               >  7:780705/1‑64 (MQ=255)
                                                                                                                                                                                                         cAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTaa                                                                                                                                                                                                                                                                               <  8:780705/64‑1 (MQ=255)
                                                                                                                                                                                                            accgaccAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACGTTTTTTCACGAACTTCACCgg                                                        >  1:36633/1‑275 (MQ=255)
                                                                                                                                                                                                                 ccAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGACCTTCAACGGGcgcg                                                   >  2:259063/1‑275 (MQ=255)
                                                                                                                                                                                                                        gTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACAt                                            <  2:515829/275‑1 (MQ=255)
                                                                                                                                                                                                                                      gtCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAg                               >  2:228650/1‑274 (MQ=255)
                                                                                                                                                                                                                                        cacaCCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCtt                                                                                                                                                           >  5:756129/1‑148 (MQ=255)
                                                                                                                                                                                                                                                cTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTgc                                                                                                                            >  2:490170/1‑171 (MQ=255)
                                                                                                                                                                                                                                                 tGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCCGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTATACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCTCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGACTTTGTTTGa                    >  2:368195/1‑274 (MQ=255)
                                                                                                                                                                                                                                                    aTTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTgcggcg                                                                                                                        <  4:490169/171‑1 (MQ=255)
                                                                                                                                                                                                                                                     ttGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTg                                                                                                                                                                                                                                           <  7:411061/55‑1 (MQ=255)
                                                                                                                                                                                                                                                     ttGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTg                                                                                                                                                                                                                                           >  8:411061/1‑55 (MQ=255)
                                                                                                                                                                                                                                                        atcatcatcGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTGTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTGATAACCaa             >  2:224827/1‑274 (MQ=255)
                                                                                                                                                                                                                                                          catcatcGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTa                                                                                                                                                                                                             <  6:40072/80‑1 (MQ=255)
                                                                                                                                                                                                                                                          catcatcGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTa                                                                                                                                                                                                             >  5:40072/1‑80 (MQ=255)
                                                                                                                                                                                                                                                                cGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCAAGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTGATAACCAACTGgcaagc    <  2:386771/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                  gcgcTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAACGTCGATATTCTTCATCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTAAATCGAGAACATCCAGAATTTGGTTGATAACAAACTGgcaagcaa  >  1:37412/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                      tAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTg                                                                                                                             <  6:756129/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                               
CTGTTCCAAAAAACGATGGAAGAATACGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAATGATGATTCAACCGTCAATTTCCTGCGCGATCTGGAAAAAGAACAGCAGCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTGATAACCAACTGGCAAGCAA  >  NC_000913/1987180‑1987712

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