Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,884,593 T→C D876G (GAT→GGT)  ygcB ← Cascade complex anti‑viral R‑loop helicase‑annealase Cas3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,884,5930TC97.7% 157.8 / NA 44D876G (GAT→GGT) ygcBCascade complex anti‑viral R‑loop helicase‑annealase Cas3
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (25/18);  total (26/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 = 1.00e+00

ATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTA  >  NC_000913/2884447‑2884736
                                                                                                                                                  |                                                                                                                                               
aTACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCACcccc                                                                                                                                             >  2:334843/1‑151 (MQ=255)
                  aagaTGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCACCCCCTGTATAGGTAATAACAAt                                                                                                                           >  1:104499/1‑151 (MQ=255)
                    gaTGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCACCCCCTGTATAGGTAATAACAATAc                                                                                                                         <  1:182321/151‑1 (MQ=255)
                        tACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCACCCCCTGTATAGGTAATAACAATACTGtt                                                                                                                     >  2:404626/1‑151 (MQ=255)
                                gCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCACCCCCTGTATAGGTAATAACAATACTGTTACCCTGc                                                                                                              <  2:31545/150‑1 (MQ=255)
                                  aCCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCACCCCCTGTATAGGTAATAACAATACTGTTACCCTGcc                                                                                                             >  2:432385/1‑149 (MQ=255)
                                  aCCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCACCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAg                                                                                                           <  1:595995/151‑1 (MQ=255)
                                       cccTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCACCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGAccc                                                                                                       <  2:408902/150‑1 (MQ=255)
                                            cTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCACCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACccatccat                                                                                                 >  2:450735/1‑151 (MQ=255)
                                               aaGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCACCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCatccatcc                                                                                               <  1:668528/150‑1 (MQ=255)
                                               aaGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCACCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCatccatcc                                                                                               <  1:610634/150‑1 (MQ=255)
                                                aGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCACCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAg                                                                                             <  1:347397/151‑1 (MQ=255)
                                                  tAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCACCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCatccatcc                                                                                               >  1:569007/1‑147 (MQ=255)
                                                                tATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCACCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTc                                                                              <  2:368606/150‑1 (MQ=255)
                                                                 aTTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCACCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCa                                                                            <  1:637097/151‑1 (MQ=255)
                                                                       cTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCACCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCaa                                                                       >  1:737470/1‑150 (MQ=255)
                                                                        ttattaTCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCACCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCaa                                                                       >  2:719342/1‑149 (MQ=255)
                                                                          attatCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCACCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAATGt                                                                    >  2:496438/1‑150 (MQ=255)
                                                                          attatCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCACCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTa                                                                   <  2:724726/151‑1 (MQ=255)
                                                                          attatCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCACCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTa                                                                   <  1:239329/151‑1 (MQ=255)
                                                                             atCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCACCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAAcc                                                                >  2:643928/1‑151 (MQ=255)
                                                                               cGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCACCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAAccc                                                               >  2:14098/1‑150 (MQ=255)
                                                                                  cATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCACCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCt                                                           >  1:105737/1‑151 (MQ=255)
                                                                                      gATCACAATCCTTCCCGAAGTTATTTGGGATTTCAAGGGCTACCTCTGGTAATCCCTTCATCCCCTGTATAGGTAATAACGATAGGGTTACCCGGCCAGACCCCTCGATCCAGATTCTGTTTACCTTCCAGCCAAGGTAACCAATCTTCAt                                                       >  1:124364/1‑151 (MQ=255)
                                                                                       aTAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCACCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATc                                                      >  2:604246/1‑151 (MQ=255)
                                                                                          aCAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCACCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATctac                                                   >  2:639411/1‑151 (MQ=255)
                                                                                          aCAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCACCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATctac                                                   <  2:598997/151‑1 (MQ=255)
                                                                                           cAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCACCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATctact                                                  >  2:297059/1‑151 (MQ=255)
                                                                                                aTTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCACCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTc                                             <  1:572926/151‑1 (MQ=255)
                                                                                                  tCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCACCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTc                                             >  1:23919/1‑149 (MQ=255)
                                                                                                     cGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCACCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGaaa                                        <  1:34439/151‑1 (MQ=255)
                                                                                                           tATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCACCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTAc                                  >  2:531717/1‑151 (MQ=255)
                                                                                                              ttGGGATTTGCAGGGATGACTCTGGTCATCCCTTCACCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCGTCATCTACTACTTCAGAAAAACTACg                                 <  1:770477/149‑1 (MQ=255)
                                                                                                               tGGGATTTGCAGGGATGACTCTGGTCATCCCTTCACCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGtt                               <  2:72275/150‑1 (MQ=255)
                                                                                                                    tttGCAGGGATGACTCTGGTCATCCCTTCACCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAgg                         >  2:432164/1‑151 (MQ=255)
                                                                                                                      tGCAGGGATGACTCTGGTCATCCCTTCACCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTg                       >  1:726979/1‑151 (MQ=255)
                                                                                                                       gCAGGGATGACTCTGGTCATCCCTTCACCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTg                       >  2:657633/1‑150 (MQ=255)
                                                                                                                         aGGGATGACTCTGGTCATCCCTTCACCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGa                      >  1:41759/1‑149 (MQ=255)
                                                                                                                              tGACTCTGGTCATCCCTTCACCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTAc               >  1:289877/1‑151 (MQ=255)
                                                                                                                                aCTCTGGTCATCCCTTCACCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTc                                             <  1:343410/119‑1 (MQ=255)
                                                                                                                                aCTCTGGTCATCCCTTCACCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTc                                             >  2:343444/1‑119 (MQ=255)
                                                                                                                                      gTCATCCCTTCACCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGAcgc       <  2:527454/151‑1 (MQ=255)
                                                                                                                                         aTCCCTTCACCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGAt    >  2:487433/1‑151 (MQ=255)
                                                                                                                                             cTTCACCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTa  >  2:426553/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
ATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTA  >  NC_000913/2884447‑2884736

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑