Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913646,7470TG46.9% 61.8 / 17.1 49T255P (ACC→CCC) citG2‑(5''‑triphosphoribosyl)‑3'‑dephosphocoenzyme‑A synthase
Reads supporting (aligned to +/- strand):  ref base T (3/23);  new base G (22/1);  total (25/24)
Fisher's exact test for biased strand distribution p-value = 1.07e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCC  >  NC_000913/646602‑646897
                                                                                                                                                   |                                                                                                                                                      
atTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGa                                                                                                                                                    <  1:612083/150‑1 (MQ=255)
       aaTTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATTgcccc                                                                                                                                            >  2:626573/1‑151 (MQ=255)
       aaTTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGGCGCAAGccccc                                                                                                                                            >  1:704133/1‑151 (MQ=255)
        aTTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCt                                                                                                                                           <  2:576547/151‑1 (MQ=255)
         tttATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGAATGccccc                                                                                                                                            >  2:1024772/1‑149 (MQ=255)
          ttATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGCATACCCCCctt                                                                                                                                         >  1:786161/1‑151 (MQ=255)
           tATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCttt                                                                                                                                         <  2:267579/150‑1 (MQ=255)
           tATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGAACAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGCGGGGTCAATGCCCCCcttt                                                                                                                                        >  1:883916/1‑151 (MQ=255)
                   gCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGGATGCCCCCTTTTTGCaata                                                                                                                                 >  1:513458/1‑150 (MQ=255)
                   gCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATGCCCCCctttt                                                                                                                                       >  2:442934/1‑144 (MQ=255)
                     ttAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGAATGCCCCCCTTTTGtaataa                                                                                                                                >  2:763650/1‑149 (MQ=255)
                            atTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGAATGCCCCCCCTTTGTAATAATtttt                                                                                                                           >  1:297725/1‑147 (MQ=255)
                               aaaTCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCt                                                                                                                     <  2:443051/150‑1 (MQ=255)
                                     gtgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGTCGCCTGGACTGAGATTTCGTTCGATCCACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGCGGTCGAAAGCCCCCTTTTTGCAATACTGTTTGCGCCTCGCGTTg                                                                                                              <  1:330024/151‑1 (MQ=255)
                                     gtgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTcgcg                                                                                                                 <  1:649636/148‑1 (MQ=255)
                                     gtgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCt                                                                                                               <  1:659908/150‑1 (MQ=255)
                                     gtgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCCTTTTTCAATAATGTTTTTGCCCCgcgc                                                                                                                >  2:346342/1‑149 (MQ=255)
                                       gCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCCCATGGCCCCCTTTTgtagataagttt                                                                                                                           >  1:371066/1‑125 (MQ=255)
                                          aaaaaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGTGCCTCGCGCTGTAGcc                                                                                                         <  2:742416/151‑1 (MQ=255)
                                             aaCCAGGTAAGGATCAGGAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGGTGAAAGACCCCCTTTTTTAATAAATTTTGTGCCCcgccgcgcagac                                                                                                         >  1:440009/1‑139 (MQ=255)
                                                 aGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGCATGCCCCCCTTTTTTAATAAAATTTTTGCCCCCCGCCGTAGCCCgcgc                                                                                                    >  2:282885/1‑149 (MQ=255)
                                                    tAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTGGGAATGCCCCCTTTTTCAAAAAAGTTTTGGGCTCCCGCCGTAGGCacccccgg                                                                                                 >  1:301596/1‑142 (MQ=255)
                                                     aaGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGcc                                                                                               <  1:763404/150‑1 (MQ=255)
                                                           cAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTc                                                                                         <  2:589854/150‑1 (MQ=255)
                                                           cAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTc                                                                                         <  1:924499/150‑1 (MQ=255)
                                                           cAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTc                                                                                         <  2:533750/150‑1 (MQ=255)
                                                            aGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGc                                                                                       <  2:741670/151‑1 (MQ=255)
                                                                       aCTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGAAAGCCCCCCTTTTTTAATAAAGTTTTTGCCCCGCCCTGTAGCCAGCGCCGGCCCCCCCCCCCACcccgc                                                                               >  2:396964/1‑148 (MQ=255)
                                                                         tGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATGCCCCCCTTTTTTAATAAATTTTTCGCCCCCCGCCGTAGCCCGCGCCGGcccccccccccccccccc                                                                               >  2:928247/1‑134 (MQ=255)
                                                                            cgccTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCCTTTTTCAATAATGTTTTTGCCCCCCGCTGTAGCCAGCGCAAGCCCCCCCCCCCCCCCCGCGCTGGaa                                                                        >  2:915974/1‑150 (MQ=255)
                                                                              ccTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGACTGCCCCCTTTTTGTAATAAAAATTGCGCCGCGCGCCGTATCCAGCGCAGCCCCCCCCCCCCACCCCGCGCTGCCAcat                                                                     >  1:1034271/1‑149 (MQ=255)
                                                                                      gagaTTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATc                                                              <  2:818376/150‑1 (MQ=255)
                                                                                      gagaTTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCg                                                             <  2:519552/151‑1 (MQ=255)
                                                                                      gagaTTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGAAGCCCCCCTTTTTCAATAATTTTTTTGCCCCGCGCCCTAGCCAGCGCAAGGCCCCCCCGCCCCCCCCCGGCGCAACGTTgg                                                                  >  2:440450/1‑146 (MQ=255)
                                                                                       agaTTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCg                                                             <  1:924931/150‑1 (MQ=255)
                                                                                         aTTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCg                                                          <  2:474997/151‑1 (MQ=255)
                                                                                              gttcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCCTTTTTCAATAAAGTTTTCGCCCCGCGCCGTAGCCCGCGCAGGCCCCCCCCCCCACCCCCCGCAGCCACGGTGGTATAGCCGCTggt                                                     >  1:114018/1‑149 (MQ=255)
                                                                                                 cgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGc                                                  <  2:388806/151‑1 (MQ=255)
                                                                                                     acacTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCAt                                               <  1:944225/150‑1 (MQ=255)
                                                                                                     acacTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCCTTTTTCAAAAATATTTTCGCCTCCCGCTGTAGCCAGCGCCGGCCCCCCCCCCCCCCCCCCGCTGCAACGTTGGGGTCGCCCGTGGTCGcc                                                 >  2:537988/1‑148 (MQ=255)
                                                                                                          cccTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAg                                         <  1:396750/151‑1 (MQ=255)
                                                                                                                     tGCCGGAGATAATCGAGATCGGCGGGGGGTTGAATGCCCCCCTTTTTTAATAATTTTTTCGCCCCCCGCTGTAGCCCGCGCAGGCCCCCCCCCCCCCCCCCCGGTGCAACGTTGGGATCCCCGTTGGTCGCCCTCCGTAGGAGGAAgg                                 >  2:575961/1‑148 (MQ=255)
                                                                                                                      gCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCTTTTTTCAATAATGTTTTCGCCCCCCGCCGTAGCCAGCGCAGGCCCCCCCCCCCCCCCCCCGCTGCAACCTTGGGATCGCCGGTGATCGCCCTCAGTAGGAGCAAGGGAt                              >  2:845977/1‑150 (MQ=255)
                                                                                                                          gagaTAATCGAGATCGGCGGGGGGTCGAAAGCCCCCCTTTTGCAATAATTTTTTCGCCCCGCGCTGTAGCCAGCGCAGGCCCCCCCCGCCCCCCCCCCGCGCAACGCTGGTATCGCCGTTTATCGGCATCAGTAGGAGCAAGGTATCGa                           >  2:101500/1‑149 (MQ=255)
                                                                                                                                  cGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCa                  <  1:442702/150‑1 (MQ=255)
                                                                                                                                     gaTCGGCGGGGGTTCGAAAGCCCCCTTTTTGTAATAAAGTTTTCGCCTCCCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCCCGCGATGCCACGTTGGTATCGCCCTTGATCGCCATCAGTAGGAGCAAGGTCTCCAGCAATGCCCGtt               >  1:796510/1‑150 (MQ=255)
                                                                                                                                        cggcggGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAgg           <  2:201558/151‑1 (MQ=255)
                                                                                                                                             gggggTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCt       <  2:303633/151‑1 (MQ=255)
                                                                                                                                              ggggTTCGAATGCCCCCTTGGTGCAACAATGTTTCCGCCTCGCGCTGTATGCGGCGCTGGCTCCCCTCGCCACCGCGCGATGCAACGTTGGTCTCGCGGTTGATCGCCCTCGGTAGGAGCTAGGTTTCGAGCAATGCCAGTTCAggtt         <  1:334007/148‑3 (MQ=255)
                                                                                                                                                   ttCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAAccc  >  2:325075/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
ATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCC  >  NC_000913/646602‑646897

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