Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 222,988 T→G 49.0% intergenic (+43/‑137) lpcA → / → prfH D‑sedoheptulose 7‑phosphate isomerase/predicted peptide chain release factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported222,9880TG49.0% 48.7 / 106.1 69intergenic (+43/‑137)lpcA/prfHD‑sedoheptulose 7‑phosphate isomerase/predicted peptide chain release factor
Reads supporting (aligned to +/- strand):  ref base T (24/11);  new base G (14/19);  total (38/31)
Fisher's exact test for biased strand distribution p-value = 4.99e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTGATTGAAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTA  >  Exported/222918‑223056
                                                                      |                                                                    
gTTGATTGAAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCAt                                                                      <  1:3064617/71‑1 (MQ=255)
gTTGATTGAAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCAt                                                                      <  1:2530753/71‑1 (MQ=255)
 ttgattgaAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCTGCAtt                                                                     <  1:41700/71‑1 (MQ=255)
 ttgattgaAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCAtt                                                                     <  1:570566/71‑1 (MQ=255)
         aaaaGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGct                                                             <  1:1090087/71‑1 (MQ=255)
         aaaaGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGct                                                             <  1:3333644/71‑1 (MQ=255)
         aaaaGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGct                                                             >  1:155804/1‑71 (MQ=255)
          aaaGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctg                                                            <  1:764669/71‑1 (MQ=255)
            agagATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgat                                                          >  1:3852053/1‑71 (MQ=255)
            agagATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgat                                                          >  1:1296964/1‑71 (MQ=255)
              agaTGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTgcg                                                               >  1:1077818/1‑64 (MQ=255)
                 tGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgatgcctg                                                     <  1:1196182/71‑1 (MQ=255)
                   gTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCAttt                                                                    >  1:3970016/1‑54 (MQ=255)
                       agtaagtCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCg                                                <  1:3505789/70‑1 (MQ=255)
                           agtCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCt                                           <  1:1655151/71‑1 (MQ=255)
                              cTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGAc                                        >  1:434380/1‑71 (MQ=255)
                                  tgTAGGCCGGACAAGACGTTTA‑‑TCGCATCCGGCAGTCGTGTGcagatgcct                                                      <  1:2467046/50‑1 (MQ=9)
                                  tgTAGGCCGGACAAGACGTTTA‑‑TCGCATCCGGCAGTCGTGTGCAGATGCCTGATGCGACGCTAACGCGTCt                                  <  1:3991590/70‑1 (MQ=9)
                                  tgTAGGCCGGACAAGACGTTTA‑‑TCGCATCCGGCAGTCGTGTGCAGATGCCTGATGCGACGCTAACGCGTCt                                  <  1:4325091/70‑1 (MQ=9)
                                      ggCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTgcgc                                                              >  1:1101933/1‑41 (MQ=255)
                                      ggCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTa                                >  1:4039942/1‑71 (MQ=255)
                                      ggCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTa                                <  1:3968490/71‑1 (MQ=255)
                                      ggCCGGACAAGACGTTTA‑‑TCGCATCCGGCAGTCGTGTGCAGATGCCTGATGCGACGCTAACGCGTCTTAt                               <  1:4582333/70‑1 (MQ=9)
                                      ggCCGGACAAGACGTTTA‑‑TCGCATCCGGCAGTCGTGTGCAGATGCCTGATGCGACGCTAACGCGTCTTAt                               <  1:3118110/70‑1 (MQ=9)
                                      ggCCGGACAAGACGTTTA‑‑TCGCATCCGGCAGTCGTGTGCAGATGCCTGATGCGACGCTAACGCGTCTTAt                               <  1:4683535/70‑1 (MQ=9)
                                      ggCCGGACAAGACGTTTA‑‑TCGCATCCGGCAGTCGTGTGCAGATGCCTGATGCGACGCTAACGCGTCTTAt                               >  1:2766014/1‑70 (MQ=9)
                                      ggCCGGACAAGACGTTTA‑‑TCGAATCCGGCAGTCGTGTGCAGATGCCTGATGCGACGCTAACGCGTCTTAt                               <  1:834604/70‑1 (MQ=2)
                                       gCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTa                                >  1:4619478/1‑70 (MQ=11)
                                         cGGATAAGGCGTTTACGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGt                                    <  1:2226868/64‑1 (MQ=14)
                                          ggATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGcagatgcctgatgc                                                 <  1:1756742/50‑1 (MQ=11)
                                          ggATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGcagatgcctgatgc                                                 <  1:422928/50‑1 (MQ=11)
                                          ggATAAGGCGTTTACGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGTCTTAt                               >  1:1063915/1‑68 (MQ=14)
                                          ggACAAGACGTTTA‑‑TCGCATCCGGCAGTCGTGTGCAGATGCCTGATGCGACGCTAACGCGTCTTATc                              <  1:3573976/67‑1 (MQ=9)
                                           gATAAGGCGTTTACGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGTCTTATCagg                           >  1:4350560/1‑69 (MQ=11)
                                           gATAAGGCGTTTACGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGTCTTATCagg                           >  1:5554197/1‑69 (MQ=11)
                                           gATAAGGCGTTTACGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGTCTTATCagg                           >  1:824945/1‑69 (MQ=11)
                                           gATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgatgcct                                                      <  1:1701772/44‑1 (MQ=255)
                                            aTAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgatgcctga                                                    >  1:2609340/1‑45 (MQ=255)
                                            aTAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgatgcctga                                                    >  1:1805980/1‑45 (MQ=255)
                                            aTAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgatgcctga                                                    >  1:3616525/1‑45 (MQ=255)
                                             tAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGcc                         <  1:3814930/71‑1 (MQ=11)
                                             tAAGACGTTTACGCCGCATCCGGCATTTGTGCGctga                                                           >  1:4290447/1‑37 (MQ=255)
                                               aggcGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTa                       <  1:1875605/68‑1 (MQ=11)
                                               aggcGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCt                        <  1:75157/67‑1 (MQ=11)
                                               aGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTa                       >  1:989136/1‑71 (MQ=255)
                                               aGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTa                       >  1:3666333/1‑71 (MQ=255)
                                                ggcGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTAc                      <  1:5436686/69‑1 (MQ=11)
                                                   gTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACaaa                   >  1:1943052/1‑71 (MQ=14)
                                                   gTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATg                           >  1:1076548/1‑63 (MQ=255)
                                                   gTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGc                          >  1:2237443/1‑64 (MQ=255)
                                                   gTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACaaa                   >  1:2742367/1‑71 (MQ=255)
                                                    tttATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATc                              <  1:3095830/59‑1 (MQ=11)
                                                    ttcaCGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGTCTTATCAGACCTACAAAt                  <  1:316021/68‑1 (MQ=11)
                                                     ttATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACaaact                 >  1:187324/1‑69 (MQ=11)
                                                     ttACGCCGCATCCGGCATTTGTGCGATGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATc                 >  1:2157744/1‑71 (MQ=255)
                                                      tACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCt                >  1:3034541/1‑71 (MQ=255)
                                                      tACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCt                >  1:5551855/1‑71 (MQ=255)
                                                      tACGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGc                                            >  1:5012285/1‑43 (MQ=14)
                                                      caCGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGTCTTATCAGGCCTACAAAtt                 >  1:3545990/2‑69 (MQ=11)
                                                      caCGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGTCTTATCAGACCTACAAAttc                >  1:4035262/2‑69 (MQ=11)
                                                      caCGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGTCTTATCAGACCTACAAAttc                >  1:2380082/2‑69 (MQ=11)
                                                        cgccgcATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGTCTTATCAGGCCTACAAAttct               <  1:3513868/70‑4 (MQ=11)
                                                          ccgcATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTAc            >  1:4984833/1‑71 (MQ=255)
                                                          ccgcATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTAc            <  1:5228399/71‑1 (MQ=255)
                                                           cgcATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTAcg           <  1:4007667/71‑1 (MQ=255)
                                                           cgcATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGt                                    >  1:557805/1‑46 (MQ=14)
                                                            gcATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTAcgc          >  1:3970778/1‑71 (MQ=255)
                                                            gcATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGTCTTATCAGGCCTACa                     <  1:1360600/60‑1 (MQ=11)
                                                             cATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTAcg                     <  1:1850183/58‑2 (MQ=14)
                                                               tCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGTCTTATCAGGCCTACAAAttc                >  1:1944324/1‑60 (MQ=11)
                                                                    cATTTGTGCGCTGATGCCTGATGCGACGCTGCCGCGTCTTATCATGCCTACAAATCTGTACGCGAACCg    >  1:1816658/1‑69 (MQ=255)
                                                                    cATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCg    >  1:2460913/1‑69 (MQ=255)
                                                                    cATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCg    >  1:374981/1‑69 (MQ=255)
                                                                    cATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTa  <  1:1381827/71‑1 (MQ=255)
                                                                      tttGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACaa                    <  1:1505484/51‑1 (MQ=255)
                                                                      |                                                                    
GTTGATTGAAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTA  >  Exported/222918‑223056

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: