Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,806,6320CT37.9% 17.9 / 12.1 29S181L (TCG→TTG) coaBCbifunctional phosphopantothenoylcysteine decarboxylase/phosphopantothenate‑‑cysteine ligase CoaBC
Reads supporting (aligned to +/- strand):  ref base C (5/13);  new base T (11/0);  total (16/13)
Fisher's exact test for biased strand distribution p-value = 1.29e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.25e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATA  >  NZ_CP009273/3806543‑3806712
                                                                                         |                                                                                
aGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTc                                                                                  <  2:287449/90‑1 (MQ=255)
     gTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTTGcccc                                                                             >  2:269336/1‑89 (MQ=255)
                         gATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTTGCCCGTCAAAGACCTGAAACATCTg                                                         >  2:151299/1‑90 (MQ=255)
                         gATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTg                                                         <  1:220464/90‑1 (MQ=255)
                                     gggCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCCTCAAAGACCTGAAAAATCTGAACATTATGAtt                                             >  2:119175/1‑90 (MQ=255)
                                     gggCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGGAGGGCATTTTTTGCCCGGCAACGACCCGGAACATCTGAAAATTATGAtt                                             >  2:425910/1‑90 (MQ=255)
                                                 gATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCg                                 <  2:50031/90‑1 (MQ=255)
                                                 gATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCg                                 <  1:441534/90‑1 (MQ=255)
                                                 gATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCg                                 <  2:134917/90‑1 (MQ=255)
                                                 cagccGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCg                                 <  1:181848/87‑1 (MQ=255)
                                                      gTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGAcgcg                            >  2:160448/1‑90 (MQ=255)
                                                        tAACCATTGTGGATATGGCGGTAGGGCCTTTTTTCCCCCTCCACGACCCGAAACAACTGAACATTATGATTTACCCCGGCCCCCcgcggg                          >  1:10573/1‑88 (MQ=255)
                                                        tAACCATTGTGGATATGGCGGTAGCGCCTTTTTTTCCCGTCCACCAACTTAAAAAACTTTACATTTTTATTTACCCCGGCCCCCcgcggg                          >  1:278158/1‑88 (MQ=255)
                                                            cATTGTGGATATGGCGGTAGCGCATTTTTTGCCCGTCAACGACCTGAAAAATCTGAACATTATGATTACCGCCCGCCCGACGCGGGAAcc                      >  2:65526/1‑90 (MQ=255)
                                                              ttGTGGATATGGCGGTAGCGCATTTTTTGCCCCTCAACGACCTGAAAAATCTGAACATTATGATTACCCCCCGCCCCACCCGGGAAcccg                    >  2:181848/1‑88 (MQ=255)
                                                              ttGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCa                                                                          <  1:121620/36‑1 (MQ=255)
                                                              ttGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCa                                                                          >  2:121620/1‑36 (MQ=255)
                                                                gtgGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTc                  <  2:406953/90‑1 (MQ=255)
                                                                  ggATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGa                <  1:151299/90‑1 (MQ=255)
                                                                     tatGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCCACGCGTGAACCGCTCGATcc             >  2:2355/1‑90 (MQ=255)
                                                                       tGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCg                                 >  2:146317/1‑68 (MQ=255)
                                                                       tGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCg                                 <  1:146317/68‑1 (MQ=255)
                                                                       tGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCgg           <  2:278158/90‑1 (MQ=255)
                                                                       tGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCgg           <  2:423877/90‑1 (MQ=255)
                                                                       tGGCGGGAGCGCATTTTTTTCCCCTCAACGACCTGAAAAAACTTAAAAATATGATTATCCCCGGCCCCACGCGGGGACCCCCCCATCCgg           >  1:423877/1‑90 (MQ=255)
                                                                        ggCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGt          <  2:84219/90‑1 (MQ=255)
                                                                        ggCGGGAGCGCCTTTTTTTCCCGTCAACGACCCGAAAAATCTGAACATTATGATTACCGCCCGCCCGACGCGGGGACCGCCCCGTCCggg          >  1:84219/1‑89 (MQ=255)
                                                                          cGGTAGCGCATTTTTTGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCCACGCGGGAACCCCTCGATCCGGTGc        >  1:128984/1‑90 (MQ=255)
                                                                          cGGTAGCGCATTTTTTGCCCCTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCCACGCGGGAACCGCCCGATCCGGTGc        >  1:329830/1‑90 (MQ=255)
                                                                                cgcATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTtata  <  1:269336/90‑1 (MQ=255)
                                                                                         |                                                                                
AGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATA  >  NZ_CP009273/3806543‑3806712

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