Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,151,5810GT40.0% 77.0 / 20.5 56D133Y (GAT→TAT) yegJuncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (10/23);  new base T (21/1);  total (31/25)
Fisher's exact test for biased strand distribution p-value = 9.84e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.04e-01

ATTCTAACGGCATGTTTATTGGCGTAATCAGCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCG  >  NC_000913/2151432‑2151727
                                                                                                                                                        |                                                                                                                                                   
aTTCTAACGGCATGTTTATTGGCGTAATCAGCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTg                                                                                                                                                     <  1:841173/150‑1 (MQ=255)
aTTCTAACGGCATGTTTATTCGCGTAATCAGCAATGACCCTGGTATTGTTACAAACGTAGAATATGACCAGGAGTGGAAAATCAATAAAGAAGATATTTCAGACTGGATGTATACGCGCGGTGATAAAATTTACGGCGGTTATACCATTGa                                                                                                                                                    <  2:425578/151‑1 (MQ=255)
       cGGCATGTTTATTGGCGTAATCAGCAATGACCCTGGTATTGTTGCAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATAACATTGATTCttt                                                                                                                                              >  1:623614/1‑150 (MQ=255)
                    ggCGTAATCAGCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATc                                                                                                                                <  1:413125/151‑1 (MQ=255)
                       gTAATCAGCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATTCCATTTATTCTTTTTTTGTTACCTATcc                                                                                                                               >  1:289733/1‑149 (MQ=255)
                         aaTCAGCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTTTATCATTTATCCTTTTTTGGTCACCt                                                                                                                                   >  1:720067/1‑143 (MQ=255)
                         aaTCAGCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATAACAATTATCCttttttttttgccaattcaaa                                                                                                                            >  2:466927/1‑135 (MQ=255)
                            cAGCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATAACATTTATTCTTTTTTTGTTAccccacccaaaaa                                                                                                                         >  1:453667/1‑139 (MQ=255)
                              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)
                                                                                                                                                      ttGACCTTTTATTGGTCCCCTACCAAAAAAAAGAAGCCGACGAGTAAAGAGAAAAATTAGTTCGCTAATCAGCCTTTATATTTATTTTTTAAACTGGCGGTTGATTTTTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAgcg  <  1:653128/149‑1 (MQ=255)
                                                                                                                                                        |                                                                                                                                                   
ATTCTAACGGCATGTTTATTGGCGTAATCAGCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCG  >  NC_000913/2151432‑2151727

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