Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi3,981,1020TA20.0% 116.7 / 27.3 50V72V (GTT→GTA) yifKputative APC family amino acid transporter
Reads supporting (aligned to +/- strand):  ref base T (22/18);  new base A (5/5);  total (27/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.52e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTGGTTTATGTGGATGGCGGTGGGGATCTCTGAAATCACCGCCATTGGCGTTTATGTCCAGTTCTGGTTCCCGGAGAT  >  NC_000913_3_bme_pgi/3980957‑3981251
                                                                                                                                                 |                                                                                                                                                     
gCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAaccggttaccg                                                                                                                                                   <  2:616026/150‑1 (MQ=255)
  aCCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAaccggttaccggt                                                                                                                                                 <  1:146182/150‑1 (MQ=255)
   ccATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAaccggttaccggtt                                                                                                                                                <  4:319476/150‑1 (MQ=255)
                 cTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTAt                                                                                                                                  >  3:346909/1‑150 (MQ=255)
                     ttATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATgcgc                                                                                                                              >  1:667381/1‑150 (MQ=255)
                      tATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCa                                                                                                                             >  2:59139/1‑150 (MQ=255)
                           gggCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGtt                                                                                                                        <  2:44725/150‑1 (MQ=255)
                                   aGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGc                                                                                                                <  4:346909/150‑1 (MQ=255)
                                        ccTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTAACCGgttcgagga                                                                                                                                          >  4:216160/1‑115 (MQ=255)
                                        ccTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTAACCGgttcgagga                                                                                                                                          <  3:216160/119‑5 (MQ=255)
                                          tGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCt                                                                                                         >  1:13227/1‑150 (MQ=255)
                                          tGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTACTCGACCCGGTTACCGGTTCGTTCGCAGTTTATGAGCATCGTTATATGAGAACGTTCt                                                                                                         <  2:531545/150‑1 (MQ=255)
                                               tGGGCCGGGCCTTCCGGATTGTTGGCCTATATCCTCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGc                                                                                                    <  2:11774/150‑1 (MQ=255)
                                                  gccgggccATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTAt                                                                                                 >  4:417853/1‑150 (MQ=255)
                                                      ggccATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCa                                                                                             >  2:506389/1‑150 (MQ=255)
                                                       gccATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCAc                                                                                            <  1:506389/150‑1 (MQ=255)
                                                             cGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTgg                                                                                                     >  1:346864/1‑135 (MQ=255)
                                                             cGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTgg                                                                                                     <  2:346864/135‑1 (MQ=255)
                                                               tATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTAACCGgttcga                                                                                                                                             <  1:376938/93‑2 (MQ=255)
                                                               tATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTAACCGgttcga                                                                                                                                             >  2:376938/1‑92 (MQ=255)
                                                                    ttgGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTAACCGgttcga                                                                                                                                             <  1:529865/88‑2 (MQ=255)
                                                                    ttgGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTAACCGgttcga                                                                                                                                             >  2:529865/1‑87 (MQ=255)
                                                                       gCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTgg                                                                            >  3:386735/1‑150 (MQ=255)
                                                                         cTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTgg                                                                                                     <  3:93751/123‑1 (MQ=255)
                                                                         cTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTgg                                                                                                     >  4:93751/1‑123 (MQ=255)
                                                                              tcatcGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTGGTTTAtgt                                                                     <  1:59139/150‑1 (MQ=255)
                                                                               catcGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTGGt                                                                           >  3:361639/1‑143 (MQ=255)
                                                                               catcGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTGGt                                                                           <  4:361639/143‑1 (MQ=255)
                                                                                atcGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTGGTTTATGTgg                                                                   >  3:365148/1‑150 (MQ=255)
                                                                                 tcGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCACCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTg                                                                                                      <  1:351757/114‑1 (MQ=255)
                                                                                 tcGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCACCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTg                                                                                                      >  2:351757/1‑114 (MQ=255)
                                                                                     cGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTGGTTTATGTGGATGGc                                                              >  3:180546/1‑150 (MQ=255)
                                                                                          tGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTGGTTTATGTGGATGGCGGTgg                                                         <  4:188162/150‑1 (MQ=255)
                                                                                          tGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTGGTTTATGTGGATGGCGGTgg                                                         >  1:257078/1‑150 (MQ=255)
                                                                                            ttCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTGGTTTAtgt                                                                     <  2:208811/136‑1 (MQ=255)
                                                                                            ttCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTGGTTTAtgt                                                                     >  1:208811/1‑136 (MQ=255)
                                                                                               gTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTa                                                                                >  4:223199/1‑122 (MQ=255)
                                                                                               gTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTa                                                                                <  3:223199/122‑1 (MQ=255)
                                                                                                tCTTTTTCATGATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTGGTTTATGTGGATGGCGGTGGGGAtct                                                   >  4:58552/1‑150 (MQ=255)
                                                                                                                              tgttgtTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTGGTTTATGTGGATGGCGGTGGGGATCTCTGAAATCACCGCCATTGGCGTTTa                          <  2:370961/145‑1 (MQ=255)
                                                                                                                              tgttgtTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTGGTTTATGTGGATGGCGGTGGGGATCTCTGAAATCACCGCCATTGGCGTTTa                          >  1:370961/1‑145 (MQ=255)
                                                                                                                                     gaacGAACCGGTAACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTGGTTTATGTGGATGGCGGTGGGGATCTCTGAAATCACCGCCATTGGCGTTTATGTCCAGTTCTg              <  3:61947/147‑1 (MQ=255)
                                                                                                                                       tCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTGGTTTATGTGGATGGCGGTg                                                          >  2:128487/1‑104 (MQ=255)
                                                                                                                                       tCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTGGTTTATGTGGATGGCGGTg                                                          <  1:128487/104‑1 (MQ=255)
                                                                                                                                       tCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTGGTTTATGTGGATGGCGGTGGGGATCTCCGAAATCACCGCCATTGGCGTTTAt                         <  2:165029/137‑1 (MQ=255)
                                                                                                                                       tCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTGGTTTATGTGGATGGCGGTGGGGATCTCCGAAATCACCGCCATTGGCGTTTAt                         >  1:165029/1‑137 (MQ=255)
                                                                                                                                       tCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGCTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTGGTTTATGTGGATGGCGGTGGGGATCTCTGAAATCACCGCCATTGGCGTTTATGTCCAGTTCTGGt            >  1:443948/1‑150 (MQ=255)
                                                                                                                                        cGAACCGGTAACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTAt                                                                                                 >  2:63588/1‑64 (MQ=255)
                                                                                                                                        cGAACCGGTAACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTAt                                                                                                 <  1:63588/64‑1 (MQ=255)
                                                                                                                                         gAACCGGTAACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTGGTTTATGTGGATGGCGGTGGGGATCTCTGAAATCACCGCCATTGGCGTTTATGTCCAGTTCTGGTTc          >  1:455495/1‑150 (MQ=255)
                                                                                                                                          aaCCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTGGTTTATGTGGATgg                                                               <  1:341609/96‑1 (MQ=255)
                                                                                                                                          aaCCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTGGTTTATGTGGATgg                                                               >  2:341609/1‑96 (MQ=255)
                                                                                                                                          aaCCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTGGTTTATGTGGATGGCGGTGGGGATCTCTGAAATCACCGCCATTGGCGTTTATGTCCAGTTCTGGTTcc         <  2:257078/150‑1 (MQ=255)
                                                                                                                                           accggttaccggttCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTGGTTTATGTGGATGGCGGTGGGGATCTCTGAAATCACCGCCATTGGCGTTTATGTCCAGTTCTGGTTccc        >  4:22353/1‑150 (MQ=255)
                                                                                                                                            ccggttaccggttCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTGGTTTATGTGGATGGCGGTGGGGATCTCTGAAATCACCGCCATTGGCGTTTATGTCCAGTTCTGGTTCCCg       <  2:222/150‑1 (MQ=255)
                                                                                                                                               gttaccggttCGTTCGCCGTTTATGCGCATCGTTATATGAGCTCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTGGTTTATGTGGATGGCGGTGGGGATCTCTGAAATCACCGCCATTGGCGCTTATGTCCAGTTCTGGTTGCCGgag    >  2:293595/1‑150 (MQ=255)
                                                                                                                                                 taccggttCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTGGTTTATGTGGATGGCGGTGGGGATCTCTGAAATCACCGCCATTGGCGTTTATGTCCAGTTCTGGTTCCCGGAGAt  <  2:13227/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
GCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTGGTTTATGTGGATGGCGGTGGGGATCTCTGAAATCACCGCCATTGGCGTTTATGTCCAGTTCTGGTTCCCGGAGAT  >  NC_000913_3_bme_pgi/3980957‑3981251

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: ‑