Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,036,373 C→T *578* (TAG→TAA recJ ← ssDNA exonuclease, 5' ‑‑> 3'‑specific

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,036,3730CT98.1% 163.6 / NA 52*578* (TAG→TAArecJssDNA exonuclease, 5' ‑‑> 3'‑specific
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (26/25);  minor base A (1/0);  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 = 1.00e+00

TGAGGTCCTGAATGCGATTATTTACCGGATTAATTTCAAACATGGTCTGATTTCTTTTATTGAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGCTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGCCGTTATCCGGCCAGAGGGCGGTATCGACATTA  >  NC_000913/3036240‑3036502
                                                                                                                                     |                                                                                                                                 
tGAGGTCCTGAATGCGATTATTTACCGGATTAATTTCAAACATGGTCTGATTTCTTTTATTGAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAAt                                                                                                                              <  2:1310580/139‑1 (MQ=255)
              cGATTATTTACCGGATTAATTTCAAACATGGTCTGATTTCTTTTATTGAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGt                                                                                                                <  1:263203/139‑1 (MQ=255)
               gATTATTTACCGGATTAATTTCAAACATGGTCTGATTTCTTTTATTGAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTgg                                                                                                                           <  2:1563036/127‑1 (MQ=255)
               gATTATTTACCGGATTAATTTCAAACATGGTCTGATTTCTTTTATTGAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTgg                                                                                                                           >  1:1563036/1‑127 (MQ=255)
                         cGGATTAATTTCAAACATGGTCTGATTTCTTTTATTGAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAAtt                                                                                                     >  2:1384728/1‑139 (MQ=255)
                           gATTAATTTCAAACATGGTCTGATTTCTTTTATTGAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTg                                                                                                   >  2:931407/1‑139 (MQ=255)
                           gATTAATTTCAAACATGGTCTGATTTCTTTTATTGAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTg                                                                                                   >  2:665172/1‑139 (MQ=255)
                                           ggTCTGATTTCTTTTATTGAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCcg                                                                                   >  2:1446477/1‑139 (MQ=255)
                                            gTCTGATTTCTTTTATTGAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCcgc                                                                                  <  2:910194/139‑1 (MQ=255)
                                            gTCTGATTTCTTTTATTGAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCcgc                                                                                  >  2:1727952/1‑139 (MQ=255)
                                            gTCTGATTTCTTTTATTGAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGATAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCcgc                                                                                  >  2:1004061/1‑139 (MQ=255)
                                              cTGATTTCTTTTATTGAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGa                                                                                <  1:429086/139‑1 (MQ=255)
                                                gATTTCTTTTATTGAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGaaa                                                                              <  1:1583867/139‑1 (MQ=255)
                                                gATTTCTTTTATTGAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGaaa                                                                              <  2:514401/139‑1 (MQ=255)
                                                                cTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTg                                                                                                   <  2:745477/102‑1 (MQ=255)
                                                                cTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTg                                                                                                   >  1:745477/1‑102 (MQ=255)
                                                                cTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCg                                                                          <  1:1777624/127‑1 (MQ=255)
                                                                cTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCg                                                                          >  2:1777624/1‑127 (MQ=255)
                                                                    tCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGTAACTCGTTGATATCGAGCTtat                                                          <  1:343955/139‑1 (MQ=255)
                                                                    tCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTtat                                                          <  1:770801/139‑1 (MQ=255)
                                                                        aaTGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGa                                                                >  1:1779236/1‑129 (MQ=255)
                                                                        aaTGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGa                                                                <  2:1779236/129‑1 (MQ=255)
                                                                             ggTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTgg                                                                                                                           <  2:1462562/65‑1 (MQ=255)
                                                                             ggTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTgg                                                                                                                           >  1:1462562/1‑65 (MQ=255)
                                                                             ggTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTg                                                                                                                            >  1:416862/1‑64 (MQ=255)
                                                                             ggTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTg                                                                                                                            <  2:416862/64‑1 (MQ=255)
                                                                             ggTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGtt                                                 >  1:1348825/1‑139 (MQ=255)
                                                                              gTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTg                                                >  1:410425/1‑139 (MQ=255)
                                                                                    agagCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTaaa                                                                                                                               <  1:333823/54‑1 (MQ=255)
                                                                                    agagCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTaaa                                                                                                                               >  2:333823/1‑54 (MQ=255)
                                                                                      agCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAGTTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTcgc                                        >  2:1241698/1‑139 (MQ=255)
                                                                                        cGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTAGTTGATATCGAGCTTATAAGCCAGTTGCACTTCgcgc                                      >  1:1147375/1‑139 (MQ=255)
                                                                                         gggATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGa                                                                >  1:968373/1‑112 (MQ=255)
                                                                                         gggATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGa                                                                <  2:968373/112‑1 (MQ=255)
                                                                                         gggATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCa                                     <  1:1727952/139‑1 (MQ=255)
                                                                                           gATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGaaa                                                                              <  1:1006141/96‑1 (MQ=255)
                                                                                           gATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGaaa                                                                              >  2:1661120/1‑96 (MQ=255)
                                                                                           gATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGaaa                                                                              >  2:1006141/1‑96 (MQ=255)
                                                                                           gATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGaaa                                                                              <  1:1661120/96‑1 (MQ=255)
                                                                                                     aaTCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCg                                                                          <  1:810523/90‑1 (MQ=255)
                                                                                                     aaTCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCg                                                                          >  2:810523/1‑90 (MQ=255)
                                                                                                     aaTCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGCCGTTATCCg                         >  1:292842/1‑139 (MQ=255)
                                                                                                         cACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGCCGTTATCCGGCCa                     <  1:1004061/139‑1 (MQ=255)
                                                                                                           cGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGCCGTTATCCGGCCaga                   <  1:1446477/139‑1 (MQ=255)
                                                                                                            gCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACt                                                                            >  2:1495061/1‑81 (MQ=255)
                                                                                                            gCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACt                                                                            <  1:1495061/81‑1 (MQ=255)
                                                                                                             ctctTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGCCGTTATCCGGCCAGAgg                 <  1:665172/139‑1 (MQ=255)
                                                                                                                ttttttATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGCCGTTATCCGGCCAGAGGGCg              >  2:190356/1‑139 (MQ=255)
                                                                                                                          agaagaTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGCCGTTATCCGGCCAGAGGGCGGTATCGACAt    >  1:271888/1‑139 (MQ=255)
                                                                                                                          agaagaTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGCCGTTATCCGGCCAGAGGGCGGTATCGACAt    >  1:1340247/1‑139 (MQ=255)
                                                                                                                           gaagaTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGCCGTTATCCGCCCAGAGGGCGGTATCGACAtt   <  1:931407/139‑1 (MQ=255)
                                                                                                                            aagaTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGCCGTTATCCGGCCAGAGGGCGGTATCGACATTa  >  2:284264/1‑139 (MQ=255)
                                                                                                                                     |                                                                                                                                 
TGAGGTCCTGAATGCGATTATTTACCGGATTAATTTCAAACATGGTCTGATTTCTTTTATTGAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGCTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGCCGTTATCCGGCCAGAGGGCGGTATCGACATTA  >  NC_000913/3036240‑3036502

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