Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,813,2720CT45.7% 24.9 / 16.3 35S181L (TCG→TTG) dfpfused 4'‑phosphopantothenoylcysteine decarboxylase/phosphopantothenoylcysteine synthetase, FMN‑binding
Reads supporting (aligned to +/- strand):  ref base C (4/15);  new base T (16/0);  total (20/15)
Fisher's exact test for biased strand distribution p-value = 1.90e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCTGCCACGCAGCATAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATCGCCGCCGCCGCTGCCCGTCGTGGCGCGA  >  NC_000913/3813123‑3813418
                                                                                                                                                     |                                                                                                                                                  
cGCTGCCACGCAGCATAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCg                                                                                                                                                   <  2:94978/151‑1 (MQ=255)
   tGCCACGCAGCATAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTTGccc                                                                                                                                                >  1:181399/1‑151 (MQ=255)
   tGCCACGCAGCATAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGccc                                                                                                                                                >  1:294556/1‑151 (MQ=255)
          cagcaTAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGGGCATTTTTTGCCCCTc                                                                                                                                             >  1:345049/1‑147 (MQ=255)
              aTAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGGGCCTTTTTTGCCCGTCAACGaac                                                                                                                                      >  2:89038/1‑148 (MQ=255)
                aaTTGAGAGGTCCTTGTTTCCCGTGGTGTGCTCAACTTGGAGCCAGACAGTTGCAGTCAGGCTTGTGGTGATAGCGGTCCTAGGCTAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGGGCATTTTTCTCCCGTCAACGACCTGa                                                                                                                                   <  2:314241/151‑1 (MQ=255)
                     agagGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACa                                                                                                                               <  2:249857/150‑1 (MQ=255)
                      gagGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATc                                                                                                                             <  1:171676/151‑1 (MQ=255)
                               gcttCCCGGGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGTGAATGCTCGATCCGTTAACCATTTTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGaa                                                                                                                         <  1:686209/146‑1 (MQ=255)
                                  tCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGGGCCTTTTTCGCCCCCGAAACACCTTAAACATCttcaaacttagt                                                                                                                 >  2:467999/1‑140 (MQ=255)
                                  tCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATg                                                                                                                  <  1:303457/150‑1 (MQ=255)
                                   cccGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCCTTTTTTGCCCGTCAACGACCTTAAACATCTGAACATtattat                                                                                                                >  1:672787/1‑151 (MQ=255)
                                        ggTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGGAAGGCCTTTTTTTCCCGTCAACGCCCCGAAAAAACTGAACAATTTGAtt                                                                                                               >  2:182560/1‑147 (MQ=255)
                                                           cAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTTGCCCGTCAAACACCTGAAACATCTGAACATTATGATTACCGCCGGCCCCACGCGGGaaac                                                                                        >  1:678330/1‑149 (MQ=255)
                                                           cAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTTGCCCCTCAAAGACCCGAAACATCTGAACATTATGATTACCGCCCGCCCGAcccggggaac                                                                                        >  1:355/1‑145 (MQ=255)
                                                               cAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGAAGCGCATTTTTTGCCCGTCAAAGACCTCAAAAATCTGGACATTATTATTTCCCCCCGGCCGCCGCGTGAAcccccc                                                                                    >  2:65791/1‑147 (MQ=255)
                                                                 gtgtCAGCCAGGTTTGGGGTGATATCGGTCCTGGGCGAATGCTCGTTCCGTGAACCATTGTGGTTATGGCGGTAGCGAATTTTTCGCCCGTCAACGACCGGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCg                                                                                   <  2:221020/150‑1 (MQ=255)
                                                                      aGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCCTTTTTTCCCCCCCAACGACCCGAAAAATCTGAACAATTTTATTACCTCCCGCCCCCCGCGTGAACCCCTcccacc                                                                               >  1:523757/1‑144 (MQ=255)
                                                                             cTTGTGGTAATATCGGCCCTGGGCAAATCCACGATCCGTTAACCATTGTGGATGTGGCGGTAGCGCATTTTTCGCCGGTCAACGACCTGAAACATCTGACCATGATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTg                                                                           <  2:257497/146‑1 (MQ=255)
                                                                                 tggtgATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATAGGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTatat                                                                   <  2:238376/150‑1 (MQ=255)
                                                                                      atatCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTTGCCCGCCAACGACCTGAAACATCTTAACATTATGATTACCGCCGGCCCGACGCGTGAAACCCCCCATCCGGGGCGTTTTAtct                                                                 >  2:341191/1‑147 (MQ=255)
                                                                                      agatCGGTCCTGGGGGAATGCTCGATCCGTTAACCATTGTGTATGTGGCGGTAGCGCATTTTGCGCCCGTCAACGACCTGAATCATCTTAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTa                                                              <  2:208711/148‑1 (MQ=255)
                                                                                       tatCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTTGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCCGCCCGACGCGTGAACCCCTCGATCCGGTGCGTTATAtctc                                                                >  1:406287/1‑147 (MQ=255)
                                                                                        atCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTTGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGGGAAACCCCCGATCCGGGGCGTTTTATTTCTaa                                                             >  1:436365/1‑149 (MQ=255)
                                                                                          cGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATca                                                          <  1:84312/150‑1 (MQ=255)
                                                                                            gTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCCTTTTTTCCCCCTCCACGACCTTAAACATCTGAAAATTATGATTTCCGCCGGCCCGACGCGTGGACCGCTCGATCCGGGGCGGTATATCTCTAATc                                                           >  2:71724/1‑147 (MQ=255)
                                                                                              ccTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTTGCCCGTCAAAGACCTGAAACATCTGAACATTTTGATTACCGCCGGCCCGACGCGGGAAACGCTCGCTCCCGGGCGGTATATCTCTAATCACAGcc                                                     >  2:552655/1‑150 (MQ=255)
                                                                                                tGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAAGGCGTTTTTTGCCCGCCCACCACCTGAAACATCTGAAAATTTTGTATACCGCCGGCCCGCCCCCCGAACCCCCCGACCCCGGGGGGTTTATCTCTAAACACAGcccc                                                   >  2:623037/1‑151 (MQ=255)
                                                                                                      aaTGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTTGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCgg                                                 >  1:446174/1‑147 (MQ=255)
                                                                                                         gCTCGATCCGTTAACCATGGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGAt                                           <  1:581114/150‑1 (MQ=255)
                                                                                                                     aaCCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATc                              <  1:498802/151‑1 (MQ=255)
                                                                                                                              ggATATGGCGGTAGCGCATTTTTCGCCCGTCAACGTCCTGAAACTTCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATcgccgccgcc                     <  2:80673/151‑1 (MQ=255)
                                                                                                                                        gTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATCGCCGCCGCCGCTGCCcgtc           >  1:559045/1‑151 (MQ=255)
                                                                                                                                               aTTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATCGCCGCCGCCGCTGCCCGTCGTGgcgc    >  1:434145/1‑151 (MQ=255)
                                                                                                                                                  tttCGCCCGTCAACGCCCTGAAACATCTGACCATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATCGCCGCCGCCGCTGCCCGTCGTGGCGCGa  <  2:249915/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CGCTGCCACGCAGCATAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATCGCCGCCGCCGCTGCCCGTCGTGGCGCGA  >  NC_000913/3813123‑3813418

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