Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,080,8370GT42.9% 62.1 / 26.1 49intergenic (+24/‑519)putP/efeUproline:sodium symporter/ferrous iron permease (pseudogene)
Reads supporting (aligned to +/- strand):  ref base G (8/20);  new base T (21/0);  total (29/20)
Fisher's exact test for biased strand distribution p-value = 1.28e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGATTGTAGTGTTTAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCTTCGGGGAACGCTTCAGCGATTTTGACGCCACCAA  >  NC_000913/1080697‑1080984
                                                                                                                                            |                                                                                                                                                    
gggATTGTAGTGTTTAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTca                                                                                                                                            <  2:837084/151‑1 (MQ=255)
   aTTGTAGTGTTTAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTTCGGCGGTTTTGTTTGGCTTcagc                                                                                                                                          >  1:413321/1‑150 (MQ=255)
      gTAGTGTTTAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCg                                                                                                                                      <  1:571643/151‑1 (MQ=255)
        aGTGTTTAGTTTGCTGGGTAAAGCGCCGTCAGCTGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCggg                                                                                                                                    <  2:701546/151‑1 (MQ=255)
            tttAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGGCTCCGCACCGGGttt                                                                                                                                 >  2:4743/1‑149 (MQ=255)
             ttAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTgc                                                                                                                                <  2:124523/150‑1 (MQ=255)
                gTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTGGGTTCAGCAGCGGGTTGCGCTc                                                                                                                            >  1:872104/1‑151 (MQ=255)
                  ttGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGGTTCAGCAGCGGGGTGCGCcccc                                                                                                                          >  1:467300/1‑151 (MQ=255)
                   tGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTGGCTTCAGCAGCGGGGTGCGCcccc                                                                                                                          >  1:489819/1‑150 (MQ=255)
                      tGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGGTTTGTTTTGCTTTTGCAGCGGGGTTTGCTccc                                                                                                                          >  2:885864/1‑147 (MQ=255)
                       gggTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTaa                                                                                                                      <  2:231947/150‑1 (MQ=255)
                        ggTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGATTTTTTTGGGTTCCGGCACGGGGTGCGGCCCCtt                                                                                                                        >  1:512819/1‑147 (MQ=255)
                           aaaGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAAtgtg                                                                                                                  <  2:53686/150‑1 (MQ=255)
                            aaGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTTGCTTCAGCAGCGGGTTGCGCTCCCTTTATTTGcc                                                                                                                >  1:717505/1‑151 (MQ=255)
                              gcgcCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTc                                                                                                              <  1:921580/151‑1 (MQ=255)
                                gcCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGCCTCAGCAGCGGGGTGCGCCCCCCTTATGTGCCTc                                                                                                              >  2:865704/1‑149 (MQ=255)
                                gcCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTTGCTTCAGCAGCGGGGTGCGCTCCCTTAATGTGCCTCGc                                                                                                            >  1:487945/1‑151 (MQ=255)
                                gcCGTCAGCGGCGATGCAAAAACGCCTTGCCGAGGCCGACGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGCTTCCAAAGCggg                                                                                                                                    >  2:721360/1‑127 (MQ=255)
                                    tCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGCTTTTTTTTGCTTACGCAGCGGGGTGCGCTCCCCTAATGTGCCTCGcc                                                                                                           >  1:416845/1‑148 (MQ=255)
                                     cAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGCTTCAGCAACGGGGTTCGCCCCCCTTAAGTGCCCCGcc                                                                                                           >  2:376482/1‑147 (MQ=255)
                                           cgATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTg                                                                                                <  2:685708/151‑1 (MQ=255)
                                            gATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGCTTTAGCAGCGGGGTGCGCCCCCCTAATGTGCCTCGCCCTATtaaat                                                                                                  >  2:124629/1‑147 (MQ=255)
                                            gATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTGGCTTCAGCAGCGGGGTGCGCTCCCCTAATGTGCCCCGCCCTATaaaat                                                                                                  >  2:576687/1‑147 (MQ=255)
                                            gATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGa                                                                                               <  2:584729/151‑1 (MQ=255)
                                             aTGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGGTGCGCTCCCTTAATGGGCCCCGCCATATAAAtt                                                                                                 >  1:260218/1‑148 (MQ=255)
                                              tGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGCTTCCGCAGCGGGGTGCGCTCCCTTAATGTGCCTCGCCATATAAAtttta                                                                                              >  2:915973/1‑147 (MQ=255)
                                                  aaaaCGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTg                                                                                         <  1:243016/151‑1 (MQ=255)
                                                    aaCGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGc                                                                                        <  2:998456/150‑1 (MQ=255)
                                                     aCGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTGGCTTCCGCAGGGGGTTTCGCTCCCCTTATGTGGCTCGCCCTATAAAATTAAt                                                                                             >  1:878894/1‑144 (MQ=255)
                                                       gCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAggg                                                                                    <  1:45950/151‑1 (MQ=255)
                                                         tttGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGGTTTGTTTTGCTTTAACAACGGGGTGCGCCCCCTTAATGAGGCTCGCCCTATAAATTTAATGGGGGAggg                                                                                    >  2:455461/1‑149 (MQ=255)
                                                           tGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAgc                                                                                 <  1:759214/150‑1 (MQ=255)
                                                            gCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAgcgc                                                                               <  1:1029244/151‑1 (MQ=255)
                                                            gCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAgcgc                                                                               <  1:463921/151‑1 (MQ=255)
                                                             ccGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAgcgcg                                                                              <  2:1027822/151‑1 (MQ=255)
                                                                aGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAg                                                                           <  1:246170/151‑1 (MQ=255)
                                                                   ccGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGCTTCCGCAGCGGGGTTCGCTCCCTTAAATTGCCTCGCCCTATAAATTGAATGGTGGCGGGAGCGGGCAggg                                                                         >  2:32976/1‑150 (MQ=255)
                                                                   ccGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGCTTCAGCAGCGGGTTGCGCTCCCTTTATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAgggg                                                                        >  1:931132/1‑151 (MQ=255)
                                                                    cGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGCTTCAGCAGCGGGGTGCGCTCCCTTAATGTGCCTCGCCCTATAAATTGAATGGGGCAGGGGGCGCGCAggggg                                                                       >  2:468123/1‑151 (MQ=255)
                                                                       tGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTGGCTTCAGCAACGGGGTGGGCCCCCCTAATGAGGCTCCCCCTATAAATTTAAAGGGTGCGGGAGCGCGCCgggggggg                                                                    >  2:287298/1‑151 (MQ=255)
                                                                          gcACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCa                                                                 <  2:348910/151‑1 (MQ=255)
                                                                             cTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGCTTCAGCAGCGGGGTGCGCCCCCCTTATGTGCCTCGCCCTATAAAATGAAATGGGCAGGGGGCGCGCAGGGGGCGGCCaa                                                                >  1:529701/1‑149 (MQ=255)
                                                                              tATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTGGCTTCAGCAGCGGGGTGCGCTCCCTTAATGTGCCTCGCCCTATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATc                                                              >  1:742173/1‑150 (MQ=255)
                                                                              tATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTGGCTTCAGCAGCGGGGTGCGCTCCCTTAATGTGCCTCGCCATATAAATTTAATGGGGCAGGGAGCGCGCAGGGGGCGGCCAATcg                                                             >  1:841045/1‑151 (MQ=255)
                                                                              tATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTTGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCCCCCCCTATTAATTGAATGGTGCCGGGGGGGCGCAggggg                                                                       >  2:951710/1‑141 (MQ=255)
                                                                                                   cGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCg                                        <  1:1017057/151‑1 (MQ=255)
                                                                                                        gCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCt                                    <  2:705176/150‑1 (MQ=255)
                                                                                                               aGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTGGGTTCAGCAGCGGGGTGTGGTCCCCTAATGGGGCTCGCCATATAAATTTAATGGTGCAGGGGGGGCGCAGGGGGGGGCCCATCCGCCCCGCCCCCCGCTGTCCCCGGCCTCgggg                              >  2:1015325/1‑149 (MQ=255)
                                                                                                                                        ttttGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCTTCGGGGAACGCTTCAGCGATTTTGACGccacc    >  1:443020/1‑150 (MQ=255)
                                                                                                                                           tgtttgGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCTTCGGGGAACGCTTCAGCGATTTTGACGCCACCaa  >  1:1012949/1‑149 (MQ=255)
                                                                                                                                            |                                                                                                                                                    
GGGATTGTAGTGTTTAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCTTCGGGGAACGCTTCAGCGATTTTGACGCCACCAA  >  NC_000913/1080697‑1080984

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