Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,445,219 C→T Q370* (CAG→TAG)  tamB → translocation and assembly module subunit TamB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,445,2190CT98.9% 307.2 / NA 95Q370* (CAG→TAG) tamBtranslocation and assembly module subunit TamB
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (31/62);  total (33/62)
Fisher's exact test for biased strand distribution p-value = 3.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

GTCAGACAACTGGCCGGTGGACATCACTCTCAACAGTACACTGAACGTGGAGCCGTTGAAAGGTGAAAAAGTGAAGCTGAAAATGGGCGGCGCGCTGCGCGAACAGCTGGAGATTGGCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGCCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTT  >  NC_000913/4444978‑4445461
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                  
gTCAGCCAACTGGCCGGTGGACATCACTCTCAACAGTACACTGAACGTGGAGCCGTTGAAAGGTGAAAAAGTGAAGCTGAAAATGGGCGGCGCGCTGCGCGAACAGCTGGAGATTGGCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGTAGTATCAg                                                                                                                                                                                                                                            <  2:347456/250‑1 (MQ=255)
 tCAGACAACTGGCCGGTGGACATCACTCTCAACAGTACACTGAACGTGGAGCCGTTGAAAGGTGAAAAAGTGAAGCTGAAAATGGGCGGCGCGCTGCGCGAACAGCTGGAGATTGGCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGTAGTATCAGGc                                                                                                                                                                                                                                          >  1:310903/1‑251 (MQ=255)
  cAGACAACTGGCCGGTGGACATCACTCTCAACAGTACACTGAACGTGGAGCCGTTGAAAGGTGAAAAAGTGAAGCTGAAAATGGGCGGCGCGCTGCGCGAACAGCTGGAGATTGGCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGTAGTATCAGGc                                                                                                                                                                                                                                          <  1:514704/250‑1 (MQ=255)
  aaGACAACTGGCCGGTGGACATCACTCTCAACAGTACACTGAACGTGGAGCCGTTGAAAGGTGAAAAAGTGAAGCCGAAAATGGGCGGCGCGCTGCGCGAACAGCTGGAGATTGGCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGTAGTATCAGGCg                                                                                                                                                                                                                                         <  1:627124/250‑1 (MQ=255)
    gACAACTGGCCGGTGGACATCACTCTCAACAGTACACTGAACGTGGAGCCGTTGAAAGGTGAAAAAGTGAAGCTGAAAATGGGCGGCGCGCTGCGCGAACAGCTGGAGATTGGCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTCCTGGCCGTTCCATGGTGAGAAGTAGTATCAggaggg                                                                                                                                                                                                                                       >  2:634936/1‑247 (MQ=255)
           ggCCGGTGGACATCACTCTCAACAGTACACTGAACGTGGAGCCGTTGAAAGGTGAAAAAGTGAAGCTGAAAATGGGCGGCGCGCTGCGCGAACAGCTGGAGATTGGCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCCACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTTACTGGTGAGAAGTAGTATCCGGCGGATGATctg                                                                                                                                                                                                                                >  2:263823/1‑251 (MQ=255)
             ccGGTGGACATCACTCTCAACAGTACACTGAACGTGGAGCCGTTGAAAGGTGAAAAAGTGAAGCTGAAAATGGGCGGCGCGCTGCGCGAACAGCTGGAGATTGGCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGTAGTATCAGGCGGATGATctgaa                                                                                                                                                                                                                              <  2:70035/251‑1 (MQ=255)
               ggttgccATCACGCACAACAGTACACTGAACGTGGAGCCGTTGAAAGGTGAAAAAGTGAAGCGGAAAATGGGCAGCGCGCTGCGCGAACAGCTGGAGATTGGCGTTAATCTTTCCGGTCCGTTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAATCCGGATTGCCGCTCAACGTGGAAGTGAACCGCAAACAGCTTTACTGGCCGTTTTCTGGTGAGAAGTAGTATCAGGCTAATCATctgaaac                                                                                                                                                                                                                            <  2:287612/245‑1 (MQ=255)
               ggTGGACATCACTCTCAACAGTACACTGAACGTGGAGCCGTTGAAAGGTGAAAAAGTGAAGCTGAAAATGGGCGGCGCGCTGCGCGAACAGCTGGAGATTGGCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGTAGTATCAGGCGGATGATctgaaac                                                                                                                                                                                                                            <  1:621018/251‑1 (MQ=255)
               ggTGGACATCACTCTCAACAGTACACTGAACGTGGAGCCGTTGAAAGGTGAAAAAGTGAAGCTGAAAACGGGCGGCGCGCTGCGCGAACAGCTGGAGATTGGCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGTAGTATCAGGCGGATGATctgaaac                                                                                                                                                                                                                            <  1:542815/251‑1 (MQ=255)
                gTGGACATCGCTCTCAACAGTACACTGACCGTGGCGCCGTCGAAAGGTGAAAGAGTGACGCTGAAAACGGGCGGCACGCTGCGCGAACAGCTGGAGACTGGCGTGAATCGTTCCGGTCCGGCGGATATGGACTGACGCGCCCATACGCGACTGGCGGCAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACATCTTTTCTTTCCTTTCACTGGTGAGAAGTAGTATCAGGCTGATGTTctgaaact                                                                                                                                                                                                                           <  2:615173/251‑1 (MQ=255)
                gTGGACATCACTCTCAACAGTACACTGAACGTGGAGCCGTTGAAAGGTGAAAAAGTGAAGCTGAAAATGGGCGGCGCGCTGCGCGAACAGCTGGAGATTGGCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGTAGTATCAGGCGGATGATctgaaact                                                                                                                                                                                                                           <  1:357333/251‑1 (MQ=255)
                        cACTCTCAACAGTACACTGAACGTGGAGCCGTTGAAAGGTGAAAAAGTGAAGCTGAAAATGGGCGGCGCGCTGCGCGAACAGCTGGAGATTGGCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGGCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGTAGTATCAGGCGGATGATCTGAAACTGAAACTTa                                                                                                                                                                                                                   <  2:485383/251‑1 (MQ=255)
                        cACTCTCAACAGTACACTGAACGTGGAGCCGTTGAAAGGTGAAAAAGTGAAGCTGAAAATGGGCGGCGCGCTGCGCGAACAGCTGGAGATTGGCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGTAGTATCAGGCGGATGATCTGAAACTGAAACtt                                                                                                                                                                                                                    <  1:577382/250‑1 (MQ=255)
                        cACTCTCAACAGTACACTGAACGTGGAGCCGTTGAAAGGTGAAAAAGTGAAGCTGAAAATGGGCGGCGCGCTGCGCGAACAGCTGGAGATTGGCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGTAGTATCAGGCGGATGATCTGAAACTGAAACTTa                                                                                                                                                                                                                   <  1:349413/251‑1 (MQ=255)
                                   gtaaacTGAACGTGGAGCCGTTGAAAGGTGAAAAAGTGAAGCTGAAAATGGGCGGCGCGCTGCGCGAACAGCTGGAGATTGGCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAGCAGCTTTACTGTCCGTTCACTGGTGAGAAGTAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGTCAAAATg                                                                                                                                                                                                        <  2:574616/247‑1 (MQ=255)
                                     acacGGAACGTGGAGCCGTTGAAAGGTGAAAAAGTGAAGCTGAAAATGGGCGGCGCGCTGCGCGAACAGCTGGAGAGTGGCGTTAATCTTTCCGGTCCGGTGGATATGGATTTTCGCGCCTAGACGCGACTGGCGGAAGGCTGATGGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGTAGTTTCAGACGGATGATCTTAAACTGAATCTTACCGGCAAAATGAc                                                                                                                                                                                                      <  2:427716/251‑1 (MQ=255)
                                      cacTGAACGTGGAGCCGTTGAAAGGTGAAAAAGTGAAGCTGAAAATGGGCGGCGCGCTGCGCGAACAGCTGGAGATTGGCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGTAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGAcc                                                                                                                                                                                                     <  1:200457/251‑1 (MQ=255)
                                      cacTGAACGTGGAGCCGTTGAAAGGTGAAAAAGTGAAGCTGAAAATGGGCGGCGCGCTGCGCGAACAGCTGGAGATTGGCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGTAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGAcc                                                                                                                                                                                                     <  1:263823/251‑1 (MQ=255)
                                            aCGTGGAGCCGTTGAAAGGTGAAAAAGTGAAGCTGAAAATGGGCGGCGCGCTGCGCGAACAGCTGGAGATTGGCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGTAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTac                                                                                                                                                                                               <  1:427256/251‑1 (MQ=255)
                                              gTGGAGCCGTTGAAAGGTGAAAAAGTGAAGCTGAAAATGGGCGGCGCGCTGCGCGAACAGCTGGAGATTGGCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGTAGTATCAGGCGGATGATc                                                                                                                                                                                                                                  <  2:136777/214‑1 (MQ=255)
                                              gTGGAGCCGTTGAAAGGTGAAAAAGTGAAGCTGAAAATGGGCGGCGCGCTGCGCGAACAGCTGGAGATTGGCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGTAGTATCAGGCGGATGATc                                                                                                                                                                                                                                  >  1:136777/1‑214 (MQ=255)
                                                   gCCGTTGAAAGGTGAAAAAGTGAAGCTGAAAATGGGCGGCGCGCTGCGCGAACAGCTGGAGATTGGCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGTAGTATCAGGCGGATGATCTGAAACTGAAACTTACCgg                                                                                                                                                                                                               <  2:505457/228‑1 (MQ=255)
                                                   gCCGTTGAAAGGTGAAAAAGTGAAGCTGAAAATGGGCGGCGCGCTGCGCGAACAGCTGGAGATTGGCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGTAGTATCAGGCGGATGATCTGAAACTGAAACTTACCgg                                                                                                                                                                                                               >  1:505457/1‑228 (MQ=255)
                                                      gTTGAAAGGTGAAAAAGTGAAGCTGAAAATGGGCGGCGCGCTGCGCGAACAGCTGGAGATTGGCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGTAGTATCAGGCGGATGATCTGAAACTGACACTTACCGGCAAAATGACCGATTAAACGctctct                                                                                                                                                                                      >  1:420264/1‑250 (MQ=255)
                                                      attGAAAGGTGAAAAAGTGAAGCTGAAAATGGGAGGCGCGCTGCGCGAACGGCTGGAGATTGGCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGGCTGGCGGAAGCCGGATTGCCGCTCAACGTGTAAGTGAACAGCAAACAGCTTGACGGGCCGTTCACTGGTGAGAATTAGTTTCAAGCAGATGATCTCAAACTCAAACTTACCTGCAAAATGACCGATTACACGctct                                                                                                                                                                                        <  2:294176/247‑1 (MQ=255)
                                                         gAAAGGTGAAAAAGTGAAGCTGAAAATGGGCGGCGCGCTGCGCGAACAGCTGGAGATTGGCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGTAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGc                                                                                                                                                                                  <  2:36986/251‑1 (MQ=255)
                                                          aaaGGTGAAAACGAGAAGCTGAAAATGGGCGGCGCGCTGCGCGAACAGCTGGAGATTGGAGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAAAAGCTTTACTGGCCGTTCACTGGTGAGAAGTAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGAAAAATGACCGATTACACGCTCTCTATGCg                                                                                                                                                                                 >  2:658030/1‑251 (MQ=255)
                                                          aaaGGTGAAAAAGTGAAGCTGAAAATGGGCGGCGCGCTGCGCGAACAGCTGGAGATTGGCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCACAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGTAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCg                                                                                                                                                                                 <  1:191509/251‑1 (MQ=255)
                                                           aaGGTGAAAAAGTGAAGCTGAAAATGGGCGGCGCGCTGCGCGAACAGCTGGAGATTGGCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGTAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGt                                                                                                                                                                                <  1:535528/251‑1 (MQ=255)
                                                           aaGGTGAAAAAGTGAAGCTGAAAATGGGCGGCGCGCTGCGCGAACAGCTGGAGATTGGCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGTAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGt                                                                                                                                                                                >  1:161514/1‑251 (MQ=255)
                                                           aaGGTGAAAAAGTGAAGCTGAAAATGGGCGGCGCGCTGCGCGAACAGCTGGAGATTGGCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGTAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGt                                                                                                                                                                                <  1:316677/251‑1 (MQ=255)
                                                             ggTGAAAAAGTGAAGCTTAAAATGGGCGGCGCGCTGCGCGAACAGCTGGAGATTGGCGTTAATCTTTCCGGTCCGGTTGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGTAGTATCAGGCGGATTATCTGAAACTTAAACTTACCGGCAAAATTACCGATTACACGCTCTCTATGCGTAc                                                                                                                                                                              <  1:278147/251‑1 (MQ=255)
                                                             ggTGAAAAAGTGAAGCTGAAAATGGGCGGCGCGCTGCGCGAACAGCTGGAGATTGGCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGTAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTAc                                                                                                                                                                              <  1:157623/251‑1 (MQ=255)
                                                                                aaaTGGGTGGCGCGCTGCGCGAACAGCTGGAGATTGGCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCAGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGATTTACTGGCCGTTCACTGGTGAGAAGTAGTATCAGGCGGATGATCTGAAACTGAAACTTACCCGCAAAATGATCGATTACACGCGCGCTATGCGTACGGCAGTGAAGGGACAgcag                                                                                                                                                           <  1:562625/251‑4 (MQ=255)
                                                                                  aTGGGCGGCGCGCTGCGCGAGCAGCTGGGGTTTTTGGTTAATCGTTCCGGTCCGGTGGTGATGGATTTACGCTCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGTAGTATCAGGCTGATGATCTGAATCTGAAACGTACCGGCAAAATGACCGATTTCACTCTCTCTATGCGTACGGCAGTGAAGGGACAGgagtt                                                                                                                                                         <  2:361785/251‑3 (MQ=255)
                                                                                   ggggCGGCGCGGTGCGCGAACAGCTGGAGATTGGCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGTAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATc                                                                                                                                                        <  2:371725/250‑1 (MQ=255)
                                                                                    gggCGGCGCGCTGCGCGAACAGCTGGAGATTGGCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGTAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATcc                                                                                                                                                       <  2:25159/251‑1 (MQ=255)
                                                                                     ggcggcGCGCTGCGCGAACAGCTGGAGATTGGCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGGGAACAGCAAACAGCTTTACTGGCCGGTCACTGGTGAGAAGGAGTATCAGGCGGATGATCTGCAACTGAAACTTACAGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGGGAAGGGAAAGGAGATccc                                                                                                                                                      >  1:28020/1‑251 (MQ=255)
                                                                                        ggcgCGCTGCGCGAACAGCTGGAGATTGGCGTTAATCCTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTGCACTGGTGAGAAGTAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATcccgcc                                                                                                                                                   <  2:650155/251‑1 (MQ=255)
                                                                                         gcgcgcTGCGCGAACAGCTGGAGATTGGCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGTAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATcccgccc                                                                                                                                                  <  1:443893/251‑1 (MQ=255)
                                                                                                  gcgNACAGCTGGAGATTGGCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGTAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGTAGTATCATGCGGATGATCTGAAACTGAAACTTACCGGCAAAAGGACCGATTACACTCTCTCTATGTGTACTGCAGTGAAGGGACAGGAGATCCCGCCTGCGACTAtt                                                                                                                                         <  2:608859/251‑1 (MQ=255)
                                                                                                     aaCAGCTGGAGATTGGCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGTAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGCCCGCGACTATTAcc                                                                                                                                      <  2:570915/251‑1 (MQ=255)
                                                                                                       cAGCTGGAGACTGGCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGCGAAGTAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGCCCGCGACTATTACCCt                                                                                                                                    >  2:80220/1‑251 (MQ=255)
                                                                                                           tGGAGATTGGCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGTAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATcccg                                                                                                                                                     <  2:544649/230‑1 (MQ=255)
                                                                                                           tGGAGATTGGCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGTAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATcccg                                                                                                                                                     >  1:544649/1‑230 (MQ=255)
                                                                                                            tgagaTTGGCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGGGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGTAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACGCGCTCTCGATGCGTACGGCAGTGAAGGGCCAGGAGATCCTGCCCGCGTCTATTTCCCTCGAc                                                                                                                                <  2:442230/249‑1 (MQ=255)
                                                                                                            ggataTTGGCGTTATGCTTTCCGGTCCGGTTGTTATGGATTTACGGGCCCATTCGGGTCTGGCGGACGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGAGAGAAGTATTACCAGGCGGAGGATCTGGAACTGCAGCTTACTGGCAAAATGGCCGATTACGCGCTCTTAATGCGTATGGCAGTGAAGGGACGGTAGATCCTGCGCGCGACTCTTACCCTc                                                                                                                                   <  2:391309/243‑1 (MQ=255)
                                                                                                            ggAGATTGGCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTGCTGGCCGTTCACGGGTGAGAAGTAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAAGGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGCCCGCGACTATTACCCTCGACg                                                                                                                               <  2:129550/251‑1 (MQ=255)
                                                                                                            ggAGATTGGCGTTAATCTTTCCGGGCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGTAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGCCCGCGACTATTACCCTCGatg                                                                                                                               <  2:26354/251‑3 (MQ=255)
                                                                                                                    gCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGTAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGCCCGCGACTATTACCCTCGACGCCAAAGGt                                                                                                                       <  2:148557/251‑1 (MQ=255)
                                                                                                                        tAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGTAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCATTGAAGGGACAGGAGATCCCGCCCGCGACTATTACCCTCGACGCCAAAGGTAAt                                                                                                                    <  1:601976/250‑1 (MQ=255)
                                                                                                                              ttCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGTAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGCCCGCGACTATTACCCTCGACGCCAAAGGTAATGAAcagc                                                                                                             >  1:530259/1‑251 (MQ=255)
                                                                                                                               tccggtccggTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGTAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGCCCGCGACTATTACCCTCGACGCCAAAGGTAATGAAcagc                                                                                                             <  1:480191/250‑1 (MQ=255)
                                                                                                                                  ggtccggtGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGTAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTCAAGGGACAGGAAATCCCCCCCGCGACTATTACCCTCGACGCCAAAGGAAATGAACAGTAGGt                                                                                                         >  2:342096/1‑251 (MQ=255)
                                                                                                                                   gtccggtGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGTAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGCCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTc                                                                                                        >  1:37195/1‑251 (MQ=255)
                                                                                                                                      cggtGGATATGGATTTACGCGCCAAGACGCGACTGTCGGAAGCCGGATTGCCGCTCAAGGTGTAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGTAGGATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACAGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGCACGCGACTAGTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAAt                                                                                                     <  2:316984/251‑1 (MQ=255)
                                                                                                                                        gtGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGTAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGCCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAAtct                                                                                                   <  1:570956/251‑1 (MQ=255)
                                                                                                                                          ggATATGGATTTACGCGCCCAGACGCGACTGCCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGTAGTATCAGGCGGATAATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGCCCGCGTCTATTACCCTAGACGACAAAGGTAATGCACAGCAGGTCAATCTCg                                                                                                 >  1:565278/1‑251 (MQ=255)
                                                                                                                                          ggATATGGATTTACGCGACCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGTAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGCCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAAtctc                                                                                                  <  2:459730/250‑1 (MQ=255)
                                                                                                                                             tatGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGTAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGCCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGAc                                                                                               >  2:639604/1‑250 (MQ=255)
                                                                                                                                                    ttACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGTAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGCCCGCGACTATTACCCTCGACGCCCAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCAc                                                                                       >  1:200612/1‑251 (MQ=255)
                                                                                                                                                    ttACGAGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCACACAGCTTTACTGGCCGTTCACTGGTGAGAAGTAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGCCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCAc                                                                                       >  2:270086/1‑251 (MQ=255)
                                                                                                                                                      aCGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGTAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACTCTCGCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGCCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCAcc                                                                                      <  1:72582/250‑1 (MQ=255)
                                                                                                                                                       cgcgCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGTAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGCCCGCGACTATTACCCTCGACGCCAAAGGGAATGAACAGCAGGTCAATCTCGACAAACTCACCGt                                                                                    >  2:486173/1‑251 (MQ=255)
                                                                                                                                                        gcgcCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTGAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGTAGTATCAGGCGGATGATCTGAAACTGAAACGTACCGGCAAAAAGACCGATTACACGCTCTCTATGCGTACAGCAGTGAAGGGACAGGAGATCCCGCCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTTAATCTCGACAAACTCACCGTc                                                                                   <  2:106144/251‑1 (MQ=255)
                                                                                                                                                             cAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGGGAGAAGTAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGCCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCgcggc                                                                              >  1:124571/1‑251 (MQ=255)
                                                                                                                                                             cAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACGGGCCGTTCACTGGTGAGAAGTAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAAGGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGCCCGCGACTATTACCCTCGACGCCAAAGGGAATGAACAGCAGGTCAATCTCGACAAACTCACCGGCgcggc                                                                              <  2:420264/251‑1 (MQ=255)
                                                                                                                                                               gACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGGGAGAAGTAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGCCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGgcgc                                                                            >  2:359346/1‑251 (MQ=255)
                                                                                                                                                                   cgACTGGCGGAAGCCGGAATGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGTAGTAACAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGCCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCCCCCTCgccgcgattg                                                                         >  2:77713/1‑246 (MQ=255)
                                                                                                                                                                   cgACTGGCGGAAACCGGATTGCCGCTTAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGTAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGCCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCGCGACAAACTCACCGTCGCGGCGCTGGa                                                                        <  1:36029/251‑1 (MQ=255)
                                                                                                                                                                    gACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGGGAGAAGTAGTATCAGGCGGATGATCTGAAACTGAAACTTACAGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGCCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGa                                                                        >  2:571712/1‑250 (MQ=255)
                                                                                                                                                                              aGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGTAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGCCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTg                                                             <  1:515215/251‑1 (MQ=255)
                                                                                                                                                                                ccGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTGCTGGCCGTTCACTGGTGAGAAGTAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGCCCGCGACTATTACCCTCGACGCCAAAGGAAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGa                                                            >  1:95856/1‑250 (MQ=255)
                                                                                                                                                                                   gATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGTAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCCGTGAAGGGACAGGAGATCCCGCCCGCGACTATTACCCCCGAGGCCAAAGGGAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAACACTGAACTc                                                        >  1:627258/1‑251 (MQ=255)
                                                                                                                                                                                            tCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGTAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCAGCCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGACAACTGAACTCAAGGCGTTg                                               >  1:562034/1‑251 (MQ=255)
                                                                                                                                                                                            tCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGAGAGAAGCAGTCTCAGGCGGATGATCTGAAACTGAAACTTACCGGAAAAATGACCGATTAAACACTCTCTATGCGTACGGCAGTGCAGGGACAGGAGATCCCACCCGCGACTATTACCCTAGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCACGGCGTTg                                               >  1:80402/1‑251 (MQ=255)
                                                                                                                                                                                             cAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGTAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATTCGTACGGCAGTGAAGGGACAGGAGATCCCGCCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTg                                               <  1:228558/250‑1 (MQ=255)
                                                                                                                                                                                               ccGTGGAAGTGAACAGCAAACAGCGTTACTGGCCGTTCACTGGTGCGAAGTAGTAGCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGTCCGATTACACGCTCTCTAGGCGTACGGCAGTGAAGGGACAGGAGATCCCGCCCGCGACTATTGCCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAACCTCACCGTCGCGGCGCTGGAAGCGAAAACTGAACTCAAGGCGTTGCTc                                            <  2:358464/250‑1 (MQ=255)
                                                                                                                                                                                               aCGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGTAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGCCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTc                                            >  1:657341/1‑251 (MQ=255)
                                                                                                                                                                                                   ggAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGTAGTATCAGGCGGATGGTCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGCCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACt                                        <  2:518854/251‑1 (MQ=255)
                                                                                                                                                                                                    gAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGTAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGCCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTg                                       <  1:159097/251‑1 (MQ=255)
                                                                                                                                                                                                    gAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGTAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGCCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTg                                       <  1:375675/251‑1 (MQ=255)
                                                                                                                                                                                                          aaCAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGTAGAATACGGCGAATGCTCTCAAACTGAAACTTACCAGCAAAATGCCCGATTACACGCTCTCTATGCGTACGGCAGTGCAGGGACAGGAGATACCGCCCGCGACTATTACCCTCGACGCCAAAGGTAATGGACAGCAGGTCAATCTCAACAAACTCACCGGCGCGGCGCTGGCAGGGAAAACTGAACTCAAGGCGTTGCTCGACCGGCAAc                                  >  1:269359/1‑250 (MQ=255)
                                                                                                                                                                                                               cAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGTAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGCCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGcc                             >  1:480353/1‑250 (MQ=255)
                                                                                                                                                                                                                aaaCAGCTTTACTGGCCGTTCACTGGTGAGAAGTAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGCCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCAt                           >  1:540544/1‑251 (MQ=255)
                                                                                                                                                                                                                 aaCAGCTTTACTGGCCGTTCACTGGTGAGAAGTAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGCCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCAtt                          >  1:189848/1‑251 (MQ=255)
                                                                                                                                                                                                                 aaCAGCTTTACTGGCCGTTCACTGGTGAGAAGTAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGCCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCAtt                          >  1:94768/1‑251 (MQ=255)
                                                                                                                                                                                                                    aGCTTTACTGGCCGTTCACTGGTGAGAAGTAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATGACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGAGCCCGCCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGt                       <  1:426595/251‑1 (MQ=255)
                                                                                                                                                                                                                        ttACTGGCCGTTCACTGGTGAGAAGTAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGCCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTg                     <  1:470805/249‑1 (MQ=255)
                                                                                                                                                                                                                              gCCGTTCACTGGTGAGAAGTAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGCCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTg             <  1:133581/251‑1 (MQ=255)
                                                                                                                                                                                                                                cGTTCACTGGTGAGAAGTAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGCCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAg           <  2:310903/251‑1 (MQ=255)
                                                                                                                                                                                                                                     acgggTGAGAAGTAGTATCAGGCGTATGATCTGAAGCTGAAACTCACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGGGAAGGGACAGGAGACCCCGCCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCTATTAGTTGTCGAGGTGatttaac      <  2:55048/248‑7 (MQ=255)
                                                                                                                                                                                                                                     aCTGGTGAGAAGTAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGCCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGTCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAAc      <  1:658030/251‑1 (MQ=255)
                                                                                                                                                                                                                                         gTGAGAAGTAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGCCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCtt  <  1:80220/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                  
GTCAGACAACTGGCCGGTGGACATCACTCTCAACAGTACACTGAACGTGGAGCCGTTGAAAGGTGAAAAAGTGAAGCTGAAAATGGGCGGCGCGCTGCGCGAACAGCTGGAGATTGGCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGCCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTT  >  NC_000913/4444978‑4445461

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