Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913646,7470TG41.4% 139.6 / 50.4 99T255P (ACC→CCC) citG2‑(5''‑triphosphoribosyl)‑3'‑dephosphocoenzyme‑A synthase
Reads supporting (aligned to +/- strand):  ref base T (8/50);  new base G (41/0);  total (49/50)
Fisher's exact test for biased strand distribution p-value = 3.80e-20
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACC  >  NC_000913/646599‑646896
                                                                                                                                                      |                                                                                                                                                           
aaTATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGtt                                                                                                                                                            <  2:816101/150‑1 (MQ=255)
 atatTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCg                                                                                                                                                         <  2:1591828/151‑1 (MQ=255)
   atTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCg                                                                                                                                                         >  1:982734/1‑149 (MQ=255)
     ttCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAAt                                                                                                                                                    <  1:1243006/150‑1 (MQ=255)
       ccAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGcc                                                                                                                                                <  2:1862060/151‑1 (MQ=255)
       ccAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGc                                                                                                                                                 <  2:194647/150‑1 (MQ=255)
       ccAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGCCTGcc                                                                                                                                                >  1:481920/1‑151 (MQ=255)
         aaaTTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGcccc                                                                                                                                              <  1:1175835/151‑1 (MQ=255)
         aaaTTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGcccc                                                                                                                                              <  2:1544892/151‑1 (MQ=255)
         aaaTTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGGAAGcccc                                                                                                                                              >  2:354873/1‑151 (MQ=255)
         aaaTTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGgccc                                                                                                                                              >  1:424843/1‑151 (MQ=255)
          aaTTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGGATGgcccc                                                                                                                                             >  2:2020631/1‑151 (MQ=255)
           aTTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCt                                                                                                                                            <  1:1774120/151‑1 (MQ=255)
            tttATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTTGAATTccccc                                                                                                                                             >  2:1718215/1‑149 (MQ=255)
             ttATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCtt                                                                                                                                          <  2:1022046/150‑1 (MQ=255)
             ttATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGGATTCCCCCctt                                                                                                                                         >  2:108777/1‑151 (MQ=255)
             ttATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCctt                                                                                                                                         >  2:1104968/1‑151 (MQ=255)
                tCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTGGGAAGGCCCCCtttt                                                                                                                                        >  1:525023/1‑150 (MQ=255)
                    gTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCaat                                                                                                                                  <  2:1531747/151‑1 (MQ=255)
                    gTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCaat                                                                                                                                  <  1:795657/151‑1 (MQ=255)
                     tGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGGAGGGCGCCCTTTTTACa                                                                                                                                    >  2:981898/1‑149 (MQ=255)
                     tGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTGGGATGGCCCCttttttcga                                                                                                                                   >  2:1251146/1‑145 (MQ=255)
                         tAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATg                                                                                                                              <  1:1907357/150‑1 (MQ=255)
                         tAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATg                                                                                                                              <  2:1550247/150‑1 (MQ=255)
                         tAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATGCCCCCCTTTTTCcataa                                                                                                                                >  2:1803130/1‑148 (MQ=255)
                          aaaTAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGtt                                                                                                                            <  2:1318007/151‑1 (MQ=255)
                          aaaTAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGAAGCCCCCCTTTTTTAATAATAtt                                                                                                                            >  1:1821036/1‑151 (MQ=255)
                            ataatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGGATGGCCCCCTTTTTtaataaatttt                                                                                                                           >  1:966928/1‑145 (MQ=255)
                              aatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGGTGAATGCCCCCCTTTTTtaataaaattt                                                                                                                           >  1:541394/1‑143 (MQ=255)
                               atTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTgcgc                                                                                                                       <  2:1998148/151‑1 (MQ=255)
                               atTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTgc                                                                                                                         <  1:640977/149‑1 (MQ=255)
                               atTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGAAAGCGCCCCCTTTTCGATAATttttt                                                                                                                          >  1:791033/1‑147 (MQ=255)
                                    aTCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCCTTTTTCaataaaatttttgcccc                                                                                                                    >  2:1014284/1‑137 (MQ=255)
                                        gtgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCgcgc                                                                                                                <  2:1381002/149‑1 (MQ=255)
                                        gtgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCgcgc                                                                                                                <  1:69083/149‑1 (MQ=255)
                                        gtgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCt                                                                                                               <  1:1937854/150‑1 (MQ=255)
                                            tAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCCTTTTTTAATAATGTTTTTGCCTCCCGCCGTAGc                                                                                                          >  2:1486204/1‑151 (MQ=255)
                                             aaaaaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTGGAATGGCCCCCTTTTTTTATAATGTTTTTGCCTCCCGCGGTAGcc                                                                                                         >  1:1139226/1‑151 (MQ=255)
                                                aaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATGCCCCCCTTTTTTTATAATTTTTTCGCCCCCCGCTGGAGcc                                                                                                         >  1:1959021/1‑148 (MQ=255)
                                                aaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAAAGCCCCCCTTTTTCAATAATGTTTTCGCCCCCCGCTGTAGcc                                                                                                         >  2:356894/1‑148 (MQ=255)
                                                  ccAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAgcgc                                                                                                    <  1:249744/151‑1 (MQ=255)
                                                    aGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGCATGCCCCCCTTTTTGAATAATTTTTTCGGCTCGCGGTGTAggc                                                                                                         >  2:1857069/1‑142 (MQ=255)
                                                     ggTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGGAAGCCCCCCTTTTTCAATAATATTTTTGCCCCCCGCCGTAGcc                                                                                                         >  2:1397220/1‑143 (MQ=255)
                                                      gTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGAGGGGAGGCGGGCCCCCTTTTTGTAATATGTTTTCGGCTCCGGCCGTAGccctcaccgg                                                                                                 >  2:264635/1‑142 (MQ=255)
                                                      gTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGc                                                                                               <  2:846569/151‑1 (MQ=255)
                                                        aaGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGccc                                                                                             <  1:890691/151‑1 (MQ=255)
                                                            aTCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCt                                                                                         <  1:1582718/151‑1 (MQ=255)
                                                            aTCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGGATAACCCCCTTTTTTAATTATGTTTTCGCCCCGCCCTGTAGGGAGAGCAGGccccccc                                                                                         >  2:361868/1‑150 (MQ=255)
                                                              cAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCTGAAGTCCCCCTTTTTCGCTAATGTTTTCGCCCCCCGCCGGAGcctccccgtgcccccc                                                                                           >  2:1651801/1‑134 (MQ=255)
                                                              cAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATGGCCCCCTTTTTCAATAATGTTTTTGCCCCCCGCCGGAGCCAGCGCAAGcccccccc                                                                                        >  2:663541/1‑148 (MQ=255)
                                                                gTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGcc                                                                                     <  2:2011547/151‑1 (MQ=255)
                                                                 tAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGcc                                                                                     <  1:1416051/150‑1 (MQ=255)
                                                                  aGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCa                                                                                    <  1:475788/150‑1 (MQ=255)
                                                                  aGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATGCCCCCCTTTTTTAATAATTTTTTCGCCTCCCGCCGTAGCCAGCGCAGGccccccccccccc                                                                                   >  2:1939688/1‑144 (MQ=255)
                                                                   ggTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGAAGCCCCCCTTTTTCCATAAAGTTTTGGCCCCCCGCCGTAGCCACCGCCGGCCCCCCCCGCcccc                                                                                  >  2:1772719/1‑151 (MQ=255)
                                                                      cAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACcgc                                                                                <  2:1330704/150‑1 (MQ=255)
                                                                      cAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATAGCCCCCTTTTTCAAAAAAGTTTTGCGCCCGCGCCGTTACCCACGCCGGcccccccccccccccc                                                                                >  2:1450744/1‑140 (MQ=255)
                                                                       aGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACcgcgc                                                                              <  1:2088622/151‑1 (MQ=255)
                                                                       aGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGAAAGGCCCCCTTTTTTAATAAAGTTTTCGCCCCCCCCTGTAGCCCGCGCCGGgccccccccccccccccc                                                                              >  1:1679201/1‑139 (MQ=255)
                                                                               cgccTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGAATTGCCCCCTTTTTCCATAATGTTTTCGCCCCCCGCCGGAGCCAGCGCAGGCCCCCCCCGCcccccccccccgcca                                                                       >  2:142287/1‑139 (MQ=255)
                                                                                gccTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGGATGCCCCCCTTTTTGAATAAAGTTTTTGCCTCCCCCCGTAGGCAGCGCCGGGCCCCCCCGCCCCCCCCCGGTGCaa                                                                       >  1:655746/1‑149 (MQ=255)
                                                                                gccTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCCTTTTTCAATAAAATTTGCGCCTCGCGCCGTTAGCAGCGCAGGCCCCCCCCCCCCCCCCCCGATGCaa                                                                       >  1:516353/1‑149 (MQ=255)
                                                                                 ccTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGt                                                                    <  1:329722/151‑1 (MQ=255)
                                                                                     gACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCTTCGCCACCGCGCGATGCAACGtt                                                                   <  2:52995/148‑1 (MQ=255)
                                                                                     gACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGAAAGCCCCCCTTTTTCAATAATTTTTTGGCCCCCCGCCGTTAGCAGCGCAGGGCCCCCCCCCCCCCCCCCGGTGCaa                                                                       >  2:1137289/1‑144 (MQ=255)
                                                                                        tGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATCGCCCCCTTTTTCCATAAAAGTTTCGCCTCGCGCCGTATCCAGCGCAGGCCCCCCCCGCCACCCCGCGaaccgacccttg                                                                 >  1:1324248/1‑136 (MQ=255)
                                                                                         gagaTTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATc                                                             <  2:76463/150‑1 (MQ=255)
                                                                                         gagaTTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATc                                                             <  2:958527/150‑1 (MQ=255)
                                                                                            aTTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATGCCCCCCTTTTTCAATAATGTTTTCGCCCCGCGCCGTAGCCAGCGCAGGCCCCCCCCCCCCCCCCGCGGTGCAACGGTGGTATTGcc                                                          >  2:665678/1‑150 (MQ=255)
                                                                                                cgttcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGa                                                     <  1:900832/151‑1 (MQ=255)
                                                                                                   tcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATc                                                   <  2:959246/150‑1 (MQ=255)
                                                                                                      aTACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGcc                                                <  1:132888/150‑1 (MQ=255)
                                                                                                      aTACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCa                                               <  1:611704/151‑1 (MQ=255)
                                                                                                      aTACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGGCGAAAAGCCCCCTTTTTCAAAAAAGTTTTTGCCCCCCGCTGTGGCCCGCGCAGGGCCCCCCCCCCCCCCCGCCGGGCAACGGTGGGATTGCCGGTGGTCGcc                                                >  2:34757/1‑150 (MQ=255)
                                                                                                         cacTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATTCCCCCTCTTTTGCATTATGTTTTTGCCCCGCGCCCTGGACCACGCAGGGCCCCCCCCCCCACCCCCGATGCAACGCTGGTATCGcccattgacgccctc                                             >  2:561101/1‑137 (MQ=255)
                                                                                                           cTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATGCCCCCCTTTTTCCATAATGTTTTCGCCCCCCGCCGTGGCCCGCGCCGGGCCCCCCCCCCCCCGCGCGGCGCAACGTTGGTTTCGCCGTTGATTGCCCACAGt                                          >  2:1830479/1‑151 (MQ=255)
                                                                                                            tCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAg                                           <  1:1801358/149‑1 (MQ=255)
                                                                                                             cccTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTa                                         <  1:1481886/150‑1 (MQ=255)
                                                                                                             cccTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAg                                        <  1:1540094/151‑1 (MQ=255)
                                                                                                                tGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGAATGTCCCCTTTTTTCAATAAAGTTTTCGCCTCCCCCCGTTGCCAGCGCGGGCCCCCCTCCCCCCCCCCCCCGGCCAAGCTGGGCTCGCCGGTTATCGCCACCAGTAgggg                                     >  2:2043803/1‑149 (MQ=255)
                                                                                                                 gTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGGATGCCCCCCTTTTTCAATAATATTTTCGCCTCGCGCCGGCGCCCGCGCAGGGCCCCCCCCCCACCCCGCGCAGCAACGGTGGTGTTGGCGGTGTTCGcccccacaatg                                       >  1:1424904/1‑139 (MQ=255)
                                                                                                                  tCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTGGGAAGCCCCCCTTTTTCAATAATGTTTTTGCCCCCCGCTGTGGCCCGCGCAGGCCCCCCCCCCCCCCCCCCGCCGCAACGCTGGGATCGCCGCTGGTTGCCCTCCGTAgg                                       >  1:954/1‑147 (MQ=255)
                                                                                                                   cGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCaa                                  <  2:1073227/151‑1 (MQ=255)
                                                                                                                      aCTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGt                               <  1:424955/151‑1 (MQ=255)
                                                                                                                          ccGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCg                           <  2:1933434/151‑1 (MQ=255)
                                                                                                                          ccGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCCTTTTTCAATAATGTTTTCGCCTCCCGCTGTAGCCCGCGCAGGCCCCCCCCCCCACCCCCCGATGCAACGGTGGTATCGCCCTTGATCGCCCTCAGTAGGGGCAAgggaacg                           >  2:298629/1‑146 (MQ=255)
                                                                                                                           cGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGa                          <  1:637595/151‑1 (MQ=255)
                                                                                                                             gagaTAATCGAGATCGGCGGGGGGTTGAATGCCCCCTTTTTTCAATAATATTTGGGCCTCCCGCTGTAGCCAGCGCAGGCCCCCCTCGCCCCCCCGCGACGCAACGTTGGTATTGCCGTTGATCGCCCTCAGTAGGGGGAAAGGATCg                           >  2:422162/1‑148 (MQ=255)
                                                                                                                               gaTAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCaa                      <  1:1680208/151‑1 (MQ=255)
                                                                                                                               gaTAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCaa                      <  2:109110/151‑1 (MQ=255)
                                                                                                                                        gaTCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTc             <  2:805809/151‑1 (MQ=255)
                                                                                                                                         aTCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTc             <  1:1684873/150‑1 (MQ=255)
                                                                                                                                         aTCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCa            <  2:1928432/151‑1 (MQ=255)
                                                                                                                                         aTCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCa            <  2:1417786/151‑1 (MQ=255)
                                                                                                                                          tCGGCGGGGGGTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGGATCGAGCAATGCCAGTTc             >  1:1964100/1‑149 (MQ=255)
                                                                                                                                            ggcggGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGGATCGAGCAATGCCAGTTCAgg          >  1:1689923/1‑150 (MQ=255)
                                                                                                                                             gcggGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGAt        <  2:990708/151‑1 (MQ=255)
                                                                                                                                                 ggggTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTa     <  2:1691077/151‑1 (MQ=255)
                                                                                                                                                    ggTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTaa    >  1:1647255/1‑150 (MQ=255)
                                                                                                                                                     gTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAAcc  >  1:1100998/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                           
AATATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACC  >  NC_000913/646599‑646896

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: