Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 494,033 C→G 6.0% intergenic (+10/‑43) dnaX → / → ybaB DNA polymerase III/DNA elongation factor III, tau and gamma subunits/DNA‑binding protein, putative nucleoid‑associated protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913494,0330CG 6.0% 188.3 / 2.3 59intergenic (+10/‑43)dnaX/ybaBDNA polymerase III/DNA elongation factor III, tau and gamma subunits/DNA‑binding protein, putative nucleoid‑associated protein
Reads supporting (aligned to +/- strand):  ref base C (24/30);  new base G (2/2);  total (26/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGACAACGCAGTATGTCTGCATTTGCGCTCCTCTCAGCGGCATTTGAACAACCGCGGTGCACAGCAAAAACTGGCTGAAGCGTTGAGCATGTTAAAAGGTTCAACGGTTGAACTGACTATCGTTGAAGATGATAATCCCGCGGTGCGTACGCCGCTGGAGTGGCGTCAGGCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCGCGCAGCTGGAAGTCACCGGCGAATCTGGCGCAGGTCTGGTAAAAGTGACCATCAACGGTGCACACAACTGCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGACGCAGCACGT  >  NC_000913/493735‑494321
                                                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                                
cGACAACGCAGTATGTCTGCATTTGCGCTCCTCTCAGCGGCATTTGAACAACCGCGGTGCACAGCAAAAACTGGCTGAAGCGTTGAGCATGTTAAAAGGTTCAACGGTTGAACTGACTATCGTTGAAGATGATAATCCCGCGGTGCGTACGCCGCTGGAGTGGCGTCAGGCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATAGTAAGcacaa                                                                                                                                                                                                                                                                                                >  2:402349/1‑300 (MQ=255)
 acactctttccctacccgactctcttccgatctcAGCGGCATTTGAACAACCGCGGTGCACAGCAAAAACTGGCTGAAGCGTTGAGCATGTTAAAAGGTTCAACGGTTGAACTGACTATCGTTGAAGATGATAATCCCGCGGTGCGTACGCCGCGGGAGTGGCGTCAGGCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGc                                                                                                                                                                                                                                                                                               <  1:561778/271‑1 (MQ=255)
          gtttgttgggaTTTGCGTGGCTGTCAGCGGCAGTTGAACAACCGCGGTGCAAAGCAAAAACTGGCTGAAGCGTTGAGCATGTTAAAAGGGTCAACGGTTGAACTGACGATCGTTGAAGATGATAATCCCGCGGTGCGTACGCCGCTGGAGTGGCGTCAGGCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATGATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCg                                                                                                                                                                                                                                                                                      <  1:169766/291‑1 (MQ=255)
          gtattgcTGCATTTGTGCGCCGTTCAGCGGAAGTTGAACAACCGGGTTGCACAGCAAAAACTGGCTGGAGCGTTGAGCATGTTAATAGGTTCAACTGTTGAACTTACTGTCGTTGAAGATGATGATCCCGCGGTGCGTACGCCGCTGGTGTGGCGTCATGTGCTATACGAAGAAAGACTTGCGCAGGCGCGTGAGTCCATTATTGCGGATAATAATATACAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCACCATTTGATCGTAAGCTGTGCTTACGATca                                                                                                                                                                                                                                                                                      <  1:642833/295‑2 (MQ=255)
            aTGTTTGCATTTGTGCTCCTCTGAGCGGAATTTGAACAAGCGCGGTGCAGAGCAAAAACGGTGTGAAGCGTTGTGCAGTGTGCAAGGTTCAACGGTTGAACTTACGAACGTTGAAGATGATAATCCCGCGGTGCGTACGCCGCTGGAGTGGCGTCAGGCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTc                                                                                                                                                                                                                                                                                    <  1:468118/301‑1 (MQ=255)
                           gtgcgttgaGCGGAGTTTGAACAACGGCGGTGAACAGAAAAAACTGGCTGAAGCGTTGAGCATGTTAAAAGGTTCAACGGTTGAACTGACTAGCGTTGAAGATGATAATCCCGCGGTGCGTACGCCGCTGGAGTGGCGTCAGGCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGAGCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGAtt                                                                                                                                                                                                                                                                     <  1:78363/293‑1 (MQ=255)
                                gcgtgtgttcttctgatctccGAGGTGCACAGCAAAAACTGGCTGAATCGTTGAGAATGTTAAAAGGTTAAACGGTTGAACTGACTATCGTTGAAGATGATAATCCCGTGGTGCGTACGCCGCTGGAGTGGCGTCAGGCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTgagag                                                                                                                                                                                                                                                                <  2:328405/282‑1 (MQ=255)
                                     cGGCATTTGAACAACCGCGGTGCACAGCAAAAACTGGCTGAAGCGTTGAGCATGTTAAAAGGTTCAACGGTTGAACTGACTATCGTTGAAGATGATAATCCCGCGGTGCGTACGCCGCTGGAGTGGCGTCAGGCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGCGaaa                                                                                                                                                                                                                                                           >  2:537568/1‑301 (MQ=255)
                                      ggCATTTGAACAACCGCGGTGCACAGCAAAAACTGGCTGAAGCGTTGAGCATGTTAAAAGGTTCAACGGTTGAACTGACTATCGTTGAAGATGATAATCCCGCGGTGCGTACGCCGCTGGAGTGGCGTCAGGCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGaaaa                                                                                                                                                                                                                                                          <  2:574246/301‑2 (MQ=255)
                                              aacaacCGCGGTGCACAGCAAAAACTGGTGGAAGCTTTGTGGATGTTAAAAGGTTCAACGGTTGAACTGACTATCGTTGAAGATGATAATCCCGCGGTGCGTACGCCGCTGGAGTGGCGTCAGGCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATGATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGttt                                                                                                                                                                                                                                                  <  1:536895/301‑1 (MQ=255)
                                               acaacCGCGGTGCACAGCAAAAACTGGCTGAAGCGTTGAGCATGTTAAAAGGTTCAACGGTTGAACTGACTATCGTTGAAGATGATAATCCCGCGGTGCGTACGCCGCTGGAGTGGCGTCAGGCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCCTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAAAGAAACCTATGGTTg                                                                                                                                                                                                                                                 >  1:450012/1‑301 (MQ=255)
                                                  acCGCGGTGCACAGCAAAAACTGGCTGAAGCGTTGAGCATGTTAAAAGGTTCAACGGTTGAACTGACTATCGTTGAAGATGATAATCCCGCGGTGCGTACGCCGCTGGAGTGGCGTCAGGCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTa                                                                                                                                                                                                                                              >  2:600330/1‑301 (MQ=255)
                                                   ccGCGGTGCACAGCAAAAACTGGCTGAAGCGTTGAGCATGTTAAAAGGTTCAACGGTTGAACTGACTATCGTTGAAGATGATAATCCCGCGGTGCGTACGCCGCTGGAGTGGCGTCCGCCGATATACGAAGAAAAACTTTCGCAGGCGCGCGAGTCCATTCTTGCGGATAATAATATTCATACCCTGCGTCGGTTCTTCGATGCGGAGCTGGCTGAAGAAAGTAACCGCCCCATTTGATCATAAGCACATCTTACGGTCGTCAACCCTACAGTTATTgagagagctcggacaaggtgcgtt                                                                                                                                                                                                                                             >  1:328001/1‑284 (MQ=255)
                                                           cacaGCAAAAACTGGCTGAAGCGTTGAGCATGTTAAAAGGTTCAACGGTTGAACTGACTATCGTTGAAGATGATAATCCCGCGGTGCGTACGCCGCTGGAGTGGCGTCAGGCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTc                                                                                                                                                                                                                                     >  2:123000/1‑301 (MQ=255)
                                                             caGCAAAAACTGGCTGAAGCGTTGAGCATGTTAAAAGGTTCAACGGTTGAACTGACTATCGTTGAAGATGATAATCCCGCGGTGCGTACGCCGCTGGAGTGGCGTCAGGCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTg                                                                                                                                                                                                                                   >  2:279816/1‑301 (MQ=255)
                                                             aaGCAAAAACTTGCTGGAGCGTTGAGCATGTTAAAAGGTTCAACGGTTGAACTGACTATCGTTGAAGATGATAATCCCGCGGGGCGTACGCCGCTGGAGTGGCGTCAGGCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTg                                                                                                                                                                                                                                   <  1:599569/300‑1 (MQ=255)
                                                               gttcttccgatcTGAAGCGTTGAGCATGTTAAAAGGTTCAACGGTTGAACTGACTATCGTTGAAGATGATAATCCCGCGGTGCGTACGCCGCTGGAGTGGCGTCAGGCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTggg                                                                                                                                                                                                                                 <  2:145304/290‑1 (MQ=255)
                                                                aaaaaaCTGGCTGAAGCGTTGAGCATGTTAAAAGGTTCAACGGTTGAACTGACTATCGTTGAAGATGATAATCCCGCGGTGCGTACGCCGCTGGAGTGGCGTCAGGCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGt                                                                                                                                                                                                                                <  2:7207/300‑1 (MQ=255)
                                                                 aaaagCTGTGTGAAGGGTTGAGCGTGTTAGAATGTTCAACGTTTGAACGGACTATTGTTGAAGTGGATAATCCCGCGGTGCGTACGCGGCTGGGTTGGCGTCAGCCGATATACGAAGAAAAACTTGCGCAGGCGCGCTAGTCCATTATTGCGGATAGTAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTa                                                                                                                                                                                                                               <  1:570788/301‑1 (MQ=255)
                                                                    aaCTGGCTGAAGCGTTGAGCATGTTAAAAGGTTCAACGGTTGAACTGACTATCGTTGAAGATGATAATCCCGCGGTGCGTACGCCGCTGGAGTGGCGTCAGGCGATATACGAAGAAAAACTAGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAAcc                                                                                                                                                                                                                            <  2:34045/301‑1 (MQ=255)
                                                                            gAAGCGTTGAGCATGTTAAAAGGTTCAACGGTTGAACTGACTATCGTTGAAGATGATAATCCCGCGGTGCGTACGCCGCTGGAGTGGCGTCAGGCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTgggagatcggaagagc                                                                                                                                                                                                                    >  1:145178/1‑288 (MQ=255)
                                                                            gAAGCGTTGAGCATGTTAAAAGGTTCAACGGTTGAACTGACTATCGTTGAAGATGATAATCCCGCGGTGCGTACGCCGCTGGAGTGGCGTCAGGCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGAGGaag                                                                                                                                                                                                                    >  1:388528/1‑301 (MQ=255)
                                                                                        tttttAAAAGGTTCAACGGTTGAACTGACTATCGTTGAAGATGATAATCCCGCGGTGCGTACGCCGCTGGAGTGGCGTCAGGCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCcagcag                                                                                                                                                                                                        <  2:50882/298‑1 (MQ=255)
                                                                                        gtGTTATAAGGTTAGACGGTTGTATTTATTATCGGTGAAGGTGATAATCCCGCGGTGCTTACGCCGCTGGAGTGGCGTTCGGCGATATACGAAGAATAACTTGCGCAGGCGAGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCcagcag                                                                                                                                                                                                        <  1:401835/300‑1 (MQ=255)
                                                                                                                   aCTATCGTTGAAGATGATAATGCCGCGGTGCGTATGCGGCTGGAGTGGCGTCAGGCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAaaaatgcagaaaatgcag                                                                                                                                                                             <  1:441286/301‑1 (MQ=255)
                                                                                                                     tATCGTTGAAGATGATAATCCCGCGGTGCGTACGCCGCTGGAGTGGCGTCAGGCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCCAGAAAAAATGAAGCAAATGCAGGa                                                                                                                                                                           >  1:340049/1‑301 (MQ=255)
                                                                                                                                 tgatAATCCCGCGGTGCGTACGCCGCTGGAGTGGCGTCAGGCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCGCGCa                                                                                                                                                               >  2:85356/1‑301 (MQ=255)
                                                                                                                                         cagcggcgtacgcaccgcGGGAGTGGCGTCAGGCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGCAGCAAGTATCCGCCCCATTTGATCGTCAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGACACCTATGTTTGGTACAGGCGGTCTGGGTACCCTGATGAAGCAAGCACAACCGATGCAAAAAAAAATAAAGAAAATGCAAGAAGAAATAGCACAACAGGaac                                                                                                                                                       >  1:292668/15‑300 (MQ=255)
                                                                                                                                               gcctgaCGCCACTCCAGCGGCGTCAGGCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGAAGAAAATGCAGGAAGAGATCGCGCAGCTGGAAGTCAcca                                                                                                                                                 >  2:537599/6‑300 (MQ=255)
                                                                                                                                                          ttGGAGTTGCGTCAGGCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCGCGCAGCTGGAAGTCACCGGCGAATCTGgc                                                                                                                                      <  2:528450/300‑1 (MQ=255)
                                                                                                                                                               gTGGCGTCAGGCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCGCGCAGCTGGAAGTCACCGGCGAATCTGGCGCAgg                                                                                                                                 <  2:389042/301‑1 (MQ=255)
                                                                                                                                                                 ggCGTCAGGCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCACTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGAAAGAAATCGCGCAACTGGAAGTCACCGGCAAAACTGGCGCAGGTc                                                                                                                               >  1:173078/1‑301 (MQ=255)
                                                                                                                                                                   ggTCAGGCGATATACGAAGAAAAACTTGGGCGTGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTAGCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCGCGCAGCTGGAAGTCACCGGCGAATCTGGCGCAGGTCTg                                                                                                                             <  1:85232/300‑1 (MQ=255)
                                                                                                                                                                    gTCATGCGAGTTACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCGCGCAGCTGGAAGTCACCGGCGAATCTGGCGCAGGTCTgg                                                                                                                            <  1:248483/301‑1 (MQ=255)
                                                                                                                                                                              acaCGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCGCGCAGCTGGAAGTCACCGGCGAATCTGGCGCAGGTCTGGTAAAAGTGAc                                                                                                                  <  2:439622/299‑1 (MQ=255)
                                                                                                                                                                                               cgcAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGCTGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCGCGCAGCTGGAAGTCACCGGCGAATCTGACGCAGGTCTGGTAAAAGTGACCAACAACGGtcaaacaa                                                                                                 >  1:132058/1‑294 (MQ=255)
                                                                                                                                                                                                     cgcgcgAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCGCGCAGCTGGAAGTCACCGGCGAATCTGGCGCAGGTCTGGTAAAAGTGACCATCAACGGTGCACACAACTGcc                                                                                           <  2:340512/301‑1 (MQ=255)
                                                                                                                                                                                                      ggtcggtgcctttattGCGGATAATAATGTTTAGACTCTGCTGCGGTTCTTCGAGGCGGAGCTGGATGAATAAAGGATCCGCCCCATTTGATCGTAAGGACAGCTTACGTTCGTCATCCTTAACGTGATTGGGAGAGAAACCGATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCGCGCAGCTGGAAGTCACCGGCGAATCTGGCGCAGGTCTGGTAAAAGTGACCATCAACGGTGCACACAACTGCCg                                                                                          <  1:609763/290‑1 (MQ=255)
                                                                                                                                                                                                              ccATGATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCGCGCAGCTGGAAGTCACCGGCGAATCTGGCGCAGGTCTGGTAAAAGTGACCATCAACGGTGCACACAACTGCCGTCGCGTag                                                                                  <  2:450576/301‑1 (MQ=255)
                                                                                                                                                                                                                   attGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCGCGCAGCTGGAAGTCACCGGCGAATCTGGCGCAGGTCTGGTAAAAGTTACCATCAACGGTGCACACAACTGCCGTCGCGTAGAGAta                                                                             >  2:528551/1‑300 (MQ=255)
                                                                                                                                                                                                                             ttaataTTCATACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCGCGCAGCTGGAAGTCACCGGCGAATCTGGCGCAGGTCTGGTAAAAGTGACCATCAACGGTGCACACAACTGCCGTCGCGTAGAGATCGACCCGAGcc                                                                   <  2:253074/300‑1 (MQ=255)
                                                                                                                                                                                                                                  atTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGGTTCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCGCGCAGCTGGAAGTCACCGGCGAATCTGGCGCAGGTCTGGTAAAAGTGACCATCAACGGTGCACACAACTGCCGTCGCGTAGAGATCGACCCGAGCctgctg                                                              <  1:279513/301‑1 (MQ=255)
                                                                                                                                                                                                                                     cAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAACAATGCAGAAAATGCAGGAAGAGATCGCGCAGCTGGAAGTCACCGGCGAATCTGGCGCAGGTCTGGTAAAAGTGACCATCAACGGTGCACACAACTGCCGTCGCGTAGAGATCGACCCCAGCCCGCTGGaa                                                           >  1:517996/1‑301 (MQ=255)
                                                                                                                                                                                                                                       gACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCGCGCAGCTGGAAGTCACCGGCGAATCTGGCGCAGGTCTGGTAAAAGTGACCATCAACGGTGCACACAACTGCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAga                                                         >  1:117583/1‑301 (MQ=255)
                                                                                                                                                                                                                                         cccTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCGCGCAGCTGGAAGTCACCGGCGAATCTGGCGCAGGTCTGGTAAAAGTGACCATCAACGGTGCACACAACTGCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAgacg                                                       >  2:470574/1‑301 (MQ=255)
                                                                                                                                                                                                                                         cccTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCGCGCAGCTGGAAGTCACCGGCGAATCTGGCGCAGGTCTGGTAAAAGTGACCATCAACGGTGCACACAACTGCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAgacg                                                       <  1:122876/301‑1 (MQ=255)
                                                                                                                                                                                                                                           cTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTACCCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCGCGCAGCTGGAAGTCACCGGCGAATCTGGCGCAGGTCTGGTAAAAGTGACCATCAACGGTGCACACAACTGCCGTCGCGTAGCGACCGACCCGAGCCTGCTGGAAgacgac                                                     >  1:46319/1‑301 (MQ=255)
                                                                                                                                                                                                                                               gTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCGCGCAGCTGGAAGTCACCGGCGAATCTGGCGCAGGTCTGGTAAAAGTGACCATCAACGGTGCACACAACTGCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACAAag                                                 <  2:518668/301‑1 (MQ=255)
                                                                                                                                                                                                                                                   gTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCGCGCAGCTGGAAGTCACCGGCGAATCTGGCGCAGGTCTGGTAAAAGTGACCATCAACGGTGCACACAACTGCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACAAAaagaa                                             >  2:285957/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                   ggggaggaAGAAAGTAGCCGCCGCATGTGATCGTAAGCAGAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGGCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCGCGCAGCTGGAAGTCACCGGCGAATCTGGCGCAGGTCTGGTAAAAGTGACCATCAACGGTCACACAACTGCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCtggtgg                            <  1:635086/299‑1 (MQ=255)
                                                                                                                                                                                                                                                                     ggATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCGCGCAGCTGGAAGTCACCGGCGAATCTGGCGCAGGTCTGGTAAAAGTGACCATCAACGGTGCACACAACTGCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAgagatcggaagagcacacgactgaactccaga                           >  2:131753/1‑271 (MQ=255)
                                                                                                                                                                                                                                                                        tGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCGCGCAGCTGGAAGTCACCGGCGAATCTGGCGCAGGTCTGGTAAAAGTGACCATCAACGGTGCACACAACTGCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTgc                        >  1:440337/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                         gaagaaAGTTTCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCGCGCAGCTGGAAGTCACCGGCGAATCTGGCGCAGGTCTGGTAAAAGTGACCATCAACGGTGCACACAACTGCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTgca                       <  2:173275/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                          aagcaAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCTTCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAACATGCAGGAAGAGATCGCGCAGCTGGAAGTCACCGGCGAATCTGGCGCAGGTCTGGTAAAAGTGACCATCAACGGTGCACACAACTGCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACAAAAAAAGGATGGAAGACATGGTGGCTgcag                      >  1:459250/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                               tttATCCGCCCCATTTGATCGTAAGCACAGCTTACTTTCGTTATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCGCGCAGCTGGAAGTCACCGGCGAATCTGGCGCAGGTCTGGTAAAAGTGACCATCAACGGTGCACACAACTGCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTc                 <  2:400945/299‑1 (MQ=255)
                                                                                                                                                                                                                                                                                 tATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCGCGCAGCTGGAAGTCACCGGCGAATCTGGCGCAGGTCTGGTAAAAGTGACAATCAACGGTGCACACAACTGCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACAAAACGATGCTGGAAGACCTGGTGGCAGCAGCATTCaa               >  1:400431/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                     cGCCCCAGGTGGTCGTTAGCACAGTGTACGGTTGTAATCCTTAACGTGATTGAGAGAGAAACCTATGGTTGGTAAAGGCGGTCTGGGTAATCTGAGGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAATCGATCGCGCAGCTGGAAGTCACCGGCGAATCTGGCGCAGGTCTGGTAAAAGTGACCATCAACGGTGCACACAACTGCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAacgac           <  1:536926/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                       ccccATTTGGTCGTAAGCATGGCTTATGTTCGTCATCCTTAACGTGATTTGGAGGGAAACCTATGTTTGGTAAAGGCTGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCGCGCAGCTGGAAGTCACCGGCGAATCTGGCGCAGGTCTGGTAAAAGTGACCATCAACGGTGCACACAACTGCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGACgc         <  1:285654/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                          atgacgaacGTAAGCTGTGCTTACGATCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCGCGCAGCTGGAAGTCACCGGCGAATCTGGCGCAGGTCTGGTAAAAGTGACCATCAACGGTGCACACAACTGCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGAGGCTGCAACATTCAacgaaccaca      >  1:445980/9‑295 (MQ=255)
                                                                                                                                                                                                                                                                                              cgaacGTAAGCTGTGCTTACGATCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCGCGCAGCTGGAAGTCACCGGCGAATCTGGCGCAGGTCTGGTAAAAGTGACCATCAACGGTGCACACAACTGCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGACGCAGCAcgg  >  2:27993/5‑300 (MQ=255)
                                                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                                
CGACAACGCAGTATGTCTGCATTTGCGCTCCTCTCAGCGGCATTTGAACAACCGCGGTGCACAGCAAAAACTGGCTGAAGCGTTGAGCATGTTAAAAGGTTCAACGGTTGAACTGACTATCGTTGAAGATGATAATCCCGCGGTGCGTACGCCGCTGGAGTGGCGTCAGGCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCGCGCAGCTGGAAGTCACCGGCGAATCTGGCGCAGGTCTGGTAAAAGTGACCATCAACGGTGCACACAACTGCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGACGCAGCACGT  >  NC_000913/493735‑494321

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