Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 72,493 T→C V145V (GTA→GTG)  thiQ ← thiamine ABC transporter ATP binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091372,4930TC97.5% 345.2 / ‑2.4 118V145V (GTA→GTG) thiQthiamine ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base T (0/3);  new base C (64/51);  total (64/54)
Fisher's exact test for biased strand distribution p-value = 9.29e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.70e-01

Only 100 of 118 total aligned reads displayed.

GCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGCGCCGCATCTTCCACGCTGTGCGACACCATCAATAGCGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAACATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGTACCAGACAACGCGCTAACGCCACTCGCTGTCGCTGACCGCCGGAAAGCTCGCCCGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCATAGTTGT  >  NC_000913/72251‑72738
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     
gctgAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCCCCAGAc                                                                                                                                                                                                                                                >  2:733126/4‑250 (MQ=255)
 ctgAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCGTCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACa                                                                                                                                                                                                                                               >  1:746210/3‑250 (MQ=255)
     aGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCTTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAAcgcg                                                                                                                                                                                                                                          >  1:606499/1‑251 (MQ=255)
       cACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCt                                                                                                                                                                                                                                        <  1:295812/251‑1 (MQ=255)
        aCTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCt                                                                                                                                                                                                                                        >  2:22540/1‑250 (MQ=255)
           cGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACg                                                                                                                                                                                                                                    <  1:747964/251‑1 (MQ=255)
           cGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCAATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTAACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGATCGCGCACCAGACAACGCGCTAACg                                                                                                                                                                                                                                    >  1:259762/1‑251 (MQ=255)
            cccTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGc                                                                                                                                                                                                                                   <  1:113022/250‑1 (MQ=255)
                tACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGAACCGCGCGGCATCTTCCACAATGTGCGCCACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCACACCAGCGTCAGCATCTCATGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCTCGCAATCAAATCGGCTGTTCGCGCACCAGAAAACGCTCTAACGccccg                                                                                                                                                                                                                               >  2:90919/1‑248 (MQ=255)
                 aCCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGAAAACGCGCTAACGCCAcgc                                                                                                                                                                                                                              >  1:533411/1‑249 (MQ=255)
                  ccGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACAAGACAACGCGCTAACGCCAcccg                                                                                                                                                                                                                             >  2:841057/1‑248 (MQ=255)
                    gCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGCCAACGCGCTAACGCCACgcgct                                                                                                                                                                                                                           >  1:977775/1‑251 (MQ=255)
                    catCAACAACGCATTTGTCTTGCCCTGCCAGGCGGTGCGCCCGGCGGCTATCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACCCCATCAATAACGTCATGTTTTGCTGCTGTCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTGCTCGCACCAGACAACGCGCTAACGCCACgcgct                                                                                                                                                                                                                           <  2:357921/249‑1 (MQ=255)
                      tCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACAACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACgcgctgt                                                                                                                                                                                                                         >  1:487185/1‑251 (MQ=255)
                            actcagtggtttgaCCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGa                                                                                                                                                                                                                   <  2:34842/238‑1 (MQ=255)
                             cTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGAc                                                                                                                                                                                                                  >  1:672346/1‑251 (MQ=255)
                             cTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGAc                                                                                                                                                                                                                  >  1:210375/1‑251 (MQ=255)
                               cATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAATTCGGCTGTTCGCGCACCCGACAACGCGCTAACGCCCCGCGCTGTCGCTGAccg                                                                                                                                                                                                                >  2:930543/1‑251 (MQ=255)
                                 ttGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGAccgcc                                                                                                                                                                                                              >  1:843879/1‑251 (MQ=255)
                                      cttACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTACCGCCACGCCCTGTCGCTGACCGCCGGaa                                                                                                                                                                                                          >  2:881143/3‑250 (MQ=255)
                                          cccTGCCAGGCGATTCGCCCGCCGGCGACCACCACCGAGCGCGTGGCGATCCGCTCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTc                                                                                                                                                                                                     <  2:288321/251‑1 (MQ=255)
                                               ccAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGTCCGCCGGAAAGCTCGccag                                                                                                                                                                                                <  1:640565/251‑3 (MQ=255)
                                                     gctGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGAcc                                                                                                                                                                                                                 <  2:710246/226‑1 (MQ=255)
                                                          gcCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATCACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCCGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGCGCAATACCATCGGCTGTTCGCACACCAGCCACCGCGCTAAAGCCACGCGCTGTCGCTGAACGCCGGaaatctaccacggtaaccaccgc                                                                                                                                                                                     >  2:910172/1‑231 (MQ=255)
                                                             cGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCAtt                                                                                                                                                                                  <  1:815783/251‑1 (MQ=255)
                                                              gTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGGGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCAGCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCAGCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTa                                                                                                                                                                                 <  2:670254/251‑1 (MQ=255)
                                                              gTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGCCGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCAtt                                                                                                                                                                                  >  1:983528/1‑250 (MQ=255)
                                                               ccggcgaccgccacCGAGCGCGTGGCGAGCCGCGCGGCATCTTCCACCCTGTGTGACACCATCAATAACGTCATTTTTGGCTGCTGGCAGCTCGTGCTCACCAGCTCCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCAGCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTaa                                                                                                                                                                                <  2:352093/239‑1 (MQ=255)
                                                               ccGGCTACTGCTACCGAGCGCGTGGTGTTCCGCGCCGCATCTTCCACGATGTGTGTCCCCACCAATATCGTCATTTTTTGCTGCTGGCATCTCGTGCTCACCAGCGTCAACAACTCCTGACGTAACGCCTGATCGTGCTCAGAGAACGGTTCATGGAGCAATAAAATCGGCTGGTCGCGTACCAGTGAACGCGCTTGTTTCTCTCGCTTTCGCTGACCTCCGGAAAGCTCGCCCGGTAACCGCGCCAtt                                                                                                                                                                                  <  1:399201/248‑1 (MQ=255)
                                                                          cacCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATcc                                                                                                                                                                     <  1:817435/250‑1 (MQ=255)
                                                                           acCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCCTCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGCGCAATAAAATCGGCGGTTCGCGCAACAGACAAAGCGCTAACGCCACGCGCTGGCGCTGACCGACCGCAAGTCCGCCCAGGAACCAGGCCATTAAATTATCAATccc                                                                                                                                                                    >  2:511756/1‑251 (MQ=255)
                                                                           acCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATccc                                                                                                                                                                    >  1:207389/1‑251 (MQ=255)
                                                                            ccGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATcccc                                                                                                                                                                   <  1:402227/251‑1 (MQ=255)
                                                                             cGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCACTAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTACCGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCCAACGCCACGAGCTGTCGCTGAACGCCGGTAAGCTCGCCAGGTAACCGCGCCATTAATATATCAAAccccc                                                                                                                                                                  >  2:215156/1‑250 (MQ=255)
                                                                                gcgcGTGGCGAATCGCGCGGCATCTGCCACACTGTTGGACACTATCATGAACGTCCTTTTTGGCTGCTGGCAGCTCGGGCTCACCAGCGTCAGCCTCCCCTGACGTAACGCCGGATCGTGCGCAGAGAACGGTTCACCGAGCAATAAAATCGGCTTTTCGCGCACCAGACAACGCGCTACCGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCt                                                                                                                                                               <  2:375959/251‑1 (MQ=255)
                                                                                     tGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTggcgg                                                                                                                                                          >  1:170745/1‑251 (MQ=255)
                                                                                      ggCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTggcggg                                                                                                                                                         >  1:262726/1‑251 (MQ=255)
                                                                                                   gcATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCAt                                                                                                                                            >  1:368970/1‑251 (MQ=255)
                                                                                                   gcATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCAt                                                                                                                                            >  1:580471/1‑251 (MQ=255)
                                                                                                   gcATCTTCCACACTGTGTGACACCATCAATAACGTCATGTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCAt                                                                                                                                            <  2:948025/251‑1 (MQ=255)
                                                                                                    cATCTTCCACACTGTGTGACACCATCAATACCGTCATTTTTTGCTGCTGCCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGAGCAGAGCACGGTTCATCGAGAAATACAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTACATTATCAATCCCCATCTGGCGGGCGATAGCGTGCAtt                                                                                                                                           >  1:134599/1‑251 (MQ=255)
                                                                                                       cTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTc                                                                                                                                        >  1:222964/1‑251 (MQ=255)
                                                                                                           cccacTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGGCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCct                                                                                                                                    <  1:4434/247‑1 (MQ=255)
                                                                                                            cctctgtttgaccccATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCTCCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCTCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCTCACCAGCCAACGCGCTAACGCCACGCGCTTTCGCTGACCGCCGTAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCctg                                                                                                                                   <  1:46372/238‑1 (MQ=255)
                                                                                                                 gtgCGACACCATCAATAGCGTCATTTTTTGCGGCTGGCAGCTCGTGCTCACCAGCGTCATCATCTCCTGACGTAACGCCGGATCGAGCGCAGACAACGGTTCATCGAGCAATAAAAGCGGCTGTTCGCGTACCAGACAACGCGCTAACGCCACTCGCTGTCGCTGACCGCCGGAAAGCTCGCCCGGGAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCCCctgctg                                                                                                                                <  2:520317/249‑1 (MQ=255)
                                                                                                                    ggACACCATCAATAACGGCATTTTTTGCTGCTGGCAGCCCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCACGCCAGGTAACCGCGCCATTAAATTATCAATCCCCACCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCcg                                                                                                                           <  1:583059/250‑1 (MQ=255)
                                                                                                                        acCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGtt                                                                                                                       <  1:500417/251‑1 (MQ=255)
                                                                                                                             cAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGATCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACTCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGtt                                                                                                                  <  1:467689/251‑1 (MQ=255)
                                                                                                                               ataaactcatttttgtctgcctgCAGCCCGTGCTCACCTGCGTCAGCATCTCCGGACGTAACGCCGGACCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCCCCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTc                                                                                                                <  2:358514/229‑1 (MQ=255)
                                                                                                                                      cATTTTTGGCGGCGGGCGGCCCGTGCTCCCCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAACCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCTCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTGCAATccc                                                                                                          <  2:681763/250‑1 (MQ=255)
                                                                                                                                        ttttttGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCCGAGAACGGTTCATCGCGCAATAAAATCGGCTGTTCGCGCACCAGCCCACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCCTTACATTATCAATCCCCCTCTGGCGGGCGTAAGCGTGCAATTTCTCCTGCTGTGCCGCGTTCCGTTTCAATCCCggg                                                                                                       >  2:405773/1‑250 (MQ=255)
                                                                                                                                           tttGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGCttt                                                                                                    >  1:139372/1‑251 (MQ=255)
                                                                                                                                                 gctgGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGcc                                                                                                <  1:248287/249‑1 (MQ=255)
                                                                                                                                                   tgtcAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGcccc                                                                                              >  2:242747/4‑249 (MQ=255)
                                                                                                                                                        gCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCCTCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCg                                                                                        >  2:50923/1‑250 (MQ=255)
                                                                                                                                                            gTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCCGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGCAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCCTCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTACATCCCGGCTTTACCCCCATCCCGATGtt                                                                                   >  2:114224/1‑251 (MQ=255)
                                                                                                                                                             tGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTc                                                                                  >  1:302217/1‑251 (MQ=255)
                                                                                                                                                              gCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCt                                                                                 >  1:331123/1‑251 (MQ=255)
                                                                                                                                                               cTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCCGAGAACGGTTCATCGCGCAATAAAATCGGCTGTTCGCGCACCCGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTATCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCtg                                                                                >  2:857216/1‑251 (MQ=255)
                                                                                                                                                                 ctccAGCTTCAGCATCACCTTACGTAACGCCGGATCGAGCGCAGATAACGGTTCAACCAGCAATAAAATCGGCTTGTCGCGCACCAGACAACCCGCGAACGCCACGCGCTGTGGCTGACCGCCGGAAAGCGCGCCAGGTAACCGCGCCATTAAATTATCAATCCTCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCTGATGTTCtg                                                                                <  2:501275/247‑1 (MQ=255)
                                                                                                                                                                 cACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCtgcg                                                                              <  1:944069/251‑3 (MQ=255)
                                                                                                                                                                      gCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAAcc                                                                                                                                                                                          >  2:293317/1‑138 (MQ=255)
                                                                                                                                                                       cGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCg                                                                        <  1:134390/251‑1 (MQ=255)
                                                                                                                                                                        gTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGt                                                                       >  1:147182/1‑251 (MQ=255)
                                                                                                                                                                             cATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGACAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGCCCGTCAGGt                                                                  >  2:543253/1‑251 (MQ=255)
                                                                                                                                                                              aTCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCACGGATTTAGACCCAGCCCGATGTTCTGCGCGACCGTCAgggg                                                                 >  1:281128/1‑249 (MQ=255)
                                                                                                                                                                                 cccTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACACCGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTgact                                                              <  1:621771/250‑4 (MQ=255)
                                                                                                                                                                                        gTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAgg                                                       >  1:391861/1‑251 (MQ=255)
                                                                                                                                                                                        gTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAgg                                                       <  1:627829/251‑1 (MQ=255)
                                                                                                                                                                                        gTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGAGACTGAACAgg                                                       >  1:599934/1‑251 (MQ=255)
                                                                                                                                                                                         tAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGgt                                                      >  1:46414/1‑251 (MQ=255)
                                                                                                                                                                                             gCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAAGCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGgttgtt                                                  <  1:918389/251‑1 (MQ=255)
                                                                                                                                                                                               cGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGgttgtttt                                                <  1:157907/251‑3 (MQ=255)
                                                                                                                                                                                                   gcGAGCGCAGAGAACGGTTCACCGAGAAATAAAATCGGCTGTTCGCGCACCAGACAACGCTCTAACGCCACGCGCTGTCGCTGTCCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTg                                            <  2:798219/250‑1 (MQ=255)
                                                                                                                                                                                                             agaACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATc                                  >  2:19534/1‑251 (MQ=255)
                                                                                                                                                                                                              gaACGGTTCACCGACCAATAAAACCGGCGGTTCGCGCACCAGACAACGCGCTAACGCCCCGCGCTGTCGCTGCCCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCg                                 <  2:42177/251‑1 (MQ=255)
                                                                                                                                                                                                               cgcacttccgatataaaTAAAAGCGGCTGGTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTAGCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATGTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTGGTTTTCCTGAAACAGCATCGa                                <  2:394046/236‑1 (MQ=255)
                                                                                                                                                                                                                 cGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGaca                              <  1:900230/251‑1 (MQ=255)
                                                                                                                                                                                                                 cGGTTCATCGAGCAATAAAATCGGATGTTCGCGCACCAGACAACGCGCTAACGCCCCGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGCACGGGTTGTTTTACTGAAACAGCATCGacc                              >  2:712461/1‑250 (MQ=255)
                                                                                                                                                                                                                   gTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGCAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGCTGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACAcc                            >  1:651348/1‑251 (MQ=255)
                                                                                                                                                                                                                    ttCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCAACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGAGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGATGGTTTCATGAAACAGCATCGACACCg                           >  1:865158/1‑251 (MQ=255)
                                                                                                                                                                                                                     tCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCgg                          >  1:242983/1‑251 (MQ=255)
                                                                                                                                                                                                                      cATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCTGGTAACAGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCgg                          >  1:671123/1‑250 (MQ=255)
                                                                                                                                                                                                                      cATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACCCCgg                          >  2:210217/1‑250 (MQ=255)
                                                                                                                                                                                                                       aTCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGAc                        >  1:554357/1‑251 (MQ=255)
                                                                                                                                                                                                                       aTCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGAc                        <  1:733126/251‑1 (MQ=255)
                                                                                                                                                                                                                         cGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCAATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTATCCTGCAACAGCATCGACACCgggagg                      >  2:663796/1‑247 (MQ=255)
                                                                                                                                                                                                                         agAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGACGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCAATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGAGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCAAGCCCGATGTTCTTCGCGACCGTCAGGTGCCTGGACAGGTTGTTTTCCTGAAACAGCATCGAACCCGGACgg                      >  2:663798/2‑251 (MQ=255)
                                                                                                                                                                                                                             cAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGACACAGCATCGa                                >  1:394046/1‑237 (MQ=255)
                                                                                                                                                                                                                                tAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGGTTTCCTGAAACAGCATCGACACCGGACGGCGTGAcg               >  1:328672/1‑251 (MQ=255)
                                                                                                                                                                                                                                 aaatccGGCTGTCCGCCCACCAGACACCGCGCTAACACCTCGCTCTGGCGCTGTCCGCCGGAAAGCTCGTCAGGTATCCGCGCAATTAAATGATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTACTCCGGCGGTGCCGCGTTCAGTTTCAATCCCGGATTTAGTCCCAGCCCGATGTTCGGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCTACGCCGGACGGCGTGAcgg              <  2:422089/246‑1 (MQ=255)
                                                                                                                                                                                                                                 aaaaTCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCTCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGAcg               <  1:15569/250‑1 (MQ=255)
                                                                                                                                                                                                                                  aaaTCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGAcggc             >  1:160748/1‑251 (MQ=255)
                                                                                                                                                                                                                                  aaaTCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGAcggc             <  1:719161/251‑1 (MQ=255)
                                                                                                                                                                                                                                   aaTCGCCTGGTCGCGCCCCAGACAACGCGCTAACGCCCCGCGCTGTCGCTGACCGCCGGAAAGCGCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTTCATTTTCTCCTGCGGTGCCGCGTTCAGTGTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACTCCGGACGGCGTGAcggcg            <  2:628311/251‑1 (MQ=255)
                                                                                                                                                                                                                                         cTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCACTCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACggcggcgtcg      >  1:754285/1‑247 (MQ=255)
                                                                                                                                                                                                                                           gTTCGCGCACCAGACAACGCGCTAACGAAACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACggcggcgtcgtg    >  1:880136/1‑245 (MQ=255)
                                                                                                                                                                                                                                             gcgcgCACCAGACAACGCGCTAACGCCACGCGCGGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCAGCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACggcggcgtcgtggt  <  2:53775/250‑9 (MQ=255)
                                                                                                                                                                                                                                              cgcgCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACggcggcgtcgtggt  <  1:881143/250‑9 (MQ=255)
                                                                                                                                                                                                                                              cgcgCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACggcggcgtcgtggt  <  1:339767/250‑9 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     
GCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGCGCCGCATCTTCCACGCTGTGCGACACCATCAATAGCGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAACATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGTACCAGACAACGCGCTAACGCCACTCGCTGTCGCTGACCGCCGGAAAGCTCGCCCGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCATAGTTGT  >  NC_000913/72251‑72738

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑