Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,984,7180GC36.5% 72.5 / 24.6 56P292R (CCG→CGG) yqeFshort chain acyltransferase
Reads supporting (aligned to +/- strand):  ref base G (7/26);  new base C (19/0);  total (26/30)
Fisher's exact test for biased strand distribution p-value = 1.72e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01

CCATTGACATTGACCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCT  >  NC_000913/2984571‑2984867
                                                                                                                                                    |                                                                                                                                                       
ccaTTGACATTGACCCGACGCTCATCCCACTCAAGCATCGTGCCAACCGAAAGCGCCTGGGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGgc                                                                                                                                                       <  2:180947/150‑1 (MQ=255)
 cattgacattgaCCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGgcgc                                                                                                                                                    <  1:1157220/151‑1 (MQ=255)
 cattgacattgaCCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGgcgc                                                                                                                                                    <  1:897247/151‑1 (MQ=255)
  attgacattgacCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCgg                                                                                                                                                        >  1:1107353/1‑147 (MQ=255)
     gacattgacCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAAtt                                                                                                                                                <  2:384739/151‑1 (MQ=255)
      acattgacCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCCTACCCCGCCGCaa                                                                                                                                                  >  1:291386/1‑148 (MQ=255)
          tgacCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCa                                                                                                                                            <  1:645994/150‑1 (MQ=255)
           gacCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCAt                                                                                                                                           <  1:1212929/150‑1 (MQ=255)
                  cGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGGCTAAACCCGCGCCATTTCCCTCAATGcc                                                                                                                                    >  1:248048/1‑150 (MQ=255)
                        tCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCAACCACCCCCGCACCTCCCACCAATACCgg                                                                                                                                  >  2:493226/1‑146 (MQ=255)
                                aaGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCATACAACCGCGCAAATTCCCCCAAAGCCGGATCTACACCa                                                                                                                        >  1:1187912/1‑148 (MQ=255)
                                aaGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCGTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCGACGGCTGGTCGCATTCTCCGGCGCAATTCGCAGCAATGCCGGATCTACACCaa                                                                                                                       <  2:1019415/149‑1 (MQ=255)
                                 aGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGc                                                                                                                    <  1:643847/151‑1 (MQ=255)
                                     tCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCgg                                                                                                                                                        >  1:872638/1‑112 (MQ=255)
                                      cTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCaa                                                                                                                                                  <  2:872086/116‑1 (MQ=255)
                                         gCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCaaa                                                                                                             <  1:829899/150‑1 (MQ=255)
                                         gCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATAAACCGGCGCAATTCCCACCAATTCCGGATCTACACCAACGCTGGCcaa                                                                                                             >  1:628892/1‑150 (MQ=255)
                                           cAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCATACACCCGCGCCAATCCCCTCAACACCCGGGCTACACCCACGCCgggcaatac                                                                                                          >  2:161779/1‑144 (MQ=255)
                                           cAACCGAAAGCCCCTGTGCAGTAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACTCTCCAGGCAACGGCGGGTCGCATACACCGTCACAATTCCCATCAATGCCAAATCTACACCAACGCTGGCa                                                                                                               <  2:109757/146‑1 (MQ=255)
                                            aaCCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATgcg                                                                                                         <  1:403390/151‑1 (MQ=255)
                                             acggAAAGCGGCTGTGCAGCAAACGATTCATTAACATCGATAACATCGACTTGCGCCACCCTGCCGCCTACACGCGCCACGCAACGGCTGGTCGCATACACCTTGACAATTCCGCTCAATCCGGGATCTACACCAACGCTGGCAAATgcgc                                                                                                        <  2:1218723/148‑1 (MQ=255)
                                                gAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATAAACCCGCGCAATTCCCCTCAATGCCGGATCTACACCAAAGCTGGCAAAAgcg                                                                                                         >  1:883028/1‑147 (MQ=255)
                                                 aaaGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCGCCATTCCCATCAATGCCGGGTCTACACCAACGCCGGCAAAAGCGCGaa                                                                                                     >  1:841959/1‑150 (MQ=255)
                                                  aaGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCAAACACCCGCGCAAATCCCCTCAAAGCCGGGTCTACACCAAACCCGGCAAAAgc                                                                                                          >  1:504365/1‑144 (MQ=255)
                                                      gcCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCCCAATTCCCCTCAATGCCGGGTCTTAACCAAACCTGGCCAAATCGCGAATGCggg                                                                                               >  1:1077751/1‑151 (MQ=255)
                                                       ccTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATAAACCCGCGCAATTCCCATCAATGCCGGGTCTACACCAACCCCGGCAAAAAAGCGAATGCggg                                                                                               >  1:48913/1‑150 (MQ=255)
                                                          gtgCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCa                                                                                            <  1:920271/150‑1 (MQ=255)
                                                               gcaAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCGCCATTCCCATCAAAGCCGGGTCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGc                                                                                         >  1:137862/1‑147 (MQ=255)
                                                                  aaCGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCCCATACAACCGCCCCATTCCCCTCACAGCCGGATCTACACCAACGCTCGCACAAGCGCGGAAGCGGGGCAGccacgg                                                                                    >  2:484827/1‑144 (MQ=255)
                                                                     gCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCACACACCCGCGCCATTCCCCTCAATTCCGGAGCTACACCAACGCTGGCCAAAGCCCGAACGCGGGCCGCCACGGGGaaa                                                                                >  1:753671/1‑150 (MQ=255)
                                                                       ttCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCGCAATTCCCATCCATGCCGGGTCTACACCAACCCTGGCAAATGCGCGAAAGCGGGCCAGCACCGGGaaa                                                                                >  1:387810/1‑148 (MQ=255)
                                                                         cATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCGCAATTCCCATCAATGCCGGATCTTCACCAACGCCGGCAAAAGCGCGAAAGCGGGCCAGCACGGGGAAATTCa                                                                            >  1:1118611/1‑150 (MQ=255)
                                                                         cATTAGCCTCGATAAGAGCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATTCACCGGCATAATACGCAGCACTGCGAGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCaa                                                                           <  2:1111702/151‑1 (MQ=255)
                                                                            tAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACg                                                                         <  2:1041450/150‑1 (MQ=255)
                                                                               ccTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCg                                                                     <  1:754842/151‑1 (MQ=255)
                                                                                  cGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCGCCAGGCAACGGCGGGTCGCATACTTCTGAGCGATACCCCGCACAGCCAGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTg                                                                   <  2:360743/150‑1 (MQ=255)
                                                                                  cGATAAGATCGACTTCAGCAAACTGCCAGCCTGCACGCTCCAGGCAACGGCGGGTCGCATTCTCCGGCGTAATACCCAGAAAAGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCt                                                                       <  2:1260448/146‑1 (MQ=255)
                                                                                   gATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCGCAATTCCCCTCAATGCCGGATCTACACCAACGCCGGCAAATGCGCGCATGAGGGCCCGCACCGGGAAATTCAACGCTCGTGc                                                                  >  1:91499/1‑150 (MQ=255)
                                                                                    aTAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGGATAAACCCGCCCCATTCCCCTCAAAGCCCGGACTACACCCACCCTGGGAAAAGCGCGAATGCGGGGCAGCACAGGGGAAATCAAAGCTCGTGCtt                                                                >  1:106471/1‑151 (MQ=255)
                                                                                       aGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTgg                                                             <  1:45018/151‑1 (MQ=255)
                                                                                         aTCGACTGCAGCCAACTGCCAGCCTTCACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCt                                                           <  2:750925/151‑1 (MQ=255)
                                                                                          tCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCGCAAATTCCATCAATGCCGGGTCTACACCCACGCCGGCAAAAGCGCGAATGCGGGCCAGCACCGGGGAATTTAAAGCCCGTGCTTTTGCtt                                                          >  1:868336/1‑151 (MQ=255)
                                                                                                 cAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCGCCATTCCCCTCCATGCCGGATCTACACCAACGCTGGCAAAAGCGCGGAAGCGGGCCAGCACGGGGAAAATTAAACCTCGTGCTTTTGCtt                                                          >  1:28935/1‑144 (MQ=255)
                                                                                                   gCCACCTGCCAGACGACACGCCCCAGGCAACGGCGGGTCGCATTCACCGGCGCGATTCCCATCCATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCtc                                                     <  2:293048/147‑1 (MQ=255)
                                                                                                     cAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGCGGGCTACACCCGCCCCACTCCCCTCCATGACCGGTCTCCACCCAAGCTGGCCAATGCGCGGATGCGGGCCAGCCCGGGGGAATTTCACGCCCGTGGttttgggtctctc                                                     >  2:1259405/1‑135 (MQ=255)
                                                                                                     cAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGCCGCCTACACCGGCGCCATTCCCCTCAATGCCGGACCCCCACCAACCCTGGGCAATCCGCGGATTCGGGCCCGGCCGGGGCAATTTAACGCACCGGCTTTGGCTTCACtcctcat                                                >  1:1153613/1‑150 (MQ=255)
                                                                                                      aaCTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCTGCGCGATTCCCATACGTGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCtcatcatca                                              <  2:1280234/151‑1 (MQ=255)
                                                                                                              aGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGc                                       <  1:901718/151‑1 (MQ=255)
                                                                                                               gCCTACACGCTCCAGGCAACGGCGGGGTGCATCAACCGGCGCAATTCCCACCCATGGCGGATCTACACCAACGACGGCCAAATCCCGGATTCGGGGCCGCCCGGGCAAAATCAACCCTCGTGCCTTGGCGTGGCTCATCAACATTAcccc                                       >  2:1026671/1‑148 (MQ=255)
                                                                                                                     acGCTCCAGGCAACGGCGGGGCGCATACACCGGCGCCAATTCCCTCAACGCCCGATCTACACCAACGCTGGGAAAAGCGGGAATTGGGGCCAGCCACGGGAAAATCAAAACCCGTGCTTTTGCTTCTCTCCTTATCATTACCCCAGCTg                                  >  2:221909/1‑149 (MQ=255)
                                                                                                                         tCCAGGCAACCGCGGGTCGCATTCTCCTTAGCAATACCCCGCCCACCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTgcg                                <  2:313491/147‑1 (MQ=255)
                                                                                                                                           gCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCa           >  2:1236766/1‑150 (MQ=255)
                                                                                                                                           gCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCAt          <  1:946367/151‑1 (MQ=255)
                                                                                                                                               acacCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTAcc      >  1:487064/1‑151 (MQ=255)
                                                                                                                                                cacCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTAccc     <  1:423978/151‑1 (MQ=255)
                                                                                                                                                cacCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTAccc     <  1:1261552/151‑1 (MQ=255)
                                                                                                                                                  ccGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCg    <  1:1030329/150‑1 (MQ=255)
                                                                                                                                                   cGGCGCAATTCCCATCAATCCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCTTTCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCt  <  2:706433/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                       
CCATTGACATTGACCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCT  >  NC_000913/2984571‑2984867

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