Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,984,7180GC40.5% 53.3 / 15.9 44P292R (CCG→CGG) yqeFshort chain acyltransferase
Reads supporting (aligned to +/- strand):  ref base G (8/17);  new base C (17/0);  total (25/19)
Fisher's exact test for biased strand distribution p-value = 4.41e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.79e-01

CCACCATTGACATTGACCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTG  >  NC_000913/2984568‑2984868
                                                                                                                                                       |                                                                                                                                                        
ccaccaTTGACATTGACCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCg                                                                                                                                                          <  2:330091/151‑1 (MQ=255)
        gacattgacCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAAt                                                                                                                                                  >  1:66595/1‑150 (MQ=255)
             tgacCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCATAAACCGGCGCaaattcc                                                                                                                                              >  1:380277/1‑144 (MQ=255)
                 ccGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTAGCCTACACCCGCGCCAATCCCCTc                                                                                                                                          >  1:316998/1‑148 (MQ=255)
                  cGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATg                                                                                                                                      <  2:39074/151‑1 (MQ=255)
                  cGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCATACACCCGCGCAATTTCCa                                                                                                                                            >  1:104467/1‑145 (MQ=255)
                      gCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCGCAAATCCCCTCCATGCCg                                                                                                                                   >  1:379420/1‑150 (MQ=255)
                         cATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCATACACACGCGCAAATcccccacatgaccggtc                                                                                                                               >  2:474049/1‑137 (MQ=255)
                                 tCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCCACCACCCGCGCCACTCCCCTCAAtatccggaaaacacc                                                                                                                         >  2:498971/1‑135 (MQ=255)
                                          ttGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCCAAAAACCGCGCAAATTCCCTCAAAGCCGGGTCTACACCAAcac                                                                                                                    >  2:425473/1‑143 (MQ=255)
                                              cAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATgc                                                                                                          <  1:140552/151‑1 (MQ=255)
                                                  cGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCCTACACCCGCGCAAATTCCCTCAATTCCGGGTCTACACCAACGCCGGCAAAAgcg                                                                                                         >  1:146605/1‑148 (MQ=255)
                                                  cGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGATGCCTACACCCGCCGCACTTCCCTCCAACCCGGATCTGCAACCACCCTCGCCAAtac                                                                                                          >  2:81975/1‑145 (MQ=255)
                                                   gAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCCCGCCATTTCCATCAAAGCCCGGTCTACCCCAACCCTCGCCAAAGcgcg                                                                                                       >  2:517903/1‑149 (MQ=255)
                                                                      aCGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCATACACCCGCGCAATTCCCATCAATGCCGGGTCTACACCAACGCTGGCAAAAGCGCGAATGCGGGCCCGCACggg                                                                                     >  1:432704/1‑148 (MQ=255)
                                                                       cGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTaaa                                                                                 <  1:480050/151‑1 (MQ=255)
                                                                           tCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATTCACCGTCACAATTCGGCTCAAAGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTc                                                                              <  2:39752/150‑1 (MQ=255)
                                                                            cATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCCTCAATGCCGGATCTACACCAACGCCGGCAAATGCGCGAATGCGGGGCCGCACGGGGAAATTCa                                                                             >  1:556984/1‑150 (MQ=255)
                                                                            cATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCa                                                                             <  1:256778/150‑1 (MQ=255)
                                                                            cATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCCCCAACACCCTCGCAAATCCGCTCCATGTCCGCGCAACACCAACACCGGCGACGGCGCCGACACGGGGCAGCACCGGGaa                                                                                  >  2:56671/1‑145 (MQ=255)
                                                                                      gATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCt                                                                  <  1:469579/151‑1 (MQ=255)
                                                                                       aTAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCtt                                                                 <  1:246651/151‑1 (MQ=255)
                                                                                        tAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCCTACCCCCGCGCCATTCCCCTCAATGCCGGGTCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGGAAATTCAACCCCTGTGCttt                                                                >  1:192059/1‑151 (MQ=255)
                                                                                         aaGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATTTTTTTTAATGATACGGCGACCACCGAGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTg                                                               <  2:303239/151‑1 (MQ=255)
                                                                                              cGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCtt                                                           <  2:366404/150‑1 (MQ=255)
                                                                                              cGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCGCAATTCCCATCAATGCCGGATCTACACCAACGCCGGCAAAAGCGCGGATGCGGGCCAGCACGGGGAAATTCAACGCTCGTGCTTTGGCtt                                                           >  1:275588/1‑150 (MQ=255)
                                                                                                    cAGCCAACTGCCAGCCTACACGCTCCAGGCGACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCtcat                                                    <  1:511968/151‑1 (MQ=255)
                                                                                                    cAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCGCAAATCCCCTCAATGCCCGATCTACACCAAAGCCGGCCAAAGCGCGAAAGCGGGCCAGCACGGGGAAATTCAAAGCCCGTGCCTTGGCCTccc                                                        >  2:276213/1‑145 (MQ=255)
                                                                                                      gCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGATTTGGCTTCGCtcatc                                                   <  2:379300/150‑1 (MQ=255)
                                                                                                       ccAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCACACACCCGCGCCAATTCCCTCAATGCCGGATCTACACCCACGCTGGCAAAAGCGCGAATGCGGGCCAGCACGGGGAAATTCAACCCCCGTGCCTTTGCTTCGCtcatca                                                  >  2:388035/1‑150 (MQ=255)
                                                                                                       ccAACTGCCAGCCTAAACGCTCCAGGCAACGGCGGGTCGCATACACCCCCGCCAATCCCCTCCATGCCGGGACTACACCAACGCCGGCCAAAGCGCGAAAGCGGGGCCGCACGGGGAAATTCAAACCCCGTGCTTTGGCTTCGCtc                                                      >  2:234015/1‑146 (MQ=255)
                                                                                                         aaCTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGGTCTACACCCACGCTGGCAAATGCGCGAATGCGGGCCCGCACGGGGAAATTCAAAGCCCGTGCTTTGGCTTCGCccatcatc                                                >  1:234696/1‑150 (MQ=255)
                                                                                                           cTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCGCAATTCCCCTCCAAGCCGGATCCACCCCCACACTGGCAAAAGGGCGAATGCGGGGCCCCCCGGGGAAAATCAACCCCCGTTCTTTTGCTTCGCccctcaac                                                >  2:336031/1‑146 (MQ=255)
                                                                                                                  ccTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCCATTCCCCTCCATGCCCGGTCTACACCCAACCCGGCCAAAACGCGAATGCGGGGCAGCACGGGGAAAATCAAAGCCCGGGCTTTTGCTTCGCCCCTCATCCtt                                             >  2:427238/1‑144 (MQ=255)
                                                                                                                  ccTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGAACAATTCCCCGCAATGCCAGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAg                                      <  2:580516/151‑1 (MQ=255)
                                                                                                                       acGCTCCAGGCAACGGCGGGTCGCATACACCCGCGCAATTTCCCCCAATGCCCGATCTACACCAACGCTGGCAAAAGCGCGAAAGCGGGCCAGCACGGGGAAATTAAACGCTCGTGCTTTTGCTTCGCTCACCATCATTACCGCAGCTgc                                  >  2:531318/1‑150 (MQ=255)
                                                                                                                        cGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTgcgc                                <  1:20429/151‑1 (MQ=255)
                                                                                                                            ccAGGCAGCGGCGGGTCGCATTCTCCTGAGCGATTCCCCGCACTGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATc                            <  2:117497/151‑1 (MQ=255)
                                                                                                                                  aaCGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTa                       <  1:468230/150‑1 (MQ=255)
                                                                                                                                            cGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCa            <  1:425612/151‑1 (MQ=255)
                                                                                                                                            cGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCa            <  1:562500/151‑1 (MQ=255)
                                                                                                                                            cGCATACACCGGCACGATTCCCAGCAATGCGGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCa            <  2:562268/151‑1 (MQ=255)
                                                                                                                                              cATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCAtt          >  2:172600/1‑151 (MQ=255)
                                                                                                                                              cATACACAGGCGCAATTCCAATCAATGACGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCAtt          >  2:202607/1‑151 (MQ=255)
                                                                                                                                               aTACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTa         <  2:12225/151‑1 (MQ=255)
                                                                                                                                                       ggCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGACGCATTACCCGCTg  <  1:364453/151‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                        
CCACCATTGACATTGACCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTG  >  NC_000913/2984568‑2984868

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: