Predicted mutation
evidence seq id position mutation annotation gene description
RA minE 1,216,702 T→C Q677Q (CAA→CAG prc ← carboxy‑terminal protease for penicillin‑binding protein 3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,216,7020TC91.7% 69.7 / ‑0.5 24Q677Q (CAA→CAGprccarboxy‑terminal protease for penicillin‑binding protein 3
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base C (13/9);  total (13/11)
Fisher's exact test for biased strand distribution p-value = 1.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATC  >  minE/1216634‑1216764
                                                                    |                                                              
aTCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGgctg                                                              <  1:1476558/71‑1 (MQ=255)
aTCTTGTCGCTGCTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGgctg                                                              <  1:1352208/71‑1 (MQ=255)
 tCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGgctgt                                                             >  1:1864299/1‑71 (MQ=255)
     gTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGCTGTTc                                                           <  1:248973/69‑1 (MQ=255)
            ttAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGCTGTTCCGCGGGt                                                    <  1:345761/69‑1 (MQ=255)
              aaaaaaTCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGCTGTTCCGCGGGTCTg                                                 >  1:455365/1‑70 (MQ=255)
               aaaaaTCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGgctg                                                              <  1:2893112/56‑1 (MQ=255)
               aaaaaTCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGCTGTTCCGCGGGTCTGGc                                               >  1:1264021/1‑71 (MQ=255)
                  aaTCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGCTGTTCCGCGGGTCTGGCttt                                            >  1:1197505/1‑71 (MQ=255)
                          acaATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGCTGTTCCGCGGGTCTGGCTTTTTCAAGCt                                    <  1:3136076/71‑1 (MQ=255)
                            aaTTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGCTGTTCCGCGGGTCTGGCTTTTTCAAGCTTc                                  <  1:2574357/71‑1 (MQ=255)
                              tttCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGCTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGc                                <  1:1886101/71‑1 (MQ=255)
                              tttCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGCTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGc                                >  1:486027/1‑71 (MQ=255)
                                       ccTGATTGATATTACTTGACGGGAGCGGGCTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAg                       <  1:147440/71‑1 (MQ=255)
                                        cTGATTGATATTACTTGACGGGAGCGGGCTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGt                      <  1:449336/71‑1 (MQ=255)
                                          gattgatATTACTTGACGGGAGCGGGCTGTTCCGCGGGTCTGGCTTTTTCAAg                                      >  1:352378/1‑53 (MQ=255)
                                           attgatATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGc                    <  1:3062998/70‑1 (MQ=255)
                                           attgatATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTAGCTTTTTCAAGCTTCGCCAGATCGAGTGc                    <  1:2696673/70‑1 (MQ=255)
                                           attgatATTACTTGACGGGAGCGGGCTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGc                    >  1:614443/1‑70 (MQ=255)
                                           attgatATTACTTGACGGGAGCGGGCTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGc                    >  1:290489/1‑70 (MQ=255)
                                             tgatATTACTTGACGGGAGCGGGCTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAg                       >  1:1046303/1‑65 (MQ=255)
                                               atatTACTTGACGGGAGCGGGCTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGata                <  1:1995158/70‑1 (MQ=255)
                                                      ttcaCGGGAGCGGGCTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCAcc          >  1:1767789/4‑69 (MQ=255)
                                                      ttGACGGGAGCGGGCTGTTCCGCGGGTCTGGCTTTTTCAAGCt                                    <  1:2460298/43‑1 (MQ=255)
                                                        gACGGGAGCGGGCTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGtct      >  1:874443/1‑71 (MQ=255)
                                                         aCGGGAGCGGGCTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGata                >  1:2613173/1‑60 (MQ=255)
                                                           gggAGCGGGCTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCAt   >  1:55188/1‑71 (MQ=255)
                                                              aGCGGGCTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATc  >  1:1401911/1‑69 (MQ=255)
                                                                    |                                                              
ATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATC  >  minE/1216634‑1216764

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