Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,984,7180GC47.6% 30.2 / 24.9 46P292R (CCG→CGG) yqeFshort chain acyltransferase
Reads supporting (aligned to +/- strand):  ref base G (5/17);  new base C (20/0);  total (26/20)
Fisher's exact test for biased strand distribution p-value = 1.08e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.77e-01

CCACCATTGACATTGACCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCC  >  NC_000913/2984568‑2984864
                                                                                                                                                        |                                                                                                                                                   
ccaccaTTGACATTGACCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCg                                                                                                                                                     <  1:139863/151‑1 (MQ=255)
    cattgacattgaCCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGTCGGGTCGCATACACCGtcgc                                                                                                                                                 <  2:223054/151‑1 (MQ=255)
              gacCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGTGCATACAACCCCGCCATTTCCac                                                                                                                                       >  2:479595/1‑149 (MQ=255)
                       cTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCCTACACCCGCGCCATTCCCATCAATGcc                                                                                                                                >  1:466838/1‑148 (MQ=255)
                        tCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGGCTACACCGGCGCCATTCCCCTc                                                                                                                                      >  2:697589/1‑141 (MQ=255)
                          aTCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTGGCCAATAACCGGCCCATTCCCATCCATgac                                                                                                                                >  1:166472/1‑143 (MQ=255)
                          aTCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCt                                                                                                                          <  2:107899/151‑1 (MQ=255)
                           tCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCCTAAACCCGCGCCATTCCCCTCAATGCCCGATCt                                                                                                                          >  2:296001/1‑150 (MQ=255)
                              cACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCGACCACCCGCGCAAATACCATCAATACCGGATCTac                                                                                                                        >  1:497309/1‑149 (MQ=255)
                               aCTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCATCTGCCAGCCTACACGCTCCAGGCAACGGCGTGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCt                                                                                                                          <  2:368672/146‑1 (MQ=255)
                               aCTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCCCATACACCGGCGCAAATCCCATTGATACCGGTACTATAcc                                                                                                                     >  1:258181/1‑151 (MQ=255)
                                     gCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCt                                                                                                               <  1:619040/151‑1 (MQ=255)
                                     gCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCCCGCCATTTCCCTCAATGCCCGATCTACACCAACGc                                                                                                                >  1:465669/1‑150 (MQ=255)
                                       aTCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCATACACCCGCGCAATTCCCATCAATACCGGAACTACACCAAcac                                                                                                                >  1:567029/1‑146 (MQ=255)
                                           tGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTgg                                                                                                             <  1:637469/147‑1 (MQ=255)
                                           tGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCa                                                                                                           <  1:538435/149‑1 (MQ=255)
                                                  cGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCGCCATTCCCCTCAAAGCCGGAGCTACACCAACCCCTGCAAATgcg                                                                                                     >  2:417334/1‑148 (MQ=255)
                                                    aaaGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCCTACACCCGCGCAATTTCCCTCAATTCCGGGTCTACACCAACCCTGGCAAAAGCGCGaa                                                                                                 >  1:285214/1‑150 (MQ=255)
                                                    aaaGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCGCAATTCCCCTCAATGCCGGATCTACACCAACGCTGGCAAAAGCGCGaa                                                                                                 >  1:439779/1‑150 (MQ=255)
                                                        cgcCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATAACACCGCCCCATTCCCCCCCATGCCCGCGCTACACCCACCCTGGCAAAAGCGCGCAAGCg                                                                                             >  2:180631/1‑150 (MQ=255)
                                                        cgcCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCGCCATTCCCCTCCATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGc                                                                                              >  2:243110/1‑149 (MQ=255)
                                                              tgCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGCCGCATACACCCGCGCAAATCCCATCCAAACCGGGACTACACCAACCCTGGCAAAAGCGCGAATGGGGGCCa                                                                                        >  1:521425/1‑149 (MQ=255)
                                                                  gcaAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCGCCATTCCCCTCAATACCGGCACTACACCAACGCTGGCAAAAGCGCGAATGCGGGCCAGCACgg                                                                                  >  1:167578/1‑151 (MQ=255)
                                                                  caaaaCGCTTCATTGGCCTCGATAAGGTCGACTTAGTCCAGCTTCCAGCCGTCACGCTCCAGTCGACGGCGGGGAGCTTTCGTTTTGACGATACCGCGCCCAACCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGc                                                                                      <  2:683345/145‑1 (MQ=255)
                                                                   caAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCACAAACCCGCGCACCTCCCCCCCATGCCCGGTCTACACCCACCCCGGCCAAATCGCGCAAGCGGGCCAGCCCggg                                                                                 >  2:310336/1‑151 (MQ=255)
                                                                    aaaCGCTGAATTAGCCGCGATAAGATCGACTTCAGCAAACTACCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCGATACCGAGCACTGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCa                                                                                     <  2:241479/146‑1 (MQ=255)
                                                                        gCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAAt                                                                            <  1:313498/151‑1 (MQ=255)
                                                                         cTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGCCGCCAATCCCATCAAAGCCCGATCCACACCAAAGCTGGCCAAAACGCGAATGCGGGCCCACACGGGGAAAtt                                                                           >  2:457364/1‑151 (MQ=255)
                                                                           tCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAAtt                                                                           <  1:246804/149‑1 (MQ=255)
                                                                              ttAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTTGCCTACACCCGCGCAATTCCCATCAAAGCCCGGTCTTCACCAACGCCGGCAAAAGCGCGAATGAGGGGCCGCCCGGGGAAATTCaa                                                                        >  1:255709/1‑149 (MQ=255)
                                                                               tAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCACACCACCGCGCCATTCCCCTCCACGGCCGGTCTACACCCAACCTGGCCAAAGCGCGAAAGCGGGCCAGCCCGGGGaaaattcactc                                                                     >  2:559901/1‑143 (MQ=255)
                                                                                     cGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCCGCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTg                                                                <  2:182194/150‑1 (MQ=255)
                                                                                     cGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTg                                                                <  1:548313/150‑1 (MQ=255)
                                                                                      gATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCCTACACCCGCGCAAATCCCATCAATGCCGGATCTACACCAACGCTGGCCAATGCGCGCATGCGGGCCAGCACGGGGAAATTCAACGCTCGTGCt                                                              >  2:217915/1‑151 (MQ=255)
                                                                                       aTAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATTATTTGTAATGATACGCCGACCACCGAGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCtt                                                             <  2:662593/151‑1 (MQ=255)
                                                                                           gATCGACTGCAGCCGACTGCCAGCCTACACGCGCCAGGCAACGGCGGGTCGTATTCTTCTTAATAATACGGCGACCTCCGAGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTg                                                           <  2:644094/149‑1 (MQ=255)
                                                                                            aTCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCGCCATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAAAGCGCGAAAGCGGGCCAGCACGGGGAAATTCAAAGCCCGTGCttt                                                            >  1:156045/1‑147 (MQ=255)
                                                                                              cGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTc                                                      <  1:616425/151‑1 (MQ=255)
                                                                                                      gCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCtcatca                                              <  1:239785/151‑1 (MQ=255)
                                                                                                             gCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCTGAGCGATACGCCGTACTGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTa                                        <  2:154050/150‑1 (MQ=255)
                                                                                                             gCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCACCAATTCCGGGTCTACACCAACCCTGGCAAATGCGCGGATGCGGGCCAGCACGGGGAAAATCAAAGCCCGTGCTTTTGCTTCGCTCATCAttatta                                        >  1:397318/1‑150 (MQ=255)
                                                                                                                aGCCTACACGCTCCAGGCAACGGCGGGGCGCATACACCCGCGCCATTTCCCTCCATACCGGGTCTCCACCAACGCCGGCAAAAGCGCGAAAGCGGGCCAGCCCGGGGAAAATCAACACTCCTGCTTTGGCTTCGCCCCTCATCATTTcc                                      >  2:198255/1‑149 (MQ=255)
                                                                                                                         gCTCCAGGCAACGGCGGGTCGCCTAAACCCGCGCAATTCCCCCCAAAGCCGGATCTCCACCCACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGGAAATACAAAGCCCGCGCCTTTGCCTCGCCCCTCCTCATTACCCCAGCCGAGcc                           >  2:533688/1‑151 (MQ=255)
                                                                                                                              aGGCAACGGCGGGGCGCATAAACCCGGGCCATTCCCCTCAATGCCCGAACTACCCCAACGATGGCAAATGCGCGAATGCGGGCCAGCCCGGGGAAATTCAACACCCGGGCTTTTGCCTCGCTCCACATCATTACCGCAACTGCGcc                           >  2:605892/1‑146 (MQ=255)
                                                                                                                                gCAACGGCGGGTCGCATACACCGGCGCAATGCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGttt                    <  2:484036/151‑1 (MQ=255)
                                                                                                                                     ggcggGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGga               <  1:456693/151‑1 (MQ=255)
                                                                                                                                                    cacCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTAccc  <  1:296074/151‑1 (MQ=255)
                                                                                                                                                        |                                                                                                                                                   
CCACCATTGACATTGACCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCC  >  NC_000913/2984568‑2984864

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