Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,001,421 2 bp→GC coding (381‑382/903 nt) rarD ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,001,4210CG98.8% 218.0 / NA 84G128R (GGA→CGA) rarDhypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base G (43/40);  minor base T (0/1);  total (43/41)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009134,001,4220GC100.0% 223.9 / NA 84D127E (GAC→GAGrarDhypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (43/41);  total (43/41)

AAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCCGTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCCGT  >  NC_000913/4001352‑4001488
                                                                     ||                                                                  
aaGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGcc                                                               >  1:4000656/1‑76 (MQ=255)
 agcccACAGCTGGACTAACCCCCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCg                                                              <  1:235918/73‑1 (MQ=255)
 aGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCg                                                              >  1:264810/1‑76 (MQ=255)
   tCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAg                                                            >  2:2218786/1‑76 (MQ=255)
   tCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAg                                                            >  2:1510768/1‑76 (MQ=255)
   cccACGGCGGGACTAACCCGCCCCATATCGCCAGAACCCCCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAg                                                            <  2:3204900/75‑1 (MQ=255)
     cacaGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGa                                                          >  2:2740782/1‑76 (MQ=255)
     cacaGCTGGACTAACACGCCACATATCGCCAGAATAACCGCCAGCCATTGCATCCGGCGCAATCGCTCGCCGacga                                                          >  1:2605635/1‑73 (MQ=255)
      acaGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAggca                                                         >  2:2174346/1‑74 (MQ=255)
      acaGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGaa                                                         >  1:2001807/1‑76 (MQ=255)
       caccTGAACAAACACCCCACATATCGCCAGAATCCCCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGaaa                                                        <  1:4789629/73‑1 (MQ=255)
        aGCTGGACTAGCCCGCCACATATCGCGAGAATCCCCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGaaaa                                                       <  1:3870085/76‑1 (MQ=255)
        aGCTGGACTAACACGCCACATAGCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGaaaa                                                       <  2:1442809/76‑1 (MQ=255)
         gCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAat                                                      >  2:4559916/1‑76 (MQ=255)
            ggACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGCGGAAAatcac                                                   >  2:65033/1‑75 (MQ=255)
            ggACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGCGGAAAatcac                                                   >  2:411702/1‑75 (MQ=255)
             gACTAACACCCCACATATCGCCAGAATCCCCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAatcatc                                                  <  2:4111181/76‑1 (MQ=255)
              acaaaCACCCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATcc                                                 <  1:3986662/73‑1 (MQ=255)
               cTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAAACATccc                                                >  2:3738697/1‑76 (MQ=255)
                tAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATcccc                                               >  2:2863194/1‑76 (MQ=255)
                 aaCACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCa                                              >  2:1588051/1‑76 (MQ=255)
                  acacGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCac                                             >  2:1557091/1‑75 (MQ=255)
                  acacGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAg                                             >  1:2671861/1‑76 (MQ=255)
                  acacGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAg                                             >  2:1548266/1‑76 (MQ=255)
                  acacGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAg                                             >  2:1322151/1‑76 (MQ=255)
                   cacGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGc                                            >  2:248541/1‑76 (MQ=255)
                     cGCCCCATACCGCCAGAACCCCCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGcac                                          <  2:36630/76‑1 (MQ=255)
                     cGCCACATATCCCCAGAACCCCCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGcac                                          <  2:3267811/76‑1 (MQ=255)
                      gCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAGATCATCCCCAGcaca                                         >  1:3060966/1‑76 (MQ=255)
                      gCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGcaca                                         >  1:1563539/1‑76 (MQ=255)
                      gCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGcaca                                         >  1:1986874/1‑76 (MQ=255)
                       ccACATATCGCCAGAATCCCCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACaa                                        <  1:3718552/76‑1 (MQ=255)
                         ccaTACCGCCAGAACCCCCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATg                                      <  1:2508509/75‑1 (MQ=255)
                          caTATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGt                                     >  1:3823837/1‑76 (MQ=255)
                           agatCGCCAGAACCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGtt                                    <  2:1402767/74‑1 (MQ=255)
                            tgccGCCAGAATCCCCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTc                                   <  2:616181/73‑1 (MQ=255)
                            tatCGCCAGAATCCCCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTc                                   <  2:3866041/76‑1 (MQ=255)
                            tatCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTc                                   >  2:3586835/1‑76 (MQ=255)
                            tatCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTc                                   >  1:3668195/1‑76 (MQ=255)
                             atCGCAAGAACCCCCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCa                                  <  1:4161808/76‑1 (MQ=255)
                             accGCCAGAATCCCCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCa                                  <  1:30771/74‑1 (MQ=255)
                              tCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCAc                                 >  1:572798/1‑76 (MQ=255)
                               cGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCAcc                                >  2:4417705/1‑76 (MQ=255)
                                gCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCa                               >  2:1301649/1‑76 (MQ=255)
                                gCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCa                               >  1:3129208/1‑76 (MQ=255)
                                cccAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCa                               <  1:3766119/75‑1 (MQ=255)
                                 ccAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAg                              <  2:3820757/76‑1 (MQ=255)
                                  cAGAATCAGCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGc                             <  2:4646451/76‑1 (MQ=255)
                                   aGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCg                            >  2:3848825/1‑76 (MQ=255)
                                   aGAACCACCCCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCg                            <  2:120969/76‑1 (MQ=255)
                                      accACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGt                         <  1:1618685/74‑1 (MQ=255)
                                        cgccGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTta                       <  1:4216574/74‑1 (MQ=255)
                                        cACCGCCAGCCATCGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTta                       >  2:1691765/1‑76 (MQ=255)
                                         tccGCGATCCATTGCATCCGGCGGACTGTCTCGCCGAGGAAAATCATGCGCAGCCGAATGTTCACGCGCGTGTtaa                      <  2:2098004/75‑1 (MQ=255)
                                         cccGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTtaa                      <  1:713260/75‑1 (MQ=255)
                                         cccGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTtaa                      <  1:25232/75‑1 (MQ=255)
                                          ccGCCAGCCATTGCATCCGGCGGAATCGCCCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTcaat                     >  2:858307/1‑76 (MQ=255)
                                           cGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTtaata                    <  1:1391535/76‑1 (MQ=255)
                                           cGCCAGCCATTCCATCCGGCGGAATCGCTAGCAGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGCtaata                    >  2:460176/1‑76 (MQ=255)
                                            gccagccaTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTtaataa                   >  1:2231124/1‑76 (MQ=255)
                                             ccagccaTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATaaa                  <  2:595459/76‑1 (MQ=255)
                                             ccagccaTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATaaa                  >  2:4357789/1‑76 (MQ=255)
                                              cagccaTTGCATCCGGCGGAATCGCTCGCCGAGGAAGATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAg                 <  2:3060966/76‑1 (MQ=255)
                                              cagccaTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAg                 <  1:1548267/76‑1 (MQ=255)
                                              cagccaTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAg                 <  1:2901594/76‑1 (MQ=255)
                                               agccaTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTACTAAAga                >  2:4751650/1‑75 (MQ=255)
                                               agccaTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGt                <  1:4794418/76‑1 (MQ=255)
                                                 ccaTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTaa              <  2:3132842/76‑1 (MQ=255)
                                                  caTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAACGTAAc             >  2:1644569/1‑76 (MQ=255)
                                                   aTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAAcc            >  2:4536660/1‑76 (MQ=255)
                                                   aTTGCATCCGGCGGAATCGCTAGCCGAGGAAAACCATCCCCAGCACAATGTTCAACAGCGCGCTAATAAAGTaaac            >  1:2924400/1‑74 (MQ=255)
                                                     tGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCaa          >  1:2498578/1‑76 (MQ=255)
                                                      gCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAg         >  2:3922777/1‑76 (MQ=255)
                                                       cATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCAACAGCGGGTTAATAAAGTAACCAAgg        <  1:3094827/76‑1 (MQ=255)
                                                         tCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCt      <  2:1180805/76‑1 (MQ=255)
                                                          ccGGCGGAATCGCTCGCCGATGAAAATCATGCCCAGCACAATGTTCACCGGCGGGTTAATATAGTAACCAAGGCTc     <  1:460176/76‑1 (MQ=255)
                                                          ccGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACGAAGGCTc     <  2:1930327/76‑1 (MQ=255)
                                                          ccGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTc     <  1:4048576/76‑1 (MQ=255)
                                                          ccGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTc     <  2:1447906/76‑1 (MQ=255)
                                                          ccGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTc     <  2:3129208/76‑1 (MQ=255)
                                                             gcggAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTcgct  <  1:4300288/76‑4 (MQ=255)
                                                             gcggAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTcgct  >  1:2423106/1‑73 (MQ=255)
                                                             gcggAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTcgct  <  1:4447477/76‑4 (MQ=255)
                                                             gcggAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGCTCACCAGCGGGTTAATAAAGTAACCAACGCTcgct  >  1:1893848/1‑73 (MQ=255)
                                                                     ||                                                                  
AAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCCGTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCCGT  >  NC_000913/4001352‑4001488

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