Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,806,6320CT40.5% 29.0 / 15.4 37S181L (TCG→TTG) coaBCbifunctional phosphopantothenoylcysteine decarboxylase/phosphopantothenate‑‑cysteine ligase CoaBC
Reads supporting (aligned to +/- strand):  ref base C (2/20);  new base T (15/0);  total (17/20)
Fisher's exact test for biased strand distribution p-value = 1.45e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTA  >  NZ_CP009273/3806549‑3806718
                                                                                   |                                                                                      
tGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGGGCCTTTTTTGcccct                                                                                  >  1:192171/1‑88 (MQ=255)
    aGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAAc                                                                              <  1:182822/90‑1 (MQ=255)
    aGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAAc                                                                              <  1:355040/90‑1 (MQ=255)
            ttGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGaa                                                                      <  2:256999/90‑1 (MQ=255)
            ttGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGaa                                                                      <  1:478957/90‑1 (MQ=255)
            ttGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGaa                                                                      <  1:422131/90‑1 (MQ=255)
                   gATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGGAGGGCATTTTTTTCCCGGCCACGACCCGAAAAATCTg                                                               >  2:29437/1‑90 (MQ=255)
                    atatCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTTAACGACCTGAAACATCTGa                                                              <  1:491964/90‑1 (MQ=255)
                          gTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCCTTTTTTGCCCGGCAACGACCTGAAACATCTGAACAttt                                                        >  1:271712/1‑89 (MQ=255)
                               gggCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTTGCCCGTCAAAGACCTGAAAAATCTGAACATTATGAtt                                                   >  1:34425/1‑90 (MQ=255)
                                    aaTGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTTGCCCGTCAACGACCTGAAACATCTGAACATTATGATTAccgc                                              >  1:24244/1‑90 (MQ=255)
                                    aaTGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTTGCCCGTCAACGACCTGAAACATCTGAACATTATGATTAccgc                                              >  1:108926/1‑90 (MQ=255)
                                          cGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTTGCCCCTCAACGACCCGAAAAAACTGAACATTATGATTACCGCCGGccc                                        >  2:245383/1‑90 (MQ=255)
                                           gATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCg                                       <  1:516906/90‑1 (MQ=255)
                                           gATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTCCGCCCGTCACCGACCTGAAACATTTGAACATTATGATTACCGCCGGCCCg                                       <  2:420552/90‑1 (MQ=255)
                                                        ttGTGGATATGGCGGTAGCGCATTTTTTGCCCGTCAACGACCTGAAAAATCTGAACATTATGATTACCGCCGGCCCGACGCGGGAACCGc                          >  1:233486/1‑90 (MQ=255)
                                                         tgtgGATATGGCGGTAGCGCATTTTTTTCCCCTCCACGACCCGAAACATCTGAACATTATGATTACCCCCCGCCCCACGCGTGGACCGcc                         >  2:205889/1‑89 (MQ=255)
                                                          gtgGATATGGCGGTAGCGCATTTTTTTCCCGGCAACCACCCGAAAAAACTGAAAATTATGATTTCCGCCGGCCCCACGCGGGAAcccccc                        >  2:486545/1‑86 (MQ=255)
                                                          gtgGATATGGCGGGAGCGCATTTTTCACCCCTCAACGACCCAAAAAAACTAAACATTATTATTAACCCCCGCCCGCCCCGTGAAACGccc                        >  1:509320/1‑88 (MQ=255)
                                                             gATATGGCGGTAGCGCATTTTTTCCCCGTCTACGACCCGAAAAATCTGAACATTATGATTACCCCCGGCCCGACGCGGGGAACGCTCGAt                     >  2:249554/1‑90 (MQ=255)
                                                               tatGGCGGTAGCGCATTTTTTGCCCGTCAACGACCTGAAAAATCTTAACATTATGATTACCGCCGGCCCGACGCGGGAACCGCCCGATcc                   >  1:376804/1‑90 (MQ=255)
                                                               tatGGCGGTAGCGCATTTTTCGCCCGTCAAGGCCCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATcc                   <  2:274072/90‑1 (MQ=255)
                                                               tatGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATcc                   <  2:221094/90‑1 (MQ=255)
                                                                 tGGCGGTAGCGCATTTTTTGCCCGTCAACGACCTGAAACATCTGAACATTATTATTACCGCCGGCCCGACGCGTGAACCCCTCGATCCgg                 >  1:407501/1‑90 (MQ=255)
                                                                 tGGCGGTAGCGCATTTTTTGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCCACGCGTGAACCGCCCGATcccg                 >  1:271322/1‑88 (MQ=255)
                                                                 tGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCgg                 <  1:359443/90‑1 (MQ=255)
                                                                   gCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTg               <  1:486545/90‑1 (MQ=255)
                                                                    cGGTAGCGCATTTTTTGCCCGTCAACGACCTGAAACATCTGAAAATTATGATTACCCCCGGCCCCACGCGTGAACCGCTCGATCCGGTgg              >  2:359443/1‑89 (MQ=255)
                                                                    cGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGc              <  1:534193/90‑1 (MQ=255)
                                                                    cGGTAGCGCATTTTTCGCCCGCCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGc              <  2:617795/90‑1 (MQ=255)
                                                                      gTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGt            <  1:417984/90‑1 (MQ=255)
                                                                         gcgcATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTa                                                        <  2:429000/43‑1 (MQ=255)
                                                                         gcgcATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTa                                                        >  1:429000/1‑43 (MQ=255)
                                                                         gcgcATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTtat         <  1:3435/90‑1 (MQ=255)
                                                                           gcATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTatat       <  2:24244/90‑1 (MQ=255)
                                                                              tttttCGCCCGTCAACGCCCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATAtctc    <  2:108926/90‑1 (MQ=255)
                                                                                tttCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTa  <  1:77073/90‑1 (MQ=255)
                                                                                   |                                                                                      
TGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTA  >  NZ_CP009273/3806549‑3806718

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: