Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 566,889 T→C 7.7% Q869Q (CAA→CAG)  uvrA ← uvrABC excinuclease protein A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947566,8890TC 7.7% 144.8 / 2.5 52Q869Q (CAA→CAG) uvrAuvrABC excinuclease protein A
Reads supporting (aligned to +/- strand):  ref base T (25/23);  new base C (2/2);  total (27/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

CGCAGGCAATGATCTGCCCACCCTTGGAACCGCCCTCAGGCCCCAGGTCAACCAGCCAGTCGGCAGTCTTGATCACATCCAGGTTATGCTCGATCACCACCACGGTATTGCCGTGGTCGCGCAGGCGATGCAATACGTCCAGCAGTTGCTGGATATCGGCAAAGTGCAGGCCGGTGGTAGGTTCGTCGAGGATGTACAAGGTCTTGCCGGTATCGCGCTTGGACAGTTCGCGCGACAGCTTTACCCGCTGCGCCTCACCGCCGGACAACGTGGTCGCTGATTGCCCCAG  >  NC_002947/566744‑567032
                                                                                                                                                 |                                                                                                                                               
cGCAGGCAATGATCTGCCCACCCTTGGAACCGCCCTCAGGCCCCAGGTCAACCAGCCAGTCGGCAGTCTTGATCACATCCAGGTTATGCTCGATCACCACCACGGTATTGCCGTGGTCGCGCAGGCGATGCAATACGTCCAGCAGTTGc                                                                                                                                              >  3:210799/1‑149 (MQ=255)
       aaTGATCTGCCCACCCTTGGAACCGCCCTCAGGCCCCAGGTCAACCAGCCAGTCGGCAGTCTTGATCACATCCATGTTATGCTCGATCACCACCACGGTATTGCCGTGGTCGCGCAGGCGATGCAATACGTCCAGCAACTGCTGGac                                                                                                                                         <  1:112101/147‑2 (MQ=255)
       aaTGATCTGCCCACCCTTGGAACCGCCCTCAGGCCCCAGGTCAACCAGCCAGTCGGCAGTCTTGATCACATCCAGGTTATGCTCGATCACCACCACGGTATTGCCGTGGTCGCGCAGGCGATGCAATACGTCCAGCAACTGCTGGac                                                                                                                                         >  2:112101/1‑146 (MQ=255)
                cccACCCTTGGAACCGCCCTCAGGCCCCAGGTCAACCAGCCAGTCGGCAGTCTTGATCACATCCAGGTTATGCTCGATCACCACCACGGTATTGCCGTGGTCGCGCAGGCGATGCAATACGTCCAGCAGTTGCTGGATATCGGCaaa                                                                                                                                <  2:143882/147‑1 (MQ=255)
                     ccTTGGAACCGCCCTCAGGCCCCAGGTCAACCAGCCAGTCGGCAGTCTTGATCACATCCAGGTTATGCTCGATCACCACCACGGTATTGCCGTGGTCGCGCAGGCGATGCAATACGTCCAGCAGTTGCTGGATATCGGCAAAGTGCAgg                                                                                                                         >  3:299748/1‑149 (MQ=255)
                       ttGGAACCGCCCTCAGGCCCCAGGTCAACCAGCCAGTCGGCAGTCTTGATCACATCCAGGTTATGCTCGATCACCACCACGGTATTGCCGTGGTCGCGCAGGCGATGCAATACGTCCAGCAACTGCTGGacgta                                                                                                                                      >  5:106001/1‑130 (MQ=255)
                       ttGGAACCGCCCTCAGGCCCCAGGTCAACCAGCCAGTCGGCAGTCTTGATCACATCCAGGTTATGCTCGATCACCACCACGGTATTGCCGTGGTCGCGCAGGCGATGCAATACGTCCAGCAACTGCTGGacgta                                                                                                                                      <  6:106001/134‑5 (MQ=255)
                         ggAACCGCCCTCAGGCCCCAGGTCAACCAGCCAGTCGGCAGTCTTGATCACATCCAGGTTATGCTCGATCACCACCACGGTATTGCCGTGGTCGCGCAGGCGATGCAATACGTCCAGCAGTTGCTGGATATCGGCAAAGTGCAGGCCgg                                                                                                                     >  8:144179/1‑149 (MQ=255)
                          gAACCGCCCTCAGGCCCCAGGTCAACCAGCCAGTCGGCAGTCTTGATCACATCCAGGTTATGCTCGATCACCACCACGGTATTGCCGTGGTCGCGCAGGCGATGCAATACGTCCAGCAGTTGCTGGATATCGGCAAAGTGCAGGCCggt                                                                                                                    <  7:144179/149‑1 (MQ=255)
                                   tCAGGCCCCAGGTCAACCAGCCAGTCGGCAGTCTTGATCACATCCAGGTTATGCTCGATCACCACCACGGTATTGCCGTGGTCGCGCAGGCGATGCAATACGTCCAGCAGTTGCTGGATATCGGCAAAGTGCAGGCCGGTGGTAGGTtc                                                                                                           >  7:119064/1‑149 (MQ=255)
                                     aGGCCCCAGGTCAACCAGCCAGTCGGCAGTCTTGATCACATCCAGGTTATGCTCGATCACCACCACGGTATTGCCGTGGTCGCGCAGGCGATGCAATACGTCCAGCAGTTGCTGGATATCGGCAAAGt                                                                                                                              <  4:38913/128‑1 (MQ=255)
                                        ccccAGGTCAACCAGCCAGTCGGCAGTCTTGATCACATCCAGGTTATGCTCGATCACCACCACGGTATTGCCGTGGTCGCGCAGGCGATGCAATACGTCCAGCAGTTGCTGGATATCGGCAAAGTGCAGGCCGGTGGTAGGTtcgtcta                                                                                                      <  4:114882/149‑3 (MQ=255)
                                         cccAGGTCAACCAGCCAGTCGGCAGTCTTGATCACATCCAGGTTATGCTCGATCACCACCACGGTATTGCCGTGGTCGCGCAGGCGATGCAATACGTCCAGCAGTTGCTGGATATCGGCAAAGTGCAGGCCGGTGGTAGGTTCGTCGAg                                                                                                     >  1:65103/1‑149 (MQ=255)
                                              gTCAACCAGCCAGTCGGCAGTCTTGATCACATCCAGGTTATGCTCGATCACCACCACGGTATTGCCGTGGTCGCGCAGGCGATGCAATACGTCCAGCAGTTGCTGGATATCGGCAAAGTGCAGGCCGGTGGTAGGt                                                                                                             >  3:215406/1‑136 (MQ=255)
                                              gTCAACCAGCCAGTCGGCAGTCTTGATCACATCCAGGTTATGCTCGATCACCACCACGGTATTGCCGTGGTCGCGCAGGCGATGCAATACGTCCAGCAGTTGCTGGATATCGGCAAAGTGCAGGCCGGTGGTAGGt                                                                                                             <  4:215406/136‑1 (MQ=255)
                                                 aaCCAGCCAGTCGGCAGTCTTGATCACATCCAGGTTATGCTCGATCACCACCACGGTATTGCCGTGGTCGCGCAGGCGATGCAATACGTCCAGCAGTTGCTGGATATCGGCAAAGTGCAGGCCGGTGGTAGGTTCGTCGAGGATGTAc                                                                                              <  4:299748/148‑1 (MQ=255)
                                                  aCCAGCCAGTCGGCAGTCTTGATCACATCCAGGTTATGCTCGATCACCACCACGGTATTGCCGTGGTCGCGCAGGCGATGCAATACGTCCAGCAGTTGCTGGATATCGGCAAAGTGCAGGCCGGTGGTAGGTTCGTCGAGGATGTACaa                                                                                            <  8:119064/149‑1 (MQ=255)
                                                   ccagccagTCGGCAGTCTTGATCACATCCAGGTTATGCTCGATCACCACCACGGTATTGCCGTGGTCGCGCAGGCGATGCAATACGTCCAGCAGTTGCTGGATATCGGCAAAGTGCAGGCCGGTGGTAgg                                                                                                              >  2:285182/1‑130 (MQ=255)
                                                   ccagccagTCGGCAGTCTTGATCACATCCAGGTTATGCTCGATCACCACCACGGTATTGCCGTGGTCGCGCAGGCGATGCAATACGTCCAGCAGTTGCTGGATATCGGCAAAGTGCAGGCCGGTGGTAgg                                                                                                              <  1:285182/130‑1 (MQ=255)
                                                   ccagccagTCGGCAGTCTTGATCACATCCAGGTTATGCTCGATCACCACCACGGTATTGCCGTGGTCGCGCAGGCGATGCAATACGTCCAGCAGTTGCTGGATATCGGCAAAGTGCAGGCCGGTGGTAGGTTCGTCGAGGATGTACAAg                                                                                           >  1:247314/1‑149 (MQ=255)
                                                       ccagTCGGCAGTCTTGATCACATCCAGGTTATGCTCGATCACCACCACGGTATTGCCGTGGTCGCGCAGGCGATGCAATACGTCCAGCAGTTGCTGGATATCGGCAAAGTGCAGGCCGGTGGTAGGTTCGTCGAGGAt                                                                                                  >  5:45492/1‑138 (MQ=255)
                                                       ccagTCGGCAGTCTTGATCACATCCAGGTTATGCTCGATCACCACCACGGTATTGCCGTGGTCGCGCAGGCGATGCAATACGTCCAGCAGTTGCTGGATATCGGCAAAGTGCAGGCCGGTGGTAGGTTCGTCGAGGAt                                                                                                  <  6:45492/138‑1 (MQ=255)
                                                          gTCGGCAGTCTTGATCACATCCAGGTTATGCTCGATCACCACCACGGTATTGCCGTGGTCGCGCAGGCGATGCAATACGTCCAGCAGTTGCTGGATATCGGCAAAGTGCAGGCCGGTGGTAGGTTCGTCGAGGATGTACAAGGTCTTGc                                                                                    >  5:46197/1‑149 (MQ=255)
                                                             ggCAGTCTTGATCACATCCAGGTTATGCTCGATCACCACCACGGTATTGCCGTGGTCGCGCAGGCGATGCAATACGTCCAGCAGTTGCTGGATATCGGCAAAGTGCAGGCCGGTGGTAGGTTCGTCGAGGATGTACAAGGTCTTGCCgg                                                                                 >  5:127728/1‑149 (MQ=255)
                                                             ggCAGTCTTGATCACATCCAGGTTATGCTCGATCACCACCACGGTATTGCCGTGGTCGCGCAGGCGATGCAATACGTCCAGCAGTTGCTGGATATCGGCAAAGTGCAGGCCGGTGGTAGGTTCGTCGAGGATGTACAAGGTCTTGCCgg                                                                                 <  6:46197/149‑1 (MQ=255)
                                                               cAGTCTTGATCACATCCAGGTTATGCTCGATCACCACCACGGTATTGCCGTGGTCGCGCAGGCGATGCAATACGTCCAGCAGTTGCTGGATATCGGCAAAGTGCAGGCCGGTGGTAGGt                                                                                                             <  6:8480/119‑1 (MQ=255)
                                                               cAGTCTTGATCACATCCAGGTTATGCTCGATCACCACCACGGTATTGCCGTGGTCGCGCAGGCGATGCAATACGTCCAGCAGTTGCTGGATATCGGCAAAGTGCAGGCCGGTGGTAGGt                                                                                                             >  5:8480/1‑119 (MQ=255)
                                                                aGTCTTGATCACATCCAGGTTATGCTCGATCACCACCACGGTATTGCCGTGGTCGCGCAGGCGATGCAATACGTCCAGCAGTTGCTGGATATCGGCAAAGTGCAGGCCGGTGGTAGGTTCGTCGAGGATGTACAAGGTCTTGCCGGTAt                                                                              >  6:100183/1‑149 (MQ=255)
                                                                        tCACATCCAGTTTATGCTCTATCACCACCACGTTATTGCCGTGGTCGCGCAGGCGATGCAATACGTCCAGCAGTTGCTg                                                                                                                                            <  4:69321/79‑1 (MQ=255)
                                                                        tCACATCCAGGTTATGCTCGATCACCACCACGGTATTGCCGTGGTCGCGCAGGCGATGCAATACGTCCAGCAGTTGCTGGATATCGGCAAAGTGCAGGCCGGTGGTAGGTTCGTCGAGGATGTACAAGGTCTTGCCGGTATCGCGCTTg                                                                      <  4:210799/149‑1 (MQ=255)
                                                                        tCACATCCAGGTTATGCTCGATCACCACCACGCTATTGCCGTGGTAACGCAGGCGATGCAATACGTCCAGCAGTTGCTg                                                                                                                                            >  3:69321/1‑79 (MQ=255)
                                                                          acaTCCAGGTTATGCTCGATCACCACCACGGTATTGCCGTGGTCGCGCAGGCGATGCAATACGTCCAGCAGTTGCTGGATATCGGCAAAGTGCAGGCCGGTGGTAGGTTCGTCGAGGATGTACAAGGTCTTGCCGGTATCGCGCTTGGa                                                                    <  2:65103/149‑1 (MQ=255)
                                                                           caTCCAGGTTATGCTCGATCACCACCACGGTATTGCCGTGGTCGCGCAGGCGATGCAATACGTCCAGCAGTTGCTGGATATCGGCAAAGTGCAGGCCGGTGGTAGGTTCGTCGAGGATGTACAAGGTCTTGCCGGTATcg                                                                            >  7:199506/1‑140 (MQ=255)
                                                                           caTCCAGGTTATGCTCGATCACCACCACGGTATTGCCGTGGTCGCGCAGGCGATGCAATACGTCCAGCAGTTGCTGGATATCGGCAAAGTGCAGGCCGGTGGTAGGTTCGTCGAGGATGTACAAGGTCTTGCCGGTATcg                                                                            <  8:199506/140‑1 (MQ=255)
                                                                            aTCCAGGTTATGCTCGATCACCACCACGGTATTGCCGTGGTCGCGCAGGCGATGCAATACGTCCAGCAGTTGCTGGATATCGGCAAAGTGCAGGCCGGTGGTAGGTTCGTCGAGGATGTACAAGGTCTTGCCGGTATCGCGCTTGGa                                                                    >  4:307973/1‑147 (MQ=255)
                                                                            aTCCAGGTTATGCTCGATCACCACCACGGTATTGCCGTGGTCGCGCAGGCGATGCAATACGTCCAGCAGTTGCTGGATATCGGCAAAGTGCAGGCCGGTGGTAGGTTCGTCGAGGATGTACAAGGTCTTGCCGGTATCGCGCTTGGa                                                                    >  3:324597/1‑147 (MQ=255)
                                                                                   ttATGCTCGATCACCACCACGGTATTGCCGTGGTCGCGCAGGCGATGCAATACGTCCAGCAGTTGCTGGATATCGGCAAAGTGCAGGCCGGTGGTAGGTTCGTCGAGGATGTACAAGGTCTTGCCGGTATc                                                                             >  7:308390/1‑131 (MQ=255)
                                                                                   ttATGCTCGATCACCACCACGGTATTGCCGTGGTCGCGCAGGCGATGCAATACGTCCAGCAGTTGCTGGATATCGGCAAAGTGCAGGCCGGTGGTAGGTTCGTCGAGGATGTACAAGGTCTTGCCGGTATc                                                                             <  8:308390/131‑1 (MQ=255)
                                                                                         tCGATCACCACCACGGTATTGCCGTGGTCGCGCAGGCGATGCAATACGTCCAGCAGTTGCTGGATATCGGCAAAGTGCAGGCCGGTGGTAGGTTCGTCGAGGATGTACAAGGTCTTGCCGGTATCGCGCTTGGACAGTTCGCGCGACAg                                                     >  1:18847/1‑149 (MQ=255)
                                                                                           gATCACCACCACGGTATTGCCGTGGTCGCGCAGGCGATGCAATACGTCCAGCAGTTGCTGGATATCGGCAAAGTGCAGGCCGGTGGTAGGTTCGTCGAGGATGTACAAGGTCTTGCCGGTATCGCGCTTGGACAGTTCGCGCGACAg                                                     <  2:18847/147‑1 (MQ=255)
                                                                                           gATCACCACCACGGTATTGCCGTGGTCGCGCAGGCGATGCAATACGTCAAGCAGTTGCTGGATATCGGCAAAGTGCAGGCCGGTGGTAGGTTCGTCGAGGATGTACAAGGTCTTGCCGGTATCGCGCTTGGACAGTTCGCGAGACAGCt                                                   >  6:270847/1‑149 (MQ=255)
                                                                                               accaccacGGTATTGCCGTGGTCGCGCAGGCGATGCAATACGTCCAGCAGTTGCTGGATATCGGCAAAGTGCAGGCCGGTGGTAGGTTCGTCGAGGATGTACAAGGTCTTGCCGGTATCGCGCTTGGACAGTTCGCGCGACAGCTTTAc                                               <  2:247314/149‑1 (MQ=255)
                                                                                                ccaccacGGTATTGCCGTGGTCGCGCAGGCGATGCAATACGTCCAGCAGTTGCTGGATATCGGCAAAGTGCAGGCCGGTGGTAGGTTCGTCGAg                                                                                                     <  7:214649/94‑1 (MQ=255)
                                                                                                ccaccacGGTATTGCCGTGGTCGCGCAGGCGATGCAATACGTCCAGCAGTTGCTGGATATCGGCAAAGTGCAGGCCGGTGGTAGGTTCGTCGAg                                                                                                     >  8:214649/1‑94 (MQ=255)
                                                                                                  accacGGCAATACCGTGGTCGCGCAGGCGATGCAATACGTCCAGCAGTTGCTGGATATCGGCAAAGTGCAGGCCGGTGGTAGGTTCGTCGAGGATGTACAAGGTCTTGCCGGTATCGCGCTTGGACAGTTCGCGCGACAGCTTTACCCg                                            >  4:211633/1‑149 (MQ=255)
                                                                                                      cGGTATTGCCGTGGTCGCGCAGGCGATGCAATACGTCCAGCAGTTGCTGGATATCGGCAAAGTGCAGGCCGGTGGTAGGTTCGTCGAGGATGTACAAGGTCTTGCCGGTATCGCGCTTGGACAGTTCGCGCGACAGCTTTACCCGCTgc                                        <  5:100183/149‑1 (MQ=255)
                                                                                                        gTATTGCCGTGGTCGCGCAGGCGATGCAATACGTCCAGCAGTTGCTGGATATCGGCAAAGTGCAGGCCGGTGGTAGGTTCGTCGAGGATGTACAAGGTCTTGCCGGTATCGCGCTTGGACAGTTCGCGCGACAGCTTTACCCGCTgcgc                                      <  3:307973/149‑1 (MQ=255)
                                                                                                                                     tACGTCCAGCAGTTGCTGGATATCGGCAAAGTGCAGGCCGGTGGTAGGTTCGTCGAGGATGTACAAGGTCTTGCCGGTATCGCGCTTGGACAGTTCGCGCGACAGCTTTACCCGCTGCGCCTCACCGCCGGACAACGTGGTCGCTGAtt         >  4:166601/1‑149 (MQ=255)
                                                                                                                                         tCCAGCAGTTGCTGGATATCGGCAAAGTGCAGGCCGGTGGTAGGTTCGTCGAGGATGTACAAGGTCTTGCCGGTATCGCGCTTGGACAGTTCGCGCGACAGCTTTACCCGCTGCGCCTCACCGCCGGACAACGTGGTCGCTGATTGccc     <  4:324597/149‑1 (MQ=255)
                                                                                                                                            agcagTTGCTGGATATCGGCAAAGTGCAGGCCGGTGGTAGGTTCGTCGAGGATGTACAAGGTCTTGCCGGTATCGCGCTTGGACAGTTCGCGCGACTGCTTTACCCGCTGCGCCTCACCGCCGGACAACGTGGTCGCTGCTTGCCCCAg  >  2:126580/1‑149 (MQ=255)
                                                                                                                                             gcagTTGCTGGATATCGGCAAAGTGCAGGCCGGTGGTAGGTTCGTCGAGGATGTACAAGGTCTTGCCGGTATCGCGCTTGGACAGTTCGCGCGACAGCTTTACCCGCTgcgc                                      >  5:91394/1‑112 (MQ=255)
                                                                                                                                             gcagTTGCTGGATATCGGCAAAGTGCAGGCCGGTGGTAGGTTCGTCGAGGATGTACAAGGTCTTGCCGGTATCGCGCTTGGACAGTTCGCGCGACAGCTTTACCCGCTgcgc                                      <  6:91394/112‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
CGCAGGCAATGATCTGCCCACCCTTGGAACCGCCCTCAGGCCCCAGGTCAACCAGCCAGTCGGCAGTCTTGATCACATCCAGGTTATGCTCGATCACCACCACGGTATTGCCGTGGTCGCGCAGGCGATGCAATACGTCCAGCAGTTGCTGGATATCGGCAAAGTGCAGGCCGGTGGTAGGTTCGTCGAGGATGTACAAGGTCTTGCCGGTATCGCGCTTGGACAGTTCGCGCGACAGCTTTACCCGCTGCGCCTCACCGCCGGACAACGTGGTCGCTGATTGCCCCAG  >  NC_002947/566744‑567032

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