Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,684,7120AC39.7% 123.8 / 16.3 74T152P (ACC→CCC) tusinhibitor of replication at Ter, DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (31/13);  new base C (0/29);  total (31/43)
Fisher's exact test for biased strand distribution p-value = 1.45e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGATAACCTTTCGCAAGCAGCGTTGGTCAGTCATATTCAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCAC  >  NC_000913/1684564‑1684862
                                                                                                                                                       |                                                                                                                                                          
cGATAACCTTTCGCAAGCAGCGTTGGTCAGTCATATTCAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCAcc                                                                                                                                                          >  2:255922/1‑151 (MQ=255)
  aTAACCTTTCGCAAGCAGCGTTGGTCAGTCATATTCAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCt                                                                                                                                                       >  2:437096/1‑151 (MQ=255)
      ccTTTCGCAAGCAGCGTTGGTCAGTCATATTCAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAAt                                                                                                                                                   <  1:1088703/151‑1 (MQ=255)
      ccTTTCGCAAGCAGCGTTGGTCAGTCATATTCAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAAt                                                                                                                                                   <  1:849138/151‑1 (MQ=255)
         ttCGCAAGCAGCGTTGGTCAGTCATATTCAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATg                                                                                                                                                  <  2:996993/149‑1 (MQ=255)
         ttCGCAAGCAGCGTTGGTCAGTCATATTCAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATg                                                                                                                                                  <  1:454275/149‑1 (MQ=255)
           cGCAAGCAGCGTTGGTCAGTCATATTCAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCCCCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTa                                                                                                                                             <  1:92038/151‑1 (MQ=255)
           cGCAAGCAGCGTTGGTCAGTCATATTCAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCtt                                                                                                                                              >  1:981167/1‑150 (MQ=255)
           cGCAAGCAGCGTTGGTCAGTCATATTCAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTa                                                                                                                                             <  1:106025/151‑1 (MQ=255)
            gCAAGCAGCGTTGGTCAGTCATATTCAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTAc                                                                                                                                            >  1:406774/1‑151 (MQ=255)
             cAAGCAGCGTTGGTCAGTCATATTCAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTAcc                                                                                                                                          >  2:1217391/1‑151 (MQ=255)
                gcagcGTTGGTCAGTCATATTCAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCa                                                                                                                                       <  1:771082/151‑1 (MQ=255)
                 cagcGTTGGTCAGTCATATTCAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCAc                                                                                                                                      >  1:408942/1‑151 (MQ=255)
                      ttGGTCAGTCATATTCAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTc                                                                                                                                  >  2:740560/1‑150 (MQ=255)
                           cagtcaTATTCAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGtt                                                                                                                            <  2:353719/151‑1 (MQ=255)
                             ctcaTATTCAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCt                                                                                                                          >  2:615630/2‑151 (MQ=255)
                                     cAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGAcc                                                                                                                 <  2:396938/151‑1 (MQ=255)
                                            tCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCAc                                                                                                          >  1:443906/1‑151 (MQ=255)
                                               aTAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACttt                                                                                                       >  2:202013/1‑151 (MQ=255)
                                                 aaccccAAGCCCCCGTCCGAGCATATCGTCACGGTTGAATCAGACCCCCCCCCCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCGGGGCCTGATCCCCTTTATTCCTTCCCGCCCGCCCCCCGTTCGGCCCGCCCCCGCCACTTTAc                                                                                                     <  1:721219/146‑1 (MQ=255)
                                                 aaaCTCAAGACCACGTTCGGGCATATCGTCACGTTTGAATCAGACCTCCCCCCCGCGGCACGTTTTGAATGGGTGCTTCGTCTTTTGCCGGGGCTGATCCCCTTTATTGCTTACCGCCCGCTCACCGTTCTGACCGCCCCCGCCACTTTAc                                                                                                     <  1:654411/151‑1 (MQ=255)
                                                       aaGCCCGCGTTCGAGCATATCGCCGGGTTTGAATCAAACCCCCCCCCCGCGGCGCGTTTTAAAGGGGTGCATCGTCATTTGCCGGGGCTGATCCCCCTTATTGCTTCCCGCCCGCTCACCGTTCTGCACGCCCCCGCCACTTTACGCTttg                                                                                               <  2:273655/151‑1 (MQ=255)
                                                       aaGCCCACGTTCGAGCATATCGTCAGGGTTAAATAAAACCCCCCCCCCGCGGCACTTTTTGAATGGGTGCATCGTCTTTTGCGGGGGCTGATCCCCTTTAAGGCTTCCCGCACGCTCCCCTTTCCGCACGCCCCCGCCACTTTACGCTttg                                                                                               <  1:1286496/151‑1 (MQ=255)
                                                           ccgctttCGAGCATTCCGTCAGGGTTGAACCAGACCCCCCCCCCGCGGCACTTTTTAAGGGGGTCCATCGTCTTTTGCGGGGGCTGCCCCCCTTTATTCCTTCCCGCCCGCACCCCTTTCGGCACGCCCCCGCCACTTTACGCTttgg                                                                                              <  1:911474/143‑1 (MQ=255)
                                                            cgcGTTGGAGCATATCGTCACGTTTGAACCAGACCCCCCCCCCGCGGCACTTTTTGAATGGGTGCTTCGTCTTTTGCCGGGGCTGATCCCCCTTAATGCTTCCCGCACCCTCCCCGTTCTGCACGCCCCCGCCACTTTACGCTttggt                                                                                             <  2:849130/146‑1 (MQ=255)
                                                             aCGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTggg                                                                                         >  1:293738/1‑151 (MQ=255)
                                                             aCGTTCGACAATGTCGTCAGGGTTGAATCAAACCCCCCCGCCGCGGAGCTTTTGGAAGGGGTGCTTCGTCTTTGCCCGGGGCTGATCACCCTTAATCCTTCCCCCCCGCTCCCGGTTCTGCCCGCCCCCGCCACTTTACGCTttggttgg                                                                                          <  2:1128101/150‑1 (MQ=255)
                                                               gTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCt                                                                                       >  2:502699/1‑151 (MQ=255)
                                                                ttCGAGCATATCGTCACGGTTGAACCAGACCCCCCCCCCGCGGCTCGTTTTGATGGGGTGCATCGTTTTTTGCCGGGGCTGATCCCCTTTAATGCTTCCCGCCCGCTCCCCTTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCt                                                                                       <  2:347945/150‑1 (MQ=255)
                                                                 tCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCta                                                                                      >  1:597917/1‑150 (MQ=255)
                                                                  cGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCtaat                                                                                    >  2:1212352/1‑151 (MQ=255)
                                                                   gAGCATTTCGTCGCGGTTGAATCAGACCCCCCCCCCGCGGCACGTTTTGAGGGGGTGCATCGTCTTTTCCCGGGGCTGATCCCCTTTAATGTTTCCCGCCCGCCCCCCGTTCGGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCt                                                                                       <  1:367577/147‑1 (MQ=255)
                                                                    aGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCtaata                                                                                   >  2:648654/1‑150 (MQ=255)
                                                                       atatCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAAc                                                                                >  1:1126705/1‑150 (MQ=255)
                                                                        tatGGTCACTGTTAAATCGGACCCCCCCCCCGCGCCATTTTTTAAATGGGGCCTTCGTCTTTTGCGGGGGCTGACCCCCCTTAATTCTTCCCGCACGCTCCCGTTCCGCCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAAc                                                                                <  1:2870/149‑1 (MQ=255)
                                                                           cGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATc                                                                           >  1:531209/1‑151 (MQ=255)
                                                                              cACGGTTGAATCAGACCCCCCCCCCGCGGCACGTTTTAAATGGGTGCATCGTCTTTTGCCGGGGCTGATCCCCCTTAATGCTTCCCCCCCGCTCCCCTTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCAt                                                                         <  2:688961/150‑1 (MQ=255)
                                                                              cACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCAtt                                                                        >  1:188743/1‑151 (MQ=255)
                                                                                   ttGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGAGCATCGTCATTTGCCGCGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTCCTGCACGACCCCGCCACTTTACGCTTTGGCTGGGCTAATAAACATATCATTAAGa                                                                    >  2:1146466/1‑150 (MQ=255)
                                                                                     gAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAAtt                                                                 >  1:756236/1‑151 (MQ=255)
                                                                                      aaTCAGACCCCCCCCCCGCGGCCCGTTTTAAATGGGTGCATCGTCTTTTCCGGGGCCTGATCCCCCTTAATGCTTCCCGCCCGCTCCCCTTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAAt                                                                  <  2:985287/149‑1 (MQ=255)
                                                                                      aaTCAGACCCCCCCCCCGCGGCACGTTTTGAATGGGTGCATCGTCTTTTGCCGGGGCTGATCCCCTTTAATGCTTCCCGCTCGCTCACCTTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGaa                                                                   <  1:789784/148‑1 (MQ=255)
                                                                                      aaTCAGACCCCCCCCCCGCGCCACGTTTTAAAGGGGGGCTCCGTTTTTGGCCGGGGCTGACCCCCTTTAATGCTTACCGCACGCTCCCCGTTCTGCCCACCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAAt                                                                  <  1:1047579/149‑1 (MQ=255)
                                                                                        tCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTAc                                                              >  1:542098/1‑151 (MQ=255)
                                                                                          aaaaCCCCCCCCCGCGGCACGTTTTGAAGGGGTCCATCGTCATTTGCGGGGGCTGATCCCCCTTAATTCTTCCCGCCCGCCCCCCGTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTa                                                               <  2:469536/146‑1 (MQ=255)
                                                                                            acccccccccccCGCCACGTTTTGAGGGGGTGCATCGTCTTTTGCCGGGGCTGACCCCCTTTAATCCTTCCCCCGCGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATc                                                           <  1:83008/146‑1 (MQ=255)
                                                                                            accTCCCCCCGGCGGCACTTTTTGATTGGTTTCACCGTCTTTGCCGGGGTCTATTCCCCCTTTATCGTTACCGAACGCGCCCTTTTGTGCTCGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACAt                                                            <  2:144384/147‑1 (MQ=255)
                                                                                              ccccccccCGGGCCACTTTTTGAAGGGGTAATTCTTTTTTTCCCGGGCCTCTCCCCCTTTATTCTTTCCCGCACGCTCACGGTCCTCAACGACACCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGtg                                                        <  1:79061/149‑1 (MQ=255)
                                                                                                ccccACCGCGCCCTGTTTGGAAGGGGTTCATCTCCATTTGCGGGGACCCATCTCCTTTATTCATTACCGCAAGCTCTCTTTTCCGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGaa                                                                   <  2:502262/138‑1 (MQ=255)
                                                                                                 ccccccGTGGCCCGTTTCGATTGGGTCCAGCGTCTTTGGCGGGGGTTGATCACCTTACAGGTTTACCGCACGCTGACCTTTCTCCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGtgat                                                      <  2:158888/150‑1 (MQ=255)
                                                                                                 cccACCGCGGCAGGTTTTAAATGGGTCCATCGTCTTTTGCGGGGGCTGACCACCCTTATGGCTTCCCGCACGCTCACCGTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGtgatg                                                     <  2:490318/151‑1 (MQ=255)
                                                                                                     cgccgGCAGGTTTTGAGTGGTTGACTGGTCATTTGCTGGGCCTGCCCCCCTTTTTGCCTACGCGCACGCCCAGTTTTGTGGACTACACCGCCTCTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGt                                                 <  1:42596/148‑1 (MQ=255)
                                                                                                     ccGGGGCACTTTTTGAAGGGGGGCTTCTTCTTTTGCGGGGGCTAATCCCCCTTAATGCTTCCCGCACGCCCCCCTTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGt                                                 <  2:925626/151‑1 (MQ=255)
                                                                                                     ccGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGt                                                 >  1:718852/1‑151 (MQ=255)
                                                                                                      cgcgGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTc                                                >  2:637139/1‑151 (MQ=255)
                                                                                                       gcgGCCCTTTTTGAATGGGGGCATGGTAATTTGCCGGGGCTGATCACCTTTAATGCTTCCCGCACGCTCCCCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTcc                                               <  2:1046022/151‑1 (MQ=255)
                                                                                                         ggCACGTTTTGAAGGGGGGCTTCGTCTTTGCCGGGGGCTGATCCCCCTAAATGCTTCCCGCACGCCCCCCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCt                                              <  2:1121157/150‑1 (MQ=255)
                                                                                                           gcatttttGAAGGGGGCGCATCGTATTTTGCCGGGGCTTACCCCCCTTAATGCTACCCGCACGCACCCCTTTCGGCCCGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTg                                             <  2:805450/146‑1 (MQ=255)
                                                                                                           cACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTgg                                            >  2:138563/1‑150 (MQ=255)
                                                                                                            aCGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGc                                           >  2:494464/1‑150 (MQ=255)
                                                                                                               ttttGAATGGGTGCATCGTCATTTGCCGGGGCTGATCCCCTTTAATGCTTCCCGCACGCTCCCCGTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGcaca                                        <  2:531224/150‑1 (MQ=255)
                                                                                                                   gAATGGGTACATGGTCATTTGCTGGTGCCAGTCCCCTTTAATGCTTACCGCTCTCTCCCCGTCCTACCCGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCt                                     <  2:337853/149‑1 (MQ=255)
                                                                                                                        ggtgccATCGTCATTTGCCGGGCCTGATCCCCCTTAGTACTTCCCGCCCGCTCCCCGTTCCCCACCACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGaaaa                               <  2:837892/146‑1 (MQ=255)
                                                                                                                        gggTCCATCGTCATTTGCGGGGGCTGACCCCCCTTATTGCTTCCCGCCCCCCCCCCTTTCTCCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGa                                  <  1:437085/148‑1 (MQ=255)
                                                                                                                         ggTGCTTCTTCTTTTGCGGGGGCTATCCCCCTTTATTTCTTCCCGCACGCTCCCCTTTCTGCCCGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTg                                    <  1:167146/145‑1 (MQ=255)
                                                                                                                                cGTCATTTGCCGGGGCTGATCCCCCTTAATGCTCCCCCCGCGCTCACCGTTCTCCAGGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCt                         <  2:6036/149‑1 (MQ=255)
                                                                                                                                cGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTg                        >  1:65529/1‑150 (MQ=255)
                                                                                                                                cGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTg                        >  1:272842/1‑150 (MQ=255)
                                                                                                                                  tCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGaa                      >  1:518272/1‑150 (MQ=255)
                                                                                                                                       tGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATcacc                >  2:1020829/1‑151 (MQ=255)
                                                                                                                                        gCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATcacc                >  1:986350/1‑150 (MQ=255)
                                                                                                                                         ccGGGGTTAACCCCCTTTAATGCTTCCCGCCCGCTCACCGTTCTCCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATcac                 <  1:460395/148‑1 (MQ=255)
                                                                                                                                         ccGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATaac                 >  1:211081/1‑148 (MQ=255)
                                                                                                                                         ccGGGGCGGACCCCCCTGAAGTCTTACCCCACGGCCACCGTCCGCCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATcac                 <  2:1030497/148‑1 (MQ=255)
                                                                                                                                                       aCCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCAc  >  2:843929/1‑151 (MQ=255)
                                                                                                                                                       |                                                                                                                                                          
CGATAACCTTTCGCAAGCAGCGTTGGTCAGTCATATTCAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCAC  >  NC_000913/1684564‑1684862

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