Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,837,495 G→A 17.4% R308R (CGC→CGT nlpD ← predicted outer membrane lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,837,4950GA17.4% 134.5 / 19.7 63R308R (CGC→CGTnlpDpredicted outer membrane lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (31/21);  new base A (4/7);  total (35/28)
Fisher's exact test for biased strand distribution p-value = 1.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

TGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTCGCGATAATTGCCTGTCCTTTGCT  >  Exported/1837426‑1837560
                                                                     |                                                                 
tgtgTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCcgc                                                                  >  1:4724543/1‑71 (MQ=255)
tgtgTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCcg                                                                   >  1:4044360/1‑70 (MQ=255)
 gtgtCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCcgcg                                                                 >  1:3937281/1‑71 (MQ=255)
 gtgtCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCcgcg                                                                 >  1:1290134/1‑71 (MQ=255)
 gtgtCATTATGGGCGTAAGCCGTAATGTAATCTTCACTGTGTTTAATCATGATGAGATTACCGTAGCcacg                                                                 >  1:865233/1‑71 (MQ=255)
 gtgtCATTATGGGCGTAAGCCGTAATGTAATCTTCACTGTGTTTAATCATGATGAGATTACCGTAGCcacg                                                                 >  1:3550347/1‑71 (MQ=255)
 gtgtCATTATGGGCGTAAGCCGTAATGTAATCTTCACTGTGTTTAATCATGATGAGATTACCGTAGCcacg                                                                 >  1:2447017/1‑71 (MQ=255)
 gtgtCATTATGGGCGTAAGCCGTAATGTAATCTTCACTGTGTTTAATCATGATGAGATTACAGTAGCcacg                                                                 >  1:5366931/1‑71 (MQ=255)
    tcattATGGGCGTAAGCCGTAATGTAATCTTCACTGTGTTTAATCATGATGAGATTACCGTAGCCACGCa                                                               <  1:4151842/67‑1 (MQ=255)
    tcattATGGGCGTAAGCCGTAATGTAATCTTCACTGTGTTTAATCATGATGAGATTACCGTAGCCACGCa                                                               <  1:147037/67‑1 (MQ=255)
    tcattATGGGCGTAAGCCGTAATGTAATCTTCACTGTGTTTAATCATGATGAGATTACCGTAGCCACGCAg                                                              >  1:4608661/4‑71 (MQ=255)
    tcattATGGGCGTAAGCCGTAATGTAATCTTCACTGTGTTTAATCATGATGAGATTACCGTAGCCACGCAg                                                              >  1:435145/4‑71 (MQ=255)
    tcattATGGGCGTAAGCCGTAATGTAATCTTCACTGTGTTTAATCATGATGAGATTACCGTAGCCACGCAg                                                              >  1:4552884/4‑71 (MQ=255)
     cattATGGGCGTAAGCCGTAATGTAATCTTCACTGTGTTTAATCATGATGAGATTACCGTAGCCACGCAg                                                              >  1:5488917/3‑70 (MQ=255)
      gTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAgcg                                                            >  1:171463/1‑71 (MQ=255)
      gTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAgcg                                                            >  1:279509/1‑71 (MQ=255)
      gTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAga                                                             >  1:2485721/1‑69 (MQ=255)
      attATGGGCGTAAGCCGTAATGTAATCTTCACTGTGTTTAATCATGATGAGATTACCGTAGCCACGCAgct                                                            <  1:4579251/70‑2 (MQ=255)
       ttATGGGCGTAAGCCGTAATGTAATCTTCACTGTGTTTAATCATGATGAGATTACCGTAGCCACGCAgctg                                                           <  1:2302439/71‑3 (MQ=255)
       ttATGGGCGTAAGCCGTAATGTAATCTTCACTGTGTTTAATCATGATGAGATTACCGTAGCCACGCAgctg                                                           <  1:4749363/71‑3 (MQ=255)
       ttATGGGCGTAAGCCGTAATGTAATCTTCACTGTGTTTAATCATGATGAGATTACCGTAGCCACGCAgctg                                                           <  1:5610259/71‑3 (MQ=255)
       taaTGGGCGTAAGCCGTAATGTAATCTTCACTGTGTTTAATCATGATGAGATTACCGTAGCCACGCa                                                               <  1:1184967/65‑1 (MQ=255)
        tATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAgcgc                                                           >  1:4100491/1‑70 (MQ=255)
          tGGGCGTAGGCACTCAGGTAATCATCATTTTGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGt                                                         >  1:1414191/1‑70 (MQ=255)
          tGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGtt                                                        <  1:2870711/71‑1 (MQ=255)
                  ggCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCa                                                 >  1:2287529/1‑70 (MQ=255)
                   gCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATa                                               <  1:1831448/71‑1 (MQ=255)
                      cTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTAcc                                                     <  1:4423282/62‑1 (MQ=255)
                      cTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAaac                                            <  1:2100593/71‑1 (MQ=255)
                        cAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAaacaa                                          >  1:4076349/1‑71 (MQ=255)
                          ggTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAAcg                                        >  1:1033989/1‑71 (MQ=255)
                           gTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATa                                               <  1:893411/63‑1 (MQ=255)
                             aaTCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAaac                                            >  1:775615/1‑64 (MQ=255)
                             aaTCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAaac                                            >  1:5385861/1‑64 (MQ=255)
                             aaTCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAaac                                            >  1:2230364/1‑64 (MQ=255)
                                catcatTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAg                                                              >  1:2394416/1‑43 (MQ=255)
                                catcatTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAg                                                              >  1:3919408/1‑43 (MQ=255)
                                catcatTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAg                                                              >  1:3969487/1‑43 (MQ=255)
                                catcatTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAaac                                            <  1:1029486/61‑1 (MQ=255)
                                 atcatTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAaacaa                                          >  1:2835029/1‑62 (MQ=255)
                                 atcatTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCAt                                 <  1:43273/71‑1 (MQ=255)
                                    atTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCt                               <  1:5096701/70‑1 (MQ=255)
                                    atTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCt                               <  1:5115252/70‑1 (MQ=255)
                                     ttATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATaaa                                             >  1:4966377/1‑55 (MQ=255)
                                        tGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAaacaa                                          >  1:4051478/1‑55 (MQ=255)
                                         gTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTc                         >  1:2657380/1‑71 (MQ=255)
                                          tttGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTa                                                       <  1:4623498/40‑1 (MQ=255)
                                            tgatgaTAATCAGATTACCGTAGCCGCGCAGCGCGTTAcc                                                     >  1:2486644/1‑40 (MQ=255)
                                            tgatgaTAATCAGATTACCGTAGCCGCGCAGCGCGTTAcc                                                     >  1:3657265/1‑40 (MQ=255)
                                            tgatgaTAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAg                                                   >  1:1107523/1‑42 (MQ=255)
                                            tgatgaTAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTcgc                       <  1:4085857/70‑1 (MQ=255)
                                              atgatAATCAGATTACCGTAGCCGCGCAGCGCGTTAc                                                      >  1:2184482/1‑37 (MQ=255)
                                              atgatAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCg                            <  1:5142663/63‑1 (MQ=255)
                                              atgatAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTCGCGAt                    <  1:1269571/71‑1 (MQ=255)
                                              atgatAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTCGCGAt                    <  1:4127558/71‑1 (MQ=255)
                                               tgatAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTCGCGAt                    <  1:2406421/70‑1 (MQ=255)
                                                gatAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTCGCGATaa                  >  1:2650116/1‑71 (MQ=255)
                                                   aaTCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTCGCGAt                    <  1:2561805/66‑1 (MQ=255)
                                                   aaTCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTCGCGATAATTg               >  1:274092/1‑71 (MQ=255)
                                                    aTCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTTCGGTCGCGATa                   <  1:1361126/66‑1 (MQ=255)
                                                    aTCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTCGCGATa                   <  1:4202495/66‑1 (MQ=255)
                                                     tCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTCGCGATAATTg               >  1:1283768/1‑69 (MQ=255)
                                                      cAGATTACCGTAGCCGCGCAGCGCGTTACCAGCAt                                                >  1:2079547/1‑35 (MQ=255)
                                                       aGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAaaca                                           >  1:2805378/1‑39 (MQ=255)
                                                       aGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAaac                                            >  1:2709289/1‑38 (MQ=255)
                                                       aGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTcg                        >  1:4388101/1‑58 (MQ=255)
                                                           tACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTCGCGATAATTGCCTGTcc        <  1:66452/70‑1 (MQ=255)
                                                            aCCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTCGCGATAATTGCCTGTCCt       <  1:1917216/70‑1 (MQ=255)
                                                             ccGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTCGCGATAATTGCCTGTCCt       <  1:4140410/69‑1 (MQ=255)
                                                                tAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTCGCGATAATTGCCTGTCCTTtgct  >  1:337372/1‑71 (MQ=255)
                                                                  gCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTCGCGATaa                  >  1:4486507/1‑53 (MQ=255)
                                                                    cgcgCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTCGCGATAATTGCCTGTc         >  1:4236742/1‑60 (MQ=255)
                                                                     gcgcAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTCGCGATAATTGCCt            >  1:2050050/1‑56 (MQ=255)
                                                                     gcgcAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTCGCGATAATTGCCTGTCCttt     >  1:1302657/1‑63 (MQ=255)
                                                                     |                                                                 
TGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTCGCGATAATTGCCTGTCCTTTGCT  >  Exported/1837426‑1837560

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: