Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,315,9360TA23.6% 131.3 / 13.0 55K34* (AAA→TAA) yciGprotein YciG
Reads supporting (aligned to +/- strand):  ref base T (37/5);  new base A (0/13);  total (37/18)
Fisher's exact test for biased strand distribution p-value = 5.90e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCCGTTGCGTTGCATTGTTTCATTAATA  >  NC_000913/1315802‑1316072
                                                                                                                                        |                                                                                                                                         
aCGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAAt                                                                                                                                     >  2:435662/1‑140 (MQ=255)
   ccccGAAGTTGAATTTATTTTTATCAGCTAATTTTAATAATTAACGGTAAATAAGGTTTTGCTTTTTTTATTCCCCCCGTTTTGTTGCCCGCTTTTTTTCCCCGCTTCAGATGCGCGTTGCGGATCTTTTTTAAAATTAc                                                                                                                                  <  2:469932/139‑1 (MQ=255)
      cGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTAccgc                                                                                                                               >  1:577776/1‑140 (MQ=255)
       gCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTAccgcc                                                                                                                              >  1:844605/1‑140 (MQ=255)
        cAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTAccgccg                                                                                                                             >  2:774269/1‑140 (MQ=255)
                   aTCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCg                                                                                                                 >  1:943110/1‑140 (MQ=255)
                     cTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGc                                                                                                                >  1:975635/1‑140 (MQ=255)
                          tCAGCTGATTTTAATATTAAACGGTAAATCAGGTTTTGCTTGTTTTATTCCCCCCGTTTTGTTGCCCGCTTTTTTTCCCCGCTTCAGATGCGCGTTGCGGATCTTTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTa                                                                                                           <  2:844605/140‑1 (MQ=255)
                              cTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCACCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGc                                                                                                       >  2:791710/1‑140 (MQ=255)
                                    ttAATATTTAGGGGTGAACAGGTTTTTCCCTTTTTTTTTCCCCCCGTTTTTTTGCCCCCTTTTTTTCCCCGCTCCAGAGGCGGGTTGCGGTTTTTTTTAAAAATTCCCGCCGCTATGCTGACCGCCTTTACGGCCTGCGt                                                                                                 <  2:975635/140‑1 (MQ=255)
                                    ttAATAATTACCGGTGAATCAGGTTTTGCCTGTTTTATTCCCACCGTTTTGTTGCCCGCTTTTTTTCCCCGCTTCAGATGCGCGTTGGGGATCTTTTTAAAAATTCCCGCCGCTATGCTGACCGCCTTTACGGCCTGCGt                                                                                                 <  2:622896/140‑1 (MQ=255)
                                        taatTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGa                                                                                             >  1:316008/1‑140 (MQ=255)
                                        taatTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGa                                                                                             >  2:522348/1‑140 (MQ=255)
                                           ttAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGc                                                                                          >  1:622896/1‑140 (MQ=255)
                                            tAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGcc                                                                                         >  2:726079/1‑140 (MQ=255)
                                                    aaTCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACg                                                                                 >  2:410855/1‑140 (MQ=255)
                                                                 ccTGTTTTATTCCCCCCGTTTTGTTGCCCGCTTTTTTTCCCCGCTTCAGATGCGCGTTGCGGATCATTTTAAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGaa                                                                     <  2:577776/140‑1 (MQ=255)
                                                                 ccTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGaa                                                                     >  1:584532/1‑140 (MQ=255)
                                                                     aTTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAAttt                                                                 >  2:340765/1‑140 (MQ=255)
                                                                     aTTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAAttt                                                                 >  1:449999/1‑140 (MQ=255)
                                                                           taCCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGaa                                                           >  2:564355/1‑140 (MQ=255)
                                                                               accGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAaccac                                                       >  1:680419/1‑140 (MQ=255)
                                                                               accGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAaccac                                                       >  1:437104/1‑140 (MQ=255)
                                                                                ccGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAaccacc                                                      >  1:518799/1‑140 (MQ=255)
                                                                                ccGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCATTCTCACGGTCTTCGGCGAAATTTCCTGAaccacc                                                      >  2:702311/1‑140 (MQ=255)
                                                                                 cGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAAccacca                                                     >  2:772648/1‑140 (MQ=255)
                                                                                   ttttGTTGCCCCCTTTTTTCCCCCCCTTCAGAGGCGCGTTGCGGATCTTTTTAAAAATTCCCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACg                                                   <  2:210517/139‑1 (MQ=255)
                                                                                   cTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACg                                                   >  1:698715/1‑140 (MQ=255)
                                                                                   cTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACg                                                   >  1:684134/1‑140 (MQ=255)
                                                                                    tttGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGa                                                  >  1:606006/1‑140 (MQ=255)
                                                                                           gcccccTTTTTTCCCCCTTTTCAGAGGCGGTTTGGGGATTTTTTTAAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCg                                            <  2:698715/136‑1 (MQ=255)
                                                                                           aCCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATgctgtc                                           >  1:961541/1‑135 (MQ=255)
                                                                                            cccctttttttCCCCGCTTCAGATGCGCGTTGCGGATCTTTTTAAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGc                                          <  1:774269/135‑1 (MQ=255)
                                                                                            cccccTTTTTTCCCCGCTTCAGAGGCGCTTGGCGGATAATTTTAAAATTTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGc                                          <  1:591439/137‑1 (MQ=255)
                                                                                             cGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGcc                                         >  2:8226/1‑140 (MQ=255)
                                                                                             cGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGcc                                         >  2:977553/1‑140 (MQ=255)
                                                                                             cGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGCCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGcc                                         >  2:583710/1‑140 (MQ=255)
                                                                                              gCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCa                                        >  1:155491/1‑140 (MQ=255)
                                                                                              gCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCa                                        >  1:56237/1‑140 (MQ=255)
                                                                                                            cTTCAGAGGCGGGTGGGGGTTCTTTTTAAAATTTCCCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCcg                          <  1:466402/140‑1 (MQ=255)
                                                                                                             ctccAGATGCGCTTTGCGGATCTTTTTAAAAATTCCCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCcg                          <  2:606006/137‑1 (MQ=255)
                                                                                                              ttCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCcgt                         >  1:923869/1‑140 (MQ=255)
                                                                                                                      cgctTTGCGAATCTTTTTAAAATTTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCCgttgcgttgc                 <  2:584532/140‑1 (MQ=255)
                                                                                                                       gcgTTGGGGATTTTTTTAAAATTTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCCGTTGCGTTGCa                <  2:437104/140‑1 (MQ=255)
                                                                                                                        cgTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCcgttgcgttg                  >  2:930164/1‑137 (MQ=255)
                                                                                                                        cgTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCCGTTGCGTTGCAt               >  2:362376/1‑140 (MQ=255)
                                                                                                                        cgTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCCGTTGCGTTGCAt               >  1:301568/1‑140 (MQ=255)
                                                                                                                        cgTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCCGTTGCGTTGCAt               >  2:878691/1‑140 (MQ=255)
                                                                                                                        cgTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCCGTTGCGTTGCAt               >  1:339924/1‑140 (MQ=255)
                                                                                                                        cgTTGCGGATCATTTTAAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCcgttgcgttg                  <  1:930164/137‑1 (MQ=255)
                                                                                                                          ttGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCCGTTGCGTTGCATTg             <  1:878691/140‑1 (MQ=255)
                                                                                                                           tGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCCGTTGCGTTGCATTGt            >  1:891581/1‑140 (MQ=255)
                                                                                                                            gCGGATCTTTTTTAAATTTCCCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCCGTTGCGTTGCATTGtt           <  2:684134/140‑1 (MQ=255)
                                                                                                                            gCGAATCTTTTTAAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCCGTTGCGTTGCATTGtt           <  2:923869/140‑1 (MQ=255)
                                                                                                                               gattttttttAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCCGTTGCGTTGCATTGTTTCa        <  1:791710/135‑1 (MQ=255)
                                                                                                                                 tCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCCGTTGCGTTGCATTGTTTCAtt      <  1:772648/140‑1 (MQ=255)
                                                                                                                                   ttttttAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGCCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCCGTTGCGTTGCATTGTTTCATTaa    <  1:583710/139‑1 (MQ=255)
                                                                                                                                    tttttAAAATTTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCCGTTGCGTTGCATTGTTTCATTaa    <  1:435662/140‑1 (MQ=255)
                                                                                                                                     ttttAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCCGTTGCGTTGCATTGTTTCATTAata  <  1:362376/140‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                         
ACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCCGTTGCGTTGCATTGTTTCATTAATA  >  NC_000913/1315802‑1316072

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