Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,174,6220GA30.6% 53.1 / 11.7 36P182S (CCC→TCC) ycfTinner membrane protein YcfT
Reads supporting (aligned to +/- strand):  ref base G (16/9);  new base A (0/11);  total (16/20)
Fisher's exact test for biased strand distribution p-value = 5.53e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCATACCGAATCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCAGGTGCCAGCCACTGGTTCAGCGCACTTAACGCCAGCCACTGCAC  >  NC_000913/1174373‑1174868
                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                        
gcATACCGAATCGGTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGa                                                                                                                                                                                                                                                         <  1:511842/251‑1 (MQ=255)
    aCCGAATCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACg                                                                                                                                                                                                                                                     <  1:146955/250‑1 (MQ=255)
                           aaaaCAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACAACGGCGTGGGAACGAAATTAACCGCCACACTcagcagg                                                                                                                                                                                                                            >  2:20795/1‑250 (MQ=255)
                                         aTCACCACAATCGATACCAGAGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTCAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGACATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGa                                                                                                                                                                                                               >  2:537359/1‑250 (MQ=255)
                                              cacaATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTg                                                                                                                                                                                                         >  1:132846/1‑251 (MQ=255)
                                                 catctacaCCAGCGAAAACAGCAGCAAGATAGAAAACACCCAGGCCCCGCCCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCACCGGCCCCTCTTAAACACAGTTCATTATTGTTGCGCAAAACCATGCGCAAAGGCTGTAAAAAGGAAAATTGCGGATCACCCTGTTCATTCCCCACCACGGGGGGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCa                                                                                                                                                                                                      <  2:118731/244‑1 (MQ=255)
                                                      aTACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCa                                                                                                                                                                                                      >  1:118731/1‑246 (MQ=255)
                                                         ccAGCGAAAACAGCAGAGAAAAATAAAACAACCCGCCCCCGCCCGCCAGAACGGTAAGAAAAGAAGCCAGCAGAGGAGGGCGGCGCACCGCCCCCTCTTAAACACAGGTCATTATTGTCGCGCCAAACCAGGCGCCAGGGCTGAAATAAGGCAAATTGCGGATCACCCTGTCCACCCCCCCCCACGGCGTGGAAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGa                                                                                                                                                                                              <  2:6043/251‑1 (MQ=255)
                                                          cAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGAc                                                                                                                                                                                             >  1:293499/1‑251 (MQ=255)
                                                               acaacgGCAGCGAGAGATAAAACAACCAGGCCCCGCCCGCAAGACGGGTAAAAAAAGAACCCATCAGCAGAGGGCGGCGAAACGGCACCTCTTAAAACGTGGTCATTATTGTGGCGCCAACCCAGGCACCAGGGCTGTAATATGAAAAATTGCGGCCCACCCGGTTCTCTCCCCCCCACGGCGGTGAAAGAAAAGTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCGGACGGCTa                                                                                                                                                                                        <  2:345808/251‑1 (MQ=255)
                                                                aaacagaaggtataTAGAAAACACCCAGGCACCGCCCGCCAAAGAGGAAAACAAAGAAGCCCTCAGCGGATGGCGGCGCCACGGCACTGCTATAACACAGGTCATTATTGTTGCGCCAACCCGCGCGCCAAGGCTGTAATAGGCCAAAGTGCGAATCCCCCGGTTAATCCCCCCCCACGGCGTGGAAACGAATTTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTaa                                                                                                                                                                                       <  2:430321/238‑1 (MQ=255)
                                                                                     caacCAGGCCCCGACCGCCAAAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCCCCTCTTAAACACAGGTTATTATTGTTGCGCCAAACCATGCGCCAAGGCTGAAAAAAGGAAATTGGCGGATCACTCTGTTCATCCCCCCCCACGGCGGGGAAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTaaaaa                                                                                                                                                                                    <  1:522522/233‑1 (MQ=255)
                                                                                     caacCAGGCACCGACCGCCAGAACGGTACGCACAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTaaaaa                                                                                                                                                                                    >  2:522522/1‑233 (MQ=255)
                                                                                          aGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGACATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAACATTTTACATACCACGAAATAgcc                                                                                                                                                              >  1:467632/1‑248 (MQ=255)
                                                                                           gcccccgcccgCCAGACCGGTAAGAAAAGAAGCAATCAGCAGATGGGGGCACACCGGCCCCTCATAAACACAGTTATTTATGGTTGCGCCAACCCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGACCACCCTTTTCATCCCCCACCCCGGGGTGGAAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACaa                                                                                                                                                            <  2:186374/243‑1 (MQ=255)
                                                                                                         gAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTAATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAgccg                                                                                                                                                                                            >  2:551694/1‑202 (MQ=255)
                                                                                                         gAACGGTAAGAAAAAAACCCATCGGCAGATGGGGGCCCAACGGCACCTCTTAACCCCAGGTAATTTTTGGTGCGCAAAACAAGGCCCAAAGGTTTTAAAAAGGAAATGGGCGTATCACACTGTTAATCCCCCCCCAGGGCGTGGAAACAAAATAAGCCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACg                                                                                                                                                                                            <  1:551694/205‑1 (MQ=255)
                                                                                                               aaaaaaaaaacaccaccagaagctGGCGGAGCAACTGCTCCTTTTTAACACTGGTCATTTTTTGTGCGCCACACCATGCGCAAAGGCTGATATAAGGAATATTGCGGATCACCTTGTTTTTCCCCCACAACGTCGTGGGAACAAAAATAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGCTACcac                                                                                                                                        <  2:51852/228‑1 (MQ=255)
                                                                                                                    aaaGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTAATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCt                                                                                                                                   >  1:370718/1‑251 (MQ=255)
                                                                                                                                   agATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCACACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCCGTAAAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATCCCACAAGCTGGTGCTGGCGGTGat                                                                                                                    >  1:346192/1‑251 (MQ=255)
                                                                                                                                             ccAACGGCCCTTCTGTAACACAGGTCATTATTTGTGCGCCAAACCATCCACCAAGTCTGTAATAAGGCATATGGCGGCTCACACTTTCCATTCCCCCCCGCGGCGTGGGAACGAATTAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGc                                                                                                          <  1:540069/250‑1 (MQ=255)
                                                                                                                                                  ggCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGa                                                                                                      >  1:319144/1‑250 (MQ=255)
                                                                                                                                                    cACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGGTCAATCCGTGCAGGaaac                                                                                                   >  2:257113/1‑249 (MQ=255)
                                                                                                                                                     aCCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTc                                                                                                  >  1:440261/1‑251 (MQ=255)
                                                                                                                                                             accacaGGTCATTATGTGTGGCCAAACCCATGCCCCAAGGCTGTAAAAAGGAAAATTGCGGTCCACCCTTTTCATCCCCCCCCCCGGTGGGGAATCAAAATAAACCTCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGTg                                                                                           <  2:358173/248‑1 (MQ=255)
                                                                                                                                                                                      cAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGGGGAATCGGCATAAGCGGCATTgg                                                                     >  1:136130/1‑249 (MQ=255)
                                                                                                                                                                                       aacccATGCGCAAAGGCTGAAAAAAGACAAATTGCGGATCACACTGTTCATCCCCCACCACGGCGTGGAAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAgg                                                                  <  1:351393/248‑1 (MQ=255)
                                                                                                                                                                                         aCCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTg                                                                                                                                  >  2:404830/1‑185 (MQ=255)
                                                                                                                                                                                          cccaggcgcAAAGGCTGTAATAAGGAAAATGGCGAATCACCCTGTTCATCCCCCACCACGGCGTGGAAACAAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACATATACCACAAGCTg                                                                                                                                  <  1:404830/180‑1 (MQ=255)
                                                                                                                                                                                                    aagccgaaaataaGGCAAAGGGCGAATCACACTGTCCACCCCCCCCCAGGGGGTGGAAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAAt                                                      <  1:257113/243‑1 (MQ=255)
                                                                                                                                                                                                             ataaGGCAAATTGCGAAACCCCCTGCTCATCCCCCCCCGCGGGGTGGAAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCAg                                             <  2:129466/251‑1 (MQ=255)
                                                                                                                                                                                                              aaaGGAAAATTGCGGATCCCACTGTTCATTCCCCACCAGGGCGGGGGAACGAATTAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCAgg                                            <  2:132846/250‑1 (MQ=255)
                                                                                                                                                                                                                     aatggCGGAACACACTGTTCATCCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCAGGTgccagc                                     <  1:33625/247‑1 (MQ=255)
                                                                                                                                                                                                                              atccgccgttttatgccccACCCGGGCGTGGAACCAAAATTTACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGTGGAATCGGCATAATCGCCATTGGAGGCATTGCTTAAATCGCGCTCAGGTGCCAGCCACTGGtt                             <  1:510772/236‑1 (MQ=255)
                                                                                                                                                                                                                                              cccACCACGGGGTGGAAAATAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATGCCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGTCATTGGAGGCATTGCTTAAATCGCGCGCAGGTGCCAGCCACTGGTTCAGCGCACTTAACgcc             <  2:474869/251‑1 (MQ=255)
                                                                                                                                                                                                                                                         gTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGGACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCCGGTGCCAGCCACTGGGTCAGCGCCCTTAAAGCCAGACACTGcac  >  1:421566/1‑251 (MQ=255)
                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                        
GCATACCGAATCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCAGGTGCCAGCCACTGGTTCAGCGCACTTAACGCCAGCCACTGCAC  >  NC_000913/1174373‑1174868

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