Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,461,1370CT42.9% 31.1 / 12.6 35S82F (TCC→TTC) paaJ3‑oxoadipyl‑CoA/3‑oxo‑5,6‑dehydrosuberyl‑CoA thiolase
Reads supporting (aligned to +/- strand):  ref base C (9/11);  new base T (15/0);  total (24/11)
Fisher's exact test for biased strand distribution p-value = 5.50e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGTGCGAAACCCGCGTCTCGATGCGGAGTGTATCGATGATGTGATCCTCGGCTGTGCTAATCAGGCGGGAGAAGATAACCGTAACGTAGCCCGGATGGCGACTTTACTGGCGGGGCTGCCGCAGAGTGTTTCC‑GGCACAACCATTAACCG‑CTTGTGTGGTTCCGGGCTGGACGCACTGGGGTTTGCCGCACGGGCGATTAAAGCGGGCGATGGCGATTTGCTGATCGCCGGTGGCGTGGAGTCAATGTCACGGGCACCGTTTGTTATGGGCAAGGCAGCCAG  >  NC_000913/1461005‑1461287
                                                                                                                                    |                                                                                                                                                        
tgGTGCGAAACCCGCGTCTCGATGCGGCGTGTATCGATGATGTGATCCTCGGCTGTGCTAATCAGGCGGGAGAAGATAACCGTAACGTAGCCCGGATGGCGACTTTACTGGCGGGGCCGCCGCAGAATGTTTTC‑GGCACAAACCCTAAcc                                                                                                                                        >  1:78596/1‑150 (MQ=255)
           ccGCGTCTCGATGCGGAGTGTATCGATGATGTGATCCTCGGCTGTGCTAATCAGGCGGGAGAAGATAACCGTAACGTAGCCCGGATGGCGACTTTACTGGCGGGGGTGCCGCAGAGGGGTTTC‑GGCACAACCATTAACCC‑CTTGTGTgg                                                                                                                             >  1:345440/1‑149 (MQ=255)
                 ctcGATGCGGAGTGTATCGATGATGTGATCCTCGGCTGTGCTAATCAGGCGGGAGAAGATAACCGTAACGTAGCCCGGATGGCGACTTTACTGGCGGGGGCGGCGCCGAGGGGTTTC‑GGGGAAAACAATAAACG‑CTTGTGTGGGTCCcggg                                                                                                                     >  1:51382/1‑150 (MQ=255)
                   cGATGCGGAGTGTATCGATGATGTGATCCTCGGCTGTGCTAATCAGGCGGGAGAAGATAACCGTAACGTAGCCCGGATGGCGACTTTACTGGCGGGGGTGCCGCCGAGTGTTTTC‑GGGGCAACCATTAAACG‑CCTGGGTGGGTCCcgg                                                                                                                      >  2:318020/1‑148 (MQ=255)
                       gCGGAGTGTATCGATGATGTGATCCTCGGCTGTGCTAATCAGGCGGGAGAAGATAACCGTAACGTAGCCCGGATGGCGACTTTACTGGCGGGGGTGGCGCCGAGGGTTTCCGGGCCCACCCATAAACC‑‑CCTGGGTGGGTCCggg                                                                                                                      >  2:229868/1‑144 (MQ=255)
                         ggggtgtATCGATGATGTGATCCTCGGCTGTGCTAATCAGGCGGGAGAAGATAACCGTAACGTAGCCCGGATGGCGACTTTACTGGCGGGGGTGCCGCAGAGGGTTTCC‑GGCCCAAACATTAACCG‑CTTGTGTGGTTTCGGGCTgg                                                                                                                  >  1:186964/4‑146 (MQ=255)
                          gAGTGTATCGATGATGTGATCCTCGGCTGTGCTAATCAGGCGGGAGAAGATAACCGTAACGTAGCCCGGATGGCGACTTTACTGGCGGGGGCGCCGCAGAGTGGTTCC‑CGCGCCACCACTAACCG‑CCTGTGTGGTTCCGGGCTGGGCg                                                                                                               >  2:513034/1‑148 (MQ=255)
                                aTCGATGATGTGATCCTCGGCTGTGCTAATCAGGCGGGAGAAGATAACCGTAACGTAGCCCGGATGGCGACTTTACTGGCGGGGGTGCCGCAGAGGGGTTTC‑GGCCCAACCATTAACCC‑CTTGTGGGGGTCCCGGGGGGGCGCACCgggg                                                                                                       >  1:1416/1‑150 (MQ=255)
                                  cGATGATGTGATCCTCGGCTGTGCTAATCAGGCGGGAGAAGATAACCGTAACGTAGCCCGGATGGCGACTTTACTGGCGGGGCTGCCGCAGAGTGTTTCC‑GGCACAACCATTAACCG‑CTTGTGTGGTTCCGGGCTGGACGCACTGGGGtt                                                                                                     <  1:20434/150‑1 (MQ=255)
                                  cGATGATGTGATCCTCGGCTGTGCTAATCAGGCGGGAGAAGATAACCGTAACGTAGCCCGGATGGCGACTTTACTGGCGGGGCTGCCGCAGAGGGTTTTC‑CGCACAACCCTTAAACG‑CTTGTGTGGGTTCGGGCTGGGCGCCCTGGGGgtt                                                                                                    >  1:157488/1‑151 (MQ=255)
                                      gatgTGATCCTCGGCTGTGCTAATCAGGCGGGAGAAGATAACCGTAACGTAGCCCGGATGGCGACTTTACTGGCGGGGCTGCCGCAGAGTGTTTCC‑GGCACAACCATTAACCG‑CTTGTGTGGTTCCGGGCTGGACGCACTGGGGTTTGCCg                                                                                                <  2:626012/151‑1 (MQ=255)
                                                   gCTGTGCTAATCAGGCGGGAGAAGATAACCGTAACGTAGCCCGGATGGCGACTTTACTGGCGGGGGCGCCGCCGAGGGGTTTC‑GGGGCAACCATTAACCG‑CTTGGTGGGGTCCGGGGCGGGCGCCCGGGGGGTTGCCGCCCgg                                                                                           >  2:91890/1‑143 (MQ=255)
                                                        gCTAATCAGGCGGGAGAAGATAACCGTAACGTAGCCCGGATGGCGACTTTACTGGCGGGGGCGCCGCCGAGTGTTTTC‑GGCCCAAACATTAAACG‑CTTGTGTGGGTCCGGGGTGGGCGCAATGGGGGTTTCCGCCCCGGGGGTTAAAACgg                                                                              >  1:380738/1‑151 (MQ=255)
                                                         cTAATCAGGCGGGAGAAGATAACCGTAACGTAGCCCGGATGGCGACTTTACTGGCGGGGGCGCCGCCGAGTGTTTTC‑TGCCCAACCCTTAACCC‑CGTGGGGGGGTCCCGGGGGGGGGGCCTTGGGGTTGCCCCACGGGTGAGTTAAGCggg                                                                             >  2:56012/1‑151 (MQ=255)
                                                           aaTCAGGCGGGAGAAGATAACCGTAACGTAGCCCGGATGGCGACTTTACTGGCGGGGCTGCCGCAGAGTGTTTCC‑GGCACAACCATTAACCG‑CTTGTGTGGTTCCGGGCTGGACGCACTGGGGTTTGCCGCACGGGCGATTAAAGCGggc                                                                            <  1:680158/150‑1 (MQ=255)
                                                            aTCAGGCGGGAGAAGATAACCGTAACGTAG‑CCGGATGGCGACTTTACTGGCGGGGGTGCCGCAGAGTGTTTTC‑GGGACAACCATTTACCGCCTGGTGGGGTCCCGGG‑TGGGCGCACTGGGGGTTGCCGCACGGGCGGTTAAAACggg                                                                             >  1:322230/1‑147 (MQ=255)
                                                                     gagaAGATAACCGTAACGTAGCCCGGATGGCGACTTTACTGGCGGGGCTGCCGCAGAGTGTTTCC‑GGCACAACCATTAACCG‑CTTGTGTGGTTCCGGGCTGGACGCACTGGGGTTTGCCGCACGGGCGATTAAAGCGGGCGATGGCGAttt                                                                 <  2:477464/151‑1 (MQ=255)
                                                                       gaagaTAACCGTAACGTAGCCCGGATGGCGACTTTACTGGCGGGGGTGCCGCAGAGGGTTTTC‑GGCACAAACATTAAACC‑CTTGTGTGGGTTCCGGCTGGGCGCCCCGGGGGTTGCCGCACGGGCGGTTTAAACGGGGGGTGGGGAttt                                                                 >  1:556657/1‑149 (MQ=255)
                                                                             aaCCGTAACGTAGCCCGGATGGCGACTTTACTGGCGGGGGTGCCGCAGAGTGTTTCC‑GGCCCAACCATTAACCG‑CCTGTGGGGGTCCGGGCTGGGCGCACTGGGGGTTGCCGCCCGGGCGGTTAAAGCGGGGGGTGGGGCGTTGCTGAt                                                           >  2:18635/1‑149 (MQ=255)
                                                                              aCCGTAACGTAGCCCGGATGGCGACTTTACTGGCGGGGGTGCCGCAGAGTGTTTTC‑GGCACCACCATTAACCG‑CCTGTGTGGTTCCGGGCTGGACGCACTGGGGTTTGCCGCACGGGGGATTAAAACGGGGGGTGGCGATTTTCTGAt                                                           >  1:777930/1‑148 (MQ=255)
                                                                               ccGTAACGTAGCCCGGATGGCGACTTTACTGGCGGGGGTGCCGCAGAGTGTTTTC‑CGCACAACCATTAACCG‑CTTGTGTGGTTCCGGGCTGGACGCACTGGGGGTTTCCGCACGGGCGGTTAAAACGGGCGATGGCGGTTTGCTGATCGc                                                        >  1:813591/1‑150 (MQ=255)
                                                                                        aGCCCGGATGGCGACTTTACTGGCGGGGCTGCCGCAGAGTGTTTCC‑GGCACAACCATTAACCG‑CTTGTGTGGTTCCGGGCTGGACGCACTGGGGTTTGCCGCACGGGCGATTAAAGCGGGCGATGGCGATTTGCTGATCGCCGGTGGCGt                                               <  1:801520/150‑1 (MQ=255)
                                                                                        aGCCCGGATGGCGACTTTACTGGCGGGGCCGCCCCCAAGGGGTTCC‑GGCCCAACCCTTAACCG‑CTTGTGTGGTTCCGCGGCGGGCGCACTGGGGGTTGCCCCACCGGGGGTTTAAGCGGGCGGCGGCGACTTGCTGGTccccccgggcggg                                              >  2:43090/1‑142 (MQ=255)
                                                                                          cccGGATGGCGACTTTACTGGCGGGGCTGCCGCAGAGTGTTTCC‑GGCACAACCATTAACCG‑CTTGTGTGGTTCCGGGCTGGACGCACTGGGGTTTGCCGCACGGGCGATTAAAGCGGGCGATGGCGATTTGCTGATCGCCGGTGGCGTGGa                                            <  1:87407/151‑1 (MQ=255)
                                                                                           ccGGATGGCGACTTTACTGGCGGGGCTGCCGCAGAGTGTTTCC‑GGCACAACCATTAACCG‑CTTGTGTGGTTCCGGGCTGGACGCACTGGGGTTTGCCGCACGGGCGATTAAAGCGGGCGATGGCGATTTGCTGATCGCCGGTGGCGTGGAg                                           <  2:133194/151‑1 (MQ=255)
                                                                                            cGGATGGCGACTTTACTGGCGGGGCTGCCGCAGAGTGTTTTC‑GGCACAACCATTAACCG‑CTTGTGTGGGTCCGGGCTGGGCGCACTGGGGGTTGCCGCACGGGCGATTAAAACGGGCGATGGCGATTTGCTGATCGCCGGGGGCGTGGAg                                           >  1:354049/1‑150 (MQ=255)
                                                                                                tGGCGACTTTACTGGCGGGGGTGCCGCAGAGGGTTTTC‑CGGACAAACATTAAACG‑CTTGTGTGGGTTCGGGCGGGACGCACTGGGGGTTGCCGCACGGGCGGTTTAAGGGGGCGGGGGCGATTTGCTGATCGCCGGGGGCGGGGGGGCaa                                       >  2:523331/1‑150 (MQ=255)
                                                                                                               cGGGGCTGCCGCAGAGTGTTTCC‑GGCACAACCATTAACCG‑CTTGTGTGGTTCCGGGCTGGACGCACTGGGGTTTGCCGCACGGGCGATTAAAGCGGGCGATGGCGATTTGCTGATCGCCGGTGGCGTGGAGTCAATGTCACGGGCACCGtt                       <  2:699829/151‑1 (MQ=255)
                                                                                                                ggggCTGCCGCAGAGTGTTTCC‑GGCACAACCATTAACCG‑CTTGTGTGGTTCCGGGCTGGACGCACTGGGGTTTGCCGCACGGGCGATTAAAGCGGGCGATGGCGATTTGCTGATCGCCGGTGGCGTGGAGTCAATGTCACGGGCACCGttt                      <  2:408143/151‑1 (MQ=255)
                                                                                                                ggggCTGCCGCAGAGTGTTTCC‑GGCACAACCATTAACCG‑CTTGTGTGGTTCCGGGCTGGACGCACTGGGGTTTGCCGCACGGGCGATTAAAGCGGGCGATGGCGATTTGCTGATCGCCGGTGGCGTGGAGTCAATGTCACGGGCACCGttt                      >  1:657166/1‑151 (MQ=255)
                                                                                                                      gCCGCAGAGTGTTTCC‑GGCACAACCATTAACCG‑CTTGTGTGGTTCCGGGCTGGACGCACTGGGGTTTGCCGCACGGGCGATTAAAGCGGGCGTTGGCGATTTGCTGATCGCCGGTGGCGTGGAGTCAATGTCACGGGCACCGTTTGTTATg                <  1:318078/151‑1 (MQ=255)
                                                                                                                       ccGCAGAGTGTTTCC‑GGCACAACCATTAACCG‑CTTGTGTGGTTCCGGGCTGGACGCACTGGGGTTTGCCGCACGGGCGATTAAAGCGGGCGATGGCGATTTGCTGATCGCCGGTGGCGTGGAGTCAATGTCACGGGCACCGTTTGTTATgg               <  1:102319/151‑1 (MQ=255)
                                                                                                                         gCAGAGTGTTTCC‑GGCACAACCATTAACCG‑CTTGTGTGGTTCCGGGCTGGACGCACTGGGGTTTGCCGCACGGGCGATTAAAGCGGGCGATGGCGATTTGCTGATCGCCGGTGGCGTGGAGTCAATGTCACGGGCACCGTTTGTTATGGGc             >  2:88643/1‑151 (MQ=255)
                                                                                                                           agagTGTTTCC‑GGCACAACCATTAACCG‑CTTGTGTGGTTCCGGGCTGGACGCACTGGGGTTTGCCGCACGGGCGATTAAAGCGGGCGATGGCGATTTGCTGATCGCCGGTGGCGTGGAGTCAATGTCACGGGCACCGTTTGTTATGGGCa            >  1:36577/1‑150 (MQ=255)
                                                                                                                                  ttCC‑GGCACAACCATTAACCG‑CTTGTGTGGTTCCGGGCTGGACGCACTGGGGTTTGCCGCACGGGCGATTAAAGCGGGCGATGGCGATTTGCTGATCGCCGGTGGCGTGGAGTCAATGTCACGGGCACCGTTTGTTATGGGCAAGGCAGcc    >  1:35484/1‑151 (MQ=255)
                                                                                                                                    cc‑GGCACAACCATTAACCG‑CTTGTGTGGTTCCGGGCTGGACGCACTGGGGTTTGCCGCACGGGCGATTAAAGCGGGCGATGGCGATTTGCTGATCGCCGGTGGCGTGGAGTCAATGTCACGGGCACCGTTTGTTATGGGCAAGGCAGCCAg  <  2:826112/151‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                                        
TGGTGCGAAACCCGCGTCTCGATGCGGAGTGTATCGATGATGTGATCCTCGGCTGTGCTAATCAGGCGGGAGAAGATAACCGTAACGTAGCCCGGATGGCGACTTTACTGGCGGGGCTGCCGCAGAGTGTTTCC‑GGCACAACCATTAACCG‑CTTGTGTGGTTCCGGGCTGGACGCACTGGGGTTTGCCGCACGGGCGATTAAAGCGGGCGATGGCGATTTGCTGATCGCCGGTGGCGTGGAGTCAATGTCACGGGCACCGTTTGTTATGGGCAAGGCAGCCAG  >  NC_000913/1461005‑1461287

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