Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,818,584 G→A *275* (TAG→TAA dinD → DNA damage‑inducible protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,818,5840GA96.1% 156.1 / ‑0.8 51coding (825/825 nt)dinDDNA damage‑inducible protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (23/26);  minor base . (1/1);  total (24/27)
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

CACACACTTTGACGTGGGTCGCAAGGTGAGGCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAGCAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACAATGGCTAAATATTGCACCTTTCTTTTCCCCCTCAGTT  >  NC_000913/3818449‑3818693
                                                                                                                                       |                                                                                                             
cacacaCTTTGACGTGGGTCGCAAGGTGAGGCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAg                                                                                                            >  1:1278276/1‑139 (MQ=255)
  cacacTTTGACGTGGGTCGCAAGGTGAGGCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGaa                                                                                                          >  2:1494413/1‑139 (MQ=255)
  cacacTTTGACGTGGGTCGCAAGGTGAGGCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGaa                                                                                                          <  1:1494413/139‑1 (MQ=255)
   acacTTTGACGTGGGTCGCAAGGTGAGGCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGaaa                                                                                                         >  2:1011174/1‑139 (MQ=255)
         tGACGTGGGTCGCAAGGTGAGGCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTaaaa                                                                                                   <  2:278874/139‑1 (MQ=255)
                   cGCAAGGTGAGGCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTttg                                                                                         <  1:10081/139‑1 (MQ=255)
                             ggCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTgaga                                                                               >  2:1225625/1‑139 (MQ=255)
                              gCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTgag                                                                                >  1:1134837/1‑137 (MQ=255)
                              gCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTgag                                                                                <  2:1134837/137‑1 (MQ=255)
                               cAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGaca                                                                             <  2:1278276/139‑1 (MQ=255)
                                  aCCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTttgtt                                                                                       <  1:908434/126‑1 (MQ=255)
                                  aCCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTttgtt                                                                                       >  2:908434/1‑126 (MQ=255)
                                        cAGGTACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCg                                                                    <  1:1665178/139‑1 (MQ=255)
                                               ttGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCa                                                             <  2:279365/139‑1 (MQ=255)
                                                           tGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTgaga                                                                               <  1:1441568/109‑1 (MQ=255)
                                                           tGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTgaga                                                                               >  2:1441568/1‑109 (MQ=255)
                                                           tGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGc                                                 >  2:250756/1‑139 (MQ=255)
                                                            gCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTATAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCa                                                <  2:888359/139‑1 (MQ=255)
                                                            gCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCa                                                <  1:1225625/139‑1 (MQ=255)
                                                                  ggagTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGt                                                   >  2:1235314/1‑130 (MQ=255)
                                                                  ggagTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGt                                                   <  1:1235314/130‑1 (MQ=255)
                                                                      ttGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCa                                                             >  2:1334896/1‑116 (MQ=255)
                                                                      ttGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCa                                                             <  1:1334896/116‑1 (MQ=255)
                                                                      ttGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACaa                                      >  2:349965/1‑139 (MQ=255)
                                                                      ttGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACaa                                      >  2:919013/1‑139 (MQ=255)
                                                                      ttGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGAGCAGTGAACACaa                                      >  2:266229/1‑139 (MQ=255)
                                                                        gCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACaa                                      >  1:392075/1‑137 (MQ=255)
                                                                        gCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACaa                                      <  2:392075/137‑1 (MQ=255)
                                                                             ccccGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTgaga                                                                               <  1:1650505/91‑1 (MQ=255)
                                                                             ccccGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTgaga                                                                               >  2:1650505/1‑91 (MQ=255)
                                                                                          cATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCg                                                                         >  2:736105/1‑84 (MQ=255)
                                                                                          cATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCg                                                                         <  1:736105/84‑1 (MQ=255)
                                                                                                gCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACAATGGCTAAATATTGCACCTTTCTTTTc            <  1:266229/139‑1 (MQ=255)
                                                                                                    ttGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCa                                                             <  2:93628/86‑1 (MQ=255)
                                                                                                    ttGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCa                                                             >  1:93628/1‑86 (MQ=255)
                                                                                                          aaCAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACAATGGCTAAATATTGCACCTTTCTTTTCCCCCTCAGtt  <  1:349965/139‑1 (MQ=255)
                                                                                                              gtgtAAAAATTACAGAGAAGAAAT‑‑‑‑‑‑AACTAAAAACTCTTTTTGTTGATTGAGACACCCg                                                                         <  1:204713/58‑1 (MQ=255)
                                                                                                              gtgtAAAAATTACAGAGAAGAAAT‑‑‑‑‑‑AACTAAAAACTCTTTTTGTTGATTGAGACACCCg                                                                         >  2:204713/1‑58 (MQ=255)
                                                                                                               tgtAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATg                                                                      >  2:1039109/1‑66 (MQ=255)
                                                                                                               tgtAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATg                                                                      <  1:1039109/66‑1 (MQ=255)
                                                                                                                  aaaaaTTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCa                                                             <  2:322740/72‑1 (MQ=255)
                                                                                                                  aaaaaTTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCa                                                             >  1:322740/1‑72 (MQ=255)
                                                                                                                         aCAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAacac                                        >  2:1211046/1‑86 (MQ=255)
                                                                                                                         aCAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAacac                                        <  1:1211046/86‑1 (MQ=255)
                                                                                                                         aCAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAAcaca                                       >  2:1142496/1‑87 (MQ=255)
                                                                                                                         aCAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAAcaca                                       <  1:1142496/87‑1 (MQ=255)
                                                                                                                           agagAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGaa                                           <  2:1671413/81‑1 (MQ=255)
                                                                                                                           agagAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGaa                                           >  1:1671413/1‑81 (MQ=255)
                                                                                                                           agagAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACAATGGCTAAATATTg                         >  2:160457/1‑99 (MQ=255)
                                                                                                                           agagAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACAATGGCTAAATATTg                         <  1:160457/99‑1 (MQ=255)
                                                                                                                                    aTAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGcc                                                              <  1:38950/53‑1 (MQ=255)
                                                                                                                                       |                                                                                                             
CACACACTTTGACGTGGGTCGCAAGGTGAGGCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAGCAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACAATGGCTAAATATTGCACCTTTCTTTTCCCCCTCAGTT  >  NC_000913/3818449‑3818693

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