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

CAGTCACTCGAGATGAAGCAGTCGATATCAACACTTATATTTTTG‑ATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTG‑CCTCGGCCTTT‑CCCGACCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGCGT  >  NC_000913/1510241‑1510531
                                                                                                                                                          |                                                                                                                                            
cAGTCACTCGAGATGAAGCAGTCGATATCAACACTTATATTTTTG‑ATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTG‑CCTCGGCCTTT‑CCCGACCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGt                                                                                                                                              >  2:722911/1‑151 (MQ=255)
 aGTCACTCGAGATGAAGCAGTCGATATCAACACTTATACTTTTG‑ATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTG‑CCTCGGCCTTT‑CCCGACCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTa                                                                                                                                             <  2:165038/151‑1 (MQ=255)
    cACTCGAGATGAAGCAGTCGATATCAACACTTATATTTTTG‑ATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTG‑CCTCGGCCTTT‑CCCGACCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAg                                                                                                                                           <  1:942770/150‑1 (MQ=255)
        cGAGATGAAGCAGTCGATATCAACACTTATATTTTTG‑ATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTG‑CCTCGGCCTTT‑CCCGACCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCaaaa                                                                                                                                      <  1:770617/151‑1 (MQ=255)
            aTGAAGCAGTCGATATCAACACTTATATTTTTG‑ATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTG‑CCTCGGCCTTT‑CCCGACCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGc                                                                                                                                  <  1:342570/151‑1 (MQ=255)
                      cGATATCAACACTTATATTTTTG‑ATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTG‑CCTCGGCCTTT‑CCCGACCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGACACAGCGAt                                                                                                                         >  2:90466/1‑150 (MQ=255)
                          atCAACACTTATATTTTTG‑ATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTG‑CCTCGGCCTTT‑CCCGACCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCg                                                                                                                    >  1:199623/1‑151 (MQ=255)
                             aaCACTTATATTTTTG‑ATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTG‑CCTCGGCCTTT‑CCCGACCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTg                                                                                                                  >  1:886834/1‑150 (MQ=255)
                               cttttATTTTTTTG‑ATATTTTGCAGCCCGGCCGCGACCCGTCGGTC‑TTCCCTTTGCCCCCGGCCTTC‑CCCGCCCGGCTATTTCCCCCCCCCC‑AAAAACAAACCCCCGCGCGGGCCGAGGAAAAAAAAACCGCCACAGCGATGAGc                                                                                                                     <  2:197792/142‑1 (MQ=255)
                                   taattttttGTATTTGTGGCAGGCCCGCCTCGTCCCGTCGGTCTCTCCTTTCC‑CGCCGCCCTTC‑CCCCCCCCGCGTCTTTCCCCCCCCC‑AACAACAAACCCGCTCGCGGGCGCGGGTAACAAAAACCGCCACAGCGATGAGCGTGAt                                                                                                                <  2:691030/143‑1 (MQ=255)
                                    atatTTTTG‑ATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTG‑CCTCGGCCTTT‑CCCGACCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGaaa                                                                                                           >  1:855880/1‑150 (MQ=255)
                                     tttttttt‑ATTTGTTGCAGCCCTGCCGCGTTCGGTCGTTCTTCCCTTTTC‑CCCCGCCCTTC‑CCCGCCCCGCGATTTTCCCCCCCCC‑AAAAAATAACCCGCTCGCTGGCGCAGGTAAGAAAAACCGCCACAGCGATGAGCGTGATTGaaa                                                                                                           <  1:816696/147‑1 (MQ=255)
                                     tatTTTTG‑ATATGTTGCGGGCCAGTCGCGCTCCGTCGGTTTTTCCGTTTG‑CCTCGGCCTTC‑CCCGCCCCGCGTCTTTCCCCCCCCC‑AACAACAAACCCGCTCGCTGGCGCAGGAAAGCAAAACCGCCACAGCGATGAGCGTGATTg                                                                                                              <  1:461461/146‑1 (MQ=255)
                                       tttttA‑ATTTGTTGCAGGCCAGCCGCGACCCGTCGGTCTTCCCGTTTC‑CCTCGGCCTTC‑CCCGACCCGCGACTTTCCCCCCCCC‑AACAACTAACCCGCTCGCTGGCGCAGATAAACAAAACCGCCACAGCGATGAGCGTGATTGaaa                                                                                                           <  1:554986/147‑1 (MQ=255)
                                           ttgttttgtggcggcccAGCCGCGACCCGTCGTTCTTCCCTTTTC‑CCCCGCCCTTC‑CCCGCCCCGCGTTTTTCCCCCCCCC‑AACAACAAACCCGCTCGCTGGCGCAGGTAAGAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTAc                                                                                                     <  1:232549/135‑1 (MQ=255)
                                              atatGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTG‑CCTCGGCCTTT‑CCCGACCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGCTTGAAAACTTACCGcc                                                                                                 >  2:628953/1‑151 (MQ=255)
                                                 tGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTG‑CCTCGGCCTTT‑CCCGACCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAgg                                                                                              >  2:121917/1‑151 (MQ=255)
                                                        ggCCGGTCGCGACCCGTCGTTCGTCCCTTTTC‑CCCCGGCCTTC‑CCCCACCCGCGATTTTCCCCCCCCA‑AAAAACAAACCCGCCCGCTGGCGCAGGAAAAAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACg                                                                                         <  1:622037/149‑1 (MQ=255)
                                                                gcgcTCCGTGGTTCTTCCCTTTCG‑CCTGGCCCTTC‑CCCGCCCCGCGTTTTTCCCCCCCCC‑AACAACAAACCCGCTCGGGGGCGCGGAAAAGAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGc                                                                                <  1:21690/150‑1 (MQ=255)
                                                                     cgggcgttttttCCTTTTT‑CCGCGGCTTTT‑CCCCCCCCGCCATTTTCCCCCCCCC‑AAAAACAAAACCGCTCCTGGGCGGGGAAAAGAAAACCCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTc                                                                             <  1:152191/138‑1 (MQ=255)
                                                                              gTTCCGTTTG‑CCTCGGCCTTT‑CCCGACCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGACGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCt                                                                  >  2:501305/1‑150 (MQ=255)
                                                                                      tttcccccgCCTTTCCCCGCCCCGCGTTTTTCCCCCCCCA‑AACAACAAACCCCCCCGGGGGCGCAGAAAACAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTAt                                                                      <  1:925827/132‑1 (MQ=255)
                                                                                      tG‑CCTCGGCCTTT‑CCCGACCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTAt                                                         >  2:14902/1‑151 (MQ=255)
                                                                                                 ttt‑CCCGACCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAgg                                               >  1:677957/1‑151 (MQ=255)
                                                                                                  tccccccccccGCGGCTTTCCCCCCCCA‑AACATCTAACCCGCCCGCTGGCGAGGGAAAGAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACa                                                 <  1:258119/141‑1 (MQ=255)
                                                                                                       ccgccccgggcTTTTCCCCCCCCAAAAAACAAACCCCCCCGGGGGCGGAGTAAAGAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAg                                                <  1:362037/136‑1 (MQ=255)
                                                                                                           ccGCGTTTTTCCCCCCCCC‑AACAACAAACCCGCCCGCGGGCGAAGGAAAAAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCa                                      <  1:154940/151‑1 (MQ=255)
                                                                                                              cgcttttcCCCCCCCC‑AACAACAAACCCGCTCGCGGGCGAAGGAAAGAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCAtt                                    <  1:217906/146‑1 (MQ=255)
                                                                                                                ccTTTTCCCCCCCCAAACAACAAACCCCCTCGCGGCCGCAGGAAAAAAAACCCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTtc                                  <  2:277863/150‑1 (MQ=255)
                                                                                                                      cccccTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCtgat                           >  2:49922/1‑151 (MQ=255)
                                                                                                                             c‑aacaaCTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTg                    >  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)
                                                                                                                                                          |                                                                                                                                            
CAGTCACTCGAGATGAAGCAGTCGATATCAACACTTATATTTTTG‑ATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTG‑CCTCGGCCTTT‑CCCGACCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGCGT  >  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: ‑