Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273170,5520AC54.5% 19.9 / 27.0 23V273G (GTA→GGA) hemLglutamate‑1‑semialdehyde 2,1‑aminomutase
Reads supporting (aligned to +/- strand):  ref base A (9/1);  new base C (0/12);  total (9/14)
Fisher's exact test for biased strand distribution p-value = 2.01e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.45e-04
Rejected as consensus: Frequency below/above cutoff threshold.

TCGGGTTACCGGAAAGCGTACCCGCCTGATAGACCGGACCCGTCGGGGCCAGCGCATCCATTACATCACGACGACCACCGAATGCGCCTACCGGCATTCCACCGCCGATGATTTTGCCGAGGCAGGTTAAATCTGGCACTACGCCGTAATAATCCTGTGCGCCAGCTAGCGCTA  >  NZ_CP009273/170463‑170636
                                                                                         |                                                                                    
tCGGGTTACCGGAAAGCGTACCCGCCTGATAGACCGGACCCGTCGGGGCCAGCGCATCCATTACATCACGACGACCACCGAATGCGCCTa                                                                                      <  2:718538/90‑1 (MQ=255)
    gTTACCGGAAAGCGTACCCGCCTGATAGACCGGACCCGTCGGGGCCAGCGCATCCATTACATCACGACGACCACCGAATGCGCCTACCgg                                                                                  >  2:459275/1‑90 (MQ=255)
            aaaGCGTCCCCGCCTGATAGCCCGGCCCCGTGGGGGCCAGCGCATCCATTACATCACGACGACCACCGAATGCGCCTCCCGGCATTCCAc                                                                          <  1:365871/90‑1 (MQ=255)
            aaaGCGTACCCGCCTGATAGACCGGACCCGTCGGGGCCAGCGCATCCATTACATCACGACGACCACCGAATGCGCCTACCGGCATTCCAc                                                                          >  1:577209/1‑90 (MQ=255)
            aaaGCGTACCCGCCTGATAGACCGGACCCGTCGGGGCCAGCGCATCCATTACATCACGACGACCACCGAATGCGCCTACCGGCATTCCAc                                                                          >  2:398713/1‑90 (MQ=255)
                 gTACCCGCCTGATAGACCGGACCCGTCGGGGCCAGCGCATCCATTACATCACGACGACCACCGAATGCGCCTACCGGCATTCCAccgccg                                                                     >  1:481291/1‑90 (MQ=255)
                                aCCGGGCCGGTCGGGCCCCCCCCATCCGTTACATCCCGGCGACCCCGGAAGGCCCCTTACTCCTCTCCACCGCCGATGATTTTGCCGAgg                                                      <  2:76236/90‑1 (MQ=255)
                                 ccGGCCCCGTCGGGCCCAGCCCACCCATTACAACACGACGCCCACCGAATGCCCCTCCCGGCATTCCACCGCCGATGATTTTGCCGAGGc                                                     <  2:98108/90‑1 (MQ=255)
                                  cGGACCCGTCGGGGCCAGCGCATCCATTACATCACGACGACCACCGAATGCGCCTACCGGCATTCCACCGCCGATGATTTTGCCGAGGCa                                                    >  1:291547/1‑90 (MQ=255)
                                       ccGTGGGGCCCGGCGCATCCATTACATCCCGACGCCCACGGAATGCCCCTCCCGGCATCCCACCGCCGATGATTTTGCCGAGGCAGGTTa                                               <  2:450723/90‑1 (MQ=255)
                                         gTCGGGGCCAGCGCATCCATTACATCACGACGACCACCGAATGCGCCTACCGGCATTCCACCGCCGATGATTTTGCCGAGGCAGGTTaaa                                             >  1:570258/1‑90 (MQ=255)
                                                 cAGCGCACCCATTACTACCCGCCGCCCCCGGAATGCCCCTCCCGGCATTCCACCGCCGATGATTTTGCCGAGGCAGGTTAAATCTGGCAc                                     <  1:685958/90‑1 (MQ=255)
                                                   gcgcATCCATTACATGACGACAACCACCGAATGCGCCTACCGGCATTCCACCGCGGATGAGTTTGCCGAGGCAGGTTAAATCTGGCACTa                                   >  2:553978/1‑90 (MQ=255)
                                                   gcgcATCCATTACATCACGACGACCACCGAATGCGCCTACCGGCATTCCACCGCCGATGATTTTGCCGAGGCAGGTTAAATCTGGCACTa                                   >  1:787110/1‑90 (MQ=255)
                                                     ccATCCTTTACATCACGACGCCCCCCGAATGCCCCTCCCGGCATCCCACCGCCGATGATTTTGCCGAGGCAGGTTAAATCTGGCACTACg                                 <  1:522674/89‑1 (MQ=255)
                                                       ctCCATTACATCACGACGCCCCCCGAATGCGCCTCCCGGCATTCCACCGCCGATGATTTTGCCGAGGCAGGTTAAATCTGGCACTACGcc                               <  2:418567/89‑1 (MQ=255)
                                                       ccccATTACACCACGCCGCCCCCCGCAGGCCCCCCCCGGCATTCCACCGCCGATGATTTTGCCGAGGCAGGTTAAATCTGGCACTACGcc                               <  2:408823/88‑1 (MQ=255)
                                                       acccATTACATCACGACGCCCCCCGAATGCCCCTCCCGGCATCCCACCGCCGATGATTTTGCCGAGGCAGGTTAAATCTGGCACTACGcc                               <  1:799967/88‑1 (MQ=255)
                                                       aTCCTTTACATCACGACGACCACCGAATGCGCCTCCCGGCTTTCCCCCGCCGATGATTTTGCCGAGGCAGGTTAAATCTGGCACTACGcc                               <  2:291547/90‑1 (MQ=255)
                                                       aTCCATTACATCACGACGACCCCCGAATGCCCCTCCCGGCATTCCACCGCCGATGATTTTGCCGAGGCAGGTTAAATCTGGCACTACGcc                               <  1:125683/90‑1 (MQ=255)
                                                       aTCCATTACACCACGCCGACCACCGAATGCGCCTCCCGGCATTCCACCGCCGATGATTTTGCCGAGGCAGGTTAAATCTGGCACTACGcc                               <  1:12976/90‑1 (MQ=255)
                                                       aTCCATTACACCACGACGACCACCGAATGCGCCTACCGGCTTTCCACCGCCGATGATTTTGCCGAGGCAGGTTAAATCTGGCACTACGcc                               <  1:830848/90‑1 (MQ=255)
                                                       aTCAATTACATCACGCCGCCCCCCGAATGCGCCTCCCGCCATTCCACCGCCGATGATTTTGCCGAGGCAGGTTAAATCTGGCACTACGcc                               <  2:748892/90‑1 (MQ=255)
                                                                                    cgcCTACCGGCATTCCACCGCCGATGATTTTGCCGAGGCAGGTTAAATCTGGCACTACGCCGTAATAATCCTGTGCGCCAGCTAGCGCTa  >  1:678721/1‑90 (MQ=255)
                                                                                         |                                                                                    
TCGGGTTACCGGAAAGCGTACCCGCCTGATAGACCGGACCCGTCGGGGCCAGCGCATCCATTACATCACGACGACCACCGAATGCGCCTACCGGCATTCCACCGCCGATGATTTTGCCGAGGCAGGTTAAATCTGGCACTACGCCGTAATAATCCTGTGCGCCAGCTAGCGCTA  >  NZ_CP009273/170463‑170636

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: