Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,984,4410AC51.2% 48.7 / 73.6 41G384G (GGT→GGG) yqeFshort chain acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (20/0);  new base C (0/21);  total (20/21)
Fisher's exact test for biased strand distribution p-value = 3.72e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.27e-09

CGAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACGC  >  NC_000913/2984291‑2984591
                                                                                                                                                      |                                                                                                                                                      
cGAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCcaca                                                                                                                                                        <  2:230666/151‑1 (MQ=255)
     atGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTg                                                                                                                                                   >  1:37267/1‑151 (MQ=255)
                     aaaaTGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACa                                                                                                                                    >  1:831172/1‑150 (MQ=255)
                       aaTGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACaa                                                                                                                                   >  2:884416/1‑149 (MQ=255)
                                 cccAAAAAATTTATGTCGATTTGCCAGCCATCATTATTCAACAAATTAATGACGGAAGATATTTATCGGTAGGGAGAAAATGCAGAGTTTTTCGTCACGTTAAAGGGTAAATGCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGt                                                                                                                       <  2:831134/151‑1 (MQ=255)
                                   cAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTcc                                                                                                                     >  2:602892/1‑151 (MQ=255)
                                    aaaaaaTGTATGTCGATTTGCCAGCCATCTTTATCCAACAGATAAATGACGGCAGATATTTATGGGTAAGGAGAAAATGCAGAGCTTTTCGTCACGTTAAATGGTCAATCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTcc                                                                                                                     <  2:160345/148‑1 (MQ=255)
                                        aaTGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGc                                                                                                                >  2:71483/1‑151 (MQ=255)
                                            tatgtCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATGGGTAAGGAGAAAATGCAGAGCTTTTCGTCACTTTAAATGGTAAATGCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCgg                                                                                                              <  2:267072/149‑1 (MQ=255)
                                                     ttGCCAGCCATAATTATTAAACAAATTAATGACGGCAGATATTTACCGGTAGGGAGAAAAGGCAGAGTTTTTCGTCACTTTAAATGGTAAATCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGc                                                                                                            <  2:44906/142‑1 (MQ=255)
                                                     ttGCCAGCCATAATTATTAAAAAAATAAAGGCGGGAAAATATTTATGGGTAAGGAGAAAAGGAGGAGTTTTTCGCCACTTTCAATGGTCAATCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACg                                                                                                      <  1:326537/148‑1 (MQ=255)
                                                         cagccaTCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCAcc                                                                                               >  2:863306/1‑151 (MQ=255)
                                                         cagccaTCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCAcc                                                                                               >  2:69890/1‑151 (MQ=255)
                                                          agctaTCATTATAAAAAAAATAAATGAGGGCATATTTTTACCGTAAAGGAGGAAATGCAGACCTTTTCGCCACGTTCATTGGTCAAGCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCAcc                                                                                               <  1:227346/150‑1 (MQ=255)
                                                             cattatTATACAACAGATTAATGACGGCAGATTTTTAGCGGAAAGGAGAAAAGGAAGAGCTTTTCGCCCCGTTCAAGGGTCATTGCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCa                                                                                                 <  1:362434/145‑1 (MQ=255)
                                                              atcatTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGCCCCCCGCCGATACAAAGCGTCGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTc                                                                                          >  1:15823/1‑151 (MQ=255)
                                                              atcatTATTAAACAGATAAAGGACGGCAGATATTTAGCGGAAAGGAGAAAATGCAGAGTTTTTCGTCACGTTAAAGGGTCAATGCCACCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTc                                                                                          <  2:82653/151‑1 (MQ=255)
                                                                catTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCAt                                                                                        >  1:416304/1‑151 (MQ=255)
                                                                  ttattCAACAGATTAATGACGGCAGATATTTAGCGGTAAGGAGGAAATGAAGAGTTTTTCGTCACTTTAAATGGTCAATGCCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATg                                                                                       <  2:106778/150‑1 (MQ=255)
                                                                  ttatcaaaCAAATAAAGGCCGGCAAATTTTTAGCGGAAGGGAGGAAATGCAGATCTTTTCGCCACGTTCAATGGTCAATCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTc                                                                                          <  2:597948/141‑1 (MQ=255)
                                                                    ataaaaCAGATAAATGACGGAAGATATTTATCGGTAGGGAGAAAATGCAGTGCTTTTCGTCTCTTACAATGGTCAATGCCACCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGaa                                                                                     <  2:267905/146‑1 (MQ=255)
                                                                        aaCAGATTAATGACGGAAAATATTTACCGGAAAGGAGAAAAGGAAGAGTTTTTCGCCCCGTTAAATGGACAATGCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCaga                                                                                <  1:809098/151‑1 (MQ=255)
                                                                         aCAAATTAAGGGCGGCAGATTTTTATGGGTAAGGAGAAAAGGCAGACCTTTTCGTCACGTTAAATGGTCAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCagag                                                                               <  1:824407/151‑1 (MQ=255)
                                                                                aagggcGGCAGATATTTCGCGGAAAGGAGAAAGCGCAGCGCTTTTCTTCACGTGAAGGGGTAAGGCCCACCCCGCTGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAg                                                                        <  1:740882/146‑1 (MQ=255)
                                                                                      ggCAGATTTTTATCGGTAGGGAGAAAAGGCAGAGCTTTTCGTCCCGTTAAAGGGTAAATCCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGAt                                                                     <  2:77202/148‑1 (MQ=255)
                                                                                       gCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCgg                                                                 >  2:660070/1‑151 (MQ=255)
                                                                                         aaatttttATCGGAAAGGAGAAAAGGCAGAGCTATTCGCCCCGTTAAATGGCCAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCg                                                                  <  2:487539/143‑1 (MQ=255)
                                                                                         aGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGc                                                                >  2:199763/1‑150 (MQ=255)
                                                                                           atatTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAAc                                                             >  1:805818/1‑151 (MQ=255)
                                                                                           atatTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAAc                                                             >  2:218426/1‑151 (MQ=255)
                                                                                              tttATCGGAAAGGAGAAAATGAAGAGTTTTTGGTCCCTTTCAATGGACAATGCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAAcc                                                            <  2:58887/149‑1 (MQ=255)
                                                                                                    ggTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAgcgc                                                    >  1:410913/1‑151 (MQ=255)
                                                                                                       aaGGAGAAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCt                                                  >  2:16922/1‑150 (MQ=255)
                                                                                                        aggaggAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTAt                                                >  1:290252/1‑151 (MQ=255)
                                                                                                              aaaTGAAGATCTTTTCGTCACGTTAAATGGTCAATCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGt                                           <  2:341144/150‑1 (MQ=255)
                                                                                                                aTGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGAc                                        >  1:539383/1‑151 (MQ=255)
                                                                                                                     gagtttttCGTCCCTTTCATTGTTCATTCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCg                                      <  1:828295/142‑1 (MQ=255)
                                                                                                                           ttCGTCACGTTAAATGGTCAACCCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGc                                 <  1:87474/147‑1 (MQ=255)
                                                                                                                            tccccccGTTCAAGGGTCAACCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATc                             <  1:100951/144‑1 (MQ=255)
                                                                                                                               tCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCgcgc                         >  2:588943/1‑151 (MQ=255)
                                                                                                                                       aaTGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattg                 >  2:530319/1‑151 (MQ=255)
                                                                                                                                         aggTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgac               >  1:219203/2‑151 (MQ=255)
                                                                                                                                                    acacCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGAc    >  2:640312/1‑151 (MQ=255)
                                                                                                                                                      acCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACGc  >  2:239869/1‑151 (MQ=255)
                                                                                                                                                      acCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACGc  >  1:796035/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                      
CGAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACGC  >  NC_000913/2984291‑2984591

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: ‑