Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,684,7490AC43.3% 157.4 / 19.1 98D164A (GAC→GCC) tusinhibitor of replication at Ter, DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (44/11);  new base C (0/42);  total (44/54)
Fisher's exact test for biased strand distribution p-value = 1.89e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCAAAGAAAACTGGAGCG  >  NC_000913/1684605‑1684899
                                                                                                                                                 |                                                                                                                                                      
acaTCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGc                                                                                                                                                  <  4:1335202/151‑1 (MQ=255)
    cAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACt                                                                                                                                              >  1:1077677/1‑151 (MQ=255)
         acccaaaGCCCACGTTCGAGCTTATCGTCACGGTTGAACCAGACCCCCCCCCCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGACCCCCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGCCCCCGCCACTTTACg                                                                                                                                         <  1:495574/146‑1 (MQ=255)
          aCTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGc                                                                                                                                        >  4:9076/1‑151 (MQ=255)
           cTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCt                                                                                                                                       >  1:416442/1‑151 (MQ=255)
           cTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCt                                                                                                                                       >  1:213146/1‑151 (MQ=255)
              aaGACCACGTTCGACCATATCGTCACGGTTGAATCAAACCCCCCCCCCGCGGCACTTTTTGAATGGGTGCATCGTCTTTTGCCGGGGCTGACCCCCTTTAGTGCTTACCGCACGCTCACCGTTCTGCACGCCCCCGCCACTTTACGCt                                                                                                                                       <  4:544092/148‑1 (MQ=255)
                  cccctttggAGCTTATCGTCACGTTTGAATCAGACCCCCCCCCCGCGCCGCTTTTTGAAGGGGTGCATCGTCTTTTGCGGGGGCTAACCCCCTTTATTGCTTCCGGCCCGCTCCCCTTTCTGCACGCCCCCGCCACTTTACGCTttggttg                                                                                                                                <  4:1425901/143‑1 (MQ=255)
                  ccccGTTCGATCATATCTTCAGGGTTGAATCAGACCCCCCCCCCGCGGCACTTTTTGATTGGGTGCTTCTTCTTTTGCGGGGGCTCATCCCCCTTAATGCTTTCCGCACGCTCCCATTTCTGCACGCCCCCGCCACTTTACGCTttggttg                                                                                                                                <  4:594937/148‑1 (MQ=255)
                  ccACTTTCGAGCATATCGTCCCGGTTAAATCAGACCCCCCCCCCGCGGCACTTTTTAAATGGGTGCATCGTCATTTGCGGGGCCTGACCCCCTTTAATGCTTCCCGCCCGCCCCCCTTTCTGCACGCCCCCGCCACTTTACGCTttggt                                                                                                                                  <  1:1140240/149‑1 (MQ=255)
                   cACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTttggttgg                                                                                                                               >  3:7077/1‑151 (MQ=255)
                     cGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGc                                                                                                                             >  1:422163/1‑151 (MQ=255)
                      gTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCt                                                                                                                            >  1:237221/1‑151 (MQ=255)
                       ttCGAGCATATCGTCAGGGTTGAATCAGACCCCCCCCCCGCGGCACGTTTTGAATGGGTGCACCGTCTTTTGCGGGGGCTGATCCCCCTTAATGCTTACCGCACGCTCCCCGTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCt                                                                                                                            <  1:543138/150‑1 (MQ=255)
                         cGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCtaat                                                                                                                         >  1:1200502/1‑151 (MQ=255)
                           aGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCtaata                                                                                                                        >  4:1379507/1‑150 (MQ=255)
                            gCATATCGTCACGGTTGAATCATAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTACTGCTTACCGCACGCTCACCGTTCTGCACGACCCAGCCACTTTACGCTTTGGTTGGGCtaata                                                                                                                        >  3:174697/1‑149 (MQ=255)
                             cATATCGTCGCGGTTAAATCAGACCCCCCCCCCGCGGCATTTTTGAAAGGGGTTCTTCGTCATTTGCCGGGGCTGATCCCCCTTAATGCTCCCCGCGCGCTCCCCGTTCGGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCtaa                                                                                                                          <  1:443198/146‑1 (MQ=255)
                             cATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAAc                                                                                                                     >  4:319896/1‑151 (MQ=255)
                              atatCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAAc                                                                                                                     >  4:531370/1‑150 (MQ=255)
                                  cgtcctggTTGAAGCAGAACCCCCCCCCGCGGCTCTTTTTGAAGGGGGGCTTCTTCTTTTGCCGGGTCTGATCCCCCTTATTGCTTACCGCCCGCCCACCGTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACatat                                                                                                                 <  4:654092/144‑1 (MQ=255)
                                   gTCACGGTTAAATCAGACCCCCCCCCCGCGGCACGTTTTGAGTGGGTGCATCTTTTTTTGCCGGGGCTCATCCCCTTTATTGCTTCCCGCACGCTCACTTTTCTGCCCGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATc                                                                                                                <  4:985626/150‑1 (MQ=255)
                                     cgcGGTTGAACCAGACCCCCCCCCCGCGGCACGTTTTAAAGGGGTGCACGGTCTTTGGCGGGGGCTGACCCCCTTTAATGTTTCCCGCACGCTCCCCGTTCCGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCa                                                                                                               <  1:1042171/147‑1 (MQ=255)
                                     cACGGTTGAATCAGACCCCCCCCCCGCGGCACTTTTTGAATGGGTGCATCGTCTTTTGCGGGGGCTGACCCCCCTTAATGCTTCCCGCCCGCTCCCCGTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACatat                                                                                                                 <  1:1125483/147‑1 (MQ=255)
                                     cACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCAtt                                                                                                             >  4:203030/1‑151 (MQ=255)
                                     cACGGTTAAATCAAACCCCCCCCCCGCGGCACGTTTTGAATGGGTGCATCTTCTTTTGCCGGGGCTGACCCCCCTTAATGCTTCCCGCCCGCTCACCTTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCAtt                                                                                                             <  1:364925/151‑1 (MQ=255)
                                       cGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTa                                                                                                            >  4:689952/1‑150 (MQ=255)
                                         tttaaaTAAGACCCCCCCCCCGCGGCTCTTTTTGAATGGTTTCGTCTTAATTTGCCGGTTATTCTCACTTTTATTCATGCCCGCGGGCTAACTGTTCCGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAg                                                                                                          <  4:944065/146‑1 (MQ=255)
                                             aaTCAGACCCCCCCCCCGCGGCACGTTTTGAATGGGTGCATCTTCTTTTGCCGGGGCTGATCCCCTTTAATGCTTCCCGCACGCTCACCTTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAAttt                                                                                                     <  4:162980/151‑1 (MQ=255)
                                              accAGACCCCCCCCCCGCGGCCCTTTTTGAAGGGGTGCCTCGTCTTTGCCCGGGGCTGATCCCCCTTAATGCTTCCCGCACCCTCCCCTTCCTCCCCGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAAt                                                                                                       <  1:277067/146‑1 (MQ=255)
                                               tCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTAc                                                                                                   >  4:977017/1‑151 (MQ=255)
                                                cAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTAc                                                                                                   >  4:1460122/1‑150 (MQ=255)
                                                 agacccccccccCGCGGCCTTTTTTGAAGGGGTGCTTCGTCTTTTGCGGGGGCTGATCCCCTTTATTGCTTACCGCCCGCTCACCTTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACAt                                                                                                 <  4:926051/145‑1 (MQ=255)
                                                 aaacccccccccCGCGGCACGTTTTGAATGGGTGCTTCTTCTTTTGCGGGGGCTGACCCCCCTTAATGCTTCCCGCCCGCCCCCCTTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTAc                                                                                                   <  1:363725/143‑1 (MQ=255)
                                                 aGAACCCCCCCCCGCGGCACTTTTTGAATGGGTGCTTCGTCATTTGCCGGGGCTGACCCCCTTTAATGCTTCCCGCACGCTCACCTTTCTGCACGACACCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTa                                                                                                    <  4:149404/148‑1 (MQ=255)
                                                   acccccccccCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCGGGGGCTGATCCCCCTTAATGCTTCCCGCACGCTCACCGTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCg                                                                                               <  1:1058096/147‑1 (MQ=255)
                                                     ccccccccccCGGCACTTTTTAAATGGGTGCTTCTCCTTTTCCGGGGCGTATTCCCCCTAAGTGCTTCCCGCCCCCCCCCCGTTCTGCCCGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGtg                                                                                             <  1:269895/149‑1 (MQ=255)
                                                     ccccccccCGCGGCCCTTTTTGAATGGGTGCATCGTCTTTTGCCGGGGCTAACCCCCCTTAATGCTTCCCGCCCGCTCACCGTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGt                                                                                              <  1:1136427/148‑1 (MQ=255)
                                                       ccccACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGtgat                                                                                           >  1:1354363/1‑151 (MQ=255)
                                                        cccACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGtgatg                                                                                          >  4:1294411/1‑151 (MQ=255)
                                                        cccACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGtgatg                                                                                          >  1:48169/1‑151 (MQ=255)
                                                          ccccgcgcccttttttAAATGGGTTCTTCCTCTTTGGCCGGGGTTCACCCCCTTTTAGCCTTACCACACGCCAACAATTCTTCCCATCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGaa                                                                                        <  3:60276/139‑1 (MQ=255)
                                                           aCCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGtgatga                                                                                         >  4:648604/1‑149 (MQ=255)
                                                            ccGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGt                                                                                      >  1:278502/1‑151 (MQ=255)
                                                                gccatgTTTTGAATGGGTGCATCGTAATTTGCCGGGGCTGACCCCCTTTATTGCTTCCCGCCCGCTCCCCGTTCTCCCCGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTg                                                                                  <  4:1452412/146‑1 (MQ=255)
                                                                  ctttttttGATGGGGTGCCTCGTCTTTTGCGGGGTCTGTCCCCCTTTAATGCTTCCCGCCCGCCCACCTTTCTTCGCACCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCt                                                                                   <  3:107416/144‑1 (MQ=255)
                                                                  cACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGc                                                                                >  4:1429389/1‑151 (MQ=255)
                                                                   aCGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGc                                                                                >  4:299423/1‑150 (MQ=255)
                                                                   aCGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGc                                                                                >  1:300429/1‑150 (MQ=255)
                                                                      ttttGAAGGGGTGCATCGTCTTTTGCCGGGGCTGATCCCCTTTAATGCTTCCCGCACCCTCACCTTTCTCCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGca                                                                               <  1:1262461/148‑1 (MQ=255)
                                                                      ttttAAAGGGGTGCATCTTCTTTTGCCGGGGCTGATCCCCTTTATTGCTTCCCGCCCGCCCACCTTTCTCCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGcaca                                                                             <  1:556570/150‑1 (MQ=255)
                                                                      tttgGAAGGGGTGCATCGTCATTGGCGGGGGCTGACCCCCCTAAATGCTTACCGCACGCTCCCCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGca                                                                               <  1:658058/148‑1 (MQ=255)
                                                                       tttGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAg                                                                            >  1:279850/1‑150 (MQ=255)
                                                                           aaTGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTgg                                                                        >  4:419287/1‑150 (MQ=255)
                                                                           aaTGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTgg                                                                        >  4:796560/1‑150 (MQ=255)
                                                                              gggTGCTTCTTCTTTTGCGGGGGCTGACCCCCTTTATTCCTTCCCGCCCGCCCCCCTTTCGGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGaaaa                                                                    <  4:1276270/151‑1 (MQ=255)
                                                                               ggTGCATCGTCATTTGCGGGGGCTGATCCCCCTTAATGCTTCCCGCCCGCTCCCCGTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGaaaaa                                                                   <  1:619991/151‑1 (MQ=255)
                                                                                    aTCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCt                                                              >  1:979997/1‑151 (MQ=255)
                                                                                      cGTCTTTTGCCGGGGCTGATCACCTTTAATGCTTCCCGCCCCCCCACCTTTCTCCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGc                                                                <  4:63072/147‑1 (MQ=255)
                                                                                       gTCATTTGCGGGGGCTGACCCCCCTTAATGCTTCCCGCGCGCTCCCCTTTCTCCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTg                                                             <  4:1172181/149‑1 (MQ=255)
                                                                                         cacttGCCGGGGCAGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACACCTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGaaa                                                          >  3:97803/4‑150 (MQ=255)
                                                                                           tttGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATc                                                        >  3:94352/1‑150 (MQ=255)
                                                                                              gCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATcacc                                                     >  1:1356940/1‑150 (MQ=255)
                                                                                               ccGGGGCTGACCCCCCTTAATGCTTCCCGCCCGCTCCCCGTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATcac                                                      <  1:1084718/148‑1 (MQ=255)
                                                                                                 ggggCTGATCCCCCTTATTGCTTCCCGCCCGCTCCCCTTTCTGCCCGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACg                                                  <  1:208225/150‑1 (MQ=255)
                                                                                                 ggggCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGc                                                 >  1:1189814/1‑151 (MQ=255)
                                                                                                     cTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAgtg                                             >  4:868705/1‑151 (MQ=255)
                                                                                                         aTCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCg                                          >  4:66220/1‑151 (MQ=255)
                                                                                                           cccccTTAATGCTTCCCGCACGCTCCCCTTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGc                                         <  4:91876/148‑1 (MQ=255)
                                                                                                             cctttAATGTTTCCCGCCCACTCCCCTTTCTGCGCGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAgtgt                                            <  4:1274786/142‑1 (MQ=255)
                                                                                                                ttAATGCTTACCGCACGCTCACCGTTCTGCCCGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGt                                    <  4:828266/150‑1 (MQ=255)
                                                                                                                ttAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTg                                   >  1:1000864/1‑151 (MQ=255)
                                                                                                                ttAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTg                                   >  4:1065359/1‑151 (MQ=255)
                                                                                                                 tattGCTCACCCCCCGCTCCCCGTTTTACACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGt                                    <  4:213132/146‑1 (MQ=255)
                                                                                                                 tAATGCTTCCCGCACGCTCACCGTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTGCATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTg                                   <  1:219704/150‑1 (MQ=255)
                                                                                                                 tAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTgg                                  >  4:805605/1‑151 (MQ=255)
                                                                                                                 tAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTgg                                  >  4:972972/1‑151 (MQ=255)
                                                                                                                  aaTGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTgg                                  >  2:187363/1‑150 (MQ=255)
                                                                                                                    tGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGAcg                               >  1:159555/1‑151 (MQ=255)
                                                                                                                       ttcccgcccgcTCCCCGTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGAc                                <  1:1294093/144‑1 (MQ=255)
                                                                                                                       ttcccgcccgcCCCCCGTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGAcgc                              <  1:839248/146‑1 (MQ=255)
                                                                                                                        tcccccccccccACCTTTCTGCCCGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGAcg                               <  4:751344/136‑1 (MQ=255)
                                                                                                                          ccgcccgcTCCCCGTTCTGCCCGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCgag                         <  4:1239234/151‑1 (MQ=255)
                                                                                                                          ccgcccgcCCCCCTTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGAcgcg                             <  1:1452211/147‑1 (MQ=255)
                                                                                                                          ccGCACGCTCCCCTTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGAcgcgcg                           <  1:505100/149‑1 (MQ=255)
                                                                                                                          ccGCACGCTCACCTTTCTGCCCGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGAcgcgcg                           <  1:567678/149‑1 (MQ=255)
                                                                                                                                  tCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGc                  >  1:545908/1‑150 (MQ=255)
                                                                                                                                    aCCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCaa                >  1:1355021/1‑150 (MQ=255)
                                                                                                                                     cctttCTGCCCGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCaa                <  4:1080973/146‑1 (MQ=255)
                                                                                                                                     ccGTTCTGCCCGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCaaa               <  1:951213/150‑1 (MQ=255)
                                                                                                                                     ccGTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCaaa               <  4:423089/150‑1 (MQ=255)
                                                                                                                                     ccGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCaa                >  4:618841/1‑149 (MQ=255)
                                                                                                                                        ttCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCAAAGaa            <  4:1414130/150‑1 (MQ=255)
                                                                                                                                          cTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCAAAGaaa           <  1:471653/149‑1 (MQ=255)
                                                                                                                                             cACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCAAAGAAAACt        <  1:702816/149‑1 (MQ=255)
                                                                                                                                             cACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCAAAGAAAACTg       <  1:824933/150‑1 (MQ=255)
                                                                                                                                               cGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCAAAGAAAACTgg      >  4:675614/1‑149 (MQ=255)
                                                                                                                                               cGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCAAAGAAAACTGGa     <  1:170840/150‑1 (MQ=255)
                                                                                                                                                 aCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCAAAGAAAACTGGAGCg  >  4:498382/1‑151 (MQ=255)
                                                                                                                                                 |                                                                                                                                                      
ACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCAAAGAAAACTGGAGCG  >  NC_000913/1684605‑1684899

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