Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,837,5790TG53.8% 27.5 / 31.5 26F86V (TTT→GTT) adeDadenine deaminase
Reads supporting (aligned to +/- strand):  ref base T (2/10);  new base G (14/0);  total (16/10)
Fisher's exact test for biased strand distribution p-value = 1.24e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.27e-04
Rejected as consensus: Frequency below/above cutoff threshold.

CATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTG  >  NZ_CP009273/3837497‑3837659
                                                                                  |                                                                                
cATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGGGCCCGGGGTTATTGa                                                                           >  2:115383/1‑90 (MQ=255)
  ttGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGGTTATTGATg                                                                         >  1:139113/1‑90 (MQ=255)
   tGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGc                                                                        <  1:381296/90‑1 (MQ=255)
    gCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGGTTATTGATGCt                                                                       >  1:466085/1‑90 (MQ=255)
       ggTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCAc                                                                    <  1:411220/90‑1 (MQ=255)
            tGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCa                                                               <  1:205478/90‑1 (MQ=255)
            tGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCa                                                               <  1:228560/90‑1 (MQ=255)
                 cAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTg                                                          <  1:242095/90‑1 (MQ=255)
                  aGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGGTTATTGATGCTCACCTGCATATTGa                                                         >  1:200351/1‑90 (MQ=255)
                    aaTACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAAt                                                       <  2:266101/90‑1 (MQ=255)
                             aTGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCatga                                              >  1:346916/1‑90 (MQ=255)
                             aTGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGGTTTTTTATGCTCCCCCGCATATTTAATCCCGCatga                                              >  2:452802/1‑90 (MQ=255)
                             aTGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGGTTTTTTATGCTCACCTGCATATTGAATCCAGCatga                                              >  2:396626/1‑90 (MQ=255)
                             aTGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGGTTTTTTATGCTCACCCGCATATTGAATCCAGCatga                                              >  2:270421/1‑90 (MQ=255)
                             aTGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGGTTATTTATGCTCACCTGCATATTGAATCCAGCatga                                              >  2:76619/1‑90 (MQ=255)
                             aTGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGGTTATTGATGCTCACCTGCATATTGAATCCCGCatga                                              >  1:305064/1‑90 (MQ=255)
                             aTGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGGGTTTTTATTCTCCCCTGCATATTGAATCCAGCatga                                              >  2:384238/1‑90 (MQ=255)
                               gCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGGTTTTTGAGGCTCCCCCGCCTATTGAATCCCGCctggtg                                            >  2:454767/1‑90 (MQ=255)
                               gCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGGTTATTGATGCTCACCTGCATATTGAATCCAGCatgatg                                            >  1:230172/1‑90 (MQ=255)
                                    ggCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGGTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGcc                                       >  1:303609/1‑90 (MQ=255)
                                               ggATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTg                            >  1:67523/1‑90 (MQ=255)
                                                 aTTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGaa                          <  1:432176/90‑1 (MQ=255)
                                                 aTTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGaa                          <  2:133013/90‑1 (MQ=255)
                                                 aTTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGaa                          <  1:375656/90‑1 (MQ=255)
                                                                gcAACGGCGGTGCCAGGGGTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGGCACTTTTTAAACCGCCACCccgccg           >  2:405112/1‑90 (MQ=255)
                                                                         gTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTg  <  1:308481/90‑1 (MQ=255)
                                                                                  |                                                                                
CATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTG  >  NZ_CP009273/3837497‑3837659

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