Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,151,6080GA43.1% 47.6 / 28.5 51E142K (GAA→AAA) yegJuncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (5/24);  new base A (19/3);  total (24/27)
Fisher's exact test for biased strand distribution p-value = 8.67e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTAC  >  NC_000913/2151462‑2151757
                                                                                                                                                     |                                                                                                                                                       
gCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTTAACCTTTTTTTGTTACCTATTCCAaagaaa                                                                                                                                                      >  1:360633/1‑149 (MQ=255)
  aaTGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAaagaagaag                                                                                                                                                   <  1:56113/151‑1 (MQ=255)
     gACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTTTAACCTTTTTTCTTTTTTTGTTACCTATCCCAaaaaaaaa                                                                                                                                                    >  1:306172/1‑147 (MQ=255)
       cccTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCcgac                                                                                                                                              <  1:57911/151‑1 (MQ=255)
       cccTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGATATTACAATTATTCTTTTTTTGTTGCCTAAACaaaaaaaaaa                                                                                                                                                    >  2:356638/1‑145 (MQ=255)
        ccTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGATATAACACTGCTCCCTTTTTTGTCACCTCTTCaaaaaaaaaaaaccac                                                                                                                                              >  2:294285/1‑144 (MQ=255)
           ggTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAg                                                                                                                                           <  1:760451/150‑1 (MQ=255)
           ggTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAg                                                                                                                                           <  2:1043906/150‑1 (MQ=255)
             tATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATAAAATTGATCCTTTTTTTGTGGCCTATACAAAAAAAAAAAcc                                                                                                                                                 >  2:759117/1‑142 (MQ=255)
               ttGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGGTATAACATTTATACTTTTTTTGTTACCTATCCCAAAAAAAAAAACGACGAGCTa                                                                                                                                        >  1:745050/1‑149 (MQ=255)
                 gTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAaga                                                                                                                                    <  1:909341/150‑1 (MQ=255)
                    aCAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCaa                                                                                                                                <  1:670284/151‑1 (MQ=255)
                       aaCGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAAt                                                                                                                             <  2:713535/151‑1 (MQ=255)
                        aCGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAAtt                                                                                                                            <  1:852877/151‑1 (MQ=255)
                           tAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATAACAATGTTCCTTTATTTGTTACCCATCCAAAAAAAGAAGACGACGAGCTAAAAGAAAAAATa                                                                                                                           >  1:796458/1‑149 (MQ=255)
                             gAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATAACCTTTATTCTTTTTTTGTGACCCATCCCAAAAAAAAAACCGACGAGGGAAGAGCaaaa                                                                                                                              >  1:244395/1‑144 (MQ=255)
                             gAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGGTTTAACATTTATCCTTTTTTTGGGGCCTATTCAAAAAAAAAAGGCGGACAGCTCAGGGGAAAAtattttt                                                                                                                       >  2:16827/1‑145 (MQ=255)
                              aaTATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGGTTTACAAATGGTTCTTTTTTGGTGGCCCTTGCAAAAAAAAAAGCCCAGGAGCTTAGAGGAAAAATaat                                                                                                                         >  1:934984/1‑146 (MQ=255)
                               atatGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTTTATCATTTATCCTTTTTTTGTCACCTATCCCAAAGAAAAAACCGACGAGCTAAGAGCAAAAATAGtt                                                                                                                        >  1:382660/1‑148 (MQ=255)
                                 atGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATATCATTTATCCTTTTTTTGTTACCTATTCAAAAAAAAAAAACCACGAACTAAAAACaaaa                                                                                                                              >  1:977147/1‑140 (MQ=255)
                                   ggCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATAACACTTGTCCTTTTTTTTGCACCTATCCCAAAGAAGGAACAGGAGAGGTAAAAACAAAATTTGtt                                                                                                                        >  2:477040/1‑144 (MQ=255)
                                    gCCAGGAGTGGAAAATCAAGAGAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATAACCATGATTCTTTTTTTGTTGCCCATTCAAAAAAAGAAGCCGACGGACCAAGAAGAAAAATTAttttccaa                                                                                                                   >  2:284964/1‑144 (MQ=255)
                                    gCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTTTACCCTTTATAATTTTTTTTTCACCCATATCAAAAAAGAAAACGGACAGAGAAGAACAAAATTAGTTCGCTa                                                                                                                   >  2:458336/1‑148 (MQ=255)
                                     ccAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCa                                                                                                               <  1:247085/151‑1 (MQ=255)
                                      cAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTTTAACATTTATTCTTTTTTTGTTACCTATCCCAAAAAAAAAAACGACGAGGCAAAAGGAAAAATAATTCGCTaa                                                                                                                  >  1:559524/1‑147 (MQ=255)
                                      cAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCa                                                                                                               <  2:559024/150‑1 (MQ=255)
                                       aagggTGGAAAATCAAAGAAGGAAATTATTTCAGCTGGATTGTTTAGCCCCGGGGTTAAATTTACCGGCGGTAATCCCTTTAACCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGTCGAGCTAAGAGCAAAATTAGTTCGCTAATCAGc                                                                                                             <  1:431839/147‑1 (MQ=255)
                                        ggAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTTTACCCTTGATCCTTTATTTGTCACCTCTCCAAAAAAAAAAACCGACCAGCTAAAAACAAAAATAATTTGCTAATCAAcc                                                                                                            >  1:125789/1‑151 (MQ=255)
                                          aGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGATATACCCTTGATCCTTTCTTTGGCGCCCAACCAAAAAAAAAAGACGGCCATCGAAGCGAAAAAATAGTTTGCtcatcaac                                                                                                             >  2:349473/1‑146 (MQ=255)
                                             ggAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTTTATCATTGAGCCTTTTTTTGTCACCTATCCCAAAAAAAAAACCGCCCAGCTAAGAACAAAAATAGTTTGCTCATTAGcc                                                                                                            >  1:761530/1‑146 (MQ=255)
                                               aaaaTCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCGTGTATTGGTCACCTATCCAAAAGAAGAAACCGACGAGCGAAGAGCAAAATTAGTTCGCTAATCAGCCGTTatat                                                                                                     <  2:697770/151‑1 (MQ=255)
                                                aaaTCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAATCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTAtatt                                                                                                    <  2:699394/151‑1 (MQ=255)
                                                 aaTCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTAtattt                                                                                                   <  2:348439/151‑1 (MQ=255)
                                                    cAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTAtatttat                                                                                                 <  1:696269/150‑1 (MQ=255)
                                                        gaaagaaGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGtt                                                                                             <  2:1052674/151‑1 (MQ=255)
                                                         aaagaaGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTTTACAATTGATTCTTTTTTTGTTACCCATCCAAAAAAAGAAAACCACGAGCCAAAAGCAAAAATAGTTTGCTAATCAAACGCTATATTTTTTGttt                                                                                            >  1:439935/1‑151 (MQ=255)
                                                         aaagaaGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATATCAATTATTCTTTTTTTGTGACCTATACAAAAAAAAAAGACGACGAGCTAAGAGCAAAAATAATTCGCTAATTAACCGCTTTAtttttttttt                                                                                            >  1:1011808/1‑151 (MQ=255)
                                                          aagaagATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAATCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTa                                                                                           <  2:786486/151‑1 (MQ=255)
                                                          aagaagATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTa                                                                                           <  1:605233/151‑1 (MQ=255)
                                                           agaagaTATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTaa                                                                                          <  2:706401/151‑1 (MQ=255)
                                                                atatTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGCTTTATCCTTTATCCCTTTTTTGTTACCCATTCCAAAAAAAAAAACGGCGGGCTAAAAGCAAAAATAGTTTGCTTATAAACCCGTATATTTTTTGtttttaatg                                                                                      >  1:580791/1‑144 (MQ=255)
                                                                         aCTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATAACAATTTTTCCTTTTTTGTTCCCCCTTCCAAAAAAAAAAACCGCGAGCCAAAAGGAAAAATAGGTTGCTAATAAACCGTTATATTTTTTTTTTTAAATGGCgggt                                                                                >  1:298284/1‑145 (MQ=255)
                                                                            ggATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAAttg                                                                         <  1:203319/150‑1 (MQ=255)
                                                                                tGTATACGCGCGGGGATAAAATTTACGGCGGGTTTACCATTTATCCTTTTATGGTCACCTATACAAAAGAAGAAGACGGCGAGATAAGAGCAAAAATAATTTGCTCATCACCCGCTATATTTTTTGTTTAAAATGGCCGGTGGATttttgt                                                                     >  2:246113/1‑151 (MQ=255)
                                                                                      cgcgcgGGGATAAAATTTACGGCGGTTTTAACATTTATCCTTTTTTTGTCACCTATCCCAAAAAAGAAACCGACGAGGTAAAAGGAAAAATAATTTGCTAATCAGCCGTTTTATTTTTTTTTTTAAATTGGGGGTGTAttgttttttt                                                                  >  1:774825/1‑143 (MQ=255)
                                                                                      cgcgcgGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCaa                                                               <  1:369052/151‑1 (MQ=255)
                                                                                       gcgcgGGGATAAAATTTACGGCGGATATAACATTTATTCTTTTTTTGGCACCTATCCCAAAAAAAAAACCGACGAGCCAAAAGGCAAAATAGTTTGCTAATCAACCGTTATATTTTTTTTTTTAAATGGCGGGTgtaattttttttt                                                                   >  1:830250/1‑135 (MQ=255)
                                                                                                   aaTTGAGGGGGGTTATACATTTATTCTTTTGTTGGTCCCCTATCCAAAAGAAGAACCCGAGAGGCTAAGAGAAAATTTGGTTCCATAACCACCCTTTATATTTATTGTTAAATCGGGCGTTTAATTTGTTGTCTTCAACCGCCATGAGaaa                                                  <  2:325227/151‑1 (MQ=255)
                                                                                                                     ctatttacccTTTATTTGTCCCCTACCCAAAAAAAGAACCGGACAGGCAAAAAGAAAATTTAGTTCCTTAATCAGCCTTTATTTTTATTTTTAAATCGGGTGTTGAAATTTTTGTCTTCACCCGCCATGAGAAACTCACTCTGCCcaca                                   <  2:976509/141‑1 (MQ=255)
                                                                                                                        ttGACCTTTTATTGGTCCCCTACCAAAAAAAAGAAGCCGACGAGTAAAGAGAAAAATTAGTTCGCTAATCAGCCTTTATATTTATTTTTTAAACTGGCGGTTGATTTTTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAgcg                                <  1:653128/149‑1 (MQ=255)
                                                                                                                            tCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAAACGCCATGAGAAACTCACTCTGCCCACAGCCCCCGTg                          >  1:884269/1‑151 (MQ=255)
                                                                                                                                ttATTGGTCACCTATCAAAAAGAAGAAGCCGACGAGCAAAGAGAAAAATTAGTTCGCTAATCAGCCTTTATATTTTTTGTTTAAACGGGCGGTTGATTTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCa                       <  1:135031/150‑1 (MQ=255)
                                                                                                                                   ttgtcccccTACCAAAAAGAAGAAGCCGACGAGCAAAGAGAAAAATTAGTTCGCAAATCAGCCTTTATATTTTTTGTTAAAACTGGGGGTTAAATTTTTGTCTTAAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAAc                     <  1:662745/142‑1 (MQ=255)
                                                                                                                                   ttggccccttacccAAAAAAGGAGCCCGACGAGTTAAGAGAAAATTTGTTTGGCAAACCGGCCGTAATTTTTTTTTTTTAACTGGGGGGTGGAATTGTTTTTTCCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCa                       <  2:968268/135‑1 (MQ=255)
                                                                                                                                   tgggTCACCTATCAAAAAGAAAAAGCCGACGAGCTAAGAGAAAAATTAGTTCGCTAATCAGCCGTTATTTTTTTTTTTAAAACGGGCGTTTGATTTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCaa                      <  1:191897/146‑1 (MQ=255)
                                                                                                                                          ccTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTTAATTGTTGTCTTTAACCGCCCTGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTat              >  1:810470/1‑149 (MQ=255)
                                                                                                                                                    agaagaagCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTAc  >  1:359128/1‑151 (MQ=255)
                                                                                                                                                    agaagaagCCGACGAGCTAAGAGAAAAATTAGTTCGCTAATCAGCCTTTATATTTTTTGTTAAAACGGGGGTTTGAATTTTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTAc  <  1:510386/151‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                       
GCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTAC  >  NC_000913/2151462‑2151757

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