Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,364,798 C→T G82D (GGC→GAC)  dsbD ← thiol‑disulfide exchange protein DsbD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,364,7980CT97.9% 319.4 / ‑4.6 95G82D (GGC→GAC) dsbDthiol‑disulfide exchange protein DsbD
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (40/52);  total (41/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCAATTGCGCGGTGGGCTGCTCTTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATG  >  NC_000913/4364552‑4365043
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
ggCAATTGCGCGGTGGGCTGCTCTTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTg                                                                                                                                                                                                                                                   <  2:287720/251‑1 (MQ=255)
ggCAATTGCGCGGTGGGCTGCTCTTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTg                                                                                                                                                                                                                                                   >  2:159348/1‑251 (MQ=255)
          cgGTGGGCTGCTCTTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGc                                                                                                                                                                                                                                         >  2:2789/1‑251 (MQ=255)
            gtcggCTGCTCTTGCTGCGGAACAGACACAGGCTTCGGCGCTGCGTTGTTGGCGACCTCTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAAACGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCCCCTGGTTGAGGGTGACTGGAAGCGTCAGCCGATCGCGGTAAAGCGCGCTTCTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGAAGCTgaaa                                                                                                                                                                                                                                       <  2:186594/248‑4 (MQ=255)
             tGGGCTGCTCTTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACg                                                                                                                                                                                                                                      >  1:546047/1‑251 (MQ=255)
                 ctgctCTTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCgg                                                                                                                                                                                                                                  >  2:719261/1‑251 (MQ=255)
                        tgctgcGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAAtttt                                                                                                                                                                                                                           >  1:361444/1‑251 (MQ=255)
                         gctgcGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAAtttt                                                                                                                                                                                                                           <  1:655396/250‑1 (MQ=255)
                                aaCAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTGCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGtt                                                                                                                                                                                                                   <  1:641720/251‑1 (MQ=255)
                                  cAGACACAGGCTGTGGCGCTGCGTTGTTGGCGCCCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGAGGAAGCTGCCCGTCGGCAATTTTCGCGTGTTcc                                                                                                                                                                                                                 >  2:647061/1‑251 (MQ=255)
                                   aGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCg                                                                                                                                                                                                                <  1:567646/251‑1 (MQ=255)
                                   aGACACAGGCTGTGGCGCTGCGTTGATGGCGACCACGTCGCATAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACACGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCTCGTCGGCAATTTTCGCGTGTTCCg                                                                                                                                                                                                                <  2:236636/251‑1 (MQ=255)
                                        caGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTa                                                                                                                                                                                                           <  2:257205/251‑1 (MQ=255)
                                          ggCTTTGGCGCTGCGTTGTGGTCGACCACGTCGCTTAACGGAACGGTTCTGGTGTCTGGCGGATAACAGAATCCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCGTTTGCCGTAAAACTCATCTTCATGCCAGACGTATTGCGGCAGCTGCACGTCGGCAATTTGCGCGTGGCCCGTCGTAAt                                                                                                                                                                                                         <  2:301021/251‑1 (MQ=255)
                                           gCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATg                                                                                                                                                                                                        >  1:25130/1‑251 (MQ=255)
                                             tgtgGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCg                                                                                                                                                                                                      >  1:559521/1‑251 (MQ=255)
                                              gtgGCGCTGCGTTGTTGTCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCgg                                                                                                                                                                                                     <  2:543725/251‑1 (MQ=255)
                                                ggCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGAt                                                                                                                                                                                                   <  1:21792/251‑1 (MQ=255)
                                                 gcgcTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTACTGCggctc                                                                                                                                                                                                  >  2:414219/1‑248 (MQ=255)
                                                   gcTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTGCGCGTGTTCCGGCGTAATGCGGATCTg                                                                                                                                                                                                <  2:73637/251‑1 (MQ=255)
                                                    cTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGt                                                                                                                                                                                               >  2:178687/1‑251 (MQ=255)
                                                      gCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCGTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGtt                                                                                                                                                                                              <  1:342198/250‑1 (MQ=255)
                                                            ttgGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGGTTACGGt                                                                                                                                                                                        >  2:599262/1‑250 (MQ=255)
                                                            ttgGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAAATCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGAGGATCTGTTTCCggc                                                                                                                                                                                        >  2:146124/1‑249 (MQ=255)
                                                                        tCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAAc                                                                                                                                                                            <  2:572000/250‑1 (MQ=255)
                                                                             tAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTc                                                                                                                                                                        <  1:439223/249‑1 (MQ=255)
                                                                                 ggaacggTTTGGTTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGAt                                                                                                                                                                  <  1:63203/251‑1 (MQ=255)
                                                                                                  tGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGGGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGGAGTAACCGTCTTTGATCTGCCAGGTCAGATTac                                                                                                                                                 >  2:349553/1‑250 (MQ=255)
                                                                                                   ggCGGATAACAGAAACCGGCATCAGCACAGGCCTGGTATGTGACAGTTAACGTCGCTCCCGCACTCTCCTGGTTGATGGTGACGGTAATCGTCAGCCGATCGCGGTAAAGCTCGCTGTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACTTCGGCAATTTTCGCGTGTTCCGGCGTAATGCTGATCTGTTTACGGTAGAGGGAGTAACCGTCTTTGATCTGCCAGGTCAGAtt                                                                                                                                                   <  2:680526/248‑1 (MQ=255)
                                                                                                   ggCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGAtt                                                                                                                                                   <  1:207117/248‑1 (MQ=255)
                                                                                                   ggCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTa                                                                                                                                                  <  2:269148/249‑1 (MQ=255)
                                                                                                     cGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGt                                                                                                                                              <  1:345851/251‑1 (MQ=255)
                                                                                                       gATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCa                                                                                                                                            <  2:579436/251‑1 (MQ=255)
                                                                                                       gATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGAGCGCTGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAAGTTTCGCGTGTTCCGGAGTAATGCGGATCTGTTTACGGTAGATGGAGTAACCGTCTTTGATCTGCCAGGTCAGAGTAAGGt                                                                                                                                              <  2:534631/249‑1 (MQ=255)
                                                                                                            cAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTg                                                                                                                                       <  1:468893/251‑1 (MQ=255)
                                                                                                               aaaCCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGtt                                                                                                                                    <  1:42866/251‑1 (MQ=255)
                                                                                                                  ccGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTg                                                                                                                                 <  1:121172/251‑1 (MQ=255)
                                                                                                                          aGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAAt                                                                                                                         <  2:163899/251‑1 (MQ=255)
                                                                                                                            cacaGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCa                                                                                                                       <  1:719261/251‑1 (MQ=255)
                                                                                                                                gCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCGGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAg                                                                                                                   <  1:437263/251‑1 (MQ=255)
                                                                                                                                  ccTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCa                                                                                                                 <  2:252946/251‑1 (MQ=255)
                                                                                                                                     ggTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAg                                                                                                              <  1:64873/251‑1 (MQ=255)
                                                                                                                                              cAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATcc                                                                                                     <  1:244810/251‑1 (MQ=255)
                                                                                                                                               aGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCcg                                                                                                    >  1:345074/1‑251 (MQ=255)
                                                                                                                                                gTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCcgc                                                                                                   >  2:185097/1‑251 (MQ=255)
                                                                                                                                                gTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCcgc                                                                                                   >  1:283279/1‑251 (MQ=255)
                                                                                                                                                      gTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGAc                                                                                             >  2:44890/1‑251 (MQ=255)
                                                                                                                                                        cGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACaa                                                                                           <  2:559521/251‑1 (MQ=255)
                                                                                                                                                         gCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACaa                                                                                           >  1:230922/1‑250 (MQ=255)
                                                                                                                                                           tCCCGCACTCGCCTGGTTGATGTTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATc                                                                                                                        >  1:571757/1‑219 (MQ=255)
                                                                                                                                                           tCCCGCACTCGCCGGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCACGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCTCGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTAAGATTAAGGTCATGTTGGTTTTGCTGAAAATc                                                                                                                        <  2:571757/219‑1 (MQ=255)
                                                                                                                                                             ccGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATtgt                                                                                      >  1:401247/1‑251 (MQ=255)
                                                                                                                                                               gCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGa                                                                                    <  1:178687/251‑1 (MQ=255)
                                                                                                                                                                  cTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACg                                                                                 >  1:249740/1‑251 (MQ=255)
                                                                                                                                                                     gCCTGGTTGCTGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCAAGTAGGCAATTTTCGCGTGTTCCGGCCTAATGCGGATCTGTTTACGGTAGAGGTAGGAACCGTATTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGTGGCCAAATTGTGACCGTcc                                                                              >  2:12259/1‑251 (MQ=255)
                                                                                                                                                                          gTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGTTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGgcg                                                                         <  1:662447/251‑1 (MQ=255)
                                                                                                                                                                             gATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTGTGCCGTAAAACTCATCTTCATGCCAGACGTCTGGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGt                                                                      <  1:414219/251‑1 (MQ=255)
                                                                                                                                                                              aTGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCgtcacgcc                                                                     >  2:121764/1‑244 (MQ=255)
                                                                                                                                                                                ggTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCg                                                                    <  1:393142/250‑1 (MQ=255)
                                                                                                                                                                                 gTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGa                                                                                              <  1:267888/223‑1 (MQ=255)
                                                                                                                                                                                 gTGACGGGAAGAGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGATTGATCCTCggggg                                                                                              >  2:267888/1‑222 (MQ=255)
                                                                                                                                                                                   gACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAAAGTCCCGGCGCGTCGact                                                                 >  2:3535/1‑248 (MQ=255)
                                                                                                                                                                                     cGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGAGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGGAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGaataat                                                              >  2:173540/1‑251 (MQ=255)
                                                                                                                                                                                       ggAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATcc                                                            >  2:273391/1‑251 (MQ=255)
                                                                                                                                                                                        gAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCg                                                           <  1:478210/251‑1 (MQ=255)
                                                                                                                                                                                         aaGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCgg                                                          >  1:162415/1‑251 (MQ=255)
                                                                                                                                                                                                   cGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACgg                                                 >  1:78633/1‑250 (MQ=255)
                                                                                                                                                                                                   cGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGa                                                <  1:44890/251‑1 (MQ=255)
                                                                                                                                                                                                       cgcgGTAAATCTCGCTTTTGCCGTCATAATCATCTTCATGCCAGACGTCTTGCGGCAGCTGCCCGTCCGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAACATCAAAATCAAAGGCTTGATCCTCGGGGCCAAATTGTGAACGTCCCGGCGCGTAGAATAATCCGGCAAAAACGGAAGtg                                            >  2:470057/1‑251 (MQ=255)
                                                                                                                                                                                                         cgGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAACCGGAAGtgct                                          >  1:357481/1‑251 (MQ=255)
                                                                                                                                                                                                          ggTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGtgctg                                         <  1:36153/251‑1 (MQ=255)
                                                                                                                                                                                                          ggTAAATCTCGCGTTTGCCGTAAAACTCATCTTCATGCCAGACGGCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGtgctg                                         <  1:146124/251‑1 (MQ=255)
                                                                                                                                                                                                                tctcGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGt                                   >  1:191549/1‑251 (MQ=255)
                                                                                                                                                                                                                     cTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGc                                >  1:505411/1‑249 (MQ=255)
                                                                                                                                                                                                                       tttGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAg                              <  1:587040/249‑1 (MQ=255)
                                                                                                                                                                                                                        ttGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGa                            <  2:235474/250‑1 (MQ=255)
                                                                                                                                                                                                                         tGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGCAGTGCTGCAAAGTAGCAGGATc                          >  2:406599/1‑251 (MQ=255)
                                                                                                                                                                                                                          gCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATc                          <  1:2789/250‑1 (MQ=255)
                                                                                                                                                                                                                           ccGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATc                          >  1:683026/1‑249 (MQ=255)
                                                                                                                                                                                                                           ccGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATc                          >  2:196964/1‑249 (MQ=255)
                                                                                                                                                                                                                           ccGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCGTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGGGCTGCAAAGTAGCAGGATCAg                        <  2:316946/251‑1 (MQ=255)
                                                                                                                                                                                                                           ccGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTATTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAg                        <  2:27796/251‑1 (MQ=255)
                                                                                                                                                                                                                               aaaaCTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTa                    <  2:361444/251‑1 (MQ=255)
                                                                                                                                                                                                                                aaaCTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGAGCAGCGTaa                   <  2:283279/251‑1 (MQ=255)
                                                                                                                                                                                                                                  aCTCATATTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTACCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTaa                   >  1:587444/1‑249 (MQ=255)
                                                                                                                                                                                                                                    tCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTaa                   >  1:616597/1‑247 (MQ=255)
                                                                                                                                                                                                                                     cATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGTATCAGCGTAAAGAt               <  1:602911/250‑1 (MQ=255)
                                                                                                                                                                                                                                     cATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATg              >  1:638001/1‑251 (MQ=255)
                                                                                                                                                                                                                                     cATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGGTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAGCGACCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAgt                <  2:424723/249‑2 (MQ=255)
                                                                                                                                                                                                                                       tCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCg            >  1:317434/1‑251 (MQ=255)
                                                                                                                                                                                                                                        cTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGt           >  2:166250/1‑251 (MQ=255)
                                                                                                                                                                                                                                           cATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTg         <  2:257595/250‑1 (MQ=255)
                                                                                                                                                                                                                                            aTGCCAGACGTCTTGCTGCAGCTGCACGTCGGCAAGTGTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTATCCGTCTTTGATCTGCCAGGTCAGATTAAGGGCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAg       <  2:324843/251‑1 (MQ=255)
                                                                                                                                                                                                                                            aTGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAg       <  2:376306/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 aGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATg  <  1:3535/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
GGCAATTGCGCGGTGGGCTGCTCTTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATG  >  NC_000913/4364552‑4365043

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