Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi3,208,2100TC20.4% 113.1 / 20.7 49S63P (TCC→CCC) ttdTL‑tartrate/succinate antiporte
Reads supporting (aligned to +/- strand):  ref base T (17/22);  new base C (5/5);  total (22/27)
Fisher's exact test for biased strand distribution p-value = 7.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.31e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CTGGCGGTCATCGCCATTATTGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAATTCGGTTATCTGGCTGATTTTCGCCGCCTTTATGTTT  >  NC_000913_3_bme_pgi/3208063‑3208359
                                                                                                                                                   |                                                                                                                                                     
cTGGCGGTCATCGCCATTATTGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATcc                                                                                                                                                     <  2:359013/150‑1 (MQ=255)
 tGGCGGTCATCGCCATTATTGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCca                                                                                                                                                    <  2:198784/150‑1 (MQ=255)
     ggTCATCGCCATTATTGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCcatcat                                                                                                                                                >  2:96648/1‑150 (MQ=255)
              cATTATTGCTCTACTTACGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATAct                                                                                                                                       >  2:362389/1‑150 (MQ=255)
                ttattGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATActct                                                                                                                                     <  2:619378/150‑1 (MQ=255)
                       tctACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTg                                                                                                                              <  1:678990/150‑1 (MQ=255)
                           cTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGgctg                                                                                                                          <  1:222144/150‑1 (MQ=255)
                           cTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGgctg                                                                                                                          <  2:568104/150‑1 (MQ=255)
                            ttCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGgctgc                                                                                                                         <  1:10552/150‑1 (MQ=255)
                            ttCCGGTTCCCGCAGGCCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGGTGGTGGGTATATCCATCCTCGCCATACTCTCTCCCTGgctgc                                                                                                                         >  3:388372/1‑150 (MQ=255)
                               cGGTTCCCGCAGGGCTGGCGAACCATACCTGGCCCTACTTTGCCGTTTTTACTGGCGTGGTCGTTGGACTGATCTTCGCACCCGTGCAGGGTGCAGTCGTGGCGATGGTGGGTAAATCCATCATCGCCCTACTCCGTCCCTGGCTGCAGt                                                                                                                      >  4:167858/1‑150 (MQ=255)
                                  ttCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATCCTCTCTCCCTGGCTGCTGTTCa                                                                                                                   <  4:312002/150‑1 (MQ=255)
                                      cGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCcat                                                                                                                                                   <  1:578909/114‑1 (MQ=255)
                                      cGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCcat                                                                                                                                                   >  2:578909/1‑114 (MQ=255)
                                      cGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGccc                                                                                                               >  1:278597/1‑150 (MQ=255)
                                       gCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCg                                                                                                              <  3:388582/150‑1 (MQ=255)
                                              tGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGagcagc                                                                                                       <  2:703258/150‑1 (MQ=255)
                                                          ccTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCTATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAg                                                                                           <  1:409185/150‑1 (MQ=255)
                                                              cctctAGTTTGCCTTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCCGGCtt                                                                                       <  4:348027/149‑1 (MQ=255)
                                                                tctACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTa                                                                                     >  3:21073/1‑150 (MQ=255)
                                                                      ttGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCa                                                                               <  4:369453/150‑1 (MQ=255)
                                                                      ttGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCa                                                                               <  2:278597/150‑1 (MQ=255)
                                                                        gCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACt                                                                             >  2:709779/1‑150 (MQ=255)
                                                                        gCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACt                                                                             <  4:349796/150‑1 (MQ=255)
                                                                                   tGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGATGGATATACCCACCATCGc                                                                                                                                             >  1:280447/1‑75 (MQ=25)
                                                                                   tGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGATGGATATACCCACCATCGc                                                                                                                                             <  2:280447/75‑1 (MQ=25)
                                                                                                  aCTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTc                                                   >  4:24805/1‑150 (MQ=255)
                                                                                                        ccTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGttt                                             <  3:24805/150‑1 (MQ=255)
                                                                                                          tCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGttttt                                           >  4:161700/1‑150 (MQ=255)
                                                                                                                ccGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCg                                                                                                              <  1:473572/77‑1 (MQ=255)
                                                                                                                ccGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCAGACTCTCTCCCTGGCTGCTGTTCAGCCCg                                                                                                              >  2:473572/1‑77 (MQ=255)
                                                                                                                ccGTGCCGGGTGCCGTGGTGGCGATGATGGATATACCcaccat                                                                                                                                                <  3:322623/43‑1 (MQ=14)
                                                                                                                ccGTGCCGGGTGCCGTGGTGGCGATGATGGATATACCcaccat                                                                                                                                                >  4:322623/1‑43 (MQ=14)
                                                                                                                ccGTGCCGGGTGCCGTGGTGGCGATGATGGATATACCcaccat                                                                                                                                                >  4:322622/1‑43 (MQ=14)
                                                                                                                ccGTGCCGGGTGCCGTGGTGGCGATGATGGATATACCcaccat                                                                                                                                                >  4:322619/1‑43 (MQ=2)
                                                                                                                ccGTGCCGGGTGCCGTGGTGGCGATGATGGATATACCcaccat                                                                                                                                                <  3:322619/43‑1 (MQ=14)
                                                                                                                ccGTGCCGGGTGCCGTGGTGGCGATGATGGATATACCcaccat                                                                                                                                                <  3:322622/43‑1 (MQ=14)
                                                                                                                      cGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAATTCGGTTa                               >  1:618869/1‑150 (MQ=255)
                                                                                                                         gtgtCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAATTCGGTTATCt                            <  4:240119/150‑1 (MQ=255)
                                                                                                                         gTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGagca                                                                                                         >  4:210616/1‑73 (MQ=255)
                                                                                                                         gTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGagca                                                                                                         <  3:210616/73‑1 (MQ=255)
                                                                                                                           gCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATCCTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAATTCGGTTATCTgg                          <  2:428525/150‑1 (MQ=255)
                                                                                                                             cGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAATTCGGTTATCTGGCt                        >  3:310025/1‑150 (MQ=255)
                                                                                                                              gtggtgGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAATTCGGTTATCTGGCTg                       >  4:161053/1‑150 (MQ=255)
                                                                                                                                 gtggCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAATTCGGTTATCTGGCTGAtt                    >  1:488786/1‑150 (MQ=255)
                                                                                                                                   ggCGATGATGGATATACCCACCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAATTCGGTTATCTGGCt                        >  1:9683/1‑144 (MQ=255)
                                                                                                                                   ggCGATGATGGATATACCCACCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAATTCGGTTATCTGGCt                        <  2:9683/144‑1 (MQ=255)
                                                                                                                                      gATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGcc                                                        >  4:207561/1‑109 (MQ=255)
                                                                                                                                      gATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGcc                                                        <  3:207561/109‑1 (MQ=255)
                                                                                                                                                   tCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAATTCGGTTATCTGGCTGATTTTCGCCGCCTTTATGttt  <  4:21073/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CTGGCGGTCATCGCCATTATTGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAATTCGGTTATCTGGCTGATTTTCGCCGCCTTTATGTTT  >  NC_000913_3_bme_pgi/3208063‑3208359

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