Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,375,091 A→G 8.9% I70T (ATA→ACA)  PP_3850 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,375,0910AG 8.9% 116.7 / 2.6 45I70T (ATA→ACA) PP_3850hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (20/21);  new base G (2/2);  total (22/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.58e-01

TTTGATTCTCTAGCTTAAATCAGCTTCTACGCCATGAGTTGCGGAT‑CCCCG‑CCGCAGCTCTTGGCTGAGCGCTGAATGAGCAGAGGTAATCCTAGCAATCGATTACCCCCTGCAGCCTAGATAACGATTCTTCCTCTTCATTTTCTATCAGAGCCCACAGCCTGGCCGCCTCGGTAAGCTCGAGCATGTCGACCAGGTCATTGTCGTCGACATCATTGCGCCGACGAGCCGCATAGGCCATCACACTGAGCACCGCAGCGCGCCCATCAGGATCGGTTACTAGGGCAAC  >  NC_002947/4374945‑4375233
                                                                                                                                                    |                                                                                                                                              
tttGATTCTCTAGCTTAAATCAGCTTCTACGCCATGAGTTGCGGATGCCCCG‑CCGCAGCTATTGGATGAGCGCTGAATGAGCAGAGGTAATCATAGAAATCGATTACCCCCTGCAGCATAAATAACGATTATACCTATTCATTTtata                                                                                                                                                >  1:79858/1‑148 (MQ=255)
tttGATTCTCTAGCTTAAATCAGCTTCTACGCCATGAGTTGCGGAT‑CCCCG‑CCGCAGCTCTTGGCTGAGCGCTGAATGAGCAGAGGTAATCCTAGCAATCGATTACCCCCTGCAGCCTAGATAACGATTCTTCCTCTTCATTTTCTAt                                                                                                                                               >  6:174369/1‑148 (MQ=255)
 ttGATTCTCTAGCTTAAATCAGCTTCTACGCCATGAGTTGCGGAT‑CCCCG‑CCGCAGCTCTTGGCTGAGCGCTGAATGAGCAGAGGTAATCCTAGCAATCGATTACCCCCTGCAGCCTAGATAACGATTCTTCCTCTTCATTTTCTAt                                                                                                                                               >  8:175471/1‑147 (MQ=255)
        tctAGCTTAAATCAGCTTCTACGCCATGAGTTGCGGAT‑CCCCG‑CCGCAGCTCTTGGCTGAGCGCTGAATGAGCAGAGGTAATCCTAGCAATCGATTACCCCCTGCAGCCTAGATAACGATTCTTCCTCTTCATTTTCTATCAGAGCCca                                                                                                                                      >  3:12743/1‑149 (MQ=255)
          tAGCTTAAATCAGCTTCTACGCCATGAGTTGCGGAT‑CCCCG‑CCGCAGCTCTTGGCTGAGCGCTGAATGAGCAGAGGTAATCCTAGCAATCGATTACCCCCTGCAGCCTAGATAACGATTCTTCCTCTTCATTTTCTATCAGAGCCca                                                                                                                                      <  2:79858/147‑1 (MQ=255)
             cTTAAATCAGCTTCTACGCCATGAGTTGCGGAT‑CCCCG‑CCGCAGCTCTTGGCTGAGCGCTGAATGAGCAGAGGTAATCCTAGCAATCGATTACCCCCTGCAGCCTAGATAACGATTCTTCCTCTTCATTTTCTATCAGAGCCCACAGcc                                                                                                                                 >  4:11303/1‑149 (MQ=255)
               tAAATCAGCTTCTACGCCATGAGTTGCGGAT‑CCCCGTCCGCAGCTCTTGGCTGAGCGCTGAATGAGCAGAGGTAATCGTAGCAATCGATTACCCCCTGCAGCCTAGATAACGATTGTTCCTCTTCATTTTCTATCAGAGCCCACAGCCt                                                                                                                                >  1:238684/1‑149 (MQ=255)
                        ttCTACGCCATGAGTTGCGGAT‑CCCCG‑CCGCAGCTCTTGGCTGAGCGCTGAATGAGCAGAGGTAATCCTAGCAATCGATTACCCCCTGCAGCCTAGATAACGATTCTTCCTCTTCATTTTCTAt                                                                                                                                               <  4:174344/124‑1 (MQ=255)
                             cGCCATGAGTTGCGGAT‑CCCCG‑CCGCAGCTCTTGGCTGAGCGCTGAATGAGCAGAGGTAATCCTAGCAATCGATTACCCCCTGCAGCCTAGATAACGATTCTTCCTCTTCATTTTCTATCAGAGCCCACAGCCTGGCCGCCTCGGTaa                                                                                                                  <  5:331023/148‑1 (MQ=255)
                                 aTGAGTTGCGGAT‑CCCCG‑CCGCAGCTCTTGGCTGAGCGCTGAATGAGCAGAGGTAATCCTAGCAATCGATTACCCCCTGCAGCCTAGATAACGATTCTTCCTCTTCATTTTCTATCAGAGcc                                                                                                                                        >  2:347509/1‑122 (MQ=255)
                                 aTGAGTTGCGGAT‑CCCCG‑CCGCAGCTCTTGGCTGAGCGCTGAATGAGCAGAGGTAATCCTAGCAATCGATTACCCCCTGCAGCCTAGATAACGATTCTTCCTCTTCATTTTCTATCAGAGcc                                                                                                                                        <  1:347509/122‑1 (MQ=255)
                                      ttGCGGAT‑CCCCG‑CCGCAGCTCTTGGCTGAGCGCTGAATGAGCAGAGGTAATCCTAGCAATCGATTACCCCCTGCAGCCTAGATAACGATTCTTCCTCTTCATTTTCTATCAGAGCCCACAGCCTGGCCGCCTCGGTAAGCTCGAGCAt                                                                                                        <  3:11303/149‑1 (MQ=255)
                                          ggAT‑CCCCG‑CCGCAGCTCTTGGCTGAGCGCTGAATGAGCAGAGGTAATCCTAGCAATCGATTACCCCCTGCAGCCTAGATAACGATTCTTCCTCTTCATTTTCTATCAGAGCCCACAGCCTGGCCGCCTCGGTAAGCTCGAGCATGTc                                                                                                     >  4:218261/1‑148 (MQ=255)
                                          ggAT‑CCCCG‑CCGCAGCTCTTGGCTGAGCGCTGAATGAGAAGAGGTAATCCTAGCAATCGATTACCCCCTGCAGCCTAGATAACGATTCTTCCTCTTCATTTTCTATCAGAGCCAACAGCCTGGCCGCCTCGGTAAGCTCGAGCATGTc                                                                                                     >  3:29349/1‑148 (MQ=255)
                                            aT‑CCCCG‑CCGCAGCTCTTGGCTGAGCGCTGAATGAGCAGAGGTAATCCTAGCAATCGATTACCCCCTGCAGCCTAGATAACGATTCTTCCTCTTCATTTTCTATCAGAGCCCACAGCCTGGCCGCCTCGGTAAGCTCGAGCATGTCGAc                                                                                                  <  7:175471/149‑1 (MQ=255)
                                             t‑CCCCG‑CCGCAGCTCTTGGCTGAGCGCTGAATGAGCAGAGGTAATCCTAGCAATCGATTACCCCCTGCAGCCTAGATAACGATTCTTCCTCTTCATTTTCTATCAGAGCCCACAGCCTGGCCGCCTCGGTAAGCTCGAGCATGTCGAcc                                                                                                 <  2:238684/149‑1 (MQ=255)
                                                       gcAGCTCTTGGGTGTGCGCTGTATGAGCAGAGGTAATCCTAGCAATCGATTTCCCCCTGCAGCCTAGATAACGATTCTTCCTCTTCATTTTCTATCAGAGCCCACAGCCTGGCCGCCTCGGTAAGCTCGAGCATGTCGACCAGGTCAtt                                                                                         <  4:29349/149‑1 (MQ=255)
                                                       gcAGCTCTTGGCTGAGCGCTGAATGAGCAGAGGTAATCCTAGCAATCGATTACCCCCTGCAGCCTAGATAACGATTCTTCCTCTTCATTTTCTATCAGAGCCCACAGCCTGGCCGCCTCGGTAAGCTCGAGCATGTCGACCAGGTCAtt                                                                                         <  7:227754/149‑1 (MQ=255)
                                                              ttGGCTGAGCGCTGAATGAGCAGAGGTAATCCTAGCAATCGATTACCCCCTGCAGCCTAGATAACGATTCTTCCTCTTCATTTTCTATCAGAGCCCACAGCCTGGCCGCCTCGGTAAGc                                                                                                                >  1:141630/1‑119 (MQ=255)
                                                                ggCTGAGCGCTGAATGAGCAGAGGTAATCCTAGCAATCGATTACCCCCTGCAGCCTAGATAACGATTCTTCCTCTTCATTTTCTATCAGAGCCCACAGCCTGGCCGCCTCGGTAAGCTCGAGCATGTCGACCAGGTCATTGTCGTCGa                                                                                 <  3:218261/148‑1 (MQ=255)
                                                                 gCTGAGCGCTGAATGAGCAGAGGTAATCCTAGCAATCGATTACCCCCTGCAGCCTAGATAACGATTCTTCCTCTTCATTTTCTATCAGAGCCCACAGCCTGGCCGCCTCGGTAAGCTCGAGCATGTCGACCAGGTCATTGTCGTCGAca                                                                               >  8:203796/1‑149 (MQ=255)
                                                                        gcTGAATGAGCAGAGGTAATCCTAGCAATCGATTACCCCCTGCAGCCTAGATAACGATTCTTCCTCTTCATTTTCTATCAGAGCCCACAGCCTGGCCGCCTCGGTAAGCTCGAGCATGTCGACCAGGTCATTGTCGTCGACATCATTg                                                                         <  5:174369/148‑1 (MQ=255)
                                                                                      ggTAATCCTAGCAATCGATTACCCCCTGCAGCCTAGATAACGATTCTTCCTCTTCATTTTCTATCAGAGCCCACAGCCTGGCCGCCTCGGTAAGCTCGAGCATGTCGACCAGGTCATTGTCGTCGACATCATTGCGCCGACGAGCCGCa                                                          <  1:139315/149‑1 (MQ=255)
                                                                                                gCAATCGATTACCCCCTGCAGCCTAGATAACGATTCTTCCTGTTCATTTTCTATCAGAGCCCACAGCCTGGCCGCCTCTTTAAGCTCGAGCATGGCGACCAGGTCATTGTCGTCGACATCATTGCGCCGACGAGCCGCATAGGCCATc                                                 <  5:352534/148‑1 (MQ=255)
                                                                                                gCAATCGATTACCCCCTGCAGCCTAGATAACGATTCTTCCTCTTCATTTTCTATCAGAGCCCACAGCCTGGCCGCCTCGGTAAGCTCGAGCATGTCGACCAGGTCATTGTCGTCGACATCATTGCGCCGACGAGCCGc                                                           >  3:284366/1‑138 (MQ=255)
                                                                                                gCAATCGATTACCCCCTGCAGCCTAGATAACGATTCTTCCTCTTCATTTTCTATCAGAGCCCACAGCCTGGCCGCCTCGGTAAGCTCGAGCATGTCGACCAGGTCATTGTCGTCGACATCATTGCGCCGACGAGCCGc                                                           <  4:284366/138‑1 (MQ=255)
                                                                                                    tCGATTACCCCCTGCAGCCTAGATAACGATTCTTCCTCTTCATTTTCTATCAGAGCCCACAGCCTGGCCGCCTCGGTAAGCTCGAGCATGTCGACCAGGTCATTGTCGTCGACATCATTGCGGCGACGAGCCGCATAGGCCATCACACt                                            <  1:192108/149‑1 (MQ=255)
                                                                                                    tCGATTACCCCCTGCAGCCTAGATAACGATTCTTCCTCTTCATTTTCTATCAGAGCCCACAGCCTGGCCGCCTCGGTAAGCTCGAGCATGTCGACCAGGTCATTGTCGTCGACATCATTGCGCCGACGAGCCGCATAGGCCATCACACt                                            <  1:186316/149‑1 (MQ=255)
                                                                                                      gATTACCCCCTGCAGCCTAGATAACGATTCTTCCTCTTCATTTTCTGTCAGAGCCCACAGCCTGGCCGCCTCGGTAAGCTCGAGCATGTCGACCAGGTCATTGTCGTCGACATCATTGCGCCGACGAGCCGCATAg                                                       >  4:214105/1‑136 (MQ=255)
                                                                                                      gATTACCCCCTGCAGCCTAGATAACGATTCTTCCTCTTCATTTTCTGTCAGAGCCCACAGCCTGGCCGCCTCGGTAAGCTCGAGCATGTCGACCAGGTCATTGTCGTCGACATCATTGCGCCGACGAGCCGCATAg                                                       >  4:218476/1‑136 (MQ=255)
                                                                                                      gATTACCCCCTGCAGCCTAGATAACGATTCTTCCTCTTCATTTTCTGTCAGAGCCCACAGCCTGGCCGCCTCGGTAAGCTCGAGCATGTCGACCAGGTCATTGTCGTCGACATCATTGCGCCGACGAGCCGCATAg                                                       <  3:218476/136‑1 (MQ=255)
                                                                                                      gATTACCCCCTGCAGCCTAGATAACGATTCTTCCTCTTCATTTTCTGTCAGAGCCCACAGCCTGGCCGCCTCGGTAAGCTCGAGCATGTCGACCAGGTCATTGTCGTCGACATCATTGCGCCGACGAGCCGCATAg                                                       <  3:214105/136‑1 (MQ=255)
                                                                                                      gATTACCCCCTGCAGCCTAGATAACGATTCTTCCTCTTCATTTTCTATCAGAGCCCACAGCCTGGCCGCCTCGGTAAGCTCGAGCATGTCGACCAGGTCATTGTCGTCGACATCATTGCGCCGACGAGCCGCATAGGCCATCACACTg                                           <  4:12743/148‑1 (MQ=255)
                                                                                                           cccccTGCAGCCTAGATAACGATTCTTCCTCTTCATTTTCTATCAGAGCCCACAGCCTGGCCGCCTCGGTAAGCTCGAGCATGTCGACCAGGTCATTGTCGTCGACATCATTGCGCCGACGAGCCGCATAGGCCATCACACTGAGCAcc                                     >  6:285918/1‑149 (MQ=255)
                                                                                                             cccTGCAGCCTAGATAACGATTCTTCCTCTTCATTTTCTATCAGAGCCCACAGCCTGGCCGCCTCGGTAAGCTCGAGCATGTCGACCAGGTCATTGTCGTCGACATCATTGCGCCGACGAGCCGCATAGGCCATCACACTGAGCACCg                                    >  6:185950/1‑148 (MQ=255)
                                                                                                              ccTGCAGCCTAGATAACGATTCTTCCTCTTCATTTTCTATCAGAGCCCACAGCCTGGCCGCCTCGGTAAGCTCGAGCATGTCGACCAGGTCATTgtc                                                                                      >  8:217546/1‑97 (MQ=255)
                                                                                                              ccTGCAGCCTAGATAACGATTCTTCCTCTTCATTTTCTATCAGAGCCCACAGCCTGGCCGCCTCGGTAAGCTCGAGCATGTCGACCAGGTCATTgtc                                                                                      <  7:217546/97‑1 (MQ=255)
                                                                                                                      cTAGATAACGATTCTTCCTCTTCATTTTCTATCAGAGCCCACAGCCTGGCCGCCTCGGTAAGCTAGAGCATGTCGACCAGGTCATTGTCGTCGACATCATTGCGCCGACGAGCCGCATAGGCCATCACACTGAGCACCGCAGCGCGcc                           >  5:191787/1‑148 (MQ=255)
                                                                                                                             aCGATTCTTCCTCTTCATTTTCTATCAGAGCCCACAGCCTGGCCGCCTCGGTAAGCTCGAGCATGTCGACCAGGTCATTGTCGTCGACATCATTGCGCCGACGAGCCGCATAGGCCATCACACTGAGCACCGCAGCGCGCCCATCAGGa                   <  7:203796/149‑1 (MQ=255)
                                                                                                                               gATTCTTCCTCTTCATTTTCTATCAGAGCCCACAGCCTGGCCGCCTCGGTAAGCTCGAGCATGTCGACCAGGTCATTGTCGTCGACATCATTGCGCCGACGAGCCGCATAGGCCATCACACTGAGCACCGCAGCGCGCCCATCAGGATc                 >  2:174541/1‑149 (MQ=255)
                                                                                                                                        tctTCATTTTCTATCAGAGCCCACAGCCTGGCCGCCTCGGTAAGCTCGAGCATGTCGACCAGGTCATTGTCGTCGACATCATTGCGCCGACGAGCCGCATAGGCCATCACACTGAGCACCGCAGCGCGCCCATCAGGATCGGTTACTAg        >  6:264548/1‑149 (MQ=255)
                                                                                                                                             aTTTTCTATCAGAGCCCACAGCCTGGCCGCCTCGGTAAGCTCGAGCATGTCGACCAGGTCATTGTCGTCGACATCATTGCGCCGACGAGCCGCATAGGCCATCACACTGa                                          <  3:183462/110‑1 (MQ=255)
                                                                                                                                             aTTTTCTATCAGAGCCCACAGCCTGGCCGCCTCGGTAAGCTCGAGCATGTCGACCAGGTCATTGTCGTCGACATCATTGCGCCGACGAGCCGCATAGGCCATCACACTGa                                          >  4:183462/1‑110 (MQ=255)
                                                                                                                                              ttttCTATCAGAGCCCACAGCCTGGCCGCCTCGGTAAGCTCGAGCATGTCGACGAGGTCATTGTCGTCGACATCATTGCGCCGACGAGCCGCATAGGCCATCACACTGAGCACCGCAGCGCGCCCATCAGGAGCGGTTACTAGGGCAAc  >  7:104050/1‑149 (MQ=255)
                                                                                                                                              ttttCTATCAGAGCCCACAGCCTGGCCGCCTCGGTAAGCTCGAGCATGTCGACCAGGTCATTGTCGTCGACATCATTGCGCCGACGAGCCGCATAGGCCATCACACTGAGCACCGCAGCGCGCCCATCAGGATCGGTTACTAgg       >  2:132498/1‑144 (MQ=255)
                                                                                                                                              ttttCTATCAGAGCCCACAGCCTGGCCGCCTCGGTAAGCTCGAGCATGTCGACCAGGTCATTGTCGTCGACATCATTGCGCCGACGAGCCGCATAGGCCATCACACTGAGCACCGCAGCGCGCCCATCAGGATCGGTTACTAgg       <  1:132498/144‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
TTTGATTCTCTAGCTTAAATCAGCTTCTACGCCATGAGTTGCGGAT‑CCCCG‑CCGCAGCTCTTGGCTGAGCGCTGAATGAGCAGAGGTAATCCTAGCAATCGATTACCCCCTGCAGCCTAGATAACGATTCTTCCTCTTCATTTTCTATCAGAGCCCACAGCCTGGCCGCCTCGGTAAGCTCGAGCATGTCGACCAGGTCATTGTCGTCGACATCATTGCGCCGACGAGCCGCATAGGCCATCACACTGAGCACCGCAGCGCGCCCATCAGGATCGGTTACTAGGGCAAC  >  NC_002947/4374945‑4375233

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