Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273642,6660AC46.3% 68.9 / 20.1 67G100G (GGT→GGGcitGtriphosphoribosyl‑dephospho‑CoA synthase CitG
Reads supporting (aligned to +/- strand):  ref base A (29/7);  new base C (0/31);  total (29/38)
Fisher's exact test for biased strand distribution p-value = 1.08e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AGCAAACGGCCAATTGCCGCACATAGCAGCCCTAAAGAAAAAATGCTGCCTTTATGCGTGTTTACGCCCGCAGTGGCGCGGAACATATCACCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCACTACAGGCACCAAATTCAAT  >  NZ_CP009273/642577‑642749
                                                                                         |                                                                                   
aGCAAACGGCCAATTGCCGCACATAGCAGCCCTAAAGAAAAAATGCTGCCTTTATGCGTGTTTACGCCCGCAGTGGCGCGGAACATATCa                                                                                     >  2:1050414/1‑90 (MQ=255)
aGCAAACGGCCAATTGCCGCACATAGCAGCCCTAAAGAAAAAATGCTGCCTTTATGCGTGTTTACGCCCGCAGTGGCGCGGAACATATCa                                                                                     >  2:457130/1‑90 (MQ=255)
aGCAAACGGCCAATTGCCGCACATAGCAGCCCTAAAGAAAAAATGCTGCCTTTATGCGTGTTTACGCCCGCAGTGGCGCGGAACATATCa                                                                                     >  2:1100158/1‑90 (MQ=255)
   aaaCGGCCAATTGCCGCACATAGCAGCCCTAAAGAAAAAATGCTGCCTTTATGCGTGTTTACGCCCGCAGTGGCGCGGAACATATCACCt                                                                                  >  2:637647/1‑90 (MQ=255)
   aaaCGGCCAATTGCCGCACATAGCAGCCCTAAAGAAAAAATGCTGCCTTTATGCGTGTTTACGCCCGCAGTGGCGCGGAACATATCACCt                                                                                  >  2:2265337/1‑90 (MQ=255)
   aaaCGGCCAATTGCCGCACATAGCAGCCCTAAAGAAAAAATGCTGCCTTTATGCGTGTTTACGCCCGCAGTGGCGCGGAACATATCACCt                                                                                  >  2:2358120/1‑90 (MQ=255)
       ggCCAATTGCCGCACATAGCACCCCAAAAAAAAAAATGCTCCCTTTATGCGTTTTTACGCCCGCAGTGGCGCGGAACATATCCCCTTCGc                                                                              <  2:2547566/90‑1 (MQ=255)
          cAATTGCCGCACATAGCAGCCCTAAAGAAAAAATGCTGCCTTTATGCGTGTTTACGCCCGCAGTGGCGCGGAACATATCCCCTTCGCAAg                                                                           <  1:1050331/90‑1 (MQ=255)
          cAATTGCCGCACATAGCAGCCCTAAAGAAAAAATGCTGCCTTTATGCGTGTTTACGCCCGCAGTGGCGCGGAACATATCACCTTCGCAAg                                                                           <  1:1307209/90‑1 (MQ=255)
            aTTGCCGCACATAGCAGCCCTAAAGAAAAAATGCTGCCTTTATGCGTGTTTACGCCCGCAGTGGCGCGGAACATATCACCTTCGCAAGcc                                                                         <  1:403005/90‑1 (MQ=255)
            aTTGCCGCACATAGCAGCCCTAAAGAAAAAATGCTGCCTTTATGCGTGTTTACGCCCGCAGTGGCGCGGAACATATCACCTTCGCAAGcc                                                                         <  1:3219110/90‑1 (MQ=255)
            aTTGCCGCACATAGCAGCCCAAAAAAAAAAATGCTGCCTTTATGCGTGTTTACGCCCGCAGTGGCGCGGAACATATCCCCTTCGCAAGcc                                                                         <  2:1645378/90‑1 (MQ=255)
            aTTGCCGCACATAGCAGCCAAAAAAAAAAAATGCTCCCTTTATGCGTGTTTCCGCCCGCAGTGGCGCGGAACACATCCCCTTCGCAAGcc                                                                         <  2:1932209/90‑1 (MQ=255)
            aTTGCCGCACATAGCACCCCAAAAAAAAAAATGCTGCCTTTATGCGTGTTTCCGCCCGCAGTGGCGCGGAACATATCCCCTTCGCAAGcc                                                                         <  2:2486729/90‑1 (MQ=255)
            aTTCCCGCACATAGCAGCCCTAAAAAAAAAATGCTGCCTTTATGCGTTTTTACGCCCGCAGTGGCGCGAAACATATCCCCTTCGCAAGcc                                                                         <  2:2489970/90‑1 (MQ=255)
             ttGCCGCACATAGCAGCCCTAAAGAAAAAATGCTGCCTTTATGCGTGTTTACGCCCGCAGTGGCGCGGAACATATCACCTTCGCAAGCCa                                                                        >  2:2553943/1‑90 (MQ=255)
                     cacagcagcCTTAAAGAAAAAACCCTCCCTTTATGCTTGTTTCCGCCGGCGGTGGCGGGAAAAATATCCCCTTCGCAAGCCATACCAAtt                                                                <  2:1059498/87‑1 (MQ=255)
                     caTAGCAGCCCTAAAGAAAAAATGCTGCCTTTATGCGTGTTTCCGCCCGCAGTGGCGCGGAACATATCCCCTTCGCAAGCCATACCAAtt                                                                <  1:1282789/90‑1 (MQ=255)
                     caTAGCACCCCTAAAAAAAAAATGCTCCCTTTATGCGTGTTTACGCCCGCAGTGGCGCGAAACATATCCCCTTCGCAAGCCATACCAAtt                                                                <  2:1830635/90‑1 (MQ=255)
                       tAGCAGCCCTAAAGAAAAAATGCTGCCTTTATGCGTGTTTACGCCCGCAGTGGCGCGGAACATATCCCCTTCGCAAGCCATACCAATTgg                                                              <  1:1113519/90‑1 (MQ=255)
                       tAGCAGCCCTAAAGAAAAAATGCTGCCTTTATGCGTGTTTACGCCCGCAGTGGCGCGGAACATATCACCTTCGCAAGCCATACCAATTgg                                                              >  1:2692406/1‑90 (MQ=255)
                       tAGCAGCCCTAAAGAAAAAATGCTGCCTTTATGCGTGTTTACGCCCGCAGTGGCGCGGAACATATCACCTTCGCAAGCCATACCAATTgg                                                              >  1:1297191/1‑90 (MQ=255)
                       tAGCAGCCCTAAAGAAAAAATGCTGCCTTTATGCGTGTTTACGCCCGCAGTGGCGCGGAACATATCACCTTCGCAAGCCATACCAATTgg                                                              >  2:1526478/1‑90 (MQ=255)
                          cagcCCTAAAGAAAAAATGCTGCCTTTATGCGTGTTTACGCCCGCAGTGGCGCGGAACATATCACCTTCGCAAGCCATACCAATTGGGCg                                                           >  1:60293/1‑90 (MQ=255)
                          cagcCCTAAAGAAAAAAGGCTGCCTTTATGCGTGTTTACGCCCGCAGTGGCGCGGAACATATCCCCTTCGCAAGCCATACCAATTGGGCg                                                           <  1:2672898/90‑1 (MQ=255)
                          cagcCCAAAAAAAAAAATGCTCCCTTTATGCGTGTTTACGCCCGCAGTGGCGCGGAACATATCCCCTTCGCAAGCCATACCAATTGGGCg                                                           <  2:167116/90‑1 (MQ=255)
                          caccccAAAAAAAAAAATGCTGCTTTTATGCGTTTTTCCGCCCGCAGTGGCGCGGAACATATCCCCTTCGCAAGCCATACCAATTGGGCg                                                           <  2:2042756/87‑1 (MQ=255)
                            gcCCTAAAGAAAAAATGCTGCCTTTATGCGTGTTTACGCCCGCAGTGGCGCGGAACATATCACCTTCGCAAGCCATACCAATTGGGCGTa                                                         >  1:2643645/1‑90 (MQ=255)
                            gcCCTAAACAAAAAATGCTCCCTTTATGCGCGTTTCCGCCCGCAGTGGGGCGGAACATATCACCTTCGCAAGCCATACCAATTGGGCGTa                                                         <  1:940356/90‑1 (MQ=255)
                            gcCCTAAAAAAAAAATGCTGCCTTTATGCGTGTTTACGCCCGCAGGGGCGCGGAACATATCCCCTTCGCAAGCCATACCAATTGGGCGTa                                                         <  2:978134/90‑1 (MQ=255)
                             cccTAAAGAAAAAATGCTGCCTTTATGCGTGTTTACGCCCGCAGTGGCGCGGAACATATCACCTTCGCAAGCCATACCAATTGGGCGTaa                                                        <  1:2358120/90‑1 (MQ=255)
                             cccTAAAGAAAAAATGCTGCCTTTATGCGTGTTTACGCCCGCAGTGGCGCGGAACATATCACCTTCGCAAGCCATACCAATTGGGCGTaa                                                        <  1:1100158/90‑1 (MQ=255)
                               caaaaGAAAAAATGCTGCCTTTATGCGTGTTTACGCCCGCAGTGGCGCGGAACATATCCCCTTCGCAAGCCATACCAATTGGGCGTAATc                                                      <  1:1940074/88‑1 (MQ=255)
                                tAAAGAAAAAATGCTGCCTTTATGCGTGTTTACGCCCGCAGTGGCGCGGAACATATCACCTTCGCAAGCCATACCAATTGGGCGTAATcc                                                     >  1:491042/1‑90 (MQ=255)
                                 aaaGAAAAAATGCTGCCTTTATGCGTGTTTACGCCCGCAGTGGCGCGGAACATATCACCTTCGCAAGCCATACCAATTGGGCGTAATCCg                                                    >  1:1437346/1‑90 (MQ=255)
                                 aaaGAAAAAATGCTGCCTTTATGCGTGTTTACGCCCGCAGTGGCGCGGAACATATCACCTTCGCAAGCCATACCAATTGGGCGTAATCCg                                                    >  2:2598950/1‑90 (MQ=255)
                                 aaaGAAAAAATGCTGCCTTTATGCGTGTTTACGCCCGCAGTGGCGCGGAACATATCACCTTCGCAAGCCATACCAATTGGGCGTAATCCg                                                    >  1:1674485/1‑90 (MQ=255)
                                 aaaGAAAAAATGCTGCCTTTATGCGTGTTTACGCCCGCAGTGGCGCGGAACATATCACCTTCGCAAGCCATACCAATTGGGCGTAATCCg                                                    >  2:160372/1‑90 (MQ=255)
                                 aaaGAAAAAATGCTGCCTTTATGCGTGTTTACGCCCGCAGTGGCGCGGAACATATCACCTTCGCAAGCCATACCAATTGGGCGTAATCCg                                                    >  2:3104364/1‑90 (MQ=255)
                                 aaaGAAAAAATGCTGCCTTTATGCGTGTTTACGCCCGCAGTGGCGCGGAACATATCACCTTCGCAAGCCATACCAATTGGGCGTAATCCg                                                    >  2:772182/1‑90 (MQ=255)
                                 aaaGAAAAAATGCTGCCTTTATGCGTGTTTACGCCCGCAGTGGCGCGGAACATATCACCTTCGCAAGCCATACCAATTGGGCGTAATCCg                                                    >  1:1875653/1‑90 (MQ=255)
                                 aaaGAAAAAATGCTGCCTTTATGCGTGTTTACGCCCGCAGTGGCGCGGAACATATCACCTTCGCAAGCCATACCAATTGGGCGTAATCCg                                                    >  1:78771/1‑90 (MQ=255)
                                    aaaaaaaGGCTTCCTTTATACGTTTTTCCGCCGGGAGTGGCGCGAAACATACCCCCCTCGCAAGCCATACCAATTGGGCGTAATCCGTgg                                                 <  2:1437346/89‑1 (MQ=255)
                                     aaaaaaTGCTGCCTTTAGGCGTTTTTCCCCCCGGGGGGGGGCGGAACATATCCCCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGa                                                <  2:1914078/90‑1 (MQ=255)
                                       aaaaTGCTGCCTTTATGCGTGTTTACGCCCGCAGTGGCGCGGAACATATCACCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGt                                              >  1:1288884/1‑90 (MQ=255)
                                           tgctgcCTTTATGCGTGTTTACGCCCGCAGTGGCGCGGAACATATCACCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCg                                          >  2:1147837/1‑90 (MQ=255)
                                           tgctgcCTTTATGCGTGTTTACGCCCGCAGTGGCGCGGAACATATCACCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCg                                          >  2:3079272/1‑90 (MQ=255)
                                                 cTTTATGCGTTTTTCCGCCCGGGGGGGCGGGGAACATATCCCCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTg                                    <  2:2814727/90‑1 (MQ=255)
                                                 cTTTATGCGTGTTTCCCCCCGCAGGGGCGCGAAACATCCCCCCTCCCCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTg                                    <  2:1288884/90‑1 (MQ=255)
                                                      tGCGTTTTTACGCCCGCAGTGGCGCGAAACATATCCCCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGcc                               <  1:2287387/90‑1 (MQ=255)
                                                      tGCGTGTTTACGCCCGCAGGGGTGCGGAACATATCACCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGcc                               <  1:1786306/90‑1 (MQ=255)
                                                         gtgtTTACGCCCGCAGTGGCGCGGAACATATCACCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCAtt                            >  1:1196757/1‑90 (MQ=255)
                                                          tgtTTACGCCCGCAGTGGCGCGGAACATATCACCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCAttt                           >  2:2401232/1‑90 (MQ=255)
                                                             ttACGCCCGCAGGGGGGCGGAACATATCACCTTCGCAAGCCATAACAATTGGGCGTAATCCGTGGAGCGCCGCTTCTGGTGCCATTTCCg                        >  2:3296631/1‑90 (MQ=255)
                                                               ccGCCCGCAGTGGCGCGGAACATATCCCCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCa                      <  1:110792/89‑1 (MQ=255)
                                                                    cGCAGTGGCGCGGAACATATCACCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCACTACa                 >  2:2634539/1‑90 (MQ=255)
                                                                     gCAGTGGCGCGGAACATATCACCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGATTCTGGTGCCATTTCCGCACTACAg                >  1:1033813/1‑90 (MQ=255)
                                                                      gggTGGCGCGGAACATCCCCCCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCACTACAgg               <  2:686824/88‑1 (MQ=255)
                                                                      cAGTGGCGCGGAACATATCCCCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCACTACAgg               <  1:400876/90‑1 (MQ=255)
                                                                      cAGTGGCGCGGAACATATCCCCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCACTACAgg               <  1:3073352/90‑1 (MQ=255)
                                                                      cAGTGGCGCGGAACATATCCCCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCACTACAgg               <  2:2509125/90‑1 (MQ=255)
                                                                        gTGGCGCGGAACATATCCCCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCACTACAGGCa             <  1:2428335/90‑1 (MQ=255)
                                                                           gcgcGGAACATATCCCCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCACTACAGGCACCa          <  1:2234960/90‑1 (MQ=255)
                                                                             gcgGAACATATCACCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCACTACAGGCACCaaa        >  2:2904437/1‑90 (MQ=255)
                                                                              ggaaaCATATCCCCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCACTACAGGCACCAAAt       <  2:1674485/87‑1 (MQ=255)
                                                                              cgGAACATATCCCCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCACTACAGGAACCAAAt       <  1:1533304/90‑1 (MQ=255)
                                                                                gAACATATCACCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCACTACAGGCACCAAATTc     >  1:344832/1‑90 (MQ=255)
                                                                                   cATATCCCCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCACTACAGGCACCAAATTCAAt  <  1:2894273/90‑1 (MQ=255)
                                                                                    atatCCCCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGc                       <  2:2771305/68‑1 (MQ=255)
                                                                                    atatCACCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGc                       >  1:2771305/1‑68 (MQ=255)
                                                                                         |                                                                                   
AGCAAACGGCCAATTGCCGCACATAGCAGCCCTAAAGAAAAAATGCTGCCTTTATGCGTGTTTACGCCCGCAGTGGCGCGGAACATATCACCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCACTACAGGCACCAAATTCAAT  >  NZ_CP009273/642577‑642749

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