Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,960,264 Δ1 bp 6.2% coding (2165/3369 nt) recC ← exonuclease V (RecBCD complex), gamma chain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,960,2640G. 6.2% 229.9 / 7.9 65coding (2165/3369 nt)recCexonuclease V (RecBCD complex), gamma chain
Reads supporting (aligned to +/- strand):  ref base G (33/28);  new base . (2/2);  total (35/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.35e-01

TGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGCCCGGTAGATAATGACTTTGCCCGATGTAGTCGATCAGTTCCTGCACCAGTACCGACGGGAAACGCTCACTGTTATCCTGAATGGAACGACCGATATAGCTGATATAGAGTTTTTGCTGCGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGGCTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTGAACGGAATTGAACGCATTGGCATCAGAGTACAAATGTTAACCGGTCCGGCGAGAAAACGCTGGCTGATACGTTCTTGATCCAGACGCTGTGCCAGTTCATCACGCAATAGTGACAGCGGCACCGCGTCGCCATACTGC  >  NC_000913/2959981‑2960551
                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                               
tgCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGCCCGGTAGATAATGACTTTGCCCGATGTAGTCGATCAGTTCCTGCACCAGTACCGACGGGAAACGCTCACTGTTATCCTGAATGGAACGACCGATATAGCTGATATAGAGTTTTTGCTGCGCGGAAAtt                                                                                                                                                                                                                                                                                            >  2:125772/1‑289 (MQ=255)
tgCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGCCCGGTAGATAATGACTTTGCCCGATGTAGTCGATCAGTTCCTGCACCAGTACCGACGGGAAACGCTCACTGTTATCCTGAATGGAACGACCGATATAGCTGATATAGAGTTTTTGCTGCGCGGAAAtt                                                                                                                                                                                                                                                                                            >  1:444833/1‑289 (MQ=255)
            ctggcttgCCGCAGGGAGGCATTACCTGGCATAGGTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGCCCGGTAGATAATGACTTTGCCCGATGTAGTCGATCAGTTCCTGCACCAGTACCGACGGGAAACGCTCACTGTTATCCTGAATGGAACGACCGATATAGCTGATATAGAGTTTTTGCTGCGCGGAAATTAACGCTTcc                                                                                                                                                                                                                                                                                   <  2:364439/286‑1 (MQ=255)
            ctggctggCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGCCCGGTAGATAATGACTTTGCCCGATGTAGTCGATCAGTTCCTGCACCAGTACCGACGGGAAACGCTCACTGTTATCCTGAATGGAACGACCGATATAGCTGATATAGAGTTTTTGCTGCGCGGAAATTAACGCTTCCa                                                                                                                                                                                                                                                                                  <  1:574188/287‑1 (MQ=255)
                                        gCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGCCCGGTAGATAATGACATTGCCCGATGTAGTCGATCAGTTCCTGCACCAGTACCGACGGGAAACGCTCACTGTTATCCTGAATGGAACGACCGATATAGCTGATATAGAGTTTTTGCTGCGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACg                                                                                                                                                                                                                                                      >  1:673725/1‑287 (MQ=255)
                                                   tCGTTCGCCTGGCCGGTGGTTTTTTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATCCGCTTTTAACATTGCCCCGCTTTCATCTCAGTTGAGCGCTTCATCTCCCTGTAGATAATGACTTTGCCCGATGTAGTCGATCCGTTCCTGCATCAGTTCCGACGGGATACGCTTTCTTTTATCCTGAATCGAACGACTGCTATCGCTTATACAGCGGTTTTGCTCCGCGGAAATTCACGATTCCAGGATCAGGTACCGGTCGTCCTTGCCATGGCTACTGTCGcc                                                                                                                                                                                                                                          >  2:457741/1‑288 (MQ=255)
                                                        cGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGCCCGGTAGATAATGACTTTGCCCGATGTAGTCGATCAGTTCCTGCACCAGTACCGACGGGAAACGCTCACTGTTATCCTGAATGGAACGACCGATATAGCTGATATAGAGTTTTTGCTGCGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGcc                                                                                                                                                                                                                                     >  2:607352/1‑287 (MQ=255)
                                                        cGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGCCCGGTAGATAATGACTTTGCCCGATGTAGTCGATCAGTTCCTGCACCAGTACCGACGGGAAACGCTCACTGTTATCCTGAATGGAACGACCGATATAGCTGATATAGAGTTTTTGCTGCGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCtt                                                                                                                                                                                                                                    <  2:139986/289‑1 (MQ=255)
                                                             ggctggTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGCCCGGTAGATAATGACTTTGCCCGATGTAGTCGATCAGTTCCTGCACCAGTACCGACGGGAAACGCTCACTGTTATCCTGAATGGAACGACCGATATAGCTGATATAGAGTTTTTGCTGCGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGtt                                                                                                                                                                                                                               <  1:499946/289‑1 (MQ=255)
                                                                        ttGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGCCCGGTAGATAATGACTTTGCCCGATGTAGTCGATCAGTTCCTGCACCAGTACCGACGGGAAACGCTCACTGTTATCCTGAATGGAACGACCGATATAGCTGATATAGAGTTTTTGCTGCGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGGCtcatc                                                                                                                                                                                                                    >  1:380251/1‑289 (MQ=255)
                                                                               tCAAACGGCATCCGGGTATGGAGGCAAGTAAGGTGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGGGCGCTTCATCGCCCGGTAGATAATGACTTTGCCCGATGTAGTCGATCAGTTCCTGCACCAGTACCGACGGGAAACGCTCACTGTTATCCTGAATGGAACGACCGATATAGCTGATATAGAGTTTTTGCTGCGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGTCTCATCAGGTCa                                                                                                                                                                                                              <  1:644644/288‑1 (MQ=255)
                                                                                                aTGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGCCCGGTAGATAATGACTTTGCCCGATGTAGTCGATCAGTTCCTGCACCAGTACCGACGGGAAACGCTCACTGTTATCCTGAATGGAACGACCGATATAGCTGATATAGAGTTTTTGCTGCGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGGCTCATCAGGTCAAAGCCCAATGGCGCAAg                                                                                                                                                                                             <  1:189796/288‑1 (MQ=255)
                                                                                                 tGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGCCCGGTAGATAATGACTTTGCCCGATGTAGTCGATCAGTTCCTGCACCAGTACCGACGGGAAACGCTCACTGTTATCCTGAATGGAACGACCGATATAGCTGATATAGAGTTTTTGCTGCGCGGAAATTAACGCTTc                                                                                                                                                                                                                                                                                    <  2:333182/200‑1 (MQ=255)
                                                                                                 tGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGCCCGGTAGATAATGACTTTGCCCGATGTAGTCGATCAGTTCCTGCACCAGTACCGACGGGAAACGCTCACTGTTATCCTGAATGGAACGACCGATATAGCTGATATAGAGTTTTTGCTGCGCGGAAATTAACGCTTc                                                                                                                                                                                                                                                                                    >  1:333182/1‑200 (MQ=255)
                                                                                                                   gcTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGCCCGGTAGATAATGACATTGCCCGATGTAGTCGATCAGTTCCTGCACCAGTACCGACGGGAAACGCTCACTGTTATCCTGAATGGAACGACCGATATAGCTGATATAGAGTTTTTGCTGCGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGGCTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCcgt                                                                                                                                                                         <  2:673725/289‑1 (MQ=255)
                                                                                                                     ttttACCCTTGCCTCGCTTTCATCCCAGTTGAGCGCTTCATCGCCCGGTAGATAATGACTTTGCCCGATGTAGTCGATCAGTTCCTGCACCAGTACCGACGGGAAACGCTCACTGTTATCCTGAATGGAACGACCGATATAGCTGATATAGAGTTTTTGCTGCGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGGCTCCTCAGGTCAAAGCCCAATGGCGCAAGCTGCCGTGGATAAAAGCcgtcg                                                                                                                                                                       >  2:174742/1‑289 (MQ=255)
                                                                                                                               gCCTCGCTTTCATCACAGTTGAGCGCTTCATCGCCCGGTAGATAATGACTTTGCCCGATGTAGTCGATCAGTTCCTGCACCAGTACCGACGGGAAACGCTCACTGTTATCCTGAATGGAACGACCGATATAGCTGATATAGAGTTTTTGCTGCGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGGCTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACg                                                                                                                                                                                       >  2:503776/1‑263 (MQ=255)
                                                                                                                               gCCTCGCTTTCATCACAGTTGAGCGCTTCATCGCCCGGTAGATAATGACTTTGCCCGATGTAGTCGATCAGTTCCTGCACCAGTACCGACGGGAAACGCTCACTGTTATCCTGAATGGAACGACCGATATAGCTGATATAGAGTTTTTGCTGCGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGGCTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACg                                                                                                                                                                                       <  1:503776/263‑1 (MQ=255)
                                                                                                                                             acaGTTGAGCGCTTCATCGCCCGGTAGATAATGACTTTGCCCGATGTAGTCGATCAGTTCCTGCACCAGTACCGACGGGAAACGCTCACTGTTATCCTGAATGGAACGACCGATATAGCTGATATAGAGTTTTTGCTGCGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGGCTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACt                                                                                                                                               <  2:84729/289‑1 (MQ=255)
                                                                                                                                                  tGAGCGCTTCATCGCCCGGTAGATAATGACTTTGCCCGATGTAGTCGATCAGTTCCTGCACCAGTACCGACGGGAAACGCTCACTGTTATCCTGAATGGAACGACCGATATAGCTGATATAGAGTTTTTGCTGCGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGGCTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTGaa                                                                                                                                          <  1:534362/289‑1 (MQ=255)
                                                                                                                                                   gAGCGCTTCATCGCCCGGTAGATAATGACTTTGCCCGATGTAGTCGATCAGTGCCTGCACCAGTACCGACGGGAAACGCTCACTGTTATCCTGAATGGAACGACCGATATAGCTGATATAGAGTTTTTGCTGCGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGGCTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTGaa                                                                                                                                          <  1:357135/288‑1 (MQ=255)
                                                                                                                                                         ttCATCGCCCGGTAGATAATGCATTTGCCCAGTGTAGTAGAGCAGTTCCTGCCCCAGAACCGACGGGAAACGCTAACTGTTATCCTGAATGGAACGACCGATATAGCTGATATAGAGTTTTTGCTGCGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGGCTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTGAACGGAAtt                                                                                                                                   <  2:446656/289‑1 (MQ=255)
                                                                                                                                                         ttCATCGCCCGGTAGATAATGACTTTTCCCGATGTAGTCGATCAGTTCCTGCACCAGTACCGACGGGACACGCTCACTGTTATCCTGAATGGACCGACCGATATAGCTGATATAGAGTTTTTGCTGCTCGGAACTTAATGCTTACAGGAACAGATAGCGGTCGTCATCGCGACGGCTACTGACGCCACGCAGTGGTTTGTGGCTCAGCAGGTACACGCCAACAGGCGCACTCTGAGGTGGATAACCCACGTCGCTCATTCCCACCCGCTATACCATTTTGATCgatgtt                                                                                                                                   >  2:127023/1‑284 (MQ=255)
                                                                                                                                                                        gtaagtaCTTTGCCCGATGTAGTCAATAAGTTCCTGCACCAGTACCGAGGGGGACAGACGACTCTTAGCCTGAATGGAACGAGCGAGATAGCGGATATAGAGTTTTTGCGGCGCGGAAATTAACGCGTCCAGGAACAGATGGCGGTCGGCATGGCGACGGCTACGGGCGCCACGCTTCGGTTTCTGGCTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTGAACGGAATTGAACGCATTGGCAt                                                                                                                     <  2:444833/282‑1 (MQ=255)
                                                                                                                                                                            tgccTTTGCCCGATGTAGTCGATCAGTTCCTGCACCAGTACCGACGGGAAACGCTCACTGTTATCCTGAATGGAACGACCGATATAGCTGATATAGAGTTTTTGCTGCGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGGCTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTGAACGGAATTGAACGCATTGGCATCagag                                                                                                                <  1:694100/286‑1 (MQ=255)
                                                                                                                                                                               cTTTGCCCGATGTAGTCGATCAGTTCCTGCACCAGTACCGACGGGAAACGCTCACTGTTATCCTGAATGGAACGACCGATATAGCTGATATAGAGTTTTTGCTGCGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGGCTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTGAACGGAATTGAACGCATTGGCATCa                                                                                                                   >  2:123149/1‑283 (MQ=255)
                                                                                                                                                                               cTTTGCCCGATGTAGTCGATCAGTTCCTGCACCAGTACCGACGGGAAACGCTCACTGTTATCCTGAATGGAACGACCGATATAGCTGATATAGAGTTTTTGCTGCGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGGCTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTGAACGGAATTGAACGCATTGGCATCa                                                                                                                   <  1:123149/283‑1 (MQ=255)
                                                                                                                                                                               cTTTGCCCGATGTAGTCGATCAGTTCCTGCACCAGTACCGACGGGAAACGCTCACTGTTATCCTGAATGGAACGACCGATATAGCTGATATAGAGTTTTTGCTGCGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGGCTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTGAACGGAATTGAACGCATTGGCATCAGAGTAc                                                                                                             >  1:54559/1‑289 (MQ=255)
                                                                                                                                                                               cTTTGCCCGATGTAGTCGATCAGTTCCTGCACCAGTACCGACGGGAAACGCTCACTGTTATCCTGAATGGAACGACCGATATAGCTGATATAGAGTTTTTGCTGCGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGGCTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTGAACGGAATTGAACGCATTGGCATCAGAGTAc                                                                                                             >  2:499503/1‑289 (MQ=255)
                                                                                                                                                                               cTTTGCCCGATGTAGTCGATCAGTTCCTGCACCAGTACCGACGGGAAACGCTCACTGTTATCCTGAATGGAACGACCGATATAGCTGATATAGAGTTTTTGCTGCGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGGCTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTGAACGGAATTGAACGCATTGGCATCAGAGTAc                                                                                                             <  1:125772/289‑1 (MQ=255)
                                                                                                                                                                                 ttGCCCGATGTAGTCGATCAGTTCCTGCACCAGTACCGACGGGAAACGCTCACTGTTATCCTGAATGGAACGACCGATATAGCTGATATAGAGTTTTTGCTGCGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGGCTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTGAACGGAATTGAACGCATTGGCATCAGAGTACaa                                                                                                           <  1:378177/289‑1 (MQ=255)
                                                                                                                                                                                  tGCCCGATGTAGTCGATCAGTTCCTGCACCAGTACCGACGGGAAACGCTCACTGTTATCCTGAATGGAACGACCGATATAGCTGATATAGAGTTTTTGCTGCGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGGCTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTccc                                                                                                                                                              >  1:6187/1‑237 (MQ=255)
                                                                                                                                                                                  tGCCCGATGTAGTCGATCAGTTCCTGCACCAGTACCGACGGGAAACGCTCACTGTTATCCTGAATGGAACGACCGATATAGCTGATATAGAGTTTTTGCTGCGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGGCTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTccc                                                                                                                                                              <  2:6187/237‑1 (MQ=255)
                                                                                                                                                                                  tGCCCGATGTAGTCGATCAGTTCCTGCACCAGTACCGACGGGAAACGCTCACTGTTATCCTGAATGGAACGACCGATATAGCTGATATAGAGTTTTTGCTGCGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGGCTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTGAACGGAATTGAACGCATTGGCATCAGAGTACaaa                                                                                                          >  1:111632/1‑289 (MQ=255)
                                                                                                                                                                                  tGCCCGATGTAGTCGATCAGTTCCTGCACCAGTACCGACGGGAAACGCTCACTGTTATCCTGAATGGAACGACCGATATAGCTGATATAGAGTTTTTGCTGCGCG‑‑AATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGGCTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTGAACGGAATTGAACGCATTGGCATCAGAGTACAAATg                                                                                                        >  2:380183/1‑289 (MQ=255)
                                                                                                                                                                                  tGCCCGATGTAGTCGATCAGTTCCTGCACCAGTACCGACGGGAAACGCTCACTGTTATCCTGAATGGAACGACCGATATAGCTGATATAGAGTTTTTGCTGCGCG‑‑AATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGGCTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTGAACGGAATTGAACGCATTGGCATCAGAGTACAAATg                                                                                                        >  2:313476/1‑289 (MQ=255)
                                                                                                                                                                                   gCCCGATGTAGTCGATCAGTTCCTGCACCAGTACCGACGGGAAACGCTCACTGTTATCCTGAATGGAACGACCGATATAGCTGATATAGAGTTTTTGCTGCACGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGTCACGCTTCGGTTTCTGGCTCATCAGGGCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTGAACGGAATTGAACGCATTGGCATCAGAGTACAAAt                                                                                                         <  1:457741/289‑1 (MQ=255)
                                                                                                                                                                                       gATGTAGACGATCAGTTCCTGCCCCAGTACCGACGGGAAACGCTCACTGTTATCCTGAATGGAACGACCGATATAGCTGATATAGAGTTTTTGCTGCGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCCACGGCTACGGTCGCCACGCTTCGGTTCCTGGCTCATCAGTTCATAGCCCAAGGGCGCAATCTGACGTGGATAAACTCCCTCGTCTCTCCCCAGCTGGCAAACCACTTTGAACGGTATTGAACGCATTAGCCTCCGCTTACAAAttcta                                                                                                     >  2:218456/1‑285 (MQ=255)
                                                                                                                                                                                        aTGTAGTCGATCAGTTCCTGCACCAGTACCGACGGGAAACGCTCACTGTTATCCTGAATGGAACGACCGATATAGCTGATATAGAGTTTTTGCTGCGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGGCTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCAAGCAGGCAAACCACTTTGAACGGAATTGAACGCATTGGCATCAGAGTACAAATGTTa                                                                                                     >  1:55706/1‑288 (MQ=255)
                                                                                                                                                                                                gATCAGTTCCTGCACCAGTACCGACGGGAAACGCTCACTGTTATCCTGAATGGAACGACCGATATAGCTGATATAGAGTTTTTGCTGCGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGGCTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTGAACGGAATTGAACGCATTGGCATCAGAGTACAAATGTTAACCGGTCCg                                                                                            >  1:469711/1‑289 (MQ=255)
                                                                                                                                                                                                 aTCAGTTCCTGCACCAGTACCGACGGGAAACGCTCACTGTTATCCTGAATGGAACGACCGATATAGCTGATATAGAGTTTTTGCTGCGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGGCTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAAcc                                                                                                                                                  <  1:279961/234‑1 (MQ=255)
                                                                                                                                                                                                 aTCAGTTCCTGCACCAGTACCGACGGGAAACGCTCACTGTTATCCTGAATGGAACGACCGATATAGCTGATATAGAGTTTTTGCTGCGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGGCTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAAcc                                                                                                                                                  >  2:279961/1‑234 (MQ=255)
                                                                                                                                                                                                      ttCCTGCACCAGTACCGACGGGAAACGCTCACTGTTATCCTGAATGGAACGACCGATATAGCTGATATAGAGTTTTTGCTGCGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGGCTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTGAACGGAATTGAACGCATTGGCATCAGAGTACAAATGTTAACCGGTCCGGCgaga                                                                                      >  1:86552/1‑289 (MQ=255)
                                                                                                                                                                                                        ccTGCACCAGTACCGACGGGAAACGCTCACTGTTATCCTGAATGGAACGACCGATATAGCTGATATAGAGTTTTTGCTGCGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGGCTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTGAACGGAATTGAACGCATTGGCATCAGAGTACAAATGTTAACCGGTCCGGCGAGaaa                                                                                    >  2:583142/1‑289 (MQ=255)
                                                                                                                                                                                                           gCACCAGTACCGACGGGAAACGCTCACTGTTATCCTGAATGGAACGACCGATATAGCTGATATAGAGTTTTTGCTGCGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGGCTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTGAACGGAATTGAACGCATTGGCATCAGAGTACAAATGTTAACCGGTCCGGCGAGAAAACg                                                                                 <  1:174742/289‑1 (MQ=255)
                                                                                                                                                                                                                   aCCGACGGGAAACGCTCACTGTTATCCTGAATGGAACGACCGATATAGCTGATATAGAGTTTTTGCTGCGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGGCTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTGAACGGAATTGAACGCATTGGCATCAGAGTACAAATGTTAACCGGTCCGGCGAGAAAACGCTGGCTGa                                                                         <  1:235270/289‑1 (MQ=255)
                                                                                                                                                                                                                      gACGGGAAACGCTCACTGTTATCCTGAATGGAACGACCGATATAGCTGATATAGAGTTTTTGCTGCGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGGCTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTGAACGGAATTGAACGCATTGGCATCAGAGTACAAATGTTAACCGGTCCGGCGAGAAAACGCTGGCTGATAc                                                                      <  1:495571/289‑1 (MQ=255)
                                                                                                                                                                                                                        cGGGAAACGCTCACTGTTATCCTGAATGGAACGACCGATATAGCTGATATAGAGTTTTTGCTGCGCG‑‑AATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGGCTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTGAACGGAATTGAACGCATTGGCATCAGAGTACAAATGTTAACCGGTCCGGCGAGAAAACGCTGGCTGATACGTTc                                                                  <  1:313476/289‑1 (MQ=255)
                                                                                                                                                                                                                        cGGGAAACGCTCACTGTTATCCTGAATGGAACGACCGATATAGCTGATATAGAGTTTTTGCTGCGCG‑‑AATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGGCTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTGAACGGAATTGAACGCATTGGCATCAGAGTACAAATGTTAACCGGTCCGGCGAGAAAACGCTGGCTGATACGTTc                                                                  <  1:380183/289‑1 (MQ=255)
                                                                                                                                                                                                                                    aCTGTTATCCTGAATGGAACGACCGATATAGCTGATATAGAGTTTTTGCTGCGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGGCTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAACCGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTGAACGGAATTGAACGCATTGGCATCCGAGTACAAATGTTAACCGGTCCGGCCAGAAAACGCGGGCTGATACGTTTCTGAACCAg                                                         >  2:150856/1‑288 (MQ=255)
                                                                                                                                                                                                                                      tttgagccTGAATGGAACGCCCTGAGTAGCGATTATAGAGTTTTTGCTGCGCTGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCGACGGTCGCCACGCTTCGGTTTCTGGCTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTGAACGGAATTGAACGCATTGGCATCAGAGTACAAATGTTAACCGGTCCGGCGAGAAAACGCTGGCTGATACGTTCTTGATCCAGa                                                        <  2:373992/281‑1 (MQ=255)
                                                                                                                                                                                                                                          aTCCTGAATGGAACGACCGATATAGCTGATATAGAGTTTTTGCTGCGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGGCTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCCGGCAAACCACTTTGAACGGAATTGAACGCATTGGCATCAGAGTACAAATGTTAACCGGTCCGGCGAGAAAAACGTGGCTGAAACGTTCTTGGTCCAGACGCtgt                                                  >  2:241830/1‑289 (MQ=255)
                                                                                                                                                                                                                                                aaTGGAACGACCGATATAGCTGATATAGAGTTTTTGCTGCGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGGCTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTGAACGGAATTGAACGCATTGGCATCAGAGTACAAATGTTAACCGGTCCGGCGAGAAAACGCTGGCTGATACGTTCTTGATCCAGACGCTGTGCCa                                              >  1:373992/1‑287 (MQ=255)
                                                                                                                                                                                                                                                 aTGGAACGACCGATATAGCGTATATAGAGTTTTTGCTGGGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGGGACGGCTACGGTCGCCAGGCTTCGGTTTCTGGCTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTGAACGGAGTTGAACGCATTGGCATCAGAGTACAAATGTTAACCGGTCCGGCGAGAAAACGCTGGCTGATACGTTCTTGATCCAGACGCTGTGCCAGtt                                           <  2:55706/289‑1 (MQ=255)
                                                                                                                                                                                                                                                               ataGCTGATATAGAGTTTTTGCTGCGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGGCTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTGAACGGAATTGAACGCATTGGCATCAGAGTACAAATGTTAACCGGTCCGGCGAGAAAACGCTGGCTGATACGTTCTTGATCCAGACGCTGTGCCAGTTCATCACGCAATAGt                             >  1:270611/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                        ataGAGTTTTTGCTGCGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGGCTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTGAACGGAATTGAACGCATTGGCCTCAGAGTACAAATGTTAACCGGTCCGGCTAGAAAACGCTGGCTGATACGTTCTTGATCCAGACGCTGTGCCAGTTCATCACGCAATAGTGACAGCgg                     >  2:659541/1‑288 (MQ=255)
                                                                                                                                                                                                                                                                            agTTTTTGCTGCGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGGCTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTGAACGGAATTGAACGCATTGGCATCAGAGTACAAATGTTAACCGGTCCGGCGAGAAAACGCTGGCTGATACGTTCTTGATCCAGACGCTGTGCCAGTtcatc                                       <  1:696940/266‑1 (MQ=255)
                                                                                                                                                                                                                                                                            agTTTTTGCTGCGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGGCTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTGAACGGAATTGAACGCATTGGCATCAGAGTACAAATGTTAACCGGTCCGGCGAGAAAACGCTGGCTGATACGTTCTTGATCCAGACGCTGTGCCAGTtcatc                                       >  2:696940/1‑266 (MQ=255)
                                                                                                                                                                                                                                                                            agTTTTTGCTGCGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGGCTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTGAACGGAATTGAACGCATTGGCATCAGAGTACAAATGTTAACCGGTCCGGCGAGAAAACGCTGGCTGATACGTTCTTGATCCAGACGCTGTGCCAGTTCATCACGCAATAGTGACAGCGGCAcc                 >  2:331206/1‑288 (MQ=255)
                                                                                                                                                                                                                                                                            agTTTTTGCTGCGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGGCTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTGAACGGAATTGAACGCATTGGCATCAGAGTACAAATGTTAACCGGTCCGGCGAGAAAACGCTGGCTGATACGTTCTTGATCCAGACGCTGTGCCAGTTCATCACGCAATAGTGACAGCGGCAcc                 >  1:153862/1‑288 (MQ=255)
                                                                                                                                                                                                                                                                             gTTTTTGCTGCGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGGCTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTGAACGGAATTGAACGCATTGGCATCAGAGTACAAATGTTAACCGGTCCGGCGAGAAAACGCTGGCTGATACGTTCTTGATCCAGACGCTGTGCCAGTTCATCACGCAATAGTGACAGCGGCACcg                >  2:175087/1‑288 (MQ=255)
                                                                                                                                                                                                                                                                             gTTTTTGCTGCGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGGCTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTGAACGGAATTGAACGCATTGGCATCAGAGTACAAATGTTAACCGGTCCGGCGAGAAAACGCTGGCTGATACGTTCTTGATCCAGACGCTGTGCCAGTTCATCACGCAATAGTGACAGCGGCACcg                <  2:86552/288‑1 (MQ=255)
                                                                                                                                                                                                                                                                                tttGCTGCGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGGCTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTGAACGGAATTGAACGCATTGGCATCAGAGTACAAATGTTAACCGGTCCGGCGAGAAAACGCTGGCTGATACGTTCTTGATCCAGACGCTGTGCCAGTTCATCACGCAATAGTGACAGCGGCACCGCGTc            >  1:520658/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                                tttGCTGCGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGGCTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTGAACGGAATTGAACGCATTGGCATCAGAGTACAAATGTTAACCGGTCCGGCGAGAAAACGCTGGCTGATACGTTCTTGATCCAGACGCTGTGCCAGTTCATCACGCAATAGTGACAGCGGCACCGCGTc            >  1:396142/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                                 ttGCTGCGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGGCTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTGAACGGAATTGAACGCATTGGCATCAGAGTACAAATGTTAACCGGTCCGGCGAGAAAACGCTGGCTGATACGTTCTTGATCCAGACGCTGTGCCAGTTCATCACGCAATAGTGACAGCGGCACCGCGTc            >  2:186870/1‑288 (MQ=255)
                                                                                                                                                                                                                                                                                           gAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGGCTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTGAACGGAATTGAACGCATTGGCATCAGAGTACAAATGTTAACCGGTCCGGCGAGAAAACGCTGGCTGATACGTTCTTGATCCAGACGCTGTGCCAGTTCATCACGCAATAGTGACAGCGGCACCGCGTCGCCATACTgc  <  2:153862/288‑1 (MQ=255)
                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                               
TGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGCCCGGTAGATAATGACTTTGCCCGATGTAGTCGATCAGTTCCTGCACCAGTACCGACGGGAAACGCTCACTGTTATCCTGAATGGAACGACCGATATAGCTGATATAGAGTTTTTGCTGCGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGGCTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTGAACGGAATTGAACGCATTGGCATCAGAGTACAAATGTTAACCGGTCCGGCGAGAAAACGCTGGCTGATACGTTCTTGATCCAGACGCTGTGCCAGTTCATCACGCAATAGTGACAGCGGCACCGCGTCGCCATACTGC  >  NC_000913/2959981‑2960551

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