Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,684,7120AC42.9% 101.5 / 10.6 63T152P (ACC→CCC) tusinhibitor of replication at Ter, DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (28/8);  new base C (0/27);  total (28/35)
Fisher's exact test for biased strand distribution p-value = 5.21e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCGATAACCTTTCGCAAGCAGCGTTGGTCAGTCATATTCAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAAT‑GGGTGCATCGTCATTTGCC‑GGGGCTGATCACCCTTAA‑TGCTTACCGCACGCTCACCGTTCTG‑CAC‑GACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTC  >  NC_000913/1684562‑1684858
                                                                                                                                                        |                                                                                                                                                     
gTCGATAACCTTTCGCAAGCAGCGTTGGTCAGTCATATTCAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAAT‑GGGTGCATCGTCATTTGCC‑GGGGCTGATCa                                                                                                                                                       <  1:117671/151‑1 (MQ=255)
    aTAACCTTTCGCAAGCAGCGTTGGTCAGTCATATTCAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAAT‑GGGTGCATCGTCATTTGCC‑GGGGCTGATCACCCt                                                                                                                                                   <  2:467286/151‑1 (MQ=255)
    aTAACCTTTCGCAAGCAGCGTTGGTCAGTCATATTCAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAAT‑GGGTGCATCGTCATTTGCC‑GGGGCTGATCACCCt                                                                                                                                                   <  1:658484/151‑1 (MQ=255)
     tAACCTTTCGCAAGCAGCGTTGGTCAGTCATATTCAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAAT‑GGGTGCATCGTCATTTGCC‑GGGGCTGATCACCCtt                                                                                                                                                  >  1:1242000/1‑151 (MQ=255)
          tttCGCAAGCAGCGTTGGTCAGTCATATTCAGCACATCAATAAACTCGAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAAT‑GGGTGCATCGTCATTTGCC‑GGGGCTGATCACCCTTAA‑TGc                                                                                                                                            >  2:1002113/1‑151 (MQ=255)
             cGCAAGCAGCGTTGGTCAGTCATATTCAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAAT‑GGGTGCATCGTCATTTGCC‑GGGGCTGATCACCCTTAA‑TGCtt                                                                                                                                          >  2:213125/1‑150 (MQ=255)
             cGCAAGCAGCGTTGGTCAGTCATATTCAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAAT‑GGGTGCATCGTCATTTGCC‑GGGGCTGATCACCCTTAA‑TGCtt                                                                                                                                          <  1:177614/150‑1 (MQ=255)
             cGCAAGCAGCGTTGGTCAGTCATATTCAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAAT‑GGGTGCATCGTCATTTGCC‑GGGGCTGATCACCCTTAA‑TGCTTa                                                                                                                                         <  2:921465/151‑1 (MQ=255)
                  gcagcGTTGGTCAGTCATATTCAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAAT‑GGGTGCATCGTCATTTGCC‑GGGGCTGATCACCCTTAA‑TGCTTACCGc                                                                                                                                     >  1:741905/1‑150 (MQ=255)
                  gcagcGTTGGTCAGTCATATTCAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAAT‑GGGTGCATCGTCATTTGCC‑GGGGCTGATCACCCTTAA‑TGCTTACCGCa                                                                                                                                    <  1:495548/151‑1 (MQ=255)
                    agcGTTGGTCAGTCATATTCAGCACATCAATAAACTCAAGACCGCGGTCGAGCATATCGTCAAGGTTGAATCAGAACTCCCCACCGCGGCAAGTTTTGAAT‑GGGTGCATCGTCATTTGCC‑GGGGCTGATCACCCTTAA‑TGCTTACCGCAc                                                                                                                                   >  1:1007248/1‑150 (MQ=255)
                         tGGTCAGTCATATTCAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAAT‑GGGTGCATCGTCATTTGCC‑GGGGCTGATCACCCTTAA‑TGCTTACCGCACGCTCa                                                                                                                              <  2:130312/150‑1 (MQ=255)
                         tGGTCAGTCATATTCAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAAT‑GGGTGCATCGTCATTTGCC‑GGGGCTGATCACCCTTAA‑TGCTTACCGCACGCTCAc                                                                                                                             >  1:780419/1‑151 (MQ=255)
                             cagtcaTATTCAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAAT‑GGGTGCATCGTCATTTGCC‑GGGGCTGATCACCCTTAA‑TGCTTACCGCACGCTCACCGtt                                                                                                                         >  1:1096231/1‑151 (MQ=255)
                              agtcaTATTCAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAAT‑GGGTGCATCGTCATTTGCC‑GGGGCTGATCACCCTTAA‑TGCTTACCGCACGCTCACCGTTc                                                                                                                        >  2:294502/1‑151 (MQ=255)
                                          cacaTCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAAT‑GGGTGCATCGTCATTTGCC‑GGGGCTGATCACCCTTAA‑TGCTTACCGCACGCTCACCGTTCTG‑CAC‑GACCCCg                                                                                                          >  1:1233152/1‑151 (MQ=255)
                                          cacaTCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAAT‑GGGTGCATCGTCATTTGCC‑GGGGCTGATCACCCTTAA‑TGCTTACCGCACGCTCACCGTTCTG‑CAC‑GACCCCg                                                                                                          <  1:97093/151‑1 (MQ=255)
                                                   aaaCTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAAT‑GGGTGCATCGTCATTTGCC‑GGGGCTGATCACCCTTAA‑TGCTTACCGCACGCTCACCGTTCTG‑CAC‑GACCCCGCCACTTTAc                                                                                                 >  2:126810/1‑151 (MQ=255)
                                                   aaaCTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAAT‑GGGTGCATCGTCATTTGCC‑GGGGCTGATCACCCTTAA‑TGCTTACCGCACGCTCACCGTTCTG‑CAC‑GACCCCGCCACTTTAc                                                                                                 >  1:639255/1‑151 (MQ=255)
                                                       tCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAAT‑GGGTGCATCGTCATTTGCC‑GGGGCTGATCACCCTTAA‑TGCTTACCGCACGCTCACCGTTCTG‑CAC‑GACCCCGCCACTTTACGCtt                                                                                             >  2:742769/1‑151 (MQ=255)
                                                         aaGCCCACTTTCGAGCATATCGTCACGGTTAAATCAAACCCCCCCCCCGCGCCACGTTTTGAAG‑GGGTGCATCGTCTTTTGCC‑GGGGCTGATCCCCCTTAA‑TGCTTCCCGCCCGCTCACCTTTCTG‑CAC‑GCCCCCGCCACTTTACGCTttg                                                                                           <  1:676996/151‑1 (MQ=255)
                                                             ccACGTTGGAGCATATCGTCCGGGTTAAACCAGAACCCCCCCCCGCGGCACTTTTTGAAT‑GGGTGCTTCGTTATTTGCC‑GGGGCTGATCCCCCTTAA‑TGCTTACCGCCCGCCCCCCTTTCTG‑CCC‑GCCCCCGCCACTTTACGCTttggt                                                                                         <  1:1278322/149‑1 (MQ=255)
                                                                cGTTCGAGCATATCGTCACGGTTAAATCAGAACCCCCCCCCGCGGCACGTTTTGAAT‑GGGTGCATCGTCTTTTGCC‑GGGGCTGATCCCCCTTAA‑TGCTTCCCGCCCCCTCCCCGTTCTG‑CAC‑GACCCCGCCACTTTACGCTTTGGTTGGGc                                                                                    <  1:1129602/151‑1 (MQ=255)
                                                                  ttCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAAT‑GGGTGCATCGTCATTTGCC‑GGGGCTGATCACCCTTAA‑TGCTTACCGCACGCTCACCGTTCTG‑CAC‑GACCCCGCCACTTTACGCTTTGGTTGGGCt                                                                                   >  2:906580/1‑150 (MQ=255)
                                                                     gAGCATTTCGTCACGTTTGAATCAGACCCCCCCCCCGCGGCACTTTTTGATT‑GGGTGCCTCTTCTTTGGCC‑GGGGCTGACCCCCTTTAT‑TTCTTCCCGCCCTCTCACCTTCCTG‑CCC‑CCCCCCGCCACTTTACGCTTTGGTTGGGCtaata                                                                               <  2:743577/151‑1 (MQ=255)
                                                                      aGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAAT‑GGGTGCATCGTCATTTGCC‑GGGGCTGATCACCCTTAA‑TGCTTACCGCACGCTCACCGTTCTG‑CAC‑GACCCCGCCACTTTACGCTTTGGTTGGGCtaat                                                                                >  2:333347/1‑149 (MQ=255)
                                                                       gCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAAT‑GGGTGCATCGTCATTTGCC‑GGGGCTGATCACCCTTAA‑TGCTTACCGCACGCTCACCGTTCTG‑CAC‑GACCCCGCCACTTTACGCTTTGGTTGGGCtaataa                                                                              >  1:648537/1‑150 (MQ=255)
                                                                         atttCGTCCCGGTTAAATCAAAACCCCCCCCCCCGGCCCTTTTTGAAT‑GGGTTCATCGTCATTGGCC‑GGGGCTGATCCCCCTTAT‑TGCTTCCCGCACGCTCACCGTTCTG‑CCC‑GACCCCGCCACTTTACGCTTTGGTTGGGCtaataa                                                                              <  2:1209490/145‑1 (MQ=255)
                                                                         atatCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAAT‑GGGTGCATCGTCATTTGCC‑GGGGCTGCTCACCCTTAA‑TGCTTACCGCACGCTCACCGTTCTG‑CAC‑GACCCCGCCACTTTACGCTTTGGTTGGGCtaataa                                                                              >  2:321214/1‑148 (MQ=255)
                                                                         atatCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAAT‑GGGTGCATCGTCATTTGCC‑GGGGCTGATCACCCTTAA‑TGCTTACCGCACGCTCACCGTTCTG‑CAC‑GACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAAc                                                                            >  1:103360/1‑150 (MQ=255)
                                                                           atCGTCGCGTTTAAACCAGCCCCCCCCCCCGCGGCTTGTTTTGAAT‑GGGTGCGTCGTTATTTCCG‑GGGGCTGATCCCCTTTAT‑TGCTTCCCGCTCGCTCCCCGTTCTG‑CAC‑GCCCCCGCCACTTTACGCTTTGGTTGGGCtaataa                                                                              <  2:466367/146‑1 (MQ=255)
                                                                             cgccACGGTTAAACAAAAACCCCCCCCCGCGGCCCGTTTAGAAG‑GGGTACTTCGTCTTTCCCT‑GGGGCTGTTCCCCTTTAT‑GGTTTCCCCCCCGCTCCCGGTCCTC‑CAA‑CCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACatat                                                                        <  1:431100/147‑1 (MQ=255)
                                                                               tCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAAT‑GGGTGCATCGTCATTTGCC‑GGGGCTGATCACCCTTAA‑TGCTTACCGCACGCTCACCGTTCTG‑CAC‑GACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCAt                                                                     >  2:782652/1‑151 (MQ=255)
                                                                                  cGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAAT‑GGGTGCATCGTCATTTGCC‑GGGGCTGATCACCCTTAA‑TGCTTACCGCACGCTCACCGTTCTG‑CAC‑GACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTa                                                                   >  2:400061/1‑150 (MQ=255)
                                                                                   ggTTGAATAAAACCCCCCCCCCGCGGCTCGTTTTGAAT‑GGTTGCATCGTCTTTTGCG‑GGGGCTGA‑CCCCCTTAATTGTTTCCCGCCCGCTCCCCTTTCTC‑CAC‑GCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTaa                                                                  <  2:918586/150‑1 (MQ=255)
                                                                                     ttaaaTCAAACCCCCCCCCCGCGTCCCTTTTTGAAT‑GGGTGCTCCGTCTTTTGCC‑GGGGCTAATCCCCCTTAT‑TCCTTCCCGCCCGCTCACGGTTCTG‑CCC‑GCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTa                                                                   <  1:831198/144‑1 (MQ=255)
                                                                                        aatcccacccccccccCGCGGCCCATTTTGACT‑GTTGGCATGCTCATTGCCC‑GGTGTTTCTCCCCCTTAA‑TGTTTTCCCCCCCCAC‑CCGTTCCA‑CCCGGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAAt                                                              <  2:1162952/139‑1 (MQ=255)
                                                                                          tCAGACCCCCCCCCCGCGGCACGTTTTGAAG‑GGGTGCATCGTCATTTGCC‑GGGGCTGACCCCCCTTAA‑TGCTTCCCGCACGCTCACCTTTCTG‑CAC‑ACCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTa                                                           <  2:13518/150‑1 (MQ=255)
                                                                                            agagCTCCCCCCCGCGCCGCTTTTTGAAG‑GGGTGCATTTCCGTTTCCC‑GGGGCAGACCCCCTTTAA‑TGCTTCCCGCACCCCCACCTTTCCC‑CCA‑GCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAAtt                                                             <  2:383085/146‑1 (MQ=255)
                                                                                                 tCCCCACCGCGGCACGTTTTGAAT‑GGGTGCATCGTCATTTGCC‑GGGGCTGATCACCCTTAA‑TGCTTACCGCACGCTCACCGTTCTG‑CAC‑GACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGtg                                                    >  2:1115925/1‑150 (MQ=255)
                                                                                                  ccccccCGCGGCACGTTTTGAAG‑GGGTGCATCTTCTTTTGCC‑GGGGCTGATCCCCCTTAA‑TCCTTCCCGCCCGCTCCCCGTTCTG‑CCC‑GCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGtg                                                    <  1:1185468/149‑1 (MQ=255)
                                                                                                       ccGCGGCACTTTTTGATG‑GGGTGCTTCTTCTTTTGCC‑GGGGCTGATCCCCCTTAA‑TGCTTCCCGCCCGCTCCCCGTTCTG‑CAC‑GACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGt                                             <  1:322710/151‑1 (MQ=255)
                                                                                                          cgGCACGTTTTGAAT‑GGGTGCATCGTCATTTGCC‑GGGGCTGATCACCCTTAA‑TGCTTACCGCACGCTCACCGTTCTG‑CAC‑GACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCt                                          >  1:1215800/1‑151 (MQ=255)
                                                                                                           ggCACTTTGTAAAG‑GGGTGCATCGTCTTTTGCG‑GGGGCTGACCACCCTTAA‑TGCTTCCCGCACCCTCCCCGTTCTG‑CAC‑GACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTc                                            <  1:1181329/148‑1 (MQ=255)
                                                                                                           ggCACGTTTTGAAG‑GGGTGCATCGTCTTTTGCC‑GGGGCTGATCCCCCTTAA‑TGCTTCCCGCCCGCTCCCCGTTCTG‑CAC‑GCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCt                                          <  1:167769/150‑1 (MQ=255)
                                                                                                           ggCACGTTTGGAAG‑GGGTGCATCTTCTTTTGCG‑GGGTCTGATCCCCCTTAA‑TGCTTCCCGCACGCTCCCCTTTCTG‑CCC‑GCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTcc                                           <  1:1048548/149‑1 (MQ=255)
                                                                                                           gccaggTTTTGAAG‑GGGTGCATCGTCTTTTGCC‑GGGGCTGATCCCCCTTAA‑TGCTTCCCGCCCCCTCCCCGTTCTC‑CAC‑GACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTcc                                           <  1:429167/144‑1 (MQ=255)
                                                                                                             cACTTTTTAAAG‑GGGTGCATCGTCTTTTGCC‑GGGGCTGATCCCCCTTAA‑TTTTTCCCGCCCGCTCCCCGTTCGG‑CAC‑GCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTgg                                        <  1:770367/150‑1 (MQ=255)
                                                                                                              cccttttGAATGGGGTGCATCGTTGTTTGCG‑GGGGCTGACCCCCTTTAT‑TGCGTACCCCCCCCTCACCTTTCTG‑CCC‑CCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTcc                                           <  2:741877/144‑1 (MQ=255)
                                                                                                                gTTTTGAAT‑GGGTGCATCGTCATTTGCC‑GGGGCTGATCACCCTTAA‑TGCTTACCGCACGCTCACCGTTCTG‑CAC‑GACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGcac                                     >  2:858320/1‑150 (MQ=255)
                                                                                                                          gggTGCTTCGTCATTTGCG‑GGGGCTGTCCCCCTTTAA‑TGTTTCCCGCCCGCTCTCCTTTCGC‑CAC‑GCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGa                              <  2:871887/148‑1 (MQ=255)
                                                                                                                              gccttcttctttGCCGGGGGCTCACCCCCCTTAT‑TGCTGCCCGCACGTTCCCCTTTCTC‑CAC‑ACCACCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGaaaaaa                         <  2:617430/141‑1 (MQ=255)
                                                                                                                                aTCGTCATTTGCC‑GGGGCTGATCACCCTTAA‑TGCTTACCGCACGCTCACCGTTCTG‑CAC‑GACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCt                     >  1:23596/1‑151 (MQ=255)
                                                                                                                                  cGTCATTTGCC‑GGGGCTGATCACCCTTAA‑TGCTTACCGCACGCTCACCGTTCTG‑CAC‑GACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTg                    >  1:1250526/1‑150 (MQ=255)
                                                                                                                                   gccATTAGCG‑GGGCCTGTTCCCCCTTAA‑TGCTACCACCTCCCTCACC‑TTCGGCCCC‑GCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGaa                  <  1:1185723/149‑1 (MQ=255)
                                                                                                                                     cttttCCC‑GGGGCTGACCCCCCTAAA‑TATTTCCCGCTCTCTCCTTCTTCTC‑CAC‑CCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTg                    <  1:1255141/145‑1 (MQ=255)
                                                                                                                                     cATTTGCC‑GGGGCTGATCACCCTTAA‑TGCTTACCGCACGCTCACCGTTCTG‑CAC‑GACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAAt                >  2:134136/1‑151 (MQ=255)
                                                                                                                                       tttGCG‑GGGGCTGACCCCCCTTAT‑TGCTTCCCGCACGCTCACCTTTCTC‑CCC‑GCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATc               <  1:941578/150‑1 (MQ=255)
                                                                                                                                        ttGCG‑GGGGCTAACCCCCTTTAA‑TGTTTCCCGCACGCTCCCCTTTCTG‑CAC‑GACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAAt                <  2:535361/148‑1 (MQ=255)
                                                                                                                                           cc‑GGGGCTGATCACCCTTAA‑TGCTTACCGCACGCTCACCGTTCTG‑CAC‑GACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATcacca           >  2:661417/1‑150 (MQ=255)
                                                                                                                                            cgggggCTGATCCCCTTCAT‑TGCTGGCTGCCCCCTCACCGTCCTC‑CCC‑GACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGGTGAAGTCCTGGCACAGCTGGAAAAAAGcc                      <  2:280660/137‑1 (MQ=255)
                                                                                                                                              ggggCTGATCCCCTTTAA‑TTCTTCCCCCACCCTCCCCTTTCTC‑CAC‑GCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATcacca           <  2:325757/148‑1 (MQ=255)
                                                                                                                                                   tGATCACCCTTAA‑TGCTTACCGCACGCTCACCGTTCTG‑CAC‑GACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAgtg    >  2:1107253/1‑150 (MQ=255)
                                                                                                                                                    gATCACCCTTAA‑TGCTTACCGCACGCTCACCGTTCTG‑CAC‑GACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTc  >  2:1200964/1‑151 (MQ=255)
                                                                                                                                                        |                                                                                                                                                     
GTCGATAACCTTTCGCAAGCAGCGTTGGTCAGTCATATTCAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAAT‑GGGTGCATCGTCATTTGCC‑GGGGCTGATCACCCTTAA‑TGCTTACCGCACGCTCACCGTTCTG‑CAC‑GACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTC  >  NC_000913/1684562‑1684858

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