Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,578,741 A→G 11.7% N18D (AAC→GAC)  PP_5589 → rhs repeat‑associated core domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,578,7410AG11.7% 86.7 / 3.3 34N18D (AAC→GAC) PP_5589rhs repeat‑associated core domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (14/16);  new base G (2/2);  total (16/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.81e-01

GCCATTGCTCCGTGTTGAAGTACAATGGTGATCATTTAGACGATCTCGTTAACGGGTATCATTTAGGGGTAGGTAAACGATTGTATCTGCCAGCTTTGATGCGCTTCGCATCGCCCGATGCCTTGAGCCCGTTTTTGAAAGGTGGTATCAACAGCTATTCTTACTGCGCCGGAGATCCTGTAAATTACTCAGACCCTACAGGGCAGGCCCCAGAGTCGCCGAAGCTAGTAAGAAAATCATCACAGACAGCCCCGTCGAAACCCTCATCTGGCTCAT  >  NC_002947/3578592‑3578867
                                                                                                                                                     |                                                                                                                              
gCCATTGCTCCGTGTTGAAGTACAATGGTGATCATTTAGACGATCTCGTTAACGGGTATCATTTAGGGGTAGGTAAACGATTGTATCTGCCAGCTTTGATGCGCTTCGCATCGCCCGATGCCTTGAGCCCGTTTTTGAAAGGTGGTATCa                                                                                                                                >  1:1088012/1‑150 (MQ=255)
             gtTGAAGTACAATGGTGATCATTTAGACGATCTCGTTAACGGGTATCATTTAGGGGTAGGTAAACGATTGTATCTGCCAGCTTTGATGCGCTTCGCATCGCCCGATGCCTTGAGCCCGTTTTTGAAAGGTGGTATCAACAGCTATTCTTa                                                                                                                   >  2:1121920/1‑150 (MQ=255)
               tGAAGTACAATGGTGATCATTTAGACGATCTCGTTAACGGGTATCATTTAGGGGTAGGTAAACGATTGTATCTGCCAGCTTTGATGCGCTTCGCATCGCCCGATGCCTTGAGCCCGTTTTTGAAAGGTGGTATCAACAGCTATTCTTACt                                                                                                                 <  1:1121920/150‑1 (MQ=255)
                     aCAATGGTGATCATTTAGACGATCTCGTTAACGGGTATCATTTAGGGGTATGTAAACGATTGTATCTGCCAGCTTTGATGCGCTTCGCATCGCCCGATGCCTTGAGCCCGTTTTTGAAAGGTGGTATCAACAGCTATTCTTACTGCGCCg                                                                                                           <  2:256196/150‑1 (MQ=255)
                       aaTGGTGATCATTTAGACGATCTCGTTAACGGGTATCATTTAGGGGTAGGTAAACGATTGTATCTGCCAGCTTTGATGCGCTTCGCATCGCCCGATGCCTTGAGCCCGTTTTTGAAAGGTGGTATCAACAGCTATTCTTACTGCGCCGga                                                                                                         <  2:484540/150‑1 (MQ=255)
                                      gACGATCTCGTTAACGGGTATCATTTAGGGGTAGGTAAACGATTGTATCTGCCAGCTTTGATGCGCTTCGCATCGCCCGATGCCTTGAGCCCGTTTTTGAAAGGTGGTATCAACAGCTATTCTTACTGCGCCGGAGATCCTGTaaa                                                                                              >  1:777712/1‑146 (MQ=255)
                                      gACGATCTCGTTAACGGGTATCATTTAGGGGTAGGTAAACGATTGTATCTGCCAGCTTTGATGCGCTTCGCATCGCCCGATGCCTTGAGCCCGTTTTTGAAAGGTGGTATCAACAGCTATTCTTACTGCGCCGGAGATCCTGTaaa                                                                                              <  2:777712/146‑1 (MQ=255)
                                      gACGATCTCGTTAACGGGTATCATTTAGGGGTAGGTAAACGATTGTATCTGCCAGCTTTGATGCGCTTCGCATCGCCCGATGCCTTGAGCCCGTTTTTGAAAGGTGGTATCAACAGCTATTCTTACTGCGCCGGAGATCCTGTAAATTAc                                                                                          >  2:1282628/1‑150 (MQ=255)
                                           tctcGTTAACGGGTATCATTTAGGGGTAGGTAAACGATTGTATCTGCCAGCTTTGATGCGCTTCGCATCGCCCGATGCCTTGAGCCCGTTTTTGAAAGGTGGTATCAACAGCTATTCTTACTGCGCCGGAGATCCTGTAAATTACTCAGa                                                                                     <  2:1088012/150‑1 (MQ=255)
                                            ctcGTTAACGGGTATCATTTAGGGGTAGGTAAACGATTGTATCTGCCAGCTTTGATGCGCTTCGCATCGCCCGATGCCTTGAGCCCGTTTTTGAAAGGTGGTATCAACAGCTATTCTTACTGCGCCg                                                                                                           <  2:686772/127‑1 (MQ=255)
                                            ctcGTTAACGGGTATCATTTAGGGGTAGGTAAACGATTGTATCTGCCAGCTTTGATGCGCTTCGCATCGCCCGATGCCTTGAGCCCGTTTTTGAAAGGTGGTATCAACAGCTATTCTTACTGCGCCg                                                                                                           >  1:686772/1‑127 (MQ=255)
                                                  aaCGGGTATCATTTAGGGGTAGGTAAACGATTGTATCTGCCAGCTTTGATGCGCTTCGCATCGCCCGATGCCTTGAGCCCGTTTTTGAAAGGTGGTATCAACAGCTATTCTTACTGCGCCGGAGATCCTGTAAATTACTCAGACCCTACa                                                                              <  1:175072/150‑1 (MQ=255)
                                                        tATCATTTAGGGGTAGGTAAACGATTGTATCTGCCAGCTTTGATGCGCTTCGCATCGCCCGATGCCTTGAGCCCGTTTTTGAAAGGTGGTATCAACAGCTATTCTTACTGCGCCGGAGATCCTGTAAATTACTCAGACCCTACAGGGCAg                                                                        >  1:728069/1‑150 (MQ=255)
                                                                  gggTAGGTAAACGATTGTATCTGCCAGCTTTGATGCGCTTCGCATCGCCCGATGCCTTGAGCCCGTTTTTGAAAGGAGGTATCAACAGCTATTCTTACTGCGCCGGAGATCCTGTAAATTACTCAGACCCTACAGGGCAGGCACCAGAGt                                                              >  1:1283044/1‑150 (MQ=255)
                                                                    gtaggtaAACGATTGTATCTGCCAGCTTTGATGCGCTTCGCATCGCCCGATGCCTTGAGCCCGTTTTTGAAAGGTGGTATCAACAGCTATTCTTACTGCGCCGGAGATCCTGTAAATTACTCAGACCCTACAGGGCAGGCCCCAGAGTCg                                                            >  2:1047683/1‑150 (MQ=255)
                                                                         taAACGATTGTATCTGCCAGCTTTGATGCGCTTCGCATCGCCCGATGCCTTGAGCCCGTTTTTGAAAGGTGGTATCAACAGCTATTCTTACTGCGCCGGAGATCCTGTAAATTACTCAGACCCTACAGGGCAGGCCCCAGAGTCGCCGaa                                                       >  2:466208/1‑150 (MQ=255)
                                                                           aaCGATTGTATCTGCCAGCTTTGATGCGCTTCGCATCGCCCGATGCCTTGAGCCCGTTTTTGAAAGGTGGTATCAACAGCTATTCTTACTGCGCCGGAGATCCTGTAAATTACTCAGACCCTACAGGGCAGGCCCCAGAGTCGCCGAAGc                                                     >  1:974576/1‑150 (MQ=255)
                                                                             cGATTGTATCTGCCAGCTTTGATGCGCTTCGCATCGCCCGATGCCTTGAGCCCGTTTTTGAAAGGTGGTATCAACAGCTATTCTTACTGCGCCGGAGATCCTGTAAATTACTCAGACCCTACAGGGCAGGCCCCAGAGTCGCCGAAGCTa                                                   <  1:1282628/150‑1 (MQ=255)
                                                                                 tGTATCTGCCAGCTTTGATGCGCTTCGCATCGCCCGATGCCTTGAGCCCGTTTTTGAAAGGTGGTATCAACAGCTATTCTTACTGCGCCGGAGATCCTGTAAATTACTCAGACCCTACAGGGCAGGCCCCAGAGTCGCCGAAGCTAGTaa                                               <  1:1047683/150‑1 (MQ=255)
                                                                                       tGCCAGCTTTGATGCGCTTCGCATCGCCCGATGCCTTGAGCCCGTTTTTGAAAGGTGGTATCAACAGCTATTCTTACTGCGCCGGAGATCCTGTAAATTACTCAGACCCTACAGGGCAGGccc                                                                    <  1:1189649/123‑1 (MQ=255)
                                                                                       tGCCAGCTTTGATGCGCTTCGCATCGCCCGATGCCTTGAGCCCGTTTTTGAAAGGTGGTATCAACAGCTATTCTTACTGCGCCGGAGATCCTGTAAATTACTCAGACCCTACAGGGCAGGccc                                                                    >  2:1189649/1‑123 (MQ=255)
                                                                                        gCCAGCTTTGATGCGCTTCGCATCGCCCGATGCCTTGAGCCCGTTTTTGAAAGGTGGTATCAACAGCTATTCTTACTGCGCCGGAGATCCTGTAAATTACTCAGACCCTACAGGGCAGGCCCCAGAGTCGc                                                           <  1:1059794/131‑1 (MQ=255)
                                                                                        gCCAGCTTTGATGCGCTTCGCATCGCCCGATGCCTTGAGCCCGTTTTTGAAAGGTGGTATCAACAGCTATTCTTACTGCGCCGGAGATCCTGTAAATTACTCAGACCCTACAGGGCAGGCCCCAGAGTCGc                                                           >  2:1059794/1‑131 (MQ=255)
                                                                                                tGATGCGCTTCGCATCGCCCGATGCCTTGAGCCCGTTTTTGAAAGGTGGTATCAACAGCTATTCTTACTGCGCCGGAGATCCTGTAAATTACTCAGACCCTACAGGGCAGGCCCCAGAGTCGCCGAAGCTAGTAAGAAAATCATCacaga                                <  2:728069/150‑1 (MQ=255)
                                                                                                        ttCGCATCGCCCGATGCCTTGAGCCCGTTTTTGAAAGGTGGTATCAACAGCTATTCTTACTGCGCCGGAGATCCTGTAAATTACTCAGACCCTACAGGGCAGGCCCCAGAGTCGCCGAAGCTAGTAAGAAAATCATCACAGACAGCCCCg                        >  1:885233/1‑150 (MQ=255)
                                                                                                             aTCGCCCGATGCCTTGAGCCCGTTTTTGAAAGGTGGTATCAACAGCTATTCTTACTGCGCCGGAGATCCTGTAAATTACTCAGACCCTACAGGGCAGGCCCCAGAGTCGCCGAAGCTAGTAAGAAAATCATCACAGACAGCCCCGTCGaa                   >  2:939459/1‑150 (MQ=255)
                                                                                                                gCCCGATGCCTTGAGCCCGTTTTTGAAAGGTGGTATCAACAGCTATTCTTACTGCGCCGGAGATCCTGTAAATTACTCAGACCCTACAGGGCAGGCCCCAGAGTCGCCGAAGCTAGTAAGAAAATCATCACAGACAGCCCCGTCGAAAcc                <  2:974576/150‑1 (MQ=255)
                                                                                                                   cGATGCCTTGAGCCCGTTTTTGAAAGGTGGTATCGACAGCTATTCTTACTGCGCCGGAGATCCTGTAAATTACTCAGACCCTACAGGGCAGGCCCCAGAGTCGc                                                           <  2:137783/104‑1 (MQ=255)
                                                                                                                   cGATGCCTTGAGCCCGTTTTTGAAAGGTGGTATCGACAGCTATTCTTACTGCGCCGGAGATCCTGTAAATTACTCAGACCCTACAGGGCAGGCCCCAGAGTCGc                                                           >  1:779323/1‑104 (MQ=255)
                                                                                                                   cGATGCCTTGAGCCCGTTTTTGAAAGGTGGTATCGACAGCTATTCTTACTGCGCCGGAGATCCTGTAAATTACTCAGACCCTACAGGGCAGGCCCCAGAGTCGc                                                           >  1:137783/1‑104 (MQ=255)
                                                                                                                   cGATGCCTTGAGCCCGTTTTTGAAAGGTGGTATCGACAGCTATTCTTACTGCGCCGGAGATCCTGTAAATTACTCAGACCCTACAGGGCAGGCCCCAGAGTCGc                                                           <  2:779323/104‑1 (MQ=255)
                                                                                                                          ttGAGCCCGTTTTTGAAAGGTGGTATCAACAGCTATTCTTACTGCGCCGGAGATCCTGTAAATTACTCAGACCCTACAGGGCAGGCCCCAGAGTCGCCGAAGCTAGTAAGAAAAtcatca                                    <  2:271973/120‑1 (MQ=255)
                                                                                                                          ttGAGCCCGTTTTTGAAAGGTGGTATCAACAGCTATTCTTACTGCGCCGGAGATCCTGTAAATTACTCAGACCCTACAGGGCAGGCCCCAGAGTCGCCGAAGCTAGTAAGAAAAtcatca                                    >  1:271973/1‑120 (MQ=255)
                                                                                                                          ttGAGCCCGTTTTTGAAAGGTGGTATCAACAGCTATTCTTACTGCGCCGGAGATCCTGTAAATTACTCAGACCCTACAGGGCAGGCCCCAGAGTCGCCGAAGCTAGTAAGAAAATCATCACAGACAGCCCCGTCGAAACCCTCATCTGGc      <  2:885233/150‑1 (MQ=255)
                                                                                                                              gCCCGTTTTTGAAAGGTGGTATCAACAGCTATTCTTACTGCGCCGGAGATCCTGTAAATTACTCAGACCCTACAGGGCAGGCCCCAGAGTCGCCGAAGCTAGTAAGAAAATCATCACAGACAGCCCCGTCGAAACCCTCATCTGGCTCAt  <  1:466208/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                              
GCCATTGCTCCGTGTTGAAGTACAATGGTGATCATTTAGACGATCTCGTTAACGGGTATCATTTAGGGGTAGGTAAACGATTGTATCTGCCAGCTTTGATGCGCTTCGCATCGCCCGATGCCTTGAGCCCGTTTTTGAAAGGTGGTATCAACAGCTATTCTTACTGCGCCGGAGATCCTGTAAATTACTCAGACCCTACAGGGCAGGCCCCAGAGTCGCCGAAGCTAGTAAGAAAATCATCACAGACAGCCCCGTCGAAACCCTCATCTGGCTCAT  >  NC_002947/3578592‑3578867

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