Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,813,2720CT45.8% 21.2 / 27.6 48S181L (TCG→TTG) dfpfused 4'‑phosphopantothenoylcysteine decarboxylase/phosphopantothenoylcysteine synthetase, FMN‑binding
Reads supporting (aligned to +/- strand):  ref base C (5/21);  new base T (22/0);  total (27/21)
Fisher's exact test for biased strand distribution p-value = 3.22e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGCCACGCAGCATAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAA‑CCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATCGCCGCCGCCGCTGCCCGTCGTGG  >  NC_000913/3813125‑3813413
                                                                                                                                                   |                                                                                                                                              
cTGCCACGCAGCATAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGcc                                                                                                                                             <  1:913936/151‑1 (MQ=255)
  gCCACGCAGCATAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGGAGCGCGATTTTCGCCCg                                                                                                                                           >  2:266540/1‑151 (MQ=255)
          gcaTAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGGAGGGGCTTTTTTGccccgccaaga                                                                                                                                    >  1:61858/1‑142 (MQ=255)
               aTTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAGCCATTGTGGATATGGCGGTTGCGCCTTTTTTGCCCGTCCACGACCTGaa                                                                                                                              >  2:894019/1‑151 (MQ=255)
               aTTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGGGCCTTTTTTGCCCCCCAACGCCCTGaa                                                                                                                              >  1:972616/1‑151 (MQ=255)
               aTTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGGGCATTTTTTGCCCCCCCACGACCCGaa                                                                                                                              >  1:978948/1‑151 (MQ=255)
               aTTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGGACTTTTCGCCCCTCCACAACCCGaa                                                                                                                              >  2:819402/1‑151 (MQ=255)
               aTTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCCTTTTTTTCCCGTCCACGACCTGaa                                                                                                                              >  1:1154651/1‑151 (MQ=255)
               aTTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCCTTTTTTTCCCCGCAACGAACttga                                                                                                                              >  2:599705/1‑148 (MQ=255)
                 ttAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAAc                                                                                                                            <  1:855999/151‑1 (MQ=255)
                 ttAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAAc                                                                                                                            <  1:819727/151‑1 (MQ=255)
                          cttgcttCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAAc                                                                                                                    <  2:479180/150‑1 (MQ=255)
                          cttgcttCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAAc                                                                                                                    <  2:939074/150‑1 (MQ=255)
                            tgcttCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGGAGCGCAATTTTTGCCCCTCCAACACCCGTAAAAACTTAACAtt                                                                                                                 >  1:404064/1‑151 (MQ=255)
                              cttCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTAt                                                                                                               <  2:688565/151‑1 (MQ=255)
                              cttCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTAt                                                                                                               <  2:477470/151‑1 (MQ=255)
                                      ggTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTTGGGCCGTTTTTGCCCCTCAACGACCTGAAACATCTGCACAATATGGttttc                                                                                                        >  2:265367/1‑147 (MQ=255)
                                                tCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTTGCGCCTTTTTTGCCCCGCAACGACCCGGAAAAACTTAAAAATTTTATTTACGCCGGcccca                                                                                              >  1:1251343/1‑148 (MQ=255)
                                                tCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGGAGCGCATTTTTTGCCCGGCAAACACCTGAAACATCTTAACATTTTGATTTCCCCCCGcccca                                                                                              >  1:1174972/1‑148 (MQ=255)
                                                tCTGGGGGCCAGACAGTGGCAGTCAGGCTTGGGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTGGCGCCTTTTTTGCCCCTCCACCACCCGAAACAACTGAAAATTATGATTAACCCCCGcccca                                                                                              >  1:1249240/1‑148 (MQ=255)
                                                    ggggCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGAcgcg                                                                                          <  1:952073/150‑1 (MQ=255)
                                                     gggCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGGCTTTTTTGCCCGGCAACGACCCGAAACAACTGAACATTTTGATTTCCGCCGGCCCGCcgcggg                                                                                        >  1:125415/1‑149 (MQ=255)
                                                      ggCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAAGGCCTTTTTTGCCCGTCCAAGAACTGAAACATCTTAACAATATGATTAACGCCGGCCCGACGCGTGa                                                                                       >  2:846448/1‑151 (MQ=255)
                                                      ggCCAGACAGTGGCAGTCAGGCTTGGGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGTGCCGTTTTTCCCCCTCCACGACCCGGAACCACTTAACATTTTGATTTACCCCGGccccccct                                                                                           >  2:654049/1‑142 (MQ=255)
                                                        ccAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGGGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAA‑c                                                                                    <  1:132242/151‑1 (MQ=255)
                                                        ccAGACAGTGGCAGTCAGGCTTGTGGGGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAA‑c                                                                                    <  1:1139886/151‑1 (MQ=255)
                                                            aCAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTTGCCCCTCAACGACCCGAAACATCTGGACATTATGATTTCCGCCGGCCCCACCCGTGaaaccc                                                                                  >  1:195130/1‑146 (MQ=255)
                                                             cAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCCTTTTTTGCCCCTCAACGACCCGAAAAATCTGAACATTATGATTACCCCCCGCCCCCCCCGTGaaacccc                                                                                 >  1:394150/1‑145 (MQ=255)
                                                             cAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGGAGCGCATTTTTCGCCCCCTAACGACCTGAAAAATTTGAACATTATTCTTACCGCCGGCCCCCCGCGCGGA‑ccccc                                                                                >  1:468165/1‑147 (MQ=255)
                                                                         ggCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGGGCATTTTTTGCCCGTCAACGACCTGAAACATCTGAACATTATGATTTCCCCCGGCCCCACCCGGGAA‑ACCCCCGATCCGGTGCgg                                                                   >  1:514488/1‑150 (MQ=255)
                                                                          gtttGTGGTGATATCGGTCCGGGGCGAATGCTCGATCGGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGCCCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAA‑CCGCTCGATCCGGTGc                                                                     <  2:610314/146‑1 (MQ=255)
                                                                             tgtgGTGATATCGGTCCTGGGCGAATGCGCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAA‑CCGCTCGATCCGGTGCGTtata                                                               <  2:602470/151‑1 (MQ=255)
                                                                               tggtgATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCACCGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAA‑CCGCTCGATCCGGTGCGTTatat                                                              <  1:633465/150‑1 (MQ=255)
                                                                                ggtgATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGGGCGTTTTTTGCCCGTCCAAGACCCGAAACAACTTAACATTATGAATAACCCCGGCCCGACGCGGGAACCCGCCCGA‑CCGGGGCGGTATAtc                                                             >  2:781911/1‑150 (MQ=255)
                                                                                 gtgATATCGGTCCTGGGGGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAA‑CCGCTCGATCCGGTGCGTTATAtct                                                            <  1:80631/150‑1 (MQ=255)
                                                                                 gtgATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAA‑CCGCTCGATCCGGTGCGTTATAtct                                                            <  1:169473/150‑1 (MQ=255)
                                                                                  tgATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAA‑CCGCTCGATCCGGTGCGTTATAtctc                                                           <  1:269434/150‑1 (MQ=255)
                                                                                                  cGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAA‑CCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCa                                          <  2:61844/151‑1 (MQ=255)
                                                                                                            aTCCGTTAACCATTGTGGATGTGGCGTTAGCGCATTTTCCGCCCGTCAACGACCTAAACCATGTGAACATTAGGATTACCGCCGGCCCGACGCGTGAA‑CCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGttt                                <  2:601277/151‑1 (MQ=255)
                                                                                                              ccGTTAACCATTGTGGATATGGCGGTAGCGCCTTTTTTGCCCCTCAACGACCTGAAAAATCTGAACATTATTATTACCGCCGGCCCGAACCCGGAA‑CCGCCCGCTCCGGCGCGTTATATTCCTAATCACAACCCCTGCAAGATGGGGTTTg                              >  2:784801/1‑151 (MQ=255)
                                                                                                                gTTAACCATTGTGGATATGGCGGTAGCGCATTTTTTGCCCGTCAAAGACCTTAAAAATCTGAACATTATGATTACCCCCGGCCCCCCGCGTGAA‑CCGCCCGATCCGGGGCGTTATTTTTCTTATCCCAGCCCCGGCAAGATGGGGTTTg                              >  1:966727/1‑149 (MQ=255)
                                                                                                                 ttAGCCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAA‑CCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCt                            <  1:894356/150‑1 (MQ=255)
                                                                                                                 ttAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCACCGGCCTGAAACATCTGAACTTTATTTTTACCGCCGGCCCGACGCGTGAA‑CCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCt                            <  2:707483/150‑1 (MQ=255)
                                                                                                                     ccATTGTGGATATGGCGGTAGGGCCTTTTTTTCCCGTCAACAACCTGAAACAACTGAACATTATGATTATCGCCCGCCCGCCGCGTGAA‑ACGCTCTCTCCCGGGCGTTTTTTCTCTTATCCCAGACCCCGCAAGATGGGGTTTTCTATcg                        >  2:777330/1‑150 (MQ=255)
                                                                                                                       aTTGTGGATATGGCGGTAGCGCATTTTTTGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGGGAA‑CCCCCCCATCCGGGGCGTTTTATCTCTAATCCCAGCCCCCGCAAGATGGGGTTTTCTATcgcc                      >  1:449070/1‑150 (MQ=255)
                                                                                                                            ggATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAAAATCTTGAAATTATGATTTCCGCCCGGCCGACGCGTGAA‑CCGCCCGATCCCGGGCGTTTTATCTCTAATCACAGCCCCCGGAAGGTGGGGTTTTCTATcgccgcccc                 >  2:147428/1‑148 (MQ=255)
                                                                                                                                atGGGGGTCGCGCTTTTTTCGCCCGCCAACGACCGGAACCCTCTAAACATTATGATTCCCGCCGGCCCGACGCGTGAA‑CCGCTCGATCCGGTTCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATCGCCGCCGCCGCTg            <  1:1144536/151‑1 (MQ=255)
                                                                                                                                           gcATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAA‑CCGGTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATCGCCGCCGCCGCTGCCCGTCGTgg  >  2:895755/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
CTGCCACGCAGCATAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAA‑CCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATCGCCGCCGCCGCTGCCCGTCGTGG  >  NC_000913/3813125‑3813413

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