Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,577,1230AC30.0% 83.8 / 15.9 50V166G (GTA→GGA) yddBputative TonB‑dependent receptor
Reads supporting (aligned to +/- strand):  ref base A (28/7);  new base C (0/15);  total (28/22)
Fisher's exact test for biased strand distribution p-value = 7.58e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TTCCTGAAGAGCCTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCACTACTGGCATCG  >  NC_000913/1576986‑1577251
                                                                                                                                          |                                                                                                                                
ttCCTGAAGAGCCTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTAcc                                                                                                                                >  2:609567/1‑140 (MQ=255)
ttCCTGAAGAGCCTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTAcc                                                                                                                                >  2:193173/1‑140 (MQ=255)
ttCCTGAAGAGCCTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTAcc                                                                                                                                >  1:410839/1‑140 (MQ=255)
  ccTGAAGAGCCTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTAcccc                                                                                                                              >  2:154289/1‑140 (MQ=255)
      aaGAGCCTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATACCCCAATTTCCCCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCa                                                                                                                          <  1:193173/140‑1 (MQ=255)
       agagCCTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCAt                                                                                                                         >  1:21585/1‑140 (MQ=255)
       agagCCTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCAt                                                                                                                         <  2:410839/140‑1 (MQ=255)
          tccTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGa                                                                                                                      >  1:440969/2‑140 (MQ=255)
          gCCTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGa                                                                                                                      >  2:587442/1‑140 (MQ=255)
                     aTGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAAt                                                                                                           >  1:358619/1‑140 (MQ=255)
                          ctctTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAAc                                                                                                      >  2:396623/1‑140 (MQ=255)
                          ctctTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAAc                                                                                                      >  1:762162/1‑140 (MQ=255)
                                ttATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGcac                                                                                                <  1:362622/140‑1 (MQ=255)
                                  aTTCTCATCGTTATGCGATGTTACCCGGTCGGAACGAGTAGTGCAATAACCAAATTTCCCCTTGCTATCACCAGCGTTGAAGCGTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACaa                                                                                              <  2:358619/140‑1 (MQ=255)
                                    tctcATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACaaaa                                                                                            <  1:561295/140‑1 (MQ=255)
                                    tctcATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACaaaa                                                                                            <  1:680463/140‑1 (MQ=255)
                                       cATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTg                                                                                         >  1:242754/1‑140 (MQ=255)
                                            atatGCGATGTTAACCAGCCCGAACGCGTAGTGCGATAACCCAATTTCCCCTTGCTATCATCAGCGTTGAAGCTTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATa                                                                                    <  2:580501/140‑1 (MQ=255)
                                                 cGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCg                                                                               <  2:191765/140‑1 (MQ=255)
                                                    tGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTca                                                                            <  2:332678/140‑1 (MQ=255)
                                                      ttAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTcaca                                                                          <  1:609567/140‑1 (MQ=255)
                                                      ttAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTGGAAGCGTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTcaca                                                                          <  2:620566/140‑1 (MQ=255)
                                                      ttAACCAGTCCGAACGCGTAGTGCGAAAACCCAATTTCACTTTGTTATCATCAGCGTTAAGCCGTTTGATCTTTGCAACAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTcaca                                                                          <  2:440969/140‑1 (MQ=255)
                                                              tCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAg                                                                  >  2:563052/1‑140 (MQ=255)
                                                              tCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAg                                                                  >  1:327927/1‑140 (MQ=255)
                                                              tCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAg                                                                  >  1:605843/1‑140 (MQ=255)
                                                                cGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTa                                                                >  2:455645/1‑140 (MQ=255)
                                                                    cgcgTAGTGGGATAACCCAATTCCCCCTTGCTATCATCGGCGTTGAAGCGTTTGATGTTTGCAGAAATTGCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCt                                                            <  2:568342/140‑1 (MQ=255)
                                                                    cgcgTAGTGCGATAACCCAATTTCCCCTTGCTATCATCAGCGTTGAGGCGTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCt                                                            <  1:729122/140‑1 (MQ=255)
                                                                         aGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTAATCTTTGCATAAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACAt                                                       <  2:290463/140‑1 (MQ=255)
                                                                          gTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATc                                                      <  1:154289/140‑1 (MQ=255)
                                                                          gTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATc                                                      >  2:286133/1‑140 (MQ=255)
                                                                           gtgcaaaaCCCAAATTGCTCCTTGCTATCACAAGCGATGAAGCGTTCAACCTTCGCATCACCTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATc                                                      <  2:368701/134‑1 (MQ=255)
                                                                            gCGATACCCCAATTTCCCCTTGCTATCATCAGCGTTGAAGCTTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCaa                                                    <  2:327927/140‑1 (MQ=255)
                                                                              gATAACCCAATTTCCCCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGa                                                  <  2:349892/140‑1 (MQ=255)
                                                                                  aCCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGataat                                              <  1:465243/140‑1 (MQ=255)
                                                                                    ccAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATAc                                            <  1:563052/140‑1 (MQ=255)
                                                                                     cAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCt                                                            >  1:49558/1‑123 (MQ=255)
                                                                                     cAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCt                                                            <  2:49558/123‑1 (MQ=255)
                                                                                     cAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATAcc                                           >  2:602965/1‑140 (MQ=255)
                                                                                        aTTTCCCCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGAc                                                         <  2:599277/124‑1 (MQ=255)
                                                                                        aTTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGAc                                                         >  1:599277/1‑124 (MQ=255)
                                                                                               cTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACa                                  >  1:706844/1‑140 (MQ=255)
                                                                                               cTTCCTATCATCAGCGTTAAAGCTTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACa                                  <  2:706844/140‑1 (MQ=255)
                                                                                                ttGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACAt                                 <  2:762162/140‑1 (MQ=255)
                                                                                                  gCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATcc                               <  1:534016/140‑1 (MQ=255)
                                                                                                           aGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAatat                      >  1:642083/1‑140 (MQ=255)
                                                                                                             cGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTg                    >  2:414913/1‑140 (MQ=255)
                                                                                                             cGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTg                    >  2:129378/1‑140 (MQ=255)
                                                                                                                 gAAGCGTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTg                <  1:587442/140‑1 (MQ=255)
                                                                                                                  aaGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGc               >  1:474968/1‑140 (MQ=255)
                                                                                                                  aaGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGc               >  2:379636/1‑140 (MQ=255)
                                                                                                                  aaGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGc               <  1:602965/140‑1 (MQ=255)
                                                                                                                       tttGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCactac          >  1:85555/1‑140 (MQ=255)
                                                                                                                       tttGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCactac          >  2:734055/1‑140 (MQ=255)
                                                                                                                               ttGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCTTAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCACTACTGGCATCg  >  2:419331/1‑140 (MQ=255)
                                                                                                                                          |                                                                                                                                
TTCCTGAAGAGCCTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCACTACTGGCATCG  >  NC_000913/1576986‑1577251

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