Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,414,756 (T)6→5 coding (38/810 nt) recT ← Rac prophage; recombination and repair protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,414,7510T.96.4% 234.4 / 1.2 55coding (43/810 nt)recTRac prophage; recombination and repair protein
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base . (27/26);  total (29/26)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGAGCTGCTGCCAGTTGCTCTTTCATTGATGGCTGGTTAATAAAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCATATTCCTTAGCCCAGCGGGGCAGTGATAATGTCTTAATAGCTGGCCA  >  NC_000913/1414656‑1414845
                                                                                               |                                                                                              
aaGAGCTGCTGCCAGTTGCTCTTTCATTGATGGCTGGTTAATAAAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGttttt                                                                                            <  2:754106/100‑1 (MQ=255)
 agagCTGCTGCCAGTTGCTCTTTCATTGATGGCTGGTTAATAAAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGtttttg                                                                                           <  1:714778/100‑2 (MQ=255)
  gagCTGCTGCCAGTTGCTATTTCATTGATGGCTGGTTCATCCCACTAATCACGTCGCTATTTTTCCCTGCTGCTGGTGCACGGTTTCCCTGAGtttttgt                                                                                          <  2:107773/100‑3 (MQ=255)
   agCTGCTGCCAGTTGCTCTTTCATTGATGGCTGGTTAATAAAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTTGc                                                                                         >  1:701249/1‑100 (MQ=255)
     ctgctgCCAGTTGCTCTTTCATTGATGGCTGGTTAATAATACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGtttttgc                                                                                          <  1:355419/97‑3 (MQ=255)
        ctgcCAGTTGCTCTTTCATTGATGGCTGGTTAATAAAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTGCAGATCg                                                                                   >  2:1034748/1‑100 (MQ=255)
          gcCAGTTGCTCTTTCATTGATGGCTGGTTAATAAAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGc                                                                                 >  1:155313/1‑100 (MQ=255)
          gcCAGTTGCTCTTTCATTGATGGCTGGTTAATAAAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGc                                                                                 >  1:766942/1‑100 (MQ=255)
          gcCAGTTGCTCTTTCATTGATGGCTGGTTAATAAAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGc                                                                                 >  1:339611/1‑100 (MQ=255)
           ccAGTTGCTCTTTCATTGATGGCTGGTTAATAAAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCt                                                                                <  2:898643/100‑1 (MQ=255)
             aGTTGCTCTTTCATTGATGGCTGGTTAATAAAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCttt                                                                              >  1:1150878/1‑100 (MQ=255)
                  ctctTTCATTGATGGCTGGTTAATAAAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGa                                                                         >  2:831113/1‑100 (MQ=255)
                   tctTTCATTGATGGCTGGTTAATAAAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGAt                                                                        >  1:191224/1‑100 (MQ=255)
                    ctTTGAGTGATGGCTGGTTAATAAATCTACTCACGTTGCTATTTATGACTGCTCCTAGTGCACGGTTTACATGCGATTTTGCAAATCGTTTTTTGGGagt                                                                       <  2:265027/100‑3 (MQ=255)
                    ctTTCATTGATGGCTGGTTAATAAAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGAtt                                                                       >  2:1102964/1‑100 (MQ=255)
                     tttCATTGATGGCTGGTTAATAAAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGATtg                                                                      <  1:521626/100‑1 (MQ=255)
                     tttCATTGATGGCTGGTTAATAAAACTAATCACGGCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGATtg                                                                      <  1:400336/100‑1 (MQ=255)
                       tCATTGATGGCTGGTTAATAAAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGATtggt                                                                    <  1:285271/100‑1 (MQ=255)
                         aTTGATGGCTGGTTAATAAAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGATtggtgg                                                                  <  2:736977/100‑1 (MQ=255)
                         aTTGATGGCTGGTTAATAAAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGATtggtgg                                                                  <  2:725087/100‑1 (MQ=255)
                         aTTGATGGCTGGTTAATAAAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGATtggtgg                                                                  >  2:626127/1‑100 (MQ=255)
                           tGATGGCTGGTTAATAAAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGATTGGTGGtt                                                                >  1:1079096/1‑100 (MQ=255)
                               ggCTGGTTAATAAAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCtt                                                            >  2:815139/1‑100 (MQ=255)
                                gCTGGTTAATAAAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTa                                                           >  1:1031308/1‑100 (MQ=255)
                                    gTTAATAAAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCa                                                       <  1:998890/100‑1 (MQ=255)
                                        ataaAACTAATCACGTCGCTATTTTTAACTGCTTCTGGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTg                                                   <  2:363356/100‑1 (MQ=255)
                                        ataaAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTg                                                   >  1:494426/1‑100 (MQ=255)
                                        ataaAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTg                                                   >  2:778736/1‑100 (MQ=255)
                                          tctaCTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCa                                                 <  1:210473/97‑1 (MQ=255)
                                          aaaaCTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGc                                                  >  2:562853/1‑99 (MQ=255)
                                          aaaaCTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCa                                                 <  1:23070/100‑1 (MQ=255)
                                          aaaaCTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCa                                                 >  2:612995/1‑100 (MQ=255)
                                           aaaCTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCat                                                >  1:1136003/1‑100 (MQ=255)
                                            aaCTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCata                                               >  2:646609/1‑100 (MQ=255)
                                             aCTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTGCATTTTGTTTTTTGCGATTGGTGGTTGCTTAGTCATTTGCaaga                                              >  2:210471/1‑97 (MQ=255)
                                                aaTCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCATATTcc                                           >  1:517620/1‑100 (MQ=255)
                                                    aCGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCATATTCCTTAg                                       <  2:682710/100‑1 (MQ=255)
                                                     cGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCATATTCCTTAg                                       <  2:1239861/99‑1 (MQ=255)
                                                      gTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCATATTCCTTAGcc                                     <  2:100527/100‑1 (MQ=255)
                                                           tATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCATATTCCTTAGCCCAGc                                 >  1:216075/1‑100 (MQ=255)
                                                               tttAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCATATTCCTTAGCCCAGCGGGGc                            <  1:1034758/100‑1 (MQ=255)
                                                                   aCTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCATATTCCTTAGCCCAGCGGGGCAGTg                        >  1:545736/1‑100 (MQ=255)
                                                                   aCTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCATATTCCTTAGCCCAGCGGGGCAGTg                        <  1:626131/100‑1 (MQ=255)
                                                                     tgctgctgGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCATATTCCTTAGCCCAGCGGGGCAGTGAt                      <  2:1082487/100‑1 (MQ=255)
                                                                        tgctgGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCATATTCCTTAGCCCAGCGGGGCAGTGATAAt                   >  2:264370/1‑100 (MQ=255)
                                                                         gctgGTGGACGGTTTCCGGGAGAGTTTGCAGATCGGCTTGTGCGATTGGTGGTTGCTTAGTAACATGCGTATCCCTTAGCACATCGGGGCAGtcctcca                   >  2:1199806/1‑93 (MQ=255)
                                                                         gctgGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCATATTCCTTAGCCCAGCGGGGCAGTGATAATg                  >  2:570923/1‑100 (MQ=255)
                                                                            ggTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCATATTCCTTAGCCCAGCGGGGCAGTGATAATGTCt               <  1:831119/100‑1 (MQ=255)
                                                                             gTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCATATTCCTTAGCCCAGCGGGGCAGTGATAATGTCtt              <  1:624934/100‑1 (MQ=255)
                                                                              tggatcgTTTCTCTGAGTTCTTGCAGATCTGCTTTTGCGATTTGTGGTTGCTTAGTCATTTGCATATTCCTTAGCCCAGCGTGGCAGTGATAATGTCTTa             <  2:776656/94‑1 (MQ=255)
                                                                                 aCGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCATATTCCTTAGCCCAGCGGGGCAGTGATAATGTCTTAATa          <  1:41504/100‑1 (MQ=255)
                                                                                  cGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCATATTCCTTAGCCCAGCGGGGCAGTGATAATGTCTTAATAg         <  1:654374/100‑1 (MQ=255)
                                                                                  cGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCATATTCCTTAGCCCAGCGGGGCAGTGATAATGTCTTAATAg         >  1:207005/1‑100 (MQ=255)
                                                                                    gTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCATATTCCTTTGCGCAGCGGGGCAGTGATAATGTCTTAATAGCt       <  2:854906/100‑1 (MQ=255)
                                                                                     tttCCCTGAGTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCATATTCCTTAGCCCAGCGGGGCAGTGATAATGTCTTAATAGCTg      <  1:1220453/100‑1 (MQ=255)
                                                                                      ttCCCTGAGTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCATATTCCTTAGCCCAGCGGGGCAGTGATAATGTCTTAATAGCTgg     <  2:766937/100‑1 (MQ=255)
                                                                                        cccTGAGTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCATATTCCTTAGCCCAGCGGGGCAGTGATAATGTCTTAATAGCTGGcc   >  2:1234486/1‑100 (MQ=255)
                                                                                        cccTGAGTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCATATTCCTTAGCCCAGCGGGGCAGTGATAATGTCTTAATAGCTGGcc   >  2:880677/1‑100 (MQ=255)
                                                                                         ccTGAGTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCATATTCCTTAGCCCAGCGGGGCAGTGATAATGTCTTAATAGCTGGCCa  <  2:155311/100‑1 (MQ=255)
                                                                                               |                                                                                              
AAGAGCTGCTGCCAGTTGCTCTTTCATTGATGGCTGGTTAATAAAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCATATTCCTTAGCCCAGCGGGGCAGTGATAATGTCTTAATAGCTGGCCA  >  NC_000913/1414656‑1414845

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