Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,798,9590GA53.6% ‑4.6 / 39.8 28D397D (GAC→GAT) icdisocitrate dehydrogenase (NADP(+))
Reads supporting (aligned to +/- strand):  ref base G (5/8);  new base A (8/7);  total (13/15)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.71e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TCATTTAAAAAGGCGAACCTAGCAATTAAAGTTGCTAGACTCACCTAATTTTATTTTGCTTATTTTAAATTTTTAATCAATTCATCTGCAAATGCTGATGTAGAAACTTCTTCAGCACCATCCATTAAACGGGCAAAGTCATAAGTAACAACTTTTGAAGCAATTGTATCTTCAATTGAATCTGTAATCTTATCTGCCGCTTCTTGCCATCCTAAATGTTCTAACATTAATACAGAACTTAAAATTACTGAAGATGGATTCAC  >  CP000730/1798822‑1799084
                                                                                                                                         |                                                                                                                             
tCATTTAAAAAGGCGAACCTAGCAATTAAAGTTGCTAGACTCACCTAATTTTATTTTGCTTATTTTAAATTTTTAATCAATTCATCTGCAAATGCTGATGTAGAAACTTCTTCAGCACCATCCATTAAACGGGCAAAATCATAAGTaa                                                                                                                     >  1:347994/1‑148 (MQ=255)
     tAAAAAGGCGAACCTAGCAATTAAAGTTGCTAGACTCACCTAATTTTATTTTGCTTATTTTAAATTTTTAATCAATTCATCTGCAAATGCTGATGTAGAAACTTCTTCAGCACCATCCATTAAACGGGCAAAGTCATAAGTAACAACt                                                                                                                >  1:90533/1‑148 (MQ=255)
            gCGAACCTAGCAATTAAAGTTGCTAGACTCACCTAATTTTATTTTGCTTATTTTAAATTTTTAATCAATTCATCTGCAAATGCTGATGTAGAAACTTCTTCAGCACCATCCATTAAACGGGCAAAATCATAAGTAACAACTTTTGAAg                                                                                                         >  2:241961/1‑148 (MQ=255)
                       aaTTAAAGTTGCTAGACTCACCTAATTTTATTTTGCTTATTTTAAATTTTTAATCAATTCATCTGCAAATGCTGATGTAGAAACTTCTTCAGCACCATCCATTAAACGGGCAAAATCATAAGTAACAACTTTTGAAGCAATTGTATCt                                                                                              <  1:226042/148‑1 (MQ=255)
                           aaaGTTGCTAGACTCACCTAATTTTATTTTGCTTATTTTAAATTTTTAATCAATTCATCTGCAAATGCTGATGTAGAAACTTCTTCAGCACCATCCATTAAACGGGCAAAGTCATAAGTAACAACTTTTGAAGCAATTGTATCTTCaa                                                                                          <  2:76664/148‑1 (MQ=255)
                           aaaGTTGCTAGACTCACCTAATTTTATTTTGCTTATTTTAAATTTTTAATCAATTCATCTGCAAATGCTGATGTAGAAACTTCTTCAGCACCATCCATTAAACGGGCAAAGTCATAAGTAACAACTTTTGAAGCAATTGTATCTTCaa                                                                                          <  2:192062/148‑1 (MQ=255)
                                  cTAGACTCACCTAATTTTATTTTGCTTATTTTAAATTTTTAATCAATTCATCTGCAAATGCTGATGTAGAAACTTCTTCAGCACCATCCATTAAACGGGCAAAATCATAAGTAACAACTTTTGAAGCAATTGTATCTTCAATTGAATc                                                                                   <  2:66067/148‑1 (MQ=255)
                                                  ttattttGCTTATTTTAAATTTTTAATCAATTCATCTGCAAATGCTGATGTAGAAACTTCTTCAGCACCATCCATTAAACGGGCAAAATCATAAGTAACAACTTTTGAAGCAATTGTATCTTCAATTGAATCTGTAATCTTATCTGcc                                                                   >  2:268564/1‑148 (MQ=255)
                                                  ttattttGCTTATTTTAAATTTTTAATCAATTCATCTGCAAATGCTGATGTAGAAACTTCTTCAGCACCATCCATTAAACGGGCAAAATCATAAGTAACAACTTTTGAAGCAATTGTATCTTCAATTGAATCTGTAATCTTATCTGcc                                                                   >  2:325531/1‑148 (MQ=255)
                                                         gCTTATTTTAAATTTTTAATCAATTCATCTGCAAATGCTGATGTAGAAACTTCTTCAGCACCATCCATTAAACGGGCAAAGTCATAAGTAACAACTTTTGAAGCAATTGTATCTTCAATTGAATCTGTAATCTTATCTGCCGcttctt                                                            <  2:90533/148‑1 (MQ=255)
                                                              ttttAAATTTTTAATCAATTCATCTGCAAATGCTGATGTAGAAACTTCTTCAGCACCATCCATTAAACGGGCAAAGTCATAAGTAACAACTTTTGAAGCAATTGTATCTTCAATTGAATCTGTAATCTTATCTGCCGCTTCTTGCCAt                                                       <  1:254887/148‑1 (MQ=255)
                                                              ttttAAATTTTTAATCAATTCATCTGCAAATGCTGATGTAGAAACTTCTTCAGCACCATCCATTAAACGGGCAAAGTCATAAGTAACAACTTTTGAAGCAATTGTATCTTCAATTGAATCTGTAATCTTATCTGCCGCTTCTTGCCAt                                                       <  1:346656/148‑1 (MQ=255)
                                                                ttAAATTTTTAATCAATTCATCTGCAAATGCTGATGTAGAAACTTCTTCAGCACCATCCATTAAACGGGCAAAATCATAAGTAACAACTTTTGAAGCAATTGTATCTTCAATTGAATCTGTAATCTTATCTGCCGCTTCTTGCCATcc                                                     >  1:50984/1‑148 (MQ=255)
                                                                        ttAATCAATTCATCTGCAAATGCTGATGTAGAAACTTCTTCAGCACCATCCATTAAACGGGCAAAGTCATAAGTAACAACTTTTGAAGCAATTGTATCTTCAATTGAATCTGTAATCTTATCTGCCGCTTCTTGCCATCCTAAATGtt                                             <  1:116796/148‑1 (MQ=255)
                                                                              aaTTCATCTGCAAATGCTGATGTAGAAACTTCTTCAGCACCATCCATTAAACGGGCAAAGTCATAAGTAACAACTTTTGAAGCAATTGTATCTTCAATTGAATCTGTAATCTTATCTGCCGCTTCTTGCCATCCTAAATGTTCTAACa                                       <  2:117749/148‑1 (MQ=255)
                                                                              aaTTCATCTGCAAATGCTGATGTAGAAACTTCTTCAGCACCATCCATTAAACGGGCAAAGTCATAAGTAACAACTTTTGAAGCAATTGTATCTTCAATTGAATCTGTAATCTTATCTGCCGCTTCTTGCCATCCTAAATGTTCTAACa                                       <  2:69692/148‑1 (MQ=255)
                                                                                   aactacAAATGCTGATGGAGAAACTTCTTCAGCACGATCCATTAAACGGGCAAAATCATAAGTAACAACTTTTGAAGCAATTGTATCTTCAATTGAATCTGTAATCTTATCTGCCGCTTCTTGCCATCCTAAATGTTCTAACATTAAt                                  <  1:35292/143‑1 (MQ=255)
                                                                                   aTCTGCAAATGCTGATGTAGAAACTTCTTCAGCACCATCCATTAAACGGGCAAAATCATAAGTAACAACTTTTGAAGCAATTGTATCTTCAATTGAATCTGTAATCTTATCTGCCGCTTCTTGCCATCCTAAATGTTCTAACATTAAt                                  <  1:137225/148‑1 (MQ=255)
                                                                                   aTCTGCAAATGCTGATGTAGAAACTTCTTCAGCACCATCCATTAAACGGGCAAAATCATAAGTAACAACTTTTGAAGCAATTGTATCTTCAATTGAATCTGTAATCTTATCTGCCGCTTCTTGCCATCCTAAATGTTCTAACATTAAt                                  <  2:382662/148‑1 (MQ=255)
                                                                                    tCTGCAAATGCTGATGTAGAAACTTCTTCAGCACCATCCATTAAACGGGCAAAATCATAAGTAACAACTTTTGAAGCAATTGTATCTTCAATTGAATCTGTAATCTTATCTGCCGCTTCTTGCCATCCTAAATGTTCTAACATTAATa                                 <  2:299896/148‑1 (MQ=255)
                                                                                      tGCAAATGCTGATGTAGAAACTTCTTCAGCACCATCCATTAAACGGGCAAAGTCATAAGTAACAACTTTTGAAGCAATTGTATCTTCAATTGAATCTGTAATCTTATCTGCCGCTTCTTGCCATCCTAAATGTTCTAACATTAATACa                               >  2:307952/1‑148 (MQ=255)
                                                                                      tGCAAATGCTGATGTAGAAACTTCTTCAGCACCATCCATTAAACGGGCAAAGTCATAAGTAACAACTTTTGAAGCAATTGTATCTTCAATTGAATCTGTAATCTTATCTGCCGCTTCTTGCCATCCTAAATGTTCTAACATTAATACa                               >  2:326389/1‑148 (MQ=255)
                                                                                                         aCTTCTTCAGCACCATCCATTAAACGGGCAAAATCATAAGTAACAACTTTTGAAGCAATTGTATCTTCAATTGAATCTGTAATCTTATCTGCCGCTTCTTGCCATCCTAAATGTTCTAACATTAATACAGAACTTAAAATTACTGAAg            >  2:276444/1‑148 (MQ=255)
                                                                                                         aCTTCTTCAGCACCATCCATTAAACGGGCAAAATCATAAGTAACAACTTTTGAAGCAATTGTATCTTCAATTGAATCTGTAATCTTATCTGCCGCTTCTTGCCATCCTAAATGTTCTAACATTAATACAGAACTTAAAATTACTGAAg            >  2:42981/1‑148 (MQ=255)
                                                                                                         aCTTCTTCAGCACCATCCATTAAACGGGCAAAATCATAAGTAACAACTTTTGAAGCAATTGTATCTTCAATTGAATCTGTAATCTTATCTGCCGCTTCTTGCCATCCTAAATGTTCTAACATTAATACAGAACTTAAAATTACTGAAg            >  2:24989/1‑148 (MQ=255)
                                                                                                             cttcAGCACCATCCATTAAACGGGCAAAGTCATAAGTAACAACTTTTGAAGCAATTGTATCTTCAATTGAATCTGTAATCTTATCTGCCGCTTCTTGCCATCCTAAATGTTCTAACATTAATACAGAACTTAAAATTACTGAAGATgg        >  1:43844/1‑148 (MQ=255)
                                                                                                             cttcAGCACCATCCATTAAACGGGCAAAGTCATAAGTAACAACTTTTGAAGCAATTGTATCTTCAATTGAATCTGTAATCTTATCTGCCGCTTCTTGCCATCCTAAATGTTCTAACATTAATACAGAACTTAAAATTACTGAAGATgg        >  1:90806/1‑148 (MQ=255)
                                                                                                                   cACCATCCATTAAACGGGCAAAATCATAAGTAACAACTTTTGAAGCAATTGTATCTTCAATTGAATCTGTAATCTTATCTGCCGCTTCTTGCCATCCTAAATGTTCTAACATTAATACAGAACTTAAAATTACTGAAGATGGATTCAc  <  1:241961/148‑1 (MQ=255)
                                                                                                                                         |                                                                                                                             
TCATTTAAAAAGGCGAACCTAGCAATTAAAGTTGCTAGACTCACCTAATTTTATTTTGCTTATTTTAAATTTTTAATCAATTCATCTGCAAATGCTGATGTAGAAACTTCTTCAGCACCATCCATTAAACGGGCAAAGTCATAAGTAACAACTTTTGAAGCAATTGTATCTTCAATTGAATCTGTAATCTTATCTGCCGCTTCTTGCCATCCTAAATGTTCTAACATTAATACAGAACTTAAAATTACTGAAGATGGATTCAC  >  CP000730/1798822‑1799084

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