Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,375,1230TG47.4% 90.4 / 11.7 57G85G (GGT→GGGdcuDputative transporter
Reads supporting (aligned to +/- strand):  ref base T (6/24);  new base G (27/0);  total (33/24)
Fisher's exact test for biased strand distribution p-value = 1.49e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTCATCGGCCAAATCATGGCGCAGTTTATTACCAGCGCCT  >  NC_000913/3374983‑3375271
                                                                                                                                            |                                                                                                                                                        
ggTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCt                                                                                                                                               <  2:423242/151‑1 (MQ=255)
ggTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCt                                                                                                                                               <  1:215444/151‑1 (MQ=255)
ggTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGGTATGcc                                                                                                                                                   >  1:434738/1‑147 (MQ=255)
ggTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGGTATGCCCGc                                                                                                                                                >  2:292249/1‑150 (MQ=255)
 gtgtCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGGGGGGGGGTa                                                                                                                                                        >  2:836410/1‑142 (MQ=255)
   gtCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTAc                                                                                                                                             <  2:100814/150‑1 (MQ=255)
   gtCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGAGGGGGGGTATGTGCGCTAc                                                                                                                                             >  1:709329/1‑150 (MQ=255)
    tCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTATGCCCggt                                                                                                                                               >  2:844446/1‑145 (MQ=255)
          tAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGCGGGGGGGTATGtccccccccttggggg                                                                                                                                     >  1:92930/1‑138 (MQ=255)
           aGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAgcg                                                                                                                                     <  1:845639/150‑1 (MQ=255)
            gggTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCa                                                                                                                                   <  1:614158/151‑1 (MQ=255)
                  ggtggtGTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAgggg                                                                                                                             <  2:615245/151‑1 (MQ=255)
                  ggtggtGTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCATGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGCGGGCGGGTTTGTCCCGCTAATGGGGCACAAAgggg                                                                                                                             >  1:379721/1‑151 (MQ=255)
                   gtggtgTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGGTATGGCCCCCACATGGAGGGCCTAgggg                                                                                                                             >  2:662561/1‑150 (MQ=255)
                       tgtACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTATGCCCCCCCCATGGGGGGGATAggggg                                                                                                                            >  1:239600/1‑146 (MQ=255)
                              gATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTATGCCCCCTACATGGAGCGCCTAGGGGCCAGTcgcg                                                                                                                    >  2:74302/1‑148 (MQ=255)
                               aTCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTcgcgcg                                                                                                                  <  2:740795/149‑1 (MQ=255)
                               aTCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGAt                                                                                                                <  2:524722/151‑1 (MQ=255)
                                 cccACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGGGGGGGGTTATTCCCGCCGCACAGAGGGGATAGGGGcaggatgcggggtgg                                                                                                              >  1:517559/1‑137 (MQ=255)
                                        aaaGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTATGGCCGCTACACGGGGCGCATAGGGGGCAGGCGCGGGGTGGGGAGcct                                                                                                      >  2:277775/1‑149 (MQ=255)
                                        aaaGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCtt                                                                                                      <  2:530197/151‑1 (MQ=255)
                                             agcGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTaa                                                                                                 <  1:686359/151‑1 (MQ=255)
                                             agcGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGAGGGGGGTTTTTCCCCCTTCATTGAGGGCATAGGGGGGAGACGCGCGGTGGTGAACTtttt                                                                                                   >  2:693221/1‑149 (MQ=255)
                                                 gCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAgccg                                                                                             <  2:98988/151‑1 (MQ=255)
                                                  cTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTTTGCCCGCCACATGGGGGGCATTGGGGGCAGTCGCGCGGTGGGGAGCTTGTTTAgcc                                                                                              >  2:111781/1‑149 (MQ=255)
                                                             tcgATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGGTTTGCCCCCCACATTGAGCGCATAGGGGGCAGGCGCGCGGTGGGGGGCTTGTTTAGCCCCCCGTTaaaa                                                                                 >  2:628099/1‑149 (MQ=255)
                                                              cgATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACt                                                                               <  2:629313/150‑1 (MQ=255)
                                                              cgATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGCGGGGGGGTATGCCCCCTCCACTGAGGGCATAAGGGGCCGGCGCGCGATGGGGCGAGTGTTTTGCCGCCCCGTTaaa                                                                                  >  1:8130/1‑148 (MQ=255)
                                                                      ataATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTATGCCCGCTACATGGGGCGCATAGGGGCCAGTCGCGCGATGGGGGGCTTGTTTAGCCGCCCGGTAAAACTCAtt                                                                           >  1:127193/1‑146 (MQ=255)
                                                                      ataATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATtcgctc                                                                      <  1:706045/151‑1 (MQ=255)
                                                                      ataATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATtcgctc                                                                      <  1:533678/151‑1 (MQ=255)
                                                                      ataATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGGTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGGGCTTGTTTAACCGCCCGTTAAAAATCAtt                                                                           >  2:330427/1‑146 (MQ=255)
                                                                       taATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTCGGCGGGTATGCCCGCTCCCAGGAGTGCAAAGAGGGCAGTCGCGCGATGGGGGGGTCGTTTTACCGCCCGGTTaaaaacaatccc                                                                         >  1:843338/1‑140 (MQ=255)
                                                                         aTGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGcc                                                                   <  1:154789/151‑1 (MQ=255)
                                                                          tGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGGTATTCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGGTGGGGAGCTTGTTTAGCCGCCCGGTAAAAATCATTCGCTCGcc                                                                   >  2:676475/1‑150 (MQ=255)
                                                                             aaTCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTTTGGCCGCCACATGGAGCGCCTAGGGGGCAGGCGCGCGGTGGGGGGGTTGTTTAGCCGCCCCGTTAAAAACATTCGCTCCCCGtat                                                               >  1:198223/1‑151 (MQ=255)
                                                                             aaTCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGGTATGCCCGCCACATGGCGCGCATAGGGGGCAGTCGCGCGGGGGGGAGCTTGTTTAGACGCCCGCTTAAACTCATTCGCTCGCCgca                                                                >  1:332450/1‑148 (MQ=255)
                                                                              aTCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGGCCGTCGCGCGATGGGGAGCTTGTTAAGCCGCCCGTTAAAAATCATTCGCTCGCCGtat                                                               >  2:942468/1‑150 (MQ=255)
                                                                                   gcgcATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTAtattat                                                          <  2:689443/150‑1 (MQ=255)
                                                                                     gcATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTATTCCCCCCCCACGGAGCGCATAGGGGGCAGTCGCGCGGGGGGGAGGTTGTTTAGACGCCCGGTTAAAACCATTCCCTCGCCGGATattttttt                                                       >  1:902075/1‑149 (MQ=255)
                                                                                       aTGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGGGGGTTTGCCCCCCACCTGGAGCGCGTAGGGGCCAGTCGCGCGGTGGTGGGCTTGTTAAGACGCCCCTTTAAAATCATTTGCTCCCCGTATTTTATTTTGt                                                     >  2:425842/1‑151 (MQ=255)
                                                                                              ccAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGGTATGCCCGCTACATGGAGCGCATAGGGGGCAGGCGCGCGATGGGGAGCTTGTTAAGCCCCCCGTTTAAAATCATTCGCCCGCCGGATATTATTTTGTCGGCaa                                               >  2:544870/1‑150 (MQ=255)
                                                                                                aaCCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACtt                                            <  1:814139/151‑1 (MQ=255)
                                                                                                aaCCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGGTATGCCCGCCACATGGAGCGCCTAGGGGCCCGGCGCGCGGTGGTGGGGTTGTTTAGCCGCCCGTTAAAACTCAATCGCTCGCCGTATATTATTTTGTCGGCCACtt                                            >  2:291210/1‑151 (MQ=255)
                                                                                                aaCCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTCGGGGGGTATTGCCCCCACATGGAGCGCCCACGGGGGAGGCGCGCGAGGGGGTGGGTTTTTAATCGCCCCGTTAAACTCATTTGCTCGCCCGATatttttttttccg                                                  >  1:521256/1‑138 (MQ=255)
                                                                                                  ccGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTAc                                          <  2:241160/151‑1 (MQ=255)
                                                                                                    gCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGt                                        <  2:113887/151‑1 (MQ=255)
                                                                                                     cATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTc                                       <  2:836825/151‑1 (MQ=255)
                                                                                                      aTTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGGTATGCCCGCTACATGGAGCGCATAGGGGGCAGTCGCGCGATGGGGAGCTTGTTAAGCCGCCCCTTAAAACTCATTCGCTCGCCGGATATTATTCTGTCGGCAACTTACGt                                        >  2:275369/1‑149 (MQ=255)
                                                                                                          ccGGATTGGGGCTGTCGATTATGGCGGTGGGCGGGGATGGCCGCCACATGGGGGGCATAGGGGGCCGGCGCGCGGTGGGGGGGTTGTTTAGCCGCCCCGTAAAAATCCTTCGCTCGCCGTATTTTATTTTGTTGGCCACTTACGTCCTc                                    >  1:781727/1‑149 (MQ=255)
                                                                                                                   ggCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTCATCGGCCAAATCAt                         <  2:239600/151‑1 (MQ=255)
                                                                                                                    gCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTCATCGGCCAAATCATg                        <  1:883316/151‑1 (MQ=255)
                                                                                                                       gTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTCATCGGCCAAATCATGgc                      <  2:254227/150‑1 (MQ=255)
                                                                                                                           aTTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTCATCGGCCAAATCATGGCGCAg                  <  2:757804/150‑1 (MQ=255)
                                                                                                                                     tGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTCATCGGCCAAATCATGGCGCAGTTTATTAcc         >  1:661875/1‑149 (MQ=255)
                                                                                                                                          ggTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTCATCGGCCAAATCATGGCGCAGTTTATTACCAGCGCCt  >  2:598188/1‑151 (MQ=255)
                                                                                                                                          ggTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTCATCGGCCAAATCATGGCGCAGTTTATTACCAGCGCCt  >  2:980250/1‑151 (MQ=255)
                                                                                                                                            |                                                                                                                                                        
GGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATATTTATAATGAAATCCTGCGCATGCTGTCCAACCGCATTGCCGGATTGGGGCTGTCGATTATGGCGGTGGGCGGTTATGCCCGCTACATGGAGCGCATAGGGGCCAGTCGCGCGATGGTGAGCTTGTTAAGCCGCCCGTTAAAACTCATTCGCTCGCCGTATATTATTCTGTCGGCAACTTACGTCATCGGCCAAATCATGGCGCAGTTTATTACCAGCGCCT  >  NC_000913/3374983‑3375271

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