Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 1,112,202 T→C 18.2% F431L (TTT→CTT)  mdoG → glucan biosynthesis protein, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090481,112,2020TC18.2% 69.4 / 11.4 33F431L (TTT→CTT) mdoGglucan biosynthesis protein, periplasmic
Reads supporting (aligned to +/- strand):  ref base T (13/14);  new base C (3/3);  total (16/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

TCAACGGGGGATGTGAAGCAGTCGAACCTGATTCGCCAGCCTGACGGTACTATCGCCTTTGTGGTCGATTTTACCGGCGCAGAGATGAAAAAACTGCCAGAGGATACCCCGGTCACAGCGCAAACCAGCATTG  >  AP009048/1112133‑1112265
                                                                     |                                                               
tCAACGGGGGATGTGAAGCAGTCGAACCTGATTCGCCAGCCTGACGGTACTATCGCCTTTGTGGTCGAttt                                                                >  1:1097231/1‑71 (MQ=255)
 cAACGGGGGATGTGAAGCAGTCGAACCTGATTCGCCAGCCTGACGGTACTATCGCCTTTGTGGTCGAttt                                                                <  1:769354/70‑1 (MQ=255)
 cAACGGGGGATGTGAAGCAGTCGAACCTGATTCGCCAGCCTGACGGTACTATCGCCTTTGTGGTCGAttt                                                                >  1:1479918/1‑70 (MQ=255)
  aaCGGGGGATGTGAAGCAGTCGAACCTGATTCGCCAGCCTGACGGTACTATCGCCTTTGTGGTCGAtctt                                                               <  1:2816408/70‑1 (MQ=255)
    cGGGGGATGTGAAGCAGTCGAACCTGATTCGCCAGCCTGACGGTACTATCGCCTTTGTGGTCGATTTTAcc                                                            <  1:2185357/71‑1 (MQ=255)
    cGGGGGATGTGAAGCAGTCGAACCTGATTCGCCAGCCTGACGGTACTATCGCCTTTGTGGTCGATCTTAcc                                                            >  1:3648251/1‑71 (MQ=255)
     gggggATGTGAAGCAGTCGAACCTGATTCGCCAGCCTGACGGTACTATCGCCTTTGTGGTCGATTTTACCg                                                           <  1:152441/71‑1 (MQ=255)
                aGCAGTCGAACCTGATTCGCCAGCCTGACGGTACTATCGCCTTTGTGGTCGATTTTACCGGCGCAGAGAt                                                 >  1:4249239/1‑70 (MQ=255)
                aGCAGTCGAACCTGATTCGCCAGCCTGACGGTACTATCGCCTTTGTGGTCGATTTTACCGGCGCAGAGATg                                                <  1:3892481/71‑1 (MQ=255)
                   aGTCGAACCTGATTCGCCAGCCTGACGGTACTATCGCCTTTGTGGTCGATTTTACCGGCGCagag                                                   >  1:1967106/1‑65 (MQ=255)
                     tCGAACCTGATTCGCCAGCCTGACGGTACTATCGCCTTTGTGGTCGAttt                                                                <  1:1889441/50‑1 (MQ=255)
                     tCGAACCTGATTCGCCAGCCTGACGGTACTATCGCCTTTGTGGTCGATCTTACCGGCGCAGAGATGaaa                                             <  1:505406/69‑1 (MQ=255)
                       gAACCTGATTCGCCAGCCTGACGGTACTATCGCCTTTTTGGTCGATCTTAc                                                             >  1:3717523/1‑51 (MQ=255)
                       gAACCTGATTCGCCAGCCTGACGGTACTATCGCCTTTGTGGTCGATTTTACCGGCGCAGAGATGaa                                              >  1:3484974/1‑66 (MQ=255)
                       gAACCTGATTCGCCAGCCTGACGGTACTATCGCCTTTATGGTCGATCTTACCGGCGCAGAGATGaaaaaa                                          <  1:2960921/70‑1 (MQ=255)
                        aaCCTGATTCGCCAGCCTGACGGTACTATCGCCTTTGTGGTCGATTTTACCGGCGCagag                                                   <  1:1698283/60‑1 (MQ=255)
                         aCCTGATTCGCCAGCCTGACGGTACTATCGCCTTTGTGGTCGATTTTAc                                                             >  1:353618/1‑49 (MQ=255)
                         aCCTGATTCGCCAGCCTGACGGTACTATCGCCTTTGTGGTCGATTTTACCGGCGCAGAGATGAAAAAACTg                                       >  1:1234099/1‑71 (MQ=255)
                          ccTGATTCGCCAGCCTGACGGTACTATCGCCTTTGTGGTCGATTTTACCGGCGCAGAGATGAAAAAACTGc                                      <  1:344828/71‑1 (MQ=255)
                           cTGATTCGCCAGCCTGACGGTACTATCGCCTTTGTGGTCGATTTTACCGGCGCAGAGATGa                                               <  1:2920584/61‑1 (MQ=255)
                            tGATTCGCCAGCCTGACGGTACTATCGCCTTTGTGGTCGATCTTa                                                              <  1:1947911/45‑1 (MQ=255)
                             gATTCGCCAGCCTGACGGTACTATCGCCTTTGTGGTCGATTTTACCGGCGCAGAGATGAAAAAACTGc                                      >  1:1968600/1‑68 (MQ=255)
                                  gCCAGCCTGACGGTACTATCGCCTTTGTGGTCGATTTTACCGGCGCAGAGATGAAAAAACTGCCAGAGGa                               <  1:1360334/70‑1 (MQ=255)
                                      gCCTGACGGTACTATCGCCTTTGTGGTCGATTTTACCGGCGCAGAGATGAAAAAACTGCCAGAGGATAccc                          >  1:3539523/1‑71 (MQ=255)
                                      gCCTGACGGTACTATCGCCTTTGTGGTCGATTTTACCGGCGCAGAGATGAAAAAACTGCCAGAGGATAc                            <  1:4443471/69‑1 (MQ=255)
                                      gCCTGACGGTACTATCGCCTTTGTGGTCGATCTTAc                                                             >  1:2699247/1‑36 (MQ=255)
                                       ccTGACGGTACTATCGCCTTTGTGGTCGATTTTACCGGCGCAGAGATGAAAAAACTGCCAGAGGATAccc                          <  1:3978331/70‑1 (MQ=255)
                                        cTGACGGTACTATCGCCTTTGTGGTCGATTTTACCGGCGCAGAGATGa                                               <  1:635508/48‑1 (MQ=255)
                                             ggTACTATCGCCTTTGTGGTCGATTTTACCGGCGCAGAGATGAAAAAACTGCCAGAGGATACCCCGGTca                    >  1:4242756/1‑70 (MQ=255)
                                               tACTATCGCCTTTGTGGTCGATTTTACCGGCGCAGAGATGAAAAAACTGCCAGAGGATACCCCGGTCACAg                 >  1:4608646/1‑71 (MQ=255)
                                                aCTATCGCCTTTGTGGTCGATTTTACCGGCGCAGAGATGAAAAAACTGCCAGAGGATACCCCGGTCACAgc                >  1:3752415/1‑71 (MQ=255)
                                                  tATCGCCTTTGTGGTCGATTTTACCGGCGCAGAGATGAAAAAACTGCCagag                                 <  1:2321914/52‑1 (MQ=255)
                                                      gCCTTTGTGGTCGATTTTACCGGCGCAGAGATGAAAAAACTGCCAGAGGATACCCCGGTCACAgcgc              >  1:2538274/1‑67 (MQ=255)
                                                       ccTTTGTGGTCGATTTTACCGGCGCAGAGATGAAAAAACTGCCa                                    <  1:470921/44‑1 (MQ=255)
                                                        cTTTGTGGTCGATTTTACCGGCGCAGAGATGAAAAAACTGCCAGAGGATACCCCGGTCACAGCGCAAAc          >  1:4435430/1‑69 (MQ=255)
                                                           tgtgGTCGATTTTACCGGCGCAGAGATGAAAAAACTGCCAGAGGa                               >  1:3193271/1‑45 (MQ=255)
                                                           tgtgGTCGATTTTACCGGCGCAGAGATGAAAAAACTGCCAGAGGATACCCCGGTCACAGCGCAAACCAGc      <  1:2676491/70‑1 (MQ=255)
                                                               gTCGATTTTACCGGCGCAGAGATGAAAAAACTGCCAGAGGATACCCCGGTCACAGCGCAAACCAGCATTg  <  1:901982/70‑1 (MQ=255)
                                                                     |                                                               
TCAACGGGGGATGTGAAGCAGTCGAACCTGATTCGCCAGCCTGACGGTACTATCGCCTTTGTGGTCGATTTTACCGGCGCAGAGATGAAAAAACTGCCAGAGGATACCCCGGTCACAGCGCAAACCAGCATTG  >  AP009048/1112133‑1112265

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: