Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,172,553 A→G 53.6% I471T (ATC→ACC)  mfd ← transcription‑repair coupling factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,172,5530AG53.6% 14.0 / 64.5 56I471T (ATC→ACC) mfdtranscription‑repair coupling factor
Reads supporting (aligned to +/- strand):  ref base A (17/9);  new base G (16/14);  total (33/23)
Fisher's exact test for biased strand distribution p-value = 4.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

CACCCGCTTCCAGCGTGGTCATTCCGGCATAACGACCGACGCCGTGCTCCAGATGGACCACCGGCTGACCAATATGCAGTTCCGCAAGGTTACGGATCAGTGTATCGGGGTTGATGGTGCGGCGAGAATCCTGACGACGACGGGCAACGCGTTCACCGAGCAGATCGCTTTCGCAAATCAGC  >  NC_000913/1172458‑1172639
                                                                                               |                                                                                      
cACCCGCTTCCAGCGTGGTCATTCCGGCATAACGACCGACGCCGTGCTCCAGATGGACCACCGGCTGACCAATATGCAGTTCCGCAAGGTTACGGGTCAg                                                                                    >  1:639432/1‑100 (MQ=255)
 aCCCGCTTCCAGCGTGGTCATTCCGGCATAACGACCGACGCCGTGCTCCAGATGGACCACCGGCTGACCAATATGCAGTTCCGCAAGGTTACGGATCAgt                                                                                   <  1:50079/100‑1 (MQ=255)
  cccGCTTCCAGCGTGGTCATTCCGGCATAACGACCGACGCCGTGCTCCAGATGGACCACCGGCTGACCAATATGCAGTTCCGCAAGGTTACGGGTCAgtg                                                                                  <  2:259216/100‑1 (MQ=255)
   ccGCTTCCAGCGTGGTGATTCCGGCATAACGACCGACGCCGTGCTCCAGATGGACCACCGGCTGACCAATATGCAGTTCCGCAAGGTTACGGGTCAgt                                                                                   <  2:2140087/98‑1 (MQ=255)
   ccGCTTCCAGCGTGGTCATTCCGGCATAACGACCGACGCCGTGCTCCAGATGGACCACCGGCTGACCAATATTCAGTTCCGCAAGGTTACGGGTCAgtg                                                                                  <  2:1541031/99‑1 (MQ=255)
     gCTTCCAGCGTGGTCATTCCGGCATAACGACCGACGCCGTGCTCCAGATGGACCACCGGCTGACCAATATGCAGTTCCGCAAGGTTACGGATCAGTGTa                                                                                >  1:639852/1‑99 (MQ=255)
      cTTCCAGCGTGGTCATTCCGGCATAACGACCGACGCCGTGCTCCAGATGGACCACCGGCTGACCAATATGCAGTTCCGCAAGGTTACGGATCAGTGTATc                                                                              >  2:1735093/1‑100 (MQ=255)
        tCCAGCGTGGTCATTCCGGCATAACGACCGACGCCGTGCTCCAGATGGACCACCGGCTGACCAATATGCAGTTCCGCAAGGTTACGGGTCAGTGTATCgg                                                                            >  1:852255/1‑100 (MQ=255)
         ccAGCGTGGTCATTCCGGCATAACGACCGACGCCGTGCTCCAGATGGACCACCGGCTGACCAATATGCAGTTCCGCAAGGTTACGGATCAGTGTATCggg                                                                           >  1:1324489/1‑100 (MQ=255)
           aGCGTGGTCATTCCGGCATAACGACCGACGCCGTGCTCCAGATGGACCACCGGCTGACCAATATGCAGTTCCGCAAGGTTACGGGTCACTGTATCGGGGt                                                                         <  2:32436/100‑1 (MQ=255)
             cGTGGTCATTCCGGCATAACGACCGACGCCGTGCTCCAGATGGACCACCGGCTGACCAATATGCAGTTCCGCAAGGTTACGGGTCAGTGTATCGGGGTTg                                                                       <  2:289515/100‑1 (MQ=255)
               tGGTGATTCCGCCATAACGACCGACGGCGGGCGCCAGATGGGTCATCGGATGACTAATATGCCGTTCAGGAAAGTTAGGGATGAgtgaaatgggg                                                                          <  2:1339192/95‑9 (MQ=255)
                ggTCATTCCGGCATAACGACCGACGCCGTGCTCCAGATGGGCCACCGGCTGACCAATATGCAGTTCCGCAAGGTTACGGGTCAGTGTATCgggg                                                                          <  2:1444632/94‑1 (MQ=255)
                 gTCATTCCGGCATAACGACCGACGCCGTGCTCCAGATGGACCACCGGCTGACCAATATGCAGTTCCGCAAGGTTACGGATCAGTGTATCGGGGTTGATgg                                                                   <  2:853395/100‑1 (MQ=255)
                  tCATTCCGGCATAACGACCGACGCCGTGCTCCAGATGGACCACCGGCTGACCAATATGCAGTTCCGCAAGGTTACGGGTCAGTGTATCGGGGTTGATGGt                                                                  >  1:289231/1‑100 (MQ=255)
                   cATTCCGGCATAACGACCGACGCCGTGCTCCAGATGGACCACCGGCTGACCAATATGCAGTTCCGCAAGGTTACGGGTCAGTGTATCGGGGTTGATGGt                                                                  >  1:793746/1‑99 (MQ=255)
                       ccGGCATAACGACCGACGCCGTGCTCCAGATGGACCACCGGCTGACCTATATGCAGTTCCGCAAGGTTACGGGTCAGTGTGTCGGGGTTGATGGTgcgg                                                              >  2:2025652/1‑99 (MQ=255)
                        cGGCATAACGACCGACGCCGTGCTCCAGATGGACCACCGGCTGACCAATATGCAGTTCCGCAAGGTTACGGGTCAGTGTATCGGGGTTGATGGTgcggcg                                                            <  2:859968/100‑1 (MQ=255)
                           cATAACGACCGACGCCGTGCTCCAGATGGACCACCGGCTGACCAATATGCAGTTCCGCAAGGTTACGGGTCAGTGTATCGGGGTTGATGGTGCGGCgaga                                                         <  1:941352/100‑1 (MQ=255)
                            aTAACGGCCGACGCCGTGCTCAAGATAGACCACCGGCTGACCAATATGCAGTTCTGCGAGGTTACGGATCAGTGTATCGGGGTTGATGGTggggcg                                                            >  2:2282411/1‑96 (MQ=255)
                             tAACGACCGACGCCGTGCTCCAGATGGACCACCGGCTGACCAATATGCAGTTCCGCAAGGTTACGGATCAGTGTATCGGGGTTGATGGTGCGGCGAGAAt                                                       >  2:461732/1‑100 (MQ=255)
                                cgaccgacGCCGTGCTCCAGATGGACCACCGGCTGACCAATATGCAGTTCCGCAAGGTTACGGATCAGTGTATCGGGGTTGATGGTGCGGCGAGAATCCt                                                    <  2:410340/100‑1 (MQ=255)
                                 gaccgacGTCGTGCTCCAGATGGACCACCGGCTGACCAATATGCAGTTCCGCAAGGTTACGGGTCAGTGTATCGGGGTTGATGGTGCGGCGAGAATcc                                                     <  1:705911/98‑1 (MQ=255)
                                     gacGCCGTGCTCCAGATGGACCACCGGCTGACCAATATGCAGTTCCGCAAGGTTACGGGTCAGTGTATCGGGGTTGATGGTGCGGCGAGAATCCTgacga                                               >  2:113399/1‑100 (MQ=255)
                                       cGCCGTGCTCCAGATGGACCACCGGCTGACCAATATGCAGTTCCGCAAGGTTACGGATCAGTGTATCGGGGTTGATGGTGCGGCGAGAATCCTgacgacg                                             >  1:1783563/1‑100 (MQ=255)
                                       cGCCGTGCTCCAGATGGACCACCCGCAGACCAATATGCAGTTGCGCCAGGTTACGGATCATTGTATCGGGGTTGATGGTGCGGCGAGAAtacaga                                                  >  2:2324504/1‑90 (MQ=255)
                                         ccGTGCTCCAGATGGACCACCGGCTGACCAATATGCAGTTCCGCAAGGTTACGGATCAGTGTATCGGGGTTGATGGTGCGGCGAGAATCCTgacgacgac                                           >  1:1752458/1‑100 (MQ=255)
                                          cGTGCTCCAGATGGACCACCGGCTGACCAATATGCAGTTCCGCAAGGTTACGGGTCAGTGTATCGGGGTTGATGGTGCGGCGAGAATCCTGacgacgacg                                          >  1:1061929/1‑100 (MQ=255)
                                            tGCTCCAGATGGACCACCGGCTGACCAATATGCAGTTCCGCAAGGTTACGGATCAGTGTATCGGGGTTGATGGTGCGGCGAGAATCCTGACGACGACggg                                        >  1:1170186/1‑100 (MQ=255)
                                             gCTCCAGATGGACCACCGGCTGACCAATATGCAGTTCCGCAAGGTTACGGATCAGTGTATCGGGGTTGATGGTGCGGCGAGAATCCTGACGACGACGGGc                                       <  1:1981993/100‑1 (MQ=255)
                                                ccAGATGGCCAACCGGTTGACCAATATGCAGTTCCGCAAGGTTACGGATCAGTGTATAGGGGTTGATGGTGCGGCGAGAATCCTGACGACGACGGTCAAc                                    >  2:2256629/1‑100 (MQ=255)
                                                ccAGATGGACCACCGGCTGACCAATATGCAGTTCCGCAAGGTTACGGGTCAGTGTATCGGGGTTGATGGTGCGGCGAGAATCCTGACGACGACGGGCAAc                                    >  1:126893/1‑100 (MQ=255)
                                                ccAGATGGACCACCGGCTGACCAATATGCAGTTCCGCAAGGTTACGGATCAGTGTATCGGGGTTGATGGTGCGGCGAGAATCCTGACGACGACGGGCAAc                                    <  1:142275/100‑1 (MQ=255)
                                                 cAGATGGACCACCGGCTGACCAATATGCAGTTCCGCAAGGTTACGGATCAGTGTATCGGGGTTGATGGTGCGGCGAGAATCCTGACGACGACGGGCAAcg                                   <  1:1483357/100‑1 (MQ=255)
                                                   gATGGACCACCGGCTGACCAATATGCAGTTCCGCAAGGTTACGGATCAGTGTATCGGGGTTGATGGTGCGGCGAGAATCCTGACGACGACGGGCAAcgcg                                 >  1:1854653/1‑100 (MQ=255)
                                                   gATGGACCACCGGCTGACCAATATGCAGTTCCGCAAGGTTACGGATCAGTGTATCGGGGTTGATGGTGCGGCGAGAATCCTGACGACGACGGGCAAcgcg                                 >  2:569035/1‑100 (MQ=255)
                                                       gACCACCGGCTGACCAATATGCAGTTCCGCAAGGTTACGGGCCAGTGTATCGGGGTTGATGGTGCGGCGAGAATCCTGCCGACGACGCGCAACGCGTTCa                             >  2:2301584/1‑100 (MQ=255)
                                                        accaccGGCTGACCAATATGCAGTTCCGCAAGGTTACGGATCAGTGTATCGGGGTTGATGGTGCGGCGAGAATCCTGACGACGACGGGCAACGCGTTCAc                            <  2:2006875/100‑1 (MQ=255)
                                                         ccaccGGCTGACCAATATGCAGTTCCGCAAGGTTACGGGTCAGTGTATCGGGGTTGATGGTGCGGCGAGAATCCTGACGACGACGGGCAACGCGTTCAcc                           >  2:1777152/1‑100 (MQ=255)
                                                           gccGGCTGACCAATATGCAGTTCCGCAAGGTTACGGGTCAGTGTATCGGGGTTGATGGTGCGGCGAGAATCCTGACGACGACGGGCAACGCGTTCACCGa                         >  1:1129/2‑100 (MQ=255)
                                                           accGGCTGACCAATATGCAGTTCCGCAAGGTTACGGGTCAGTGTATCGGGGTTGATGGTGCGGCGAGAATCCTGACGACGACGGGCAACGCGTTCACCGa                         >  1:1378590/1‑100 (MQ=255)
                                                           accGGCTGACCAATATGCAGTTCCGCAAGGTTACGGATCAGTGTATCGGGGTTGATGGTGCGGCGAGAATCCTGACGACGACGGGCAACGCGTTCACCGa                         >  2:905173/1‑100 (MQ=255)
                                                           accGGCTGACCAATATGCAGTTCCGCAAGGTTACGGATCAGTGTATCGGGGTTGATGGTGCGGCGAGAATCATGACGACGACGGGCAACGCGTTCACCGa                         >  2:570175/1‑100 (MQ=255)
                                                            ccGGCTGACCAATATGCAGTTCCGCAAGGTTACGGATCAGTGTATCGGGGTTGATGGTGCGGCGAGAATCCTGACGACGACGGGCAACGCGTTCACCGAg                        >  2:1310340/1‑100 (MQ=255)
                                                             cGGCTGACCAATATGCAGTTCCGCAAGGTTACGGGTCAGTGTATCGGGGTTGATGGTGCGGCGAGAATCCTGACGACGACGGGCAACGCGTTCACCGAGc                       <  1:1324640/100‑1 (MQ=255)
                                                              ggCTGACCAATATGCAGTTCCGCAAGGTTACGGGTCAGTGTATCGGGGTTGATGGTGCGGCGAGAATCCTGACGACGACGGGCAACGCGTTCACCGAGCa                      <  2:1129/100‑1 (MQ=255)
                                                                 tGACCAATATGCAGTTCCGCAAGGTTACGGGTCAGTGTATCGGGGTTGATGGTGCGGCGAGAATCCTGACGACGACGGGCAACGCGTTCACCGAGCAGAt                   <  2:1228105/100‑1 (MQ=255)
                                                                  gACCAATATGCAGTTCCGCAAGGTTACGGGTCAGTGTATCGGGGTTGATGGTGCGGCGAGAATCCTGACGACGACGGGCAACGCGTTCACCGAGCAGATc                  >  2:433096/1‑100 (MQ=255)
                                                                  gACCAATATGCAGTTCCGCAAGGTTACGGATCAGTGTATCGGGGTTGATGGTGCGGCGAGAATCCTGACGACGACGGGCAACGCGTTCACCGAGCAGATc                  >  2:1521499/1‑100 (MQ=255)
                                                                   aCCAATATGCAGTTCCGCAAGGTTACGGGTCAGTGTATCGGGGTTGATGGTGCGGCGAGAATCCTGACGACGACGGGCAACGCGTTCACCGAGCAGATc                  <  2:1470771/99‑1 (MQ=255)
                                                                           gCAGTTCCGCAAGGTTACGGATCAGTGTATCGGGGTTGATGGTGCGGCGAGAATCCTGACGACGACGGGCATCGCGTTCACCGAGCAGATCGCTTTCg           <  2:276607/98‑1 (MQ=255)
                                                                              gTTCCGCAAGGTTACGGGTCAGTGTATCGGGGTTGATGGTGCGGCGAGAATCCTGACGACGACGGGCAACGCGTTCACCGAGCAGATCGCTTTCGCaaa       >  1:103990/1‑99 (MQ=255)
                                                                              gTTCCGCAAGGTTACGGGTCAGTGTATCGGGGTTGATGGTGCGGCGAGAATCCTGACGACGACGGGCAACGCGTTCACCGAGCAGATCGCTTTCGCAAAt      >  2:593860/1‑100 (MQ=255)
                                                                               ttCCGCAAGGTTACGGGTCAGTGTATCGGGGTTGATGGTGCGGCGAGAATCCTGACGACGACGGGCAACGCGTTCACCGAGCAGATCGCTTTCGCAAATc     <  2:1168684/100‑1 (MQ=255)
                                                                                 ccGCAAGGTTACGGATCAGTGTATCGGGGTTGATGGTGCGGCGAGAATCCTGACGACGACGGGCAACGCGTTCACCGAGCAGATCGCTTTCGCAAATCAg   >  1:1124592/1‑100 (MQ=255)
                                                                                  cGCAAGGTTACGGGTCAGTGTATCGGGGTTGATGGTGCGGCGAGAATCCTGACGACGACGGGCAACGCGTTCACCGAGCAGATCGCTTTCGCAAATCAgc  >  1:1130858/1‑100 (MQ=255)
                                                                                               |                                                                                      
CACCCGCTTCCAGCGTGGTCATTCCGGCATAACGACCGACGCCGTGCTCCAGATGGACCACCGGCTGACCAATATGCAGTTCCGCAAGGTTACGGATCAGTGTATCGGGGTTGATGGTGCGGCGAGAATCCTGACGACGACGGGCAACGCGTTCACCGAGCAGATCGCTTTCGCAAATCAGC  >  NC_000913/1172458‑1172639

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