Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,921,549 T→A 20.5% intergenic (+24/+37) PP_1720 → / ← PP_1721 alcohol dehydrogenase/phosphoserine phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,921,5490TA20.5% 68.5 / 14.4 39intergenic (+24/+37)PP_1720/PP_1721alcohol dehydrogenase/phosphoserine phosphatase
Reads supporting (aligned to +/- strand):  ref base T (13/17);  new base A (4/4);  total (17/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.22e-01

TTGTCGCCGCAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACGGCAA  >  NC_002947/1921406‑1921690
                                                                                                                                               |                                                                                                                                             
ttGTCGCCGCAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTTccc                                                                                                                                          >  7:129092/1‑149 (MQ=255)
ttGTCGCCGCAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCCGGCCCCTTCGCGGGTAAAccc                                                                                                                                          <  8:97885/149‑1 (MQ=255)
     gccgcAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCcgcg                                                                                                                                       >  4:95569/1‑147 (MQ=255)
     gccgcAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGaa                                                                                                                                     <  4:89559/149‑1 (MQ=255)
                      ccTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCa                                                                                                                    <  4:278160/149‑1 (MQ=255)
                            cGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATaaaac                                                                                                              >  7:293020/1‑148 (MQ=255)
                                 gAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCTGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGGGTACGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCtt                                                                                                          <  4:260914/148‑1 (MQ=255)
                                 gAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCtt                                                                                                          <  3:95569/148‑1 (MQ=255)
                                 caGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGAAAGCTGGTTATGGTTATCTACCCTAGCTTGGTGTATGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCTAAGGGGCTGGCACAGCCCATAAAAACCTTc                                                                                                         <  7:45801/148‑1 (MQ=255)
                                         ccAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTcagctccag                                                                                                 >  8:36044/1‑149 (MQ=255)
                                         ccAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTcagctccag                                                                                                 >  5:141103/1‑149 (MQ=255)
                                                       ccGAGCTGGAGACCAACCAGGTGTCGTGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGcc                                                                                    <  3:70711/148‑1 (MQ=255)
                                                         gAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACACGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCAcc                                                                                 >  5:139818/1‑149 (MQ=255)
                                                          aGCTGGAGACCATGCAGGTGTCGGGCAAGCTTGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGACCCTTCGCGGGTTTACCCGCGAAggt                                                                                                                                  <  2:65744/99‑2 (MQ=255)
                                                          aGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCATGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAggg                                                                                                                                  >  1:65744/1‑99 (MQ=255)
                                                               gagaCCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGa                                                                           >  8:142284/1‑149 (MQ=255)
                                                               gagaCCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCATAGCGGGGTTAACCCCGAAAGGGCCCGCACCGCCCAAAAAAACCATAAACACAACCAAAGAGTCAGCCACCCaaat                                                                            >  6:190266/1‑145 (MQ=255)
                                                                                gggCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCgcgc                                                          <  7:142284/149‑1 (MQ=255)
                                                                                             atggtgatCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACATTCAGCTCCAGCTCAGAGTCGGCGCCCCATTGATGTCGGCATTCTGGCGCAACACTGGGTcct                                             >  7:272861/1‑148 (MQ=255)
                                                                                                gtgatCGACCCTAGCCTGGTGTACGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCggg                                           <  7:210402/148‑1 (MQ=255)
                                                                                                   ttccaCCCTTGACTGGTTTAAGCGGGGCTGGACCCTAGGCGGTTATACCCGCGAAGGGGCCGGCACAGCCCATAAAAAGCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCTCGCAACACTGCGTCCGGGTTc                                        <  8:272861/144‑1 (MQ=255)
                                                                                                       aCCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTcacca                                    >  4:57296/1‑148 (MQ=255)
                                                                                                          cTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATg                                <  3:57296/149‑1 (MQ=255)
                                                                                                               gtGGTGTAAGGTGGGCGGGCCGGTGCGGGGGTTTAGCCGGTAAGGGGCCGGCACAGCCCATAAAATGCTTCAGCTCCAGCTCAGAGTCGGCCACCCTTTGATCTCGGCATGCTGGCGGAACACTGCGTCCGGGTTCACCACATGCGGGt                           <  5:190266/148‑1 (MQ=255)
                                                                                                                  gtgtAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGa                         <  3:277377/148‑1 (MQ=255)
                                                                                                                   tgtAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATc                       <  7:36044/149‑1 (MQ=255)
                                                                                                                   tgtAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATc                       >  2:242531/1‑149 (MQ=255)
                                                                                                                       aGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCAct                   <  6:141103/149‑1 (MQ=255)
                                                                                                                            ggcgtgCCCCTGCGCGGGTGTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTcaa              <  8:293020/144‑1 (MQ=255)
                                                                                                                               tgcccccTTCGCGGGTTTACCGGCGAAGGGGCCGGCACAGCCCATAAAAACCATCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTcaacaa           <  8:129092/146‑1 (MQ=255)
                                                                                                                                   cccTTCGCGGGTAAACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCaacaacaac       >  2:157425/1‑149 (MQ=255)
                                                                                                                                     cTTCGCGGGTAAACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCaacaacaac       >  4:193200/1‑147 (MQ=255)
                                                                                                                                     cTTCGCGGGTAAACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCaacaacaac       <  3:193200/147‑1 (MQ=255)
                                                                                                                                     cTTCGCGGGTAAACCCGCGAAGGGGCCAGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCa                                                                   <  7:291152/87‑1 (MQ=38)
                                                                                                                                     cTTCGCGGGTAAACCCGCGAAGGGGCCAGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCa                                                                   >  8:291152/1‑87 (MQ=38)
                                                                                                                                      ttCGCGGGTAAACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGAGTCCGGGTTCACCACATGCGGGTGATACACTATCAACAACAACgg     >  2:112364/1‑148 (MQ=255)
                                                                                                                                         gcgGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACGGCaa  <  6:139818/148‑1 (MQ=255)
                                                                                                                                             gTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACtctc                <  8:200592/130‑1 (MQ=255)
                                                                                                                                             gTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACtctc                >  7:200592/1‑130 (MQ=255)
                                                                                                                                               |                                                                                                                                             
TTGTCGCCGCAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACGGCAA  >  NC_002947/1921406‑1921690

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