Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,042,6040AC40.7% 78.9 / 14.8 54V182G (GTG→GGG) yqfAhemolysin III family HylIII inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (21/11);  new base C (0/22);  total (21/33)
Fisher's exact test for biased strand distribution p-value = 3.19e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTG  >  NC_000913/3042457‑3042746
                                                                                                                                                   |                                                                                                                                              
cgcgTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCAccc                                                                                                                                             <  2:964214/151‑1 (MQ=255)
       tATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGc                                                                                                                                        <  2:750209/149‑1 (MQ=255)
       tATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGc                                                                                                                                        <  2:82020/149‑1 (MQ=255)
        aTCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCTCCAGAAAGTGGCACGCACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCg                                                                                                                                       <  2:76/149‑1 (MQ=255)
        aTCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGa                                                                                                                                      <  2:238714/150‑1 (MQ=255)
          cAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAAt                                                                                                                                    >  2:883971/1‑150 (MQ=255)
             gCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAAc                                                                                                                                <  2:1150536/151‑1 (MQ=255)
              ccTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAAcc                                                                                                                               >  2:336053/1‑151 (MQ=255)
              ccTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAAcc                                                                                                                               <  1:494353/151‑1 (MQ=255)
                 aCAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACCCACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCAcg                                                                                                                            <  1:1170097/151‑1 (MQ=255)
                   aaaGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCAcgc                                                                                                                           >  1:479455/1‑150 (MQ=255)
                     aGGCATACCCATTACGCCTGCCCAATATACAAATAGCTCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGGGCCAGATGGCATGGTTGTATGGCATGCGTTAGCAGACGTAGAAAATCACCACGAGCGAATAAACCAcgccgc                                                                                                                        >  2:288435/1‑151 (MQ=255)
                        cATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTa                                                                                                                     <  1:212676/151‑1 (MQ=255)
                         aTACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTAc                                                                                                                    >  1:223871/1‑151 (MQ=255)
                             ccATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGc                                                                                                                 <  2:186270/150‑1 (MQ=255)
                                    gCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAgg                                                                                                         >  1:523141/1‑151 (MQ=255)
                                            aaaataCAAAAAGACCCCCAGAAGGTGCCCCACACTCCCGCCGAGCCCGAAGCCGTGCCAGATGGCATGTTTGTATGGAATGCTTTTGCGGACGTAAAAAACCCCCCCGAGCGAATAAACCCCGCCGCCTACCGCCAGTAAGGTAACGCt                                                                                                  <  2:700230/147‑1 (MQ=255)
                                            aaTATACAAATAGACCCCAAGAAAGTGGCCCACCCTCCCGCCGAGCACAAAGCCGTGCCAGATGGCAGGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCCCCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCt                                                                                                  <  2:29071/150‑1 (MQ=255)
                                            aaTATAAAAATAGATCGCCAGAAAGGGGCACACACTCCCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCCCCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCt                                                                                                  <  1:37286/150‑1 (MQ=255)
                                             atatACAAATAGGCCCCCAGAAAGTGCCACACACTACCGCCGAGCACGAGGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAAAAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGc                                                                                                <  1:42737/151‑1 (MQ=255)
                                             atatACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGc                                                                                                >  2:1172668/1‑151 (MQ=255)
                                                 aCAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCccgc                                                                                            >  2:592002/1‑151 (MQ=255)
                                                   aaaaAGACCGCCAGAAAGGGGCACACACTCCCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCCCCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGcc                                                                                               <  1:23096/146‑1 (MQ=255)
                                                   aaaTAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCccgccg                                                                                          >  1:741925/1‑151 (MQ=255)
                                                         aTCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCt                                                                                    >  1:613118/1‑151 (MQ=255)
                                                           cGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCtt                                                                                   >  2:294003/1‑150 (MQ=255)
                                                               aaaagggggCACACCCCCCCGCGGAGCACGAAGCCGGGCCAGGGGGCGTGTTTGTAGGAAAGGGTTTTGCAGACGAAAAAAACCCCCCCGAGCAAAAAACCCCCCCCCCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACt                                                                               <  2:544033/143‑1 (MQ=255)
                                                               aGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTg                                                                              >  1:634488/1‑151 (MQ=255)
                                                                   aGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGGTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCa                                                                           >  2:253977/1‑150 (MQ=255)
                                                                            cacTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAAt                                                                 >  2:790208/1‑151 (MQ=255)
                                                                                 ccgcggAGCACAAGGCGGTGCCAGAGGGCACGTTTTTAGGAAATGCTTTTGCAGACGAAAAAAATCCCCCCGAGCAAAAAACCCCCGCCCCTTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTacca                                                            <  1:660596/151‑1 (MQ=255)
                                                                                 ccgccgAGCACGAGCCCGTGCCAGGTGGCATGGTTGTTGGGAATGCTTTTCCAGACGTAAAAAACCCCCCCGGGCAAAAAACCCCCGCCGCCTCCCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTac                                                              <  1:138545/149‑1 (MQ=255)
                                                                                    cggggCACGAGGCCGTGCCAGGGGGCGTGTTTTTAGGGAGTGCGTTTCCAGGCGAAAAAACCCCCCCCGAGCAAAAAACCCACCCCCCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTacca                                                            <  2:602844/144‑1 (MQ=255)
                                                                                         cacaAAGCCGGCCAAGAGGCCATGTTTGTAGGGAATGGTTTTGCAGACGAAAAAAACCCCCCGGAGCGAAAAACCCCCGCCCCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGc                                                       <  1:495098/148‑1 (MQ=255)
                                                                                         cACGAAGCCGTGCAAGATGGCTTGTTTGTATGAAAGGCTTTTGCAGACGTAAAAAACCCCCCCGGGCAAATAACCCCCGCCCCCTCCCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGAc                                                    <  1:48190/151‑1 (MQ=255)
                                                                                            gAAGCCGTGCCAGATTGCCTGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGAcag                                                  >  2:996869/1‑150 (MQ=255)
                                                                                              aGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGAcagcc                                                >  1:627099/1‑150 (MQ=255)
                                                                                                ccgtccaggtgggCGTGTTTTTAGGGAATGCTTTGGCAGCGGAAAAAAACCCCCCGGAGAGAAAAACCCCCCCCCCCTCCCGCAAGTAAGGTACCGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGAcagcc                                                <  2:246227/137‑1 (MQ=255)
                                                                                                ccGTGCCAGATGGCATGGTTGAATGGAAGGGTTTTGCAGCCGAAGAAAATCCCCCCGAGCGAAAAACCCGCGCCCCCTCCCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGa                                                     <  1:879902/143‑1 (MQ=255)
                                                                                                     cgagaggccATGGTTTTATGAAATGCGTTTGCAAAGATAAAAAACCCCCCCGGGCAAAAAACCCCCGCCCCCTCCCGCCAGTAAGGTACCGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGccc                                           <  2:991373/140‑1 (MQ=255)
                                                                                                        gagggCAGGGTTTTAGGGAAGGGGTTTCCAGGCGTAAAAAATCCCCCCGGGAAAAAAAACCCCCCCGCTTCCCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGccc                                           <  1:1205812/142‑1 (MQ=255)
                                                                                                        gATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGcc                                     >  1:401809/1‑151 (MQ=255)
                                                                                                           gacgtgtttGTAGGGAATGGTTTTGCAGGGGAAAAAAACCCCCCGGGGCAAAAAACCCACCCCCCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCAt                                         <  2:276541/137‑1 (MQ=255)
                                                                                                             cATGTTTGTAGGAAAGGCTTTTGCAACCGTAAAAAATCCCCCCGAGAGAAAAACCCCCGCCCCCTCCCCCCAGTAAGGTAACGCGGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGAt                                 <  1:240674/150‑1 (MQ=255)
                                                                                                             cATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGAt                                 >  2:510155/1‑150 (MQ=255)
                                                                                                                gtttGTAGGAAATGCGTTTGCACACGTAAAAAACCCCCCCCAGAAAAAAACCCCCCCCCCCCACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGccc                                           <  1:178738/135‑1 (MQ=255)
                                                                                                                      aTGGGATGCGTTTGCAGACGTAGAAAATCACCCCGCGCGAATAAACAACGCAGCATACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCAcc                         >  2:1192401/1‑149 (MQ=255)
                                                                                                                          aaTGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCAcc                         >  2:932421/1‑145 (MQ=255)
                                                                                                                               gTTTGCAGACGAAAAAAACCCCCCCGAGCAAAAAACCCCCCCCCCCTCCCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATaa                 <  1:72125/148‑1 (MQ=255)
                                                                                                                                ttacgaaatggaaaaaaCCCCCCGGCGCGAAAAACCCCCCCCCCCTCCCGCCAGTAAGGTAACGCGGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCTCCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGAt                   <  2:142682/132‑1 (MQ=255)
                                                                                                                                        aaaaaaaaaCCCCCCCGACCAAAAAACCCCCGCCGCCCACCGCCAGAAGGGTAGCGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGa         <  2:383101/142‑1 (MQ=255)
                                                                                                                                           aGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTg  >  1:1064736/1‑151 (MQ=255)
                                                                                                                                           aGAAAACCCCCCAGAGCGAAAAACCCCCCCCCCCTCCCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCg     <  1:985068/148‑1 (MQ=255)
                                                                                                                                             aaaaCCCCCCCGGGCAAAAAAACCACCCCCCACCCCGCCAATAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGt   <  2:988998/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
CGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTG  >  NC_000913/3042457‑3042746

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