Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,485,3870CT53.1% 7.2 / 43.0 33Q90Q (CAG→CAAppiApeptidylprolyl isomerase A
Reads supporting (aligned to +/- strand):  ref base C (0/15);  new base T (17/0);  total (18/15)
Fisher's exact test for biased strand distribution p-value = 1.77e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.33e-06
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TCAGCGGTACGTGCCATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGCGGGTTTGGTTTTTTCTGCTGCATCTGCTCGGTGAAACCGCCGCCCTGAATCATAAAGCCAGGAATGACGCGGTGAAAGGTAGTGTTGTTATAAAAACCGC  >  NZ_CP009273/3485301‑3485472
                                                                                       |                                                                                      
tCAGCGGTACGTGCCATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGCGGGTTTGGTTTTTTctgc                                                                                     <  1:35799/90‑1 (MQ=255)
 cAGCGGTACGTGCCATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGCGGGTTTGGTTTTTTctgct                                                                                    <  1:336380/90‑1 (MQ=255)
   gCGGTACGTGCCATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGCGGGGTTTGTTTTTTTttctgc                                                                                  >  2:610871/1‑90 (MQ=255)
   gCGGTACGTGCCATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGCGGGGTTTGTTTTTTTttctgc                                                                                  >  2:336380/1‑90 (MQ=255)
   gCGGTACGTGCCATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGCGGGGTTGGGTTTTTTttctgc                                                                                  >  2:358202/1‑90 (MQ=255)
              cATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGCGGGTTTTGTTTTTTTTGCTGCCTCTGCTCggg                                                                       >  2:143831/1‑89 (MQ=255)
              cATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGCGGGGTTGGGTTTTTTTTCTGCCTCTTCTCggg                                                                       >  2:740618/1‑89 (MQ=255)
                tCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGCGGGGTTGGTTTTTTTTGTTGCCTCTTCTCGGTgg                                                                    >  1:34904/1‑89 (MQ=255)
                           ccACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGCGGGTTTGGTTTTTTCTGCTGCATCTGCTCGGTGAAACCGCCGccc                                                         <  1:143831/90‑1 (MQ=255)
                             aCGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGCGGGGTTTGTTTTTTTTGCTGCATCTTCTCGGGGAAAACGcccccccg                                                       >  2:381148/1‑88 (MQ=255)
                             aCGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGCGGGGTTGGTTTTTTTTTCTGCTTCTGCTCGGTGGAAACccccccccc                                                        >  1:778305/1‑89 (MQ=255)
                             aCGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGCGGGGTTGGTTTTTTTTGCTGCATCTGCTCGGGGAAAACGCCGCCCTg                                                       >  2:260605/1‑90 (MQ=255)
                               gcgTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGCGGGTTTGGTTTTTTTTTTTGTATCTTCTCTTTGGAAACGcccccccgga                                                     >  1:599189/1‑86 (MQ=255)
                               gcgTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGCGGGTTTGGTTTTTTCTGCTGCATCTGCTCGGTGAAACCGCCGCCCTGaa                                                     <  2:545749/90‑1 (MQ=255)
                                        gcAGGCCGTTATCGGCTTCATTTTTGATTGGCGGGGTTTGGTTTTTTTGCTGGCCCTGGTCGGGGTAAACGCCCCCCCCAATAATaaaac                                            >  1:46708/1‑88 (MQ=255)
                                        gcAGGCCGTTATCGGCTTCATTTTTGATTGGCGGGTTTGGTTTTTTTTGTTGCATTTGCTCGGGGAAACCGCCCCCCTGAAACATaaaac                                            >  2:90183/1‑88 (MQ=255)
                                                    cGGCTTCATTTTTGATTGGCGGGTTTGGTTTTTTCTGCTGCATCTGCTCGGTGAAACCGCCGCCCTGAATCATAAAGCCAGGAATGAcgc                                <  1:88280/90‑1 (MQ=255)
                                                    cGGCTTCATTTTTGATTGGCGGGTTTGGTTTTTTCTGCTGCATCTGCTCGGTGAAACCGCCGCCCTGAATCATAAAGCCAGGAATGAcgc                                <  2:232461/90‑1 (MQ=255)
                                                    cGGCTTCATTTTTGATTGGCGGGTTTGGTTTTTTCTGCTGCATCTGCTCGGTGAAACCGCCGCCCTGAATCATAAAGCCAGGAATGAcgc                                <  2:321704/90‑1 (MQ=255)
                                                    cGGCTTCATTTTTGATTGGCGGGTTTGGTTTTTTCTGCTGCATCTGCTCGGTGAAACCGCCGCCCTGAATCATAAAGCCAGGAATGAcgc                                <  1:430507/90‑1 (MQ=255)
                                                    cGGCTTCATTTTTGATTGGCGGGTTTGGTTTTTTCTGCTGCATCTGCTCGGTGAAACCGCCGCCCTGAATCATAAAGCCAGGAATGAcgc                                <  1:272567/90‑1 (MQ=255)
                                                    cGGCTTCATTTTTGATTGGCGGGTTTGGTTTTTTCTGCTGCATCTGCTCGGTGAAACCGCCGCCCTGAATCATAAAGCCAGGAATGAcgc                                <  1:358202/90‑1 (MQ=255)
                                                    cGGCTTCATTTTTGATTGGCGGGTTTGGTTTTTTCTGCTGCATCTGCTCGGTGAAACCGCCGCCCTGAATCATAAAGCCAGGAATGAcgc                                <  2:450305/90‑1 (MQ=255)
                                                     ggCTTCATTTTTGATTGGCGGGTTTTGTTTTTTTTTCTGCATCTTCTCGGGGAAAACGCCGCCCCGGATCATAAAACCAGGAAAGAcgcg                               >  2:699014/1‑90 (MQ=255)
                                                     ggCTTCATTTTTGATTGGCGGGTTTGGTTTTTTTTGCTGCATCTGCTCGGGGAAACCGCCGCCCTGAATCATAAAGCCAGGGATGAcgcg                               >  2:25775/1‑90 (MQ=39)
                                                     ggCTTCATTTTTGATTGGCGGGGTTGGGTTTTTTTTCTGCATCTTCTCGGTGGAAACCCCCCCCTGAATCATAAAAACAGGAATGAggcg                               >  2:88280/1‑90 (MQ=255)
                                                      gCTTCATTTTTGATTGGCGGGTTTGGTTTTTTCTGCTGCATCTGCTCGGTGAAACCGCCGCCCTGAATCATAAAGCCAGGAATGACGCgg                              <  1:806361/90‑1 (MQ=255)
                                                      gCTTCATTTTTGATTGGCGGGTTTGGTTTTTTCTGCTGCATCTGCTCGGTGAAACCGCCGCCCTGAATCATAAAGCCAGGAATGACGCgg                              <  1:27219/90‑1 (MQ=255)
                                                       cTTCATTTTTGATTGGCGGGGTTGGGTTTTTTTTCTGCCTCTTCTCGGTGGAAACGCCCCCCCGAAACATAAAACCCGGAAAGACGCggg                             >  1:292890/1‑89 (MQ=255)
                                                       cTTCATTTTTGATTGGCGGGGTTGGGTTTTTTTTCTGCCTCTGCTCGGTGGAAACCCCCCCCCGAATCATAAAACCAGGAATGACGCggg                             >  1:59129/1‑89 (MQ=255)
                                                        ttCATTTTTGATTGGCGGGTTTGGTTTTTTCTGCTGCATCTGCTCGGTGAAACCGCCGCCCTGAATCATAAAGCCAGGAATGACGCGGTg                            <  2:448772/90‑1 (MQ=255)
                                                         tCATTTTTGATTGGCGGGTTTGGTTTTTTCTGCTGCATCTGCTCGGTGAAACCGCCGCCCTGAATCATAAAGCCAGGAATGACGCGGTGa                           <  1:381148/90‑1 (MQ=255)
                                                         tCATTTTTGATTGGCGGGTTTGGTTTTTTCTGCTGCATCTGCTCGGTGAAACCGCCGCCCTGAATCATAAAGCCAGGAATGACGCGGTGa                           <  1:409726/90‑1 (MQ=255)
                                                         tCATTTTTGATTGGCGGGGTTGGTTTTTTTTTCTGCATCTGCTCGGGGAAAACGCCGCCCTGAATCATAAAAACAGGGATGACGCggggg                           >  2:409726/1‑87 (MQ=255)
                                                                                   ttttCTGCTGCATCTGCTCGGTGAAACCGCCGCCCTGAATCATAAAGCCAGGAATGACGCGGTGAAAGGTAGTGTTGTTATAAAAACCGc  <  2:599189/90‑1 (MQ=255)
                                                                                       |                                                                                      
TCAGCGGTACGTGCCATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGCGGGTTTGGTTTTTTCTGCTGCATCTGCTCGGTGAAACCGCCGCCCTGAATCATAAAGCCAGGAATGACGCGGTGAAAGGTAGTGTTGTTATAAAAACCGC  >  NZ_CP009273/3485301‑3485472

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