Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,764,7850TA49.1% 50.0 / 38.7 55intergenic (‑399/+149)sufA/ydiHFe‑S cluster assembly protein/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (25/3);  new base A (0/27);  total (25/30)
Fisher's exact test for biased strand distribution p-value = 1.18e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTAA  >  NC_000913/1764635‑1764926
                                                                                                                                                         |                                                                                                                                             
tttAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCAt                                                                                                                                               >  1:462748/1‑151 (MQ=255)
   aaGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATaa                                                                                                                                             >  1:845176/1‑150 (MQ=255)
                          ctttCAGAAATAAAAAACCCCTTAAGTCTGTGCGACACAGGCTTAAGGTTTTCTCCCCCACCCGGCGCTTATCCCCGGCCCTCTCAGTGGCTTAGCTCTTGAAGGGGCGAAAAAAATAATCTAATAAAGCTACCCCGCCGTTTTAaca                                                                                                                        <  1:335187/146‑1 (MQ=255)
                           gTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAAcaca                                                                                                                      >  1:458964/1‑149 (MQ=255)
                             tCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACt                                                                                                                  >  2:55032/1‑151 (MQ=255)
                               aGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGc                                                                                                                >  2:303789/1‑151 (MQ=255)
                                  aaTAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGAt                                                                                                             >  2:439676/1‑151 (MQ=255)
                                   aaaaaaaacCCCTTAATTCTGTGCGCCACAGGTTAAGGGTTTTCTCCCCCACCCGGCGTTTATCCCCGGCACTCTCAGTGGTTTTGCTTTTGAGGGGGAGAAAAAAATAATCTCAAAAAGCTACCCCGCCGTTTTAACACAAACTGCGAtt                                                                                                            <  1:486109/146‑1 (MQ=255)
                                     aaGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGAtt                                                                                                            >  2:939653/1‑149 (MQ=255)
                                      aGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGt                                                                                                         >  2:968120/1‑151 (MQ=255)
                                         aaccccTTAAGTCTGCGCGACACAGGCTTAAGGTTTTCTCCCCCACCGGGCTTTTATCCCCGCCACTCTCAGTGGTTTAGCTTTTGAGGGGGGAAAAAAAATAATTTCAAAAACATAACCCGCCGTTTTAACACAAACTGCGATTAGtatt                                                                                                      <  1:251561/148‑1 (MQ=255)
                                             ccTTAAGTCTGCGCAACAGGGGTTAAGGGTTTTCTCCCCCACCCGGCCCTTATCCCCGCCACTCCGGGTGGTTTAGCTCTTGAGGGGGAAAAAAAAATAATCAAAAAAACCTACCCCGCCGTTTTAACACAAACTGCGATTAGTATTAttt                                                                                                  <  1:926092/151‑1 (MQ=255)
                                                     cTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACa                                                                                           >  1:452834/1‑150 (MQ=255)
                                                      tgtgCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAat                                                                                         >  2:604054/1‑151 (MQ=255)
                                                       gtgCGACACAGGTTTAAGGTTTTCTCCCCCACCGGGCGCTTATCCCCGCCACTCTAGGTGGTTTAGCTTTTGAAGGGGGAATAAAAATAATCTCAAAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACa                                                                                           <  1:231880/148‑1 (MQ=255)
                                                       gtgCGACACAGGCTTAAGGGTTTCTCCCCCACCCGGCGCTTATCTCCGGCACTCTCAGTGGTTTACCTCTTGAAGGGGCGATAAAAATAATCTAAAAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACaa                                                                                          <  2:935879/149‑1 (MQ=255)
                                                       gtgCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAat                                                                                         >  1:659340/1‑150 (MQ=255)
                                                           gACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATc                                                                                      >  1:431171/1‑149 (MQ=255)
                                                            acacGGGTTTAGGGTTTTCTCCCCCACCCGGCGTTTATCCCCGGCCCTCTCGGTGGTTTAGCTTTTGAAGGGGCAATAAAAATAATCTAAAAAAGCTACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGaa                                                                                             <  1:678334/141‑1 (MQ=255)
                                                            acacAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAgg                                                                                   >  2:377505/1‑151 (MQ=255)
                                                              acaGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCg                                                                                 >  1:222112/1‑151 (MQ=255)
                                                                 gtcTTAAGGGTTTCTCCCCCACCCGGCGCTTACCCCCGCCACTCTCAGTGGTTTAGCTCTTGAGGGGGCGATAAAAATATCCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGt                                                                               <  2:53133/148‑1 (MQ=255)
                                                                 ggCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGt                                                                               >  1:161722/1‑150 (MQ=255)
                                                                  gttgaggggTTTCCCCCCCACCCGGCGCTTATCCCCGGCACTGTAGGGGGTTTAGTTGTTGAGGGGGCAATAAAAATAATTAAAAAAAGCAAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGc                                                                                  <  1:375952/140‑1 (MQ=255)
                                                                  ggttAAGGGTTTTCCCCCCCACCCGGCGCTTATCCCCGGCACTCTCAGTGGTTTAGCTTTTGAAGGGGCAATAAAAATAATCTCAAAAATCAACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAgg                                                                                   <  1:132267/144‑1 (MQ=255)
                                                                  gCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGt                                                                               >  1:805251/1‑149 (MQ=255)
                                                                     tAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGAt                                                                           >  2:315381/1‑150 (MQ=255)
                                                                      aaGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATa                                                                          >  1:36994/1‑150 (MQ=255)
                                                                       agggtttccccccccACCCGGCTTTTACCCCCGCCACTCTCAGTGGTTTAGCTTTGAAGGGGGCAAAAAAAATAATCTAAAAAAGCTACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAg                                                                             <  1:609236/136‑1 (MQ=255)
                                                                         ggttttCTACCCCACCGGGCGTTTATCCCCGCCCCTCTCAGTGGTTTAGCTTTGAAAGGGGAAAAAAAAATAATCTAAAAAAGCAACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATaa                                                                         <  1:899207/146‑1 (MQ=255)
                                                                                 cccccATCCGGCGCTTATCCCCGCCCCTCTCGGTGGCTTAGCTTTTGAAGGGGAAAGAAAAATAATCAAAAAAAGCAAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGAt                                                                           <  1:680418/138‑1 (MQ=255)
                                                                                 cccccACCCGGCGTTATTCCCCGCCACTCTCAGGGGCTTAGTTTTTGAGGGGGAAATAAAAATAATCAAAAAAAACAAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTa                                                                <  1:498162/149‑1 (MQ=255)
                                                                                    ccAACCGGCGCTTTGCCCAGACACGCGAAGTGGTTTACCTCTGGGAGGGGATAAAGAAATAAACTAAAAATGCTACGCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTaagaa                                                            <  1:554783/151‑1 (MQ=255)
                                                                                        ccGGCGCTTATCCCCGGCACTCTCAGTGGTTTAGCTTTTGAGGGGGCGATAAAAATAATCTAAAAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTc                                                        <  2:240706/151‑1 (MQ=255)
                                                                                        ccGGCCCTTATCCCCGCCACTCTCAGTGTTTTAGCTCTTAAAGGGGCAAAAAAAATAATCTAATAAAGCTACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGt                                                         <  2:97591/150‑1 (MQ=255)
                                                                                         cGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTc                                                        >  1:319473/1‑150 (MQ=255)
                                                                                                 aTCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCt                                               >  2:675785/1‑151 (MQ=255)
                                                                                                  tctcCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTg                                              >  2:69830/1‑151 (MQ=255)
                                                                                                   ccccGGCACTCTGAGTGGCTTCGCTCTTGAAGGGGCAATAAGAATAATCTAAAAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTg                                              <  1:262136/148‑1 (MQ=255)
                                                                                                     ccGGCATTCTAAGTGGTTAAGTTTTTAAAGTGGTAAAAAGAATAATTAAAAAAAACTACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGAc                                           <  1:303776/151‑1 (MQ=255)
                                                                                                      cgcccctctAAGTGGTTTAGCTCTTGAAGGGGCGATAAAAATAATCTAAAAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACg                                          <  2:659362/146‑1 (MQ=255)
                                                                                                       ggCCCTCGCGGTGGTTTAGCTTTAGAGGGGGCAAAAAAAATAATCTCAAAAAGCAACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGAc                                           <  1:324078/149‑1 (MQ=255)
                                                                                                            tctcAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGtt                                    >  1:866338/1‑151 (MQ=255)
                                                                                                             ctcAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGtt                                    >  1:107590/1‑150 (MQ=255)
                                                                                                               cAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAAt                                 >  2:111245/1‑151 (MQ=255)
                                                                                                                      ttAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTc                          >  1:510615/1‑151 (MQ=255)
                                                                                                                         acctttGGAGGGGCAGAGAAAAATACTCTCAAAAGCCTACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGa                        <  2:969243/147‑1 (MQ=255)
                                                                                                                                aaGGGGCTAAAAAAATAATCTAAAAAAACAAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGGCATCGAACCTGGACGGAGGTTAATCCAGGTCGAttt                     <  2:334485/146‑1 (MQ=255)
                                                                                                                                   ggggCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAAc              >  2:781176/1‑151 (MQ=255)
                                                                                                                                    gggCAAAAAGAATAATCTAAAAAAGCAACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACt             <  2:363413/151‑1 (MQ=255)
                                                                                                                                     ggCGATAAAATTAATCTAAAAAAGCTACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACtt            <  1:793880/151‑1 (MQ=255)
                                                                                                                                      ggagaaaaaaataatCTCAAAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACtt            <  2:805282/142‑1 (MQ=255)
                                                                                                                                            aaaaataatCTAAAAAAGCAAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCgg        <  1:702480/146‑1 (MQ=255)
                                                                                                                                            aaaaataatCCCAAAAAGCAACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCgg        <  2:452847/146‑1 (MQ=255)
                                                                                                                                               aataatCTCAAAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCAtt    <  2:918528/150‑1 (MQ=255)
                                                                                                                                                ataatCTCAAAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTaa  <  2:462762/151‑1 (MQ=255)
                                                                                                                                                         |                                                                                                                                             
TTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTAA  >  NC_000913/1764635‑1764926

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