Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,747,1570CG20.0% 99.5 / 11.0 43A9A (GCC→GCG) ydhSputative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (16/16);  new base G (0/8);  total (18/25)
Fisher's exact test for biased strand distribution p-value = 1.33e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATCTACACCTTATTCTCGCTTCTACAGCAACAAACTCCACTTTCTATTTCTATCTTCGAGCAGGCTGACGAGGCCGGTGTCGGGATGCCATACAGTGATGAGGAAAACTCAAAAATGATGCTGGCAAAT  >  NC_000913/1747009‑1747298
                                                                                                                                                    |                                                                                                                                             
aaTTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGcc                                                                                                                                               >  8:206951/1‑149 (MQ=255)
       tCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGCAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTa                                                                                                                                        >  1:50180/1‑149 (MQ=255)
       tCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTa                                                                                                                                        >  4:356848/1‑149 (MQ=255)
         aaGGGCGGTTAGCGGCATTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATATCAACCAAACTAATCTGATACCACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAAATCCTATGGTGGGTGCCGGGCCAACg                                                                                                                                      >  6:172118/1‑149 (MQ=255)
                   aGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATCTac                                                                                                                             <  3:316153/148‑1 (MQ=255)
                          ttttCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAAAACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTCCTATTCAGCGTGCAGGGCCTAAGGGCATATACACCTTAt                                                                                                                      >  1:287501/1‑148 (MQ=255)
                               aTCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAATTTTTTTTTTTTTGTTTCGGGGACTACGGTGATCTACACCTTATTCTCg                                                                                                                 <  2:148399/148‑1 (MQ=255)
                                  ccTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATCTACACCTTATTCTCGCTTc                                                                                                             <  7:138013/149‑1 (MQ=255)
                                      tCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATCTACACCTTATTCTCGCTTCTACa                                                                                                         <  7:195907/149‑1 (MQ=255)
                                        tGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGAt                                                                                                                                 <  6:1899/123‑1 (MQ=255)
                                        tGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGAt                                                                                                                                 >  5:1899/1‑123 (MQ=255)
                                        tGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATCTACACCTTATTCTCGCTTCTACAg                                                                                                        >  8:132864/1‑148 (MQ=255)
                                            aTTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACCGGTAAGTAACACTGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATCTACACCTTATTCTCGCTTCTACAGcaaca                                                                                                   <  5:141761/149‑1 (MQ=255)
                                                   aaaTTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATCTACACCTTATTCTCGCTTCTACAGCAACAAACTcc                                                                                             >  1:388679/1‑148 (MQ=255)
                                                           aaaaTCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATCTACACCTTATTCTCGCTTCTACAGCAACAAACTCCACtttctat                                                                                    >  6:75784/1‑149 (MQ=255)
                                                                 tcAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATCTACACCTTATTCTCGCTTCTACAGCAACAAACTCCACtttctatttcta                                                                               >  1:425532/1‑148 (MQ=255)
                                                                   aaCCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTCCAGGTAAGTAACAATGAAAAAAATTTTTTTTTTGGTTTCGGGGACTACGGTGATCTACACCTTATTCTCGCTTCTACAGCAACAAACTCCACtttctatttct                                                                                <  2:65235/145‑1 (MQ=255)
                                                                   aaCCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATCTACACCTTATTCTCGCTTCTACAGCAACAAACTCCACtttctatttcta                                                                               <  7:324190/146‑1 (MQ=255)
                                                                    aCCAAACTTATCTGATAACACTAAATTCTATATTATGCGTACTGTTAAGTAACTATGTAATATTTTTTTTTTTTTGTTTCGGGGACTACGGTGATCTACACCTTATtctc                                                                                                                  <  4:361994/110‑1 (MQ=255)
                                                                    aCCAAACTTATCTGATAACAATAAATTAGAAAGAATGCGTACATGTAACTAACAATGAAAAAAATTGATATTGAGGGTGAAGGGCCTACGGGGATCTACACATTATtctc                                                                                                                  >  3:361994/1‑110 (MQ=255)
                                                                        aaCTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATCTACACCTTATTCTCGCTTCTACAGCAACAAACTCCACTTTCTATTTCTATCTTCGaa                                                                       >  2:19222/1‑148 (MQ=255)
                                                                               cTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGa                                                                                                                                  >  7:239448/1‑83 (MQ=255)
                                                                               cTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGa                                                                                                                                  <  8:239448/83‑1 (MQ=255)
                                                                                  aTAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATCTACACCTTATTCTCGCTTCTACAGCAACAAACTCCACTTTCTATTTCTATCTTCGAGCAGGCTGAc                                                              <  5:75784/148‑1 (MQ=255)
                                                                                     acacTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATCTACACCTTATTCTCGCTTCTACAGCAACAAACTCCACTTTCTATTTCTATCTTCGAGCAGGCTGACGAg                                                           <  7:372028/148‑1 (MQ=255)
                                                                                      cacTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAATTTTTTTTTTTTTTGTTTCGGGGACTACGGTGATCTACACCTTATTCTCGCTTCTACAGCAACAAACTCCACTTTCTATTTCTATCTTCGAGCAGGCTGACGAGGc                                                         <  2:50180/149‑1 (MQ=255)
                                                                                      cacTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGTGTGCCGGGCCTACGGGGATCTACACCTTATTCTCGCTTCTACAGCAACAAACTCCACTTTCTATTTCTATCTTCGAGAAGGCTGACGAGGc                                                         >  3:12149/1‑149 (MQ=255)
                                                                                        cTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCAGGGCCTACTGGGATCTACACCTTATTCTCGCTTCTACAGCAACAAACTCCACTTTCTATTTCTATCTTCGAGCAGGCTGCCGAGGCCg                                                       <  3:356848/149‑1 (MQ=255)
                                                                                                   aGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATCTACACCTTATTCTCGCTTCTACAGCAACAAACTCCACTTTCTATTTCTATCTTCGAGCAGGCTGACGAGGCCGGTGTCGGGATg                                            <  7:132864/149‑1 (MQ=255)
                                                                                                         cGTCCAGGAAAGTAACAATGAAAAATATTTATTTTGTGGGTGCGGGGCCTACGGGGATCTACACCTTATTCTCGCTTCTACAGCAACAAACTCCACTTTCTATTTCTATCTTCGAGCAGGCTGACGAGGCCGGTGTCGGGATGCCATa                                       <  2:376793/148‑1 (MQ=255)
                                                                                                              aGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATCTACACCTTATTCTCGCTTCTACAGCAACAAACTCCACTTTCTATTTCTATCTTCGAGCAGGCTGACGAGGCCGGTGTCGGGATGCCATACAGtga                                 >  8:323176/1‑149 (MQ=255)
                                                                                                               ggTAAGTAAAAATGTAAAATATTGTTTTTTTGGTTTCCGGGCCTAGGGTGATCTACACCTTATTCTCGCTTCTACAGCAACAAACTCCACTTTCTATTTCTATCTTCGAGCAGGCTGACGAGGCCGGTGTCGGGATGCCATACAGtgat                                <  2:287501/149‑1 (MQ=255)
                                                                                                                 taagtaaCAATGAAAAAATTTTTTTTTTTTGGTTCGGGGACTACGGTGATCTACACCTTATTCTCGCTTCTACAGCAACAAACTCCACTTTCTATTTCTATCTTCGAGCAGGCTGACGAGGCCGGTGTCGGGATGCCATACAGtgatga                              <  2:388679/149‑1 (MQ=255)
                                                                                                                 aaagtaaCAAAGAATATTATTTTTATTGTTGCTGCGGGGACTACGGTGATCTACACCTTATTCTCGCTTCTACAGCAACAAACTCCACTTTCTATTTCTATCTTCGAGCAGGCTGACGAGGCCGGTGTCGGGATGCCATACAGtgatga                              <  2:419213/148‑1 (MQ=255)
                                                                                                                   agtaaCAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATCTACACCTTATTCTCGCTTCTACAGCAACAAACTCCACTTTCTATTTCTATCTTCGAGCAGGCTGACGAGGCCGGTGTCGGGATGCCATACAGTGATGAg                             <  7:206951/148‑1 (MQ=255)
                                                                                                                     taaCAATTAAAAAAATTGCTATTTTTGGTTACGTTAATAATGTGATATACACCTTATTATCGCTTCTACAGCAACAAAATCCACTTACTATTTCTATCTTAGAGCATGCTGACGAGGCCGGTGTCCGCATGCCATACAGTGATGAGGa                           >  4:51718/1‑148 (MQ=255)
                                                                                                                        cAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATCTACACCTTATTCTCGCTTCTACAGCAACAAACTCCACTTTCTATTTCTATCTTCGAGCAGGCTGa                                                               <  4:344158/109‑1 (MQ=255)
                                                                                                                        cAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATCTACACCTTATTCTCGCTTCTACAGCAACAAACTCCACTTTCTATTTCTATCTTCGAGCAGGCTGa                                                               >  3:344158/1‑109 (MQ=255)
                                                                                                                           tGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATCTACACCTTATTCTCGCTTCTACAGCAACAAACTCCACTTTCTATTTCTATCTTCGAGCAGGCTGACGAGGCCGGTGTCGGGATGCCATACAGTGATGAGGAAAACTc                     >  3:311816/1‑148 (MQ=255)
                                                                                                                              aaaaaaTTGCTATTGTGGGTGCCGGGCCTACGGGGATCTACACCTTATTCTCGCTTCTACAGCAACAAACTCCACTTTCTATTTCTATCTTCGAGCAGGCTGACGAGGCCGGTGTCGGGATGCCATACAGTGATGAGGAAAACTCaaaa                 <  4:12149/149‑1 (MQ=255)
                                                                                                                               aaaaaTTGCTATTGTGGGTGCCGTGCCTACGTGGATCTACACCTTATTCTCGCTTCTACAGCAACAAACTCCACTTTCTATTTCTATCTTCGAGCAGGCTGACGAGGCCGGTGTCGGGATGCCATACAGTGATGAGGAAAACTCaaaaa                <  3:51718/149‑1 (MQ=255)
                                                                                                                                   aTTTCTTTTTTGGGTTCGGGGTCTACGGTGATCTACACATTATTCTCGCTTCTACAGCAACAAACTCCACTTTCTATTTCTATCTTCGAGCAGGCTGACGAGGCCGGTGTCGGGATGCCATACAGTGATGAGGAAAACTCAAAAatgat            <  8:176125/149‑1 (MQ=255)
                                                                                                                                   aTTGCTATTGTGGGTGCCGGGCCTACGGGGATCTACACCTTATTCTCGCTTCTACAGCAACAAACTCCACTTTCTATTTCTATCTTCGAGCAGGCTGACGAGGCCGGTGTCGGGATGCCATACAGTGATGAGGAAAACTCAAAAatgat            >  2:21759/1‑149 (MQ=255)
                                                                                                                                             tgGGTGCCGGGCCTACGGGGATCTACACCTTATTCTCGCTTCTACAGCAACAAACTCCACTTTCTATTTCTATCTTCGAACAGGCTGACGAGGCCGGTGTCGGGATGCCATACAGTGATGAGGAAAACTCAAAAATGATGCTGGCAAat  <  1:19222/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                             
AATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATCTACACCTTATTCTCGCTTCTACAGCAACAAACTCCACTTTCTATTTCTATCTTCGAGCAGGCTGACGAGGCCGGTGTCGGGATGCCATACAGTGATGAGGAAAACTCAAAAATGATGCTGGCAAAT  >  NC_000913/1747009‑1747298

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