Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,510,3910TA41.5% 48.1 / 16.8 41V130E (GTA→GAA) ydcRputative DNA‑binding transcriptional regulator and putative aminotransferase
Reads supporting (aligned to +/- strand):  ref base T (15/9);  new base A (0/17);  total (15/26)
Fisher's exact test for biased strand distribution p-value = 3.12e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGTCACTCGAGATGAAGCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCCAACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGCGT  >  NC_000913/1510241‑1510531
                                                                                                                                                          |                                                                                                                                            
cAGTCACTCGAGATGAAGCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCCAACAACTAAACCGCTCGCTGGCGCAGGt                                                                                                                                              >  2:722911/1‑151 (MQ=255)
 aGTCACTCGAGATGAAGCAGTCGATATCAACACTTATACTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCCAACAACTAAACCGCTCGCTGGCGCAGGTa                                                                                                                                             <  2:165038/151‑1 (MQ=255)
    cACTCGAGATGAAGCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCCAACAACTAAACCGCTCGCTGGCGCAGGTAAg                                                                                                                                           <  1:942770/150‑1 (MQ=255)
        cGAGATGAAGCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCCAACAACTAAACCGCTCGCTGGCGCAGGTAAGCaaaa                                                                                                                                      <  1:770617/151‑1 (MQ=255)
            aTGAAGCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCCAACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGc                                                                                                                                  <  1:342570/151‑1 (MQ=255)
                      cGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCCAACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGACACAGCGAt                                                                                                                         >  2:90466/1‑150 (MQ=255)
                          atCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCCAACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCg                                                                                                                    >  1:199623/1‑151 (MQ=255)
                             aaCACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCCAACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTg                                                                                                                  >  1:886834/1‑150 (MQ=255)
                               cttttATTTTTTTGATATTTTGCAGCCCGGCCGCGACCCGTCGGTCTTCCCTTTGCCCCCGGCCTTCCCCGCCCGGCTATTTCCCCCCCCCCAAAAACAAACCCCCGCGCGGGCCGAGGAAAAAAAAACCGCCACAGCGATGAGc                                                                                                                     <  2:197792/142‑1 (MQ=255)
                                   taattttttGTATTTGTGGCAGGCCCGCCTCGTCCCGTCGGTCTCTCCTTTCCCGCCGCCCTTCCCCCCCCCGCGTCTTTCCCCCCCCCAACAACAAACCCGCTCGCGGGCGCGGGTAACAAAAACCGCCACAGCGATGAGCGTGAt                                                                                                                <  2:691030/143‑1 (MQ=255)
                                    atatTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCCAACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGaaa                                                                                                           >  1:855880/1‑150 (MQ=255)
                                     ttttttttATTTGTTGCAGCCCTGCCGCGTTCGGTCGTTCTTCCCTTTTCCCCCGCCCTTCCCCGCCCCGCGATTTTCCCCCCCCCAAAAAATAACCCGCTCGCTGGCGCAGGTAAGAAAAACCGCCACAGCGATGAGCGTGATTGaaa                                                                                                           <  1:816696/147‑1 (MQ=255)
                                     tatTTTTGATATGTTGCGGGCCAGTCGCGCTCCGTCGGTTTTTCCGTTTGCCTCGGCCTTCCCCGCCCCGCGTCTTTCCCCCCCCCAACAACAAACCCGCTCGCTGGCGCAGGAAAGCAAAACCGCCACAGCGATGAGCGTGATTg                                                                                                              <  1:461461/146‑1 (MQ=255)
                                       tttttAATTTGTTGCAGGCCAGCCGCGACCCGTCGGTCTTCCCGTTTCCCTCGGCCTTCCCCGACCCGCGACTTTCCCCCCCCCAACAACTAACCCGCTCGCTGGCGCAGATAAACAAAACCGCCACAGCGATGAGCGTGATTGaaa                                                                                                           <  1:554986/147‑1 (MQ=255)
                                           ttgttttgtggcggcccAGCCGCGACCCGTCGTTCTTCCCTTTTCCCCCGCCCTTCCCCGCCCCGCGTTTTTCCCCCCCCCAACAACAAACCCGCTCGCTGGCGCAGGTAAGAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTAc                                                                                                     <  1:232549/135‑1 (MQ=255)
                                              atatGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCCAACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGCTTGAAAACTTACCGcc                                                                                                 >  2:628953/1‑151 (MQ=255)
                                                 tGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCCAACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAgg                                                                                              >  2:121917/1‑151 (MQ=255)
                                                        ggCCGGTCGCGACCCGTCGTTCGTCCCTTTTCCCCCGGCCTTCCCCCACCCGCGATTTTCCCCCCCCAAAAAACAAACCCGCCCGCTGGCGCAGGAAAAAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACg                                                                                         <  1:622037/149‑1 (MQ=255)
                                                                gcgcTCCGTGGTTCTTCCCTTTCGCCTGGCCCTTCCCCGCCCCGCGTTTTTCCCCCCCCCAACAACAAACCCGCTCGGGGGCGCGGAAAAGAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGc                                                                                <  1:21690/150‑1 (MQ=255)
                                                                     cgggcgttttttCCTTTTTCCGCGGCTTTTCCCCCCCCGCCATTTTCCCCCCCCCAAAAACAAAACCGCTCCTGGGCGGGGAAAAGAAAACCCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTc                                                                             <  1:152191/138‑1 (MQ=255)
                                                                              gTTCCGTTTGCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCCAACAACTAAACCGCTCGCTGGCGACGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCt                                                                  >  2:501305/1‑150 (MQ=255)
                                                                                      tttcccccgCCTTTCCCCGCCCCGCGTTTTTCCCCCCCCAAACAACAAACCCCCCCGGGGGCGCAGAAAACAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTAt                                                                      <  1:925827/132‑1 (MQ=255)
                                                                                      tGCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCCAACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTAt                                                         >  2:14902/1‑151 (MQ=255)
                                                                                                 tttCCCGACCCGCGACTTTTCCCCCTCCAACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAgg                                               >  1:677957/1‑151 (MQ=255)
                                                                                                  tccccccccccGCGGCTTTCCCCCCCCAAACATCTAACCCGCCCGCTGGCGAGGGAAAGAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACa                                                 <  1:258119/141‑1 (MQ=255)
                                                                                                       ccgccccgggcTTTTCCCCCCCCAAAAAACAAACCCCCCCGGGGGCGGAGTAAAGAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAg                                                <  1:362037/136‑1 (MQ=255)
                                                                                                           ccGCGTTTTTCCCCCCCCCAACAACAAACCCGCCCGCGGGCGAAGGAAAAAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCa                                      <  1:154940/151‑1 (MQ=255)
                                                                                                              cgcttttcCCCCCCCCAACAACAAACCCGCTCGCGGGCGAAGGAAAGAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCAtt                                    <  1:217906/146‑1 (MQ=255)
                                                                                                                ccTTTTCCCCCCCCAAACAACAAACCCCCTCGCGGCCGCAGGAAAAAAAACCCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTtc                                  <  2:277863/150‑1 (MQ=255)
                                                                                                                      cccccTCCAACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCtgat                           >  2:49922/1‑151 (MQ=255)
                                                                                                                             caacaaCTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTg                    >  1:297717/1‑151 (MQ=255)
                                                                                                                               aacaacTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTg                    >  2:276127/1‑150 (MQ=255)
                                                                                                                               aacaaaAAACCCGCTCGCTGGCGAAGAAAAAAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAAt                      <  1:821320/148‑1 (MQ=255)
                                                                                                                                  aataaacccccccgctGGCGCAGTAAAGAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCa                 <  2:476225/138‑1 (MQ=255)
                                                                                                                                   acTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCAtt               >  1:736923/1‑151 (MQ=255)
                                                                                                                                      accccGCTCTGTGGCGCAGTTAAAAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACt            <  2:679352/148‑1 (MQ=255)
                                                                                                                                      aacccGCGCGCGGGCGCAGGAAAAAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTa              <  1:326225/146‑1 (MQ=255)
                                                                                                                                          cgcccgcgGGGGCAGGAAACAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCg        <  1:142496/151‑1 (MQ=255)
                                                                                                                                           gctcgctGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCgg       >  2:514705/1‑151 (MQ=255)
                                                                                                                                                ctGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGCGt  >  1:479297/1‑151 (MQ=255)
                                                                                                                                                  ggCGCAGTAAAAAAAACCCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGc    <  2:522546/147‑1 (MQ=255)
                                                                                                                                                      cAGGAAAACAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTg           <  2:361978/136‑1 (MQ=255)
                                                                                                                                                          |                                                                                                                                            
CAGTCACTCGAGATGAAGCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCCAACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGCGT  >  NC_000913/1510241‑1510531

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