Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,887,352 T→A intergenic (‑133/+226) ygcB ← / ← cysH Cascade complex anti‑viral R‑loop helicase‑annealase Cas3/phosphoadenosine phosphosulfate reductase; PAPS reductase, thioredoxin dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,887,3520TA95.7% 228.7 / 0.4 70intergenic (‑133/+226)ygcB/cysHCascade complex anti‑viral R‑loop helicase‑annealase Cas3/phosphoadenosine phosphosulfate reductase; PAPS reductase, thioredoxin dependent
Reads supporting (aligned to +/- strand):  ref base T (3/0);  new base A (35/32);  total (38/32)
Fisher's exact test for biased strand distribution p-value = 2.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.02e-01

CAGCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTC  >  NC_000913/2887258‑2887447
                                                                                              |                                                                                               
cAGCGGGCTATTTTTTTTTTTTTTTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTg                                                                                            >  2:1410753/1‑100 (MQ=255)
cAGCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTg                                                                                            <  1:590963/100‑1 (MQ=255)
cAGCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTg                                                                                            >  1:460856/1‑100 (MQ=255)
cAGCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTg                                                                                            <  2:818252/100‑1 (MQ=255)
cAGCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTg                                                                                            <  2:346776/100‑1 (MQ=255)
cAGCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTg                                                                                            >  1:213572/1‑100 (MQ=255)
cAGCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTg                                                                                            <  1:1550592/100‑1 (MQ=255)
 aGCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTgg                                                                                           <  1:827155/100‑1 (MQ=255)
   cGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGt                                                                                         >  1:1035279/1‑100 (MQ=255)
      gCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCg                                                                                      <  2:594302/100‑1 (MQ=255)
             ttCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACTTTCGAATAGTTTGGGTTCGAACGCTg                                                                               <  1:688205/100‑1 (MQ=255)
             ttCGATGCTGCTCAGAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTg                                                                               <  2:1075763/100‑1 (MQ=255)
                gATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGcc                                                                            >  1:1487581/1‑100 (MQ=255)
                gATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGcc                                                                            >  2:244513/1‑100 (MQ=255)
                gATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGcc                                                                            >  1:912460/1‑100 (MQ=255)
                gATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGcc                                                                            <  1:340481/100‑1 (MQ=255)
                gATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGc                                                                             >  1:1276655/1‑99 (MQ=255)
                    ctgctCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAg                                                                        <  1:1562047/100‑1 (MQ=255)
                       ctCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGtt                                                                     >  2:1115122/1‑100 (MQ=255)
                        tCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTg                                                                    >  2:1129234/1‑100 (MQ=255)
                         cATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGa                                                                   >  2:1255128/1‑100 (MQ=255)
                          aTAACTATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGa                                                                   >  2:1145578/1‑99 (MQ=255)
                          aTAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGAt                                                                  <  2:237751/100‑1 (MQ=255)
                            aaaTATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGAtcg                                                                >  2:962007/1‑98 (MQ=255)
                                atCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAat                                                            >  2:807350/1‑100 (MQ=255)
                                   gAATTACAGTGATATACACATGTCTTGCTTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTTGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATc                                                         >  2:232485/1‑100 (MQ=255)
                                       tACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCt                                                     >  2:870961/1‑100 (MQ=255)
                                        aCAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTg                                                    <  1:1171477/100‑1 (MQ=255)
                                         cAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTgg                                                   <  1:226454/100‑1 (MQ=255)
                                          aGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTggg                                                  <  2:1487572/100‑1 (MQ=255)
                                          aGTGATATAAACATGTCTTGCCTCACTGCAAAGAAGATACTAACATGCGAATAGTATAGGTTCGAGCGCTGGCCTCAGGTTGATAGAAATATCGCCTggg                                                  <  2:110638/100‑1 (MQ=255)
                                           gTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTgggg                                                 >  1:1373324/1‑100 (MQ=255)
                                            tGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGc                                                <  1:69742/100‑1 (MQ=255)
                                            tGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGc                                                >  2:511192/1‑100 (MQ=255)
                                             gATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCt                                               >  2:1433327/1‑100 (MQ=255)
                                               tataCACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCttt                                             <  1:1427624/100‑1 (MQ=255)
                                               tataCACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCttt                                             <  1:1958/100‑1 (MQ=255)
                                                 taCACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTAGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTg                                           >  1:1108846/1‑100 (MQ=255)
                                                 taCACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACACTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTg                                           <  2:851903/100‑1 (MQ=255)
                                                  acacATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGt                                          <  1:34212/100‑1 (MQ=255)
                                                    acaTGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTcc                                        >  1:990009/1‑100 (MQ=255)
                                                     caTGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCa                                       <  1:1452654/100‑1 (MQ=255)
                                                       cttCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATc                                     <  2:757489/98‑1 (MQ=255)
                                                        gTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCt                                    <  1:694491/100‑1 (MQ=255)
                                                        gTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCt                                    <  2:1373317/100‑1 (MQ=255)
                                                         tCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATcag                                   >  1:757491/1‑98 (MQ=255)
                                                         tCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTg                                   <  1:1551218/100‑1 (MQ=255)
                                                          cTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTgg                                  >  1:753271/1‑100 (MQ=255)
                                                              ncTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTACGAACGCTGGCCTCAGGTTAATAGCAATATCGACTGGGGCTTTTGTCAATCTGGa                                 >  2:862355/2‑97 (MQ=255)
                                                               cTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAAcc                              <  2:1043157/99‑1 (MQ=255)
                                                                   cTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTcgcg                         >  2:164076/1‑100 (MQ=255)
                                                                   cTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTcgcg                         >  1:766257/1‑100 (MQ=255)
                                                                   cTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTcgcg                         >  1:630210/1‑100 (MQ=255)
                                                                   cTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTcgc                          >  2:800425/1‑99 (MQ=255)
                                                                   cTGCAAAGAAGAGACTAACATACGAATAAGTTGGGTTCGAACGCTAGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTcacg                         >  2:151740/1‑100 (MQ=255)
                                                                      cAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCATCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACATCGCGAAt                      >  2:639221/1‑100 (MQ=255)
                                                                        aagaagAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGc                    <  1:557374/100‑1 (MQ=255)
                                                                         agaagaGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCt                   <  1:1203383/100‑1 (MQ=255)
                                                                          gaagaGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGTAACCTCGCGAATGc                    >  2:880175/1‑98 (MQ=255)
                                                                            agagACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTaa                >  2:273633/1‑100 (MQ=255)
                                                                             gagaCTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAAc               <  1:244513/100‑1 (MQ=255)
                                                                               gaCTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGc             >  1:1279627/1‑100 (MQ=255)
                                                                               gaCTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGc             >  1:440640/1‑100 (MQ=255)
                                                                                aCTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGcc            <  2:1222742/100‑1 (MQ=255)
                                                                                 cTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCa           >  1:585583/1‑100 (MQ=255)
                                                                                   aaCATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGa         >  2:115976/1‑100 (MQ=255)
                                                                                     cATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACa       <  1:1410761/100‑1 (MQ=255)
                                                                                     cATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACa       >  2:985705/1‑100 (MQ=255)
                                                                                      aTTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGGATGCTTAACGCCAGAc        >  2:1338117/1‑98 (MQ=255)
                                                                                          gAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTc  <  2:1561405/100‑1 (MQ=255)
                                                                                              |                                                                                               
CAGCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTC  >  NC_000913/2887258‑2887447

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