Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,377,1400AC24.0% 101.8 / 11.6 50K103N (AAA→AAC) ycjSputative oxidoreductase YcjS
Reads supporting (aligned to +/- strand):  ref base A (26/12);  new base C (0/12);  total (26/24)
Fisher's exact test for biased strand distribution p-value = 2.23e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTACGAACATACCCTGATGGCACTGGAAGCGGGCTGCCATGTGATGTGCGAAAAACCGCCCGCCATGACGCCAGAACAGGCGCGGGAAATGTGCGATACCGCGCGCAAACTGGGCAAGGTGCTGGCCTACGACTTTCACCATCGTTTTGCGCTCGATACGCAACAGCTGCGTGAACAGGTGACCAACGGCGT  >  NC_000913/1377012‑1377277
                                                                                                                                |                                                                                                                                         
ggggACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTACGAACATCCCCTGATGGCACTGGAAGCGGGCTGCCATGTGATGTGCGAAAACCCGCCCGCCAt                                                                                                                                <  2:322073/139‑1 (MQ=255)
ggggACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTACGAACATCCCCTGATGGCACTGGAAGCGGGCTGCCATGTGATGTGCGAAAAACCGCCCGCCAt                                                                                                                                <  1:493911/139‑1 (MQ=255)
  ggACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTACGAACATACCCTGATGGCACTGGAAGCGGGCTGCCATGTGATGTGCGAAAAACCGCCCGCCATGa                                                                                                                              >  2:1221634/1‑140 (MQ=255)
  ggACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTACGAACATACCCTGATGGCACTGGAAGCGGGCTGCCATGTGATGTGCGAAAAACCGCCCGCCATGa                                                                                                                              >  1:1219897/1‑140 (MQ=255)
   gACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCTTTTTCATTACGAACATCCCCTGATGGCCCGGGAAGCGGGCTGCCATGTGATGTGCGAAAACCCGCCCGCCATGAc                                                                                                                             <  2:763574/140‑1 (MQ=255)
   gACGATCCGCAGGCCATGCTGCGGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCCCCTACCCGTTTTCATTACGAACATCCCCTGATGGCACTGGAAGCGGGCTGCCATGTGATGTGCGAAAACCCGCCCGCCATGAc                                                                                                                             <  2:895665/140‑1 (MQ=255)
     cGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCCCACCTAACCGTTTTCATTACGAACATCCCCTGATGGCACTGGAAGCGGGCTGCCATGTGATGTGCGAAAAACCGCCCGCCATGACGc                                                                                                                           <  1:81853/140‑1 (MQ=255)
        tCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGGGGTTAGCGTCTGCTCACCTAACCGTTTTCATTACGAACATACCCTGATGGCACTGGAAGCGGGCTGCCATGTGATGTGCGAAAAACCGCCCGCCATGACGCCAg                                                                                                                        >  2:991513/1‑140 (MQ=255)
              gcccATGCTGCGGGGGGTGAACCCTGATGTGGTTAGCGTCTGCTCACCTACCCTTTTTCATTACGAACATCCCCTGAGGGCACGGGAAGCGGGCTGCCATGTGATGTGCAAAAACCCGCCCGCCATGACGCCAGAACAgg                                                                                                                  <  2:366944/138‑1 (MQ=255)
                  aTGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTACGAACATCCCCTGATGGCACTGGAAGCGGGCTGCCATGTGATGTGCGAAAACCCGCCCGCCATGACGCCAGAACAGGcgcg                                                                                                              <  2:121043/140‑1 (MQ=255)
                  aTGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTACGAACATACCCTGATGGCACGGGAAGCGGGCTGCCATGTGATGTGCGAAAAACCGCCCGCCATGACGCCAGAACAGGcgcg                                                                                                              <  1:99697/140‑1 (MQ=255)
                      tgctgGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTACGAACATACCCTGATGGCACTGGAAGCGGGCTGCCATGTGATGTGCGAAAAACCGCCCACCATGACGCCAGAACAGGCGCGGGaa                                                                                                          <  1:1238205/140‑1 (MQ=255)
                      tgctgGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTACGAACATACCCTGATGGCACTGGAAGCGGGCTCCCATGTGATGTGCGAAAAACCGCCCGCCATGACGCCAGAACAGGCGCGGGaa                                                                                                          <  1:1335742/140‑1 (MQ=255)
                      ggcgggCGGTGAACCCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTACGAACATCCCCTGATGGCACTGGAGGCGGGCTGCCATGTGATGTGCGAAAAACCGCCCGCCATGACGCCAGAACAGGCGCGGGaa                                                                                                          <  2:940025/136‑1 (MQ=255)
                         tgGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTACGAACATACCCCGATTGCACTGGGAACGGGCTGCCATGTGATGTGCGAAAAACCGCCCGCCATGACGCCAGAACAGGCGCCGGGAAtg                                                                                                       >  2:969524/1‑140 (MQ=255)
                                         tgtgGTTAGCGTCTGCTCACCTAACCGTTTTCATTACGAACATACCCTGATGGCACTGGAAGCGGGCTGCCATGTGATGTGCGAAAAACCGCCCGCCATGACGCCAGAACAGGCGCGGGAAATGTGCGATACCGCGCGCa                                                                                       >  2:635161/1‑140 (MQ=255)
                                             gTTAGCGTCTGCTCACCTAACCGTTTTCATTACGAACATACCCTGATGGCACTGGAAGCGGGCTGCCATGTGATGTGCGAAAAACCGCCCGCCATGACGCCAGAACAGGCGCGGGAAATGTGCGATACCGCGCGCAAACt                                                                                   >  2:899954/1‑140 (MQ=255)
                                                aGCGTCTGCTCACCTAACCGTTTTCATTACGAACATACCCTGATGGCACTGGAAGCGGGCTGCCATGTGATGTGCGAAAAACCGCCCGCCATGACGCCAGAACAGGCGCGGGAAATGTGCGATACCGCGCGCAAACTggg                                                                                >  1:782309/1‑140 (MQ=255)
                                                aGCGTCTGCTCACCTAACCGTTTTCATTACGAACATACCCTGATGGCACTGGAAGCGGGCTGCCATGTGATGTGCGAAAAACCGCCCGCCATGACGCCAGAACAGGCGCGGGAAATGTGCGATACCGCGCGCAAACTggg                                                                                >  2:122133/1‑140 (MQ=255)
                                                  ggTCTGCGCACGTAACAGTTTTCATGACGAACACACCCTGATGCCACTGGAAGGGGCCTGCGATGTGGTGTGGAAAAACCCGCCCGCCATGACGCCAGAACAGGCGCGGGAAATGTGCGATACCGCGCGCAAACTGGGCa                                                                              <  2:507638/139‑1 (MQ=255)
                                                  cggtggCTCCCCTAACCGTTTTAATTACGAACATCCCCTGATGGCACTGGAAGGGGGCTGCCATGTGAGGTGCGAAAACCCGCCCGCCATGACGCCAGAACAGGCGCGGGAAATGTGCGATACCGCGCGCAAACTGGGCa                                                                              <  2:789752/135‑1 (MQ=255)
                                                  cGTCTGCTCACCTAACCGTTTTCATTACGAACATCCCCTGATGGCACTGGAAGCGGGCTGCCATGTGATGTGCGAAAAACCGCCCGCCATGACGCCAGAACAGGCGCGGGAAATGTGCGATACCGCGCGCAAACTGGGCa                                                                              <  1:1221634/140‑1 (MQ=255)
                                                      tGCTCACCTAACCGTTTTCATTACGAACATACCCTGATGGCACTGGAAGCGGGCTGCCATGTGATGTGCGAAAAACCGCCCGCCATGACGCCAGAACAGGCGCGGGAAATGTGCGATACCGCGCGCAAACTGGGCAAGGt                                                                          >  1:1485975/1‑140 (MQ=255)
                                                       gCTCACCTAACCGTTTTCATTACGAACATACCCTGATGGCACTGGAAGCGGCCTGCCATGTGATGTGCGAAAAACCGCCCGCCATGACGCCAGAACAGGCGCGGGAAATGTGCGATACCGCGCGCAAACTGGGCAAGGTg                                                                         >  1:1063397/1‑140 (MQ=255)
                                                        cTCACCTAACCGTTTTCATTACGAACATACCCTGATGGCACTGGAAGCGGGCTGCCATGTGATGTGCGAAAAACCGCCCGCCATGACGCCAGAACAGGCGCGGGAAATGTGCGATACCGCGCGCAAACTGGGCAAGGTGc                                                                        >  2:1052547/1‑140 (MQ=255)
                                                              tAACCGTTTTCATTACGAACATACCCTGATGGCACTGGAAGCGGGCTGCCATGTGATGTGCGAAAAACCGCCCGCCATGACGCCAGAACAGGCGCGGGAAATGTGCGATACCGCGCGCAAACTGGGCAAGGTGCTGGCCt                                                                  >  2:1450484/1‑140 (MQ=255)
                                                                aCCGTTTTCATTACGAACATACCCTGATGGCACTGGAAGCGGGCTGCCATGTGATGTGCGAAAAACCGCCCGCCATGACGCCAGAACAGGCGCGGGAAATGTGCGATACCGCGCGCAAACTGGGCAAGGTGCTGGCCTAc                                                                >  2:1119425/1‑140 (MQ=255)
                                                                   gTTTTTATTACGAACATCCCCTGATGCCACTGGAGGCGGGCTGCCATGGGATGGGCGAAAACCCGCCCGCCATGACGCCAGAACAGGCGCGGGAAATGTGCGATACCGCGCGCAAACTGGGCAAGGTGCTGGCCTACGAc                                                             <  1:991513/140‑1 (MQ=255)
                                                                       tCATTACGAACATACCCTGATGGCACTGGAAGCGGGCTGCCATGTGATGTGCGAAAAACCGCCCGCCATGACGCCAGAACAGGCGCGGGAAATGTGCGATACCGCGCGCAAACTGGGCAAGGTGCTGGCCTACGACTTTc                                                         >  2:875082/1‑140 (MQ=255)
                                                                        cATTACGAACATCCCCTGATGGCACTGGAAGCGGGCTGCCATGTGATGTGCGAAAACCCGCCCGCCATGACGCCAGAACAGGCGCGGGAAATGTGCGATACCGCGCGCAAACTGGGCAAGGTGCTGGCCTACGACTTTCa                                                        <  2:146327/140‑1 (MQ=255)
                                                                         aTTACGAACATACCCTGATGGCACTGGAAGCGGGCTGCCATGTGATGTGCGAAAAACCGCCCGCCATGACGCCAGAACAGGCGCGGGAAATGTGCGATACCGCGCGCAAACTGGGCAAGGTGCTGGCCTACGACTTTCAc                                                       >  2:207422/1‑140 (MQ=255)
                                                                                 aatccccTGATGGCCCTGAAAGCGGGCTGCCAGGTGATGTGCGAAGACCCGCCCGCCATGACGCCAGAACAGGCGCGGGAAATGTGCGATACCGCGCGCAAACTGGGCAAGGTGCTGGCCTACGACTTTCACCATCGttt                                               <  2:75797/136‑1 (MQ=255)
                                                                                   tACCCTGATGGCCCTGGAGGCGGGCTGCCAGGTGATGTGCAAAAAACCGCCCGCCATGACGCCAGAACAGGCGCGGGAAATGTGCGATACCGCGCGCAAACTGGGCAAGGTGCTGGCCTACGACTTTCACCATCGTTTTg                                             <  1:969524/140‑1 (MQ=255)
                                                                                         gATGGCACTGGAAGCGGGCTGCCATGTGATGTGCGAAAAACCGCCCGCCATGACGCCAGAACAGGCGCGGGAAATGTGCGATACCGCGCGCAAACTGGGCAAGGTGCTGGCCTACGACTTTCACCATCGTTTTGCGCTCg                                       <  1:1119425/140‑1 (MQ=255)
                                                                                          gtgccACTGGAGGGGGGCTGCCATGTGATGTGCAAAAACCCGCCCGCCATGACGCCAGAACAGGCGCGGGAAATGTGCGATACCGCGCGCAAACTGGGCAAGGTGCTGGCCTACGACTTTCACCATCGTTTTGCGCTCGa                                      <  2:1485975/136‑1 (MQ=255)
                                                                                          aTGGCACTGGAAGGGGGCTGCCATGTGATGTGCAAAAACCCGCCCGCCATGACGCCAGAACAGGCGCGGGAAATGTGCGATACCGCGCGCAAACTGGGCAAGGTGCTGGCCTACGACTTTCACCATCGTTTTGCGCTCGa                                      <  2:43807/140‑1 (MQ=255)
                                                                                          aTGGCACTGGAAGCGGGCTGCCATGTGATGTGCGAAAAACCGCCCGCCATGACGCCAGAACAGGCGCGGGAAATGTGCGATACCGCGCGCAAACTGGGCAAGGTGCTGGCCTACGACTTTCACCATCGTTTTGCGCTCGa                                      >  2:1052548/1‑140 (MQ=255)
                                                                                          aTGGCACTGGAAGCGGGCTGCCATGTGATGTGCGAAAAACCGCCCGCCATGACGCCAGAACAGGCGCGGGAAATGTGCGATACCGCGCGCAAACTGGGCAAGGTGCTGGCCTACGACTTTCACCATCGTTTTGCGCTCGa                                      >  1:1147004/1‑140 (MQ=255)
                                                                                            ggCACTGGAAGCGGGCTGCCATGTGATGTGCGAAAAACCGCCCGCCATGACGCCAGAACAGGCGCGGGAAATGTGCGATACCGCGCGCAAACTGGGCAAGGTGCTGGCCTACGACTTTCACCATCGTTTTGCGCTCGATa                                    >  1:1146908/1‑140 (MQ=255)
                                                                                               aCTGGAAGCGGGCTGCCATGTGATGTGCGAAAAACCGCCCGCCATGACGCCAGAACAGGCGCGGGAAATGTGCGATACCGCGCGCAAACTGGGCAAGGTGCTGGCCTACGACTTTCACCATCGTTTTGCGCTCGATACGc                                 >  2:597456/1‑140 (MQ=255)
                                                                                                   gAAGCGGGCTGCCATGTGATGTGCGAAAAACCGCCCGCCATGACGCCAGAACAGGCGCGGGAAATGTGCGATACCGCGCGCAAACTGGGCAAGGTGCTGGCCTACGACTTTCACCATCGTTTTGCGCTCGATACGCAACa                             <  1:875082/140‑1 (MQ=255)
                                                                                                    aaGCGGGCTGCCATGTGATGTGCGAAAAACCGCCCGCCATGACGCCAGAACAGGCGCGGGAAATGTGCGATACCGCGCGCAAACTGGGCAAGGTGCTGGCCTACGACTTTCACCATCGTTTTGCGCTCGATACGCAACAg                            >  1:1144748/1‑140 (MQ=255)
                                                                                                     aGCGGGCTGCCATGTGATGTGCGAAAAACCGCCCGCCATGACGCCAGAACAGGCGCGGGAAATGTGCGATACCGCGCGCAAACTGGGCAAGGTGCTGGCCTACGACTTTCACCATCGTTTTGCGCTCGATACGCAACAGc                           <  1:1052548/140‑1 (MQ=255)
                                                                                                      gggggCTGCCATGTGATGGGCAAAAACCCGCCCGCCATGACGCCAGAACAGGCGCGGGAAATGTGCGATACCGCGCGCAAACTGGGCAAGGTGCTGGCCTACGACTTTCACCATCGTTTTGCGCTCGATACGCAACAGCt                          <  2:1219897/138‑1 (MQ=255)
                                                                                                      gCGGGCTGCCATGTGATGTGCGAAAAACCGCCCGCCATGACGCCAGAACAGGCGCGGGAAATGTGCGATACCGCGCGCAAACTGGGCAAGGTGCTGGCCTACGACTTTCACCATCGTTTTGCGCTCGATACGCAACAGCt                          >  2:95975/1‑140 (MQ=255)
                                                                                                           cTGCCATGTGATGTGCGAAAAACCGCCCGCCATGACGCCAGAACAGGCGCGGGAAATGTGCGATACCGCGCGCAAACTGGGCAAGGTGCTGGCCTACGACTTTCACCATCGTTTTGCGCTCGATACGCAACAGCTGCGTg                     >  2:549854/1‑140 (MQ=255)
                                                                                                                      tgtgCGAAAAATCGCCCGCCATGACGCCAGAACAGGCGCGGGAAATGTGCGATACCGCGCGCAAACTGGGCAAGGTGCTGGCCTACGACTTTCACCATCGTTTTGCGCTCGATACGCAACAGCTGCGTGAACAGGTGAcc          >  1:138007/1‑140 (MQ=255)
                                                                                                                        tgCGAAAAACCGCCCGCCATGACGCCAGAACAGGCGCGGGAAATGTGCGATACCGCGCGCAAACTGGGCAAGGTGCTGGCCTACGACTTTCACCATCGTTTTGCGCTCGATACGCAACAGCTGCGTGAACAGGTGACCaa        >  1:845592/1‑140 (MQ=255)
                                                                                                                        tgCGAAAAACCGCCCGCCATGACGCCAGAACAGGCGCGGGAAATGTGCGATACCGCGCGCAAACTGGGCAAGGTGCTGGCCTACGACTTTCACCATCGTTTTGCGCTCGATACGCAACAGCTGCGTGAACAGGTGACCaa        >  1:867663/1‑140 (MQ=255)
                                                                                                                        tgCGAAAAACCGCCCGCCATGACGCCAGAACAGGCGCGGGAAATGTGCGATACCGCGCGCAAACTGGGCAAGGTGCTGGCCTACGACTTTCACCATCGTTTTGCGCTCGATACGCAACAGCTGCGTGAACAGGGGACCaa        >  2:708324/1‑140 (MQ=255)
                                                                                                                         gCGAAAACCCGCCCGCCATGACGCCAGAACAGGCGCGGGAAATGTGCGATACCGCGCGCAAACTGGGCAAGGTGCTGGCCTACGACTTTCACCATCGTTTTGCGCTCGATACGCAACAGCTGCGTGAACAGGTGACCAAc       <  1:1083159/140‑1 (MQ=255)
                                                                                                                           gAAAAACCGCCCGCCATGACGCCAGAACAGGCGCGGGAAATGTGCGATACCGCGCGCAAACTGGGCAAGGTGCTGGCCTACGACTTTCACCATCGTTTTGCGCTCGATACGCAACAGCTGCGTGAACAGGTGACCAACgg     <  1:95975/140‑1 (MQ=255)
                                                                                                                           gAAAAACCGCCCGCCATGACGCCAGAACAGGCGCGGGAAATGTGCGATACCGCGCGCAAACTGGGCAAGGTGCTGGCCTACGACTTTCACCATCGTTTTGCGCTCGATACGCAACAGCTGCGTGAACAGGTGACCAACgg     >  1:1065873/1‑140 (MQ=255)
                                                                                                                           aaaaaaCCGCCCGCCATGACGCCAGAACAGGCGCGGGAAATGTGCGATACCGCGCGCAAACTGGGCAAGGTGCTGGCCTACGACTTTCACCATCGTTTTGCGCTCGATACGCAACAGCTGCGTGAACAGGTGACCAACgg     <  2:1146908/139‑1 (MQ=255)
                                                                                                                              aaaCCGCCCGCCATGACGCCAGAACAGGCGCGGGAAATGTGCGATACCGCGCGCAAACTGGGCAAGGTGCTGGCCTACGACTTTCACCATCGTTTTGCGCTCGATACGCAACAGCTGCGTGAACAGGTGACCAACGGCGt  <  1:540598/140‑1 (MQ=255)
                                                                                                                                |                                                                                                                                         
TGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTACGAACATACCCTGATGGCACTGGAAGCGGGCTGCCATGTGATGTGCGAAAAACCGCCCGCCATGACGCCAGAACAGGCGCGGGAAATGTGCGATACCGCGCGCAAACTGGGCAAGGTGCTGGCCTACGACTTTCACCATCGTTTTGCGCTCGATACGCAACAGCTGCGTGAACAGGTGACCAACGGCGT  >  NC_000913/1377012‑1377277

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