Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,083,832 (A)7→8 intergenic (‑35/+103) speB ← / ← speA agmatinase/arginine decarboxylase, biosynthetic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,083,8251.A93.8% 107.1 / 0.4 32intergenic (‑28/+110)speB/speAagmatinase/arginine decarboxylase, biosynthetic
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base A (19/11);  total (20/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

CCAAAGGCATTGGAAACCAGTGAGTTATCGTATTGATGACCTAAGGTGCTCATGGATGGACTCCTATTACAAAGTCGATATAAAAAAACCCTTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCA  >  NC_000913/3083745‑3083900
                                                                                 |                                                                           
ccgaaGGCATTGGAAACCAGTAAGTTATCGTATTGATGACCTAAGGTGCTCATGGATGGACTCCTATTACAAAGTCGatat                                                                              <  1:567141/78‑1 (MQ=255)
ccAAAGGCATTGGAAACCAGTGAGTTATCGTATTGATGACCTAAGGTGCTCATGGATGGACTCCTATTACAAAGTCGatat                                                                              <  1:567140/81‑1 (MQ=255)
   aaGGCATTGGAAACCAGTGAGTTATCGTATTGATGACCTAAGGTGCTCATGGATGGACTCCTATTACAAAGTCGATATaaa                                                                          <  2:567143/81‑1 (MQ=255)
   aaGGCATTGGAAACCAGTGAGTTATCGTATTGATGACCTAAGGTGCTCATGGATGGACTCCTATTACAAAGTCGATATaaa                                                                          <  2:567142/81‑1 (MQ=255)
       cATTGGAAACCAGTGAGTTATCGTATTGATGACCTAAGGTGCTCATGGATGGACTCCTATTACAAAGTCGATATaaaaaa                                                                       >  1:567201/1‑80 (MQ=255)
        aTTGGAAACCAGTGAGTTATCGTATTGATGACCTAAGGTGCTCATGGATGGACTCCTATTACAAAGTCGATATaaaaaaaa                                                                     >  1:567178/1‑80 (MQ=255)
           ggAAACCAGTGAGTTATCGTATTGATGACCTAAGGTGCTCATGGATGGACTCCTATTACAAAGTCGATATAAAAAAAaccc                                                                  >  2:567180/1‑80 (MQ=255)
            gAAACCAGTGAGTTATCGTATTGATGACCTAAGGTGCTCATGGATGGACTCCTATTACAAAGTCGATATAAAAAAAaccct                                                                 <  1:567144/81‑3 (MQ=255)
             aaaCCAGTGAGTTATCGTATTGATGACCTAAGGTGCTCATGGATGGACTCCTATTACAAAGTCGATATAAAAAAAACCCtt                                                                 >  2:567179/1‑81 (MQ=255)
              aaCCAGTGAGTTATCGTATTGATGACCTAAGGTGCTCATGGATGGACTCCTATTACAAAGTCGATATAAAAAAAACCCTTc                                                                <  1:567145/81‑1 (MQ=255)
              aaCCAGTGAGTTATCGTATTGATGACCTAAGGTGCTCATGGATGGACTCCTATTACAAAGTCGATATAAAAAAAACCCTTc                                                                >  2:567198/1‑81 (MQ=255)
              aaCCAGTGAGTTATCGTATTGATGACCTAAGGTGCTCATGGATGGACTCCTATTACAAAGTCGATATAAAAAAAACCCTTc                                                                >  2:567191/1‑81 (MQ=255)
              aaCCAGTGAGTTATCGTATTGATGACCTAAGGTGCTCATGGATGGACTCCTATTACAAAGTCGATATAAAAAAAACCCTTc                                                                >  2:567181/1‑81 (MQ=255)
                ccAGTGAGTTATCGTATTGATGACCTAAGGTGCTCATGGATGGACTCCTATTACAAAGTCGATATAAAAAAAACCCTTCcg                                                              >  1:567205/1‑81 (MQ=255)
                 cAGTGAGTTATCGTATTGATGACCTAAGGTGCTCATGGATGGACTCCTATTACAAAGTCGATATAAAAAAAACCCTTCcgc                                                             >  2:567147/1‑81 (MQ=255)
                 cAGTGAGTTATCGTATTGATGACCTAAGGTGCTCATGGATGGACTCCTATTACAAAGTCGATATAAAAAAAACCCTTCcgc                                                             <  1:567146/81‑1 (MQ=255)
                  aGTGAGTTATCGTATTGATGACCTAAGGTGCTCATGGATGGACTCCTATTACAAAGTCGATATAAAAAAAACCCTTCcgcg                                                            >  2:567184/1‑81 (MQ=255)
                      aGTTATCGTATTGATGACCTAAGGTGCTCATGGATGGACTCCTATTACAAAGTCGATATAAAAAAAACCCTTCCGcgtcgt                                                        <  2:567148/81‑1 (MQ=255)
                       gTTATCGTATTGATGACCTAAGGTGCTCATGGATGGACTCCTATTACAAAGTCGATATAAAAAAAACCCTTCCGCGTCGtt                                                       <  1:567149/81‑1 (MQ=255)
                          aTCGTATTGATGACCTAAGGTGCTCATGGATGGACTCCTATTACAAAGTCGATATAAAAAAAACCCTTCCGCGTCGTTAgg                                                    >  2:567195/1‑81 (MQ=255)
                          aTCGTATTGATGACCTAAGGTGCTCATGGATGGACTCCTATTACAAAGTCGATATAAAAAAAACCCTTCCGCGTCGTTAgg                                                    >  2:567182/1‑81 (MQ=255)
                                ttGATGACCTAAGGTGCTCATGGATGGACTCCTATTACAAAGTCGATATAAAAAAAACCCTTCCGCGTCGTTAGGCCcgac                                              <  2:567150/81‑1 (MQ=255)
                                ttGATGACCTAAGGTGCTCATGGATGGACTCCTATTACAAAGTCGATATAAAAAAAACCCTTCCGCGTCGTTAGGCCcga                                               <  2:567151/80‑1 (MQ=255)
                                      aCCTAAGGTGCTCATGGATGGACTCCTATTACAAAGTCGATATAAAAAAAACCCTTCCGCGTCGTTAGGCCCGACGAgg                                          >  1:567185/1‑79 (MQ=255)
                                          aaGGTGCTCATGGATGGACTCCTATTACAAAGTCGATATAAAAAAAACCCTTCCGCGTCGTTAGGCCCGACGAGGAAGGGt                                    >  2:567199/1‑81 (MQ=255)
                                           aGGTGCTCATGGATGGGCTCCTATTAAAAAGTCGATGTAAAAAAAACCCTTCCGCGTCGTTAGGACCGACGAGGAAGGGt                                    <  2:567152/80‑1 (MQ=255)
                                               gCTCATGGATGGACTCCTATTACAAAGTCGATATAAAAAAAACCCTTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGAt                               >  2:567197/1‑81 (MQ=255)
                                                cTCATGGATGGACTCCTATTACAAAGTCGATATAAAAAAAACCCTTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGAtt                              >  1:567210/1‑81 (MQ=255)
                                                       atggaCTCCTATTACAAAGTCGATATAAAAAAAACCCTTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTcaca                       <  2:567153/81‑1 (MQ=255)
                                                       atggaCTCCTATTACAAAGTCGATATAAAAAAAACCCTTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTcac                        >  2:567173/1‑80 (MQ=255)
                                                               cTATTACAAAGTCGATATAAAAAAAACCCTTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAAt                 <  2:567154/79‑1 (MQ=255)
                                                               cTATTACAAAGTCGATATAAAAAAAACCCTTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATtg               >  1:567207/1‑81 (MQ=255)
                                                                tATTACAAAGTCGATATAAAAAAAACCCTTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATtgt              >  1:567189/1‑81 (MQ=255)
                                                                 aTTACAAAGTCGATATAAAAAAAACCCTTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATtgtg             >  1:567202/1‑81 (MQ=255)
                                                                    aCAAAGTCGATATAAAAAAAACCCTTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCg          >  2:567196/1‑81 (MQ=255)
                                                                     cAAAGTCGATATAAAAAAAACCCTTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCg          >  1:567193/1‑80 (MQ=255)
                                                                        aGTCGATATAAAAAAAACCCTTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGAtt      <  1:567155/81‑1 (MQ=255)
                                                                            gATATAAAAAAAACCCTTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCa  <  1:567156/81‑1 (MQ=255)
                                                                                 |                                                                           
CCAAAGGCATTGGAAACCAGTGAGTTATCGTATTGATGACCTAAGGTGCTCATGGATGGACTCCTATTACAAAGTCGATATAAAAAAACCCTTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCA  >  NC_000913/3083745‑3083900

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