Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,813,2720CT37.5% 31.0 / 10.5 32S181L (TCG→TTG) dfpfused 4'‑phosphopantothenoylcysteine decarboxylase and phosphopantothenoylcysteine synthetase
Reads supporting (aligned to +/- strand):  ref base C (3/17);  new base T (12/0);  total (15/17)
Fisher's exact test for biased strand distribution p-value = 2.02e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GAATGGCGAATGACCTGGTATCGACGATTTGTCTGGCTACACCTGCGCCTGTAGCCGTGCTCCCCGCCATGAACCAGCAGATGTACCGTGCCGCTGCCACGCAGCATAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATCGCCGCCGCCGCTGCCCGTCGTGGCGCGAACGTCACGCTGGTATCAGGTCCGGTTTCACTACCGACGCCACCGTTTGTTAAACGTGTTGATGTGATGACCGCGCTGGAAATGGAAGCCGCCGTGAATGCT  >  NC_000913/3813032‑3813519
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                         
gAATGGCAAATGCCCTGGTATCGACGTTTTGTCTGGCTACCCCTGCGCCTGTAGCCGTGCCCCCCGCCATGAACCAGCAGATGTACCGTGCCGCTGCCACGCAGCATAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAAc                                                                                                                                                                                                                                                <  1:271196/251‑1 (MQ=255)
                  tATCGACGATTTGTCTGGCTACACCTGCGCCTGTAGCCGTGCTCCCCGCCATGAACCAGCAGATGTACCGTGCCGCTGCCACGCAGCATAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCCTTTTTTTCCCGTCAACGACCTTAAAAATCTGaa                                                                                                                                                                                                                              >  1:311195/1‑250 (MQ=255)
                       acgaTTTGTCTGGCTACACCTGCGCCTGTAGCCGTGCTCCCCGCCATGAACCAGCAGATGTACCGTGCCGCTGCCACGCAGCATAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGGCCCTTTTTTTCCCCGCAACAACCTGAAAAATATGAAAATTAt                                                                                                                                                                                                                        >  1:278250/1‑251 (MQ=255)
                                         ccTGCGCCTGTAGCCGTGCTCCCCGCCATGAACCAGCAGATGTACCGTGCCGCTGCCACGCAGCATAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGGAGGGGCTTTTTTGCCCGTCAAAGACCTGAAACATCTGAACATTATGATTACCGCCGGGCCgcc                                                                                                                                                                                                      >  1:185586/1‑249 (MQ=255)
                                               ccTGTAGCCGTGCTCCCCGCCATGAACCAGCAGATGTACCGTGCCGCTGCCACGCAGCATAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTTTCCCGTCAACGACCTGAAACATCTGAACAATATGATTACCGCCGGCCCGACCCGTg                                                                                                                                                                                                 >  1:221905/1‑250 (MQ=255)
                                                  gTAGCCGTGCTCCCCGCCATGAACCAGCAGATGTAACGTGCCGCTGCCACGCAGCATAATTTAGAGGTGATTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGGAGAGGATTTTTTTCCCGTCAACCAACTGAAAAAACTAAACCATATAATTAccgcccgacccacgccggaacc                                                                                                                                                                                             >  2:271014/1‑234 (MQ=255)
                                                         ggccccccGCCATGAACCAGCAGATGTACCGTGCCGCTGCCACGCAGCATAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGAt                                                                                                                                                                                      <  1:222720/247‑1 (MQ=255)
                                                                 gCCATGAACCAGCAGATGTACCGTGCCGCTGCCACGCAGCATAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCCTTTTTTGCCCCGCAAAGACCTGAAACATCTGAACAATATTATTACCCCCGGCCCCACGCGTGAACCGCTCCATCCGGTGCg                                                                                                                                                                              >  1:85074/1‑251 (MQ=255)
                                                                  ccATGAACCAGCAGATGTACCGTGCCGCTGCCACGCAGCATAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGt                                                                                                                                                                             <  1:30168/251‑1 (MQ=255)
                                                                   cATGAACCAGCAGATGTACCGTGCCGCTGCCACGCAGCATAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCCTTTTTTGCCCCTCAACGACCCGAAACATCTGAACATTTTGATTACCGCCGGCCCGACCCGTGAACCCCTCGAACCCGTGCGtt                                                                                                                                                                            >  1:373633/1‑251 (MQ=255)
                                                                       aaCCAGCAGATGTACCGTGCCGTTTCCACGCAGCAAAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCTTCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACTTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTatat                                                                                                                                                                        <  1:151456/251‑1 (MQ=255)
                                                                         ccAGCAGATGTACCGTGCCGCTGCCACGCAGCATAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTTCCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATAtct                                                                                                                                                                      <  1:396530/251‑1 (MQ=255)
                                                                                    ccgggtgcgctgacccGCAGCAGAATGTAGAGTGGCTTGCTTCCGGTGTTTTGCTCATCGGGGGGCCAGACAGGGGAAGTGAGGCTGGTGGTGATATCGGTCCTGGGCGATGGCTCGATCCGTTAACAATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGGAACCGACCTGAAACACTTGAACTTGAGGATTCCCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATcacat                                                                                                                                                            <  2:383058/235‑2 (MQ=255)
                                                                                                              ttAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTTGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCCCGGACCGGTGCGGTATATCTCTTATCACAGCTCCGGCAAGATGGGGTTTGCTATcg                                                                                                                                  >  1:348869/1‑250 (MQ=255)
                                                                                                              ttAGAGGGGCTGGCTTCCCGTGTTTTGCTCATCTGGGGGCAAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATcg                                                                                                                                  <  2:348869/250‑1 (MQ=255)
                                                                                                                  agGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCCTTTTTTGCCCCTCAACGAACTGAAACATCTGAACATTATGATTACCGCCCGCCCGAAGCGGGAAACGCGCGATCCGGGGCGTTATATCTCTAATCCCAGCTCCGCCCAGATGGGTTTTTCTTTtgccgcc                                                                                                                             >  1:84612/1‑251 (MQ=255)
                                                                                                                      gcttgcttCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCCGAAACATCTGAACATTATGATTACCGCCCGCCCCACGCGTGAACCGCCCCATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGGTTTGCTTTctccgccgccc                                                                                                                         >  1:91432/1‑250 (MQ=255)
                                                                                                                                            aTCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCCTTTTTTGCCCGGCAACGACCTGAAACATCTGAACATTTTGATTACCGCCGGCCCGACGCGGGAACCGCTCGATCCGGGGCGTTTTATCTCTAATCACAGCTCCCGCCAGAAGGGGTTTGCTATTGCCCCCGCCCCTGCCCCTCGTGGGGCGAACGt                                                                                                   >  1:48813/1‑251 (MQ=255)
                                                                                                                                                   ggcaaGACGGTGGCAGTCAGGTTGGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATCGCCGCCGCCGCTGCCCGTCGTGGCGCGAACGTCACGCTg                                                                                            <  2:311195/247‑1 (MQ=255)
                                                                                                                                                   ggCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATCGCCGCCGCCGCTGCCCGTCGTGGCGCGAACGTCACGCTg                                                                                            <  1:159260/251‑1 (MQ=255)
                                                                                                                                                      cAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTTGCCCGTCAACGACCTGAAAAATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGGTATATCTCTAATAACAGCTCCCGCAAGATGGGGTTTGCTAt                                                                                                                                    >  1:268635/1‑208 (MQ=255)
                                                                                                                                                      cAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTAt                                                                                                                                    <  2:268635/208‑1 (MQ=255)
                                                                                                                                                              ggCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATCGCCGCCGCCGCTGCCCGTCGTGGCGCGAACGTCACGCTGGTATCAGGTcc                                                                                 <  1:79556/251‑1 (MQ=255)
                                                                                                                                                              ggCAGTCAGGCTTGGGGGTATTTCGGCCCGGGGCGAATGCTCGACCCGTTAACCATTGGGGATATGGCGGTAGCGCATTTTTCCCCCGTCAACGACCTGAAACATCGGAACATTATGATTCCCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATCGCCGCCGCCGCTGCCCGTCGTGGCGCGAACGTCACGCTGGTATCAGGTcc                                                                                 <  2:48813/251‑1 (MQ=255)
                                                                                                                                                                       gCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATc                                                                                                                                                                                                                                   <  2:244381/97‑1 (MQ=255)
                                                                                                                                                                       gCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAAACCTTGTGGCTACGCCGGTAGCGCAATTTTTGCCCCGCAACGACCCCAAACATCCGCAAAATTTGATTACCGCCGGCACGACGCGTGAAACGCCCGATCCGGGGGGTTATTTCTCTAATCACAGCTCCGGCAAAATGGGGTTTGCTCTCGCCCCCGCCGGTGCCCGGCGGGGCGCGAACGGCACGCTGGtgtacggtccggctttcct                                                                        >  1:400298/1‑233 (MQ=255)
                                                                                                                                                                                            gggCGAATGCTCTATCCGTTAACCATGTTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATCGCCGCCGCCGCTGCCCGTCGTGGCGCGAACGTCACGCTGGTAGCAGGTCCGGTTTCACTACCGACGCCACCGTTTGtt                                                     <  1:106009/249‑1 (MQ=255)
                                                                                                                                                                                                        gATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATCGCCGCCGCCGCTGCCCGTCGTGGCGCGAACGTCACGCTGGTATCAGGTCCGGTTTCACTACCGACGCCACCGTTTGTTAAACGTGTtgatgt                                       <  1:79247/251‑1 (MQ=255)
                                                                                                                                                                                                                         ggATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATCGCCGCCGCCGCTGCCCGTCGTGGCGCGAACGTCACGCTGGTATCAGGTCCGGTTTCACTACCGACGCCACCGTTTGTTAAACGTGTTGATGTGATGACCGCGCTGGaaa                      <  1:224986/251‑1 (MQ=255)
                                                                                                                                                                                                                                    ggtggcgctttttCGCCCCGTCAACGCCCTGAAACTTCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATCGCCGCCGCCGCTGCCCGTCGTGGCGCGAACGTCACGCTGGTATCAGGTCCGGTTTCACTACCGACGCCACCGTTTGTTAAACGTGTTGATGTGATGACCGCGCTGGAAATGGAAgccg             <  2:78685/243‑1 (MQ=255)
                                                                                                                                                                                                                                      gcgcATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATCGCCGCCGCCGCTGCCCGTCGTGGCGCGAACGTCACGCTGGTATCAGGTCCGGTTTCACTACCGACGCCACCGTTTGTTAAACGTGTTGATGTGATGACCGCGCTGGAAATGGAAGCCGCCGt         >  2:405827/1‑251 (MQ=255)
                                                                                                                                                                                                                                              tttCGCCCGTCAACGACCTGAAACATATGCACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATCGCCGCCGCCGCTGCCCGTCGTGGCGCGAACGTCAAGCTGGTATCAGGTCCGGTTTCACTACCGACGCCACCGTTTGTTAAACGTGTTGATGTGATGACCGCGCTGGAAATGGAAGGCGCCGTGAATGCt  >  2:217180/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                         
GAATGGCGAATGACCTGGTATCGACGATTTGTCTGGCTACACCTGCGCCTGTAGCCGTGCTCCCCGCCATGAACCAGCAGATGTACCGTGCCGCTGCCACGCAGCATAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATCGCCGCCGCCGCTGCCCGTCGTGGCGCGAACGTCACGCTGGTATCAGGTCCGGTTTCACTACCGACGCCACCGTTTGTTAAACGTGTTGATGTGATGACCGCGCTGGAAATGGAAGCCGCCGTGAATGCT  >  NC_000913/3813032‑3813519

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: