Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,806,4320AC70.8% 6.8 / 18.4 24T114T (ACA→ACCcoaBCbifunctional phosphopantothenoylcysteine decarboxylase/phosphopantothenate‑‑cysteine ligase CoaBC
Reads supporting (aligned to +/- strand):  ref base A (4/3);  new base C (0/17);  total (4/20)
Fisher's exact test for biased strand distribution p-value = 3.29e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.93e-05
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GCCCCTGCCACGGCAGATTTGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCTGGTATCGACGATTTGTCTGGCTACACCTGCGCCTGTAGCCGTGCTCCCCGCCATGAACCAGCAGATGTACCGTGCCGCTGCCACGCAGCATAATT  >  NZ_CP009273/3806352‑3806502
                                                                                |                                                                      
gCCCCTGCCACGGCAGATTTGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCTGGTATCGACGATTTGTCTGGCTACACCTGCGCCt                                                               <  2:20917/90‑1 (MQ=255)
  cccTGCCACGGCAGATTTGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCTGGTATCGACGATTTGTCTGGCTACACCTGCGCCTGt                                                             <  2:814054/90‑1 (MQ=255)
    cTGCCACGGCAGATTTGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCTGGTATCGACGATTTGTCTGGCTACACCTGCGCCTGTAg                                                           >  1:521523/1‑90 (MQ=255)
                   tGATTGCCCGTGTTGCTGCCGGAATGGCGAATGCCCTGGTATCGACGATTTGTCTGGCTACCCCTGCGCCTGTAGCCGTGCCCCCCGCCa                                            <  2:278720/90‑1 (MQ=255)
                       tGCCTGTGTGGCTGCCGGACGGGCGAAGGCCCGGGTCCCGAGGTTTTGTCGGGCTACCCCTCCGCCTGTACCGGTGCTCCCCGCCATGaa                                        <  1:653222/90‑1 (MQ=255)
                          cggtgtTGCTGCCGAAATGGCGAATGCCCTGGTATCGACTTTTTGTCTGGCTACCCCTGCGCCTGTAGCCGGCCCCCCCGCCATGAACca                                     <  1:476716/88‑1 (MQ=255)
                          ccGTTTTGCTGCGGGAATGGCGAATGCCCTGGTATCGCCGATTTGTCTGGCTACCCCTGCGCCTGTAGCCGTGCCCCCCGCCATGAACca                                     <  2:504750/90‑1 (MQ=255)
                          ccGTGTTGCTGCGGGAAGGGCGAATGCCCTGGTATCGACGTTTTGTCTGGCTACCCCTGCGCCTGTAGCCGTGCTCCCCGCCATGAACca                                     <  2:521524/90‑1 (MQ=255)
                          ccGTGTTGCTGCCGGAATGGCGAATGCCCTGTTATCGACGATTTGTCTGGCTCCCCCTGCGCCTGTAGCCGTGCCCCCCGCCATGAACca                                     <  2:478561/90‑1 (MQ=255)
                          ccGTGTTGCTGCCGGAAGGGCGAATGCCCTGGTATCGACGATTTGTCTGGCTACCCCTGCGCCTGTAGCCGTGCTCCCCGCCATGAACca                                     <  1:306085/90‑1 (MQ=255)
                          ccGTGTTGCGGCGGGAATGGCGAATGACCTGGTATCGACGTTTTGTCTGGCTACCCCTGCGCCTGTAGCCGTGCCCCCCGCCATGAACca                                     <  2:514842/90‑1 (MQ=255)
                            gtgtTGCTGCCGGAAGGGCAAATGCCCTGGTATCGACGTTTTGTCTGGCTACCCCTGCGCCTGTAGCCGTGCCCCCCGCCATGAACcagc                                   <  1:274181/90‑1 (MQ=255)
                            gggtTGCTGCGGGAAGGGGAAAGGCCCTGGTACCGACGTTTTGTCGGGCTCCCCCTGCGCTTGTGGCCGTGCCCCCCGCCATGAACcagc                                   <  1:735941/88‑1 (MQ=255)
                                                  tGACCTGGTATCGACGATTTGTCTGGCTACACCTGCGCCTGTAGCCGTGCTCCCCGCCATGAACCAGCAGATGTACCGTGCCGCTGCCAc             >  1:428066/1‑90 (MQ=255)
                                                    cccTGGTATCGACGTTTTGTCGGGCTCCCCCTGCGCCTGTAGCCGTGCTCCCCGCCATGAACCAGCAGATGTACCGTGCCGCTGCCACgc           <  1:38857/89‑1 (MQ=255)
                                                      cTGGTATCGACGATTTGTCTGGCTACACCTGCGCCTGTAGCCGTGCTCCCCGCCATGAACCAGCAGATGTACCGTGCCGCTGCCACgcag         >  2:53124/1‑90 (MQ=255)
                                                         gTATCGACGTTTTGTCTGGCTCCCCCTGCCCCTGTAGCCGTGCCCCCCGCCATGAACCAGCAGATGTACCGTGCCGCTGCCACGCAGCAt      <  1:278341/90‑1 (MQ=255)
                                                         gTATCGACGTTTTGTCTGGCTACCCCTGCGCCTGTAGCCGTGCTCCCCGCCATGAACCAGCAGATGTACCGTGCCGCTGCCACGCAGCAt      <  1:413111/90‑1 (MQ=255)
                                                         gTATCGACGTTTTGTCTGGCTACCCCTGCGCCTGTAGCCGTGCCCCCCGCCATGAACCAGCAGATGTACCGTGCCGCTGCCACGCAGCAt      <  1:291930/90‑1 (MQ=255)
                                                         gTATCGACGTTTTGTCGGGCTACCCCTGCGCCTGTAGCCGTGCCCCCCGCCATGAACCAGCAGATGTACCGTGCCGCTGCCACGCAGCAt      <  2:510905/90‑1 (MQ=255)
                                                         gTATCGACGATTTGTCTGGCTCCCCCTGCGCCTGTAGCCGTCCTCCCCGCCATGAACCAGCAGATGTACCGTGCCGCTGCCACGCAGCAt      <  1:660565/90‑1 (MQ=255)
                                                         gTATCGACGATTTGTCTGGCTACACCTGCGCCTGTAGCCGTGCTCCCCGCCATGAACCAGCAGATGTACCGTGCCGCTGCCACGCAGCAt      <  2:428066/90‑1 (MQ=255)
                                                         gTATCGACGATTTGTCGGGCTACCCCTGCGCCTGTAGCCGTGCCCCCCGCCATGAACCAGCAGATGTACCGTGCCGCTGCCACGCAGCAt      <  1:790849/90‑1 (MQ=255)
                                                             cgacgaTTTGTCTGGCTACACCTGCGCCTGTAGCCGTGCTCCCCGCCATGAACCAGCAGATGTACCGTGCCGCTGCCACGCAGCATAAtt  >  1:78067/1‑90 (MQ=255)
                                                                                |                                                                      
GCCCCTGCCACGGCAGATTTGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCTGGTATCGACGATTTGTCTGGCTACACCTGCGCCTGTAGCCGTGCTCCCCGCCATGAACCAGCAGATGTACCGTGCCGCTGCCACGCAGCATAATT  >  NZ_CP009273/3806352‑3806502

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