Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 965,844 A→T 8.5% V17D (GTT→GAT)  ydeM ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported965,8440AT 8.5% 165.6 / 6.9 59V17D (GTT→GAT) ydeMconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (25/29);  new base T (3/2);  total (28/31)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.25e-01

CAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAAACACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAGGC  >  Exported/965777‑965914
                                                                   |                                                                      
cAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAAaca                                                                      >  1:751341/1‑70 (MQ=255)
cAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAAaca                                                                      >  1:2178624/1‑70 (MQ=255)
cAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAAaca                                                                      >  1:1745797/1‑70 (MQ=255)
 agagTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAAacac                                                                     <  1:882199/70‑1 (MQ=255)
   agTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAAACACAtt                                                                  <  1:3197065/71‑1 (MQ=255)
            atccatccATTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTa                                                          >  1:1935879/1‑70 (MQ=255)
            atccatccATTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAAACACATTTGAGATTa                                                          >  1:369230/1‑70 (MQ=255)
            atccatccATTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAAACACATTTGAGATTa                                                          >  1:2648572/1‑70 (MQ=255)
            atccatccATTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAAACACATTTGAGATTa                                                          >  1:2167239/1‑70 (MQ=255)
            atccatccATTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAAACACATTTGAGATTAc                                                         <  1:674098/71‑1 (MQ=255)
                atccatTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAACCAGTAAACACATTTGAGATTACAtt                                                      >  1:3287565/1‑70 (MQ=255)
                atccatTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAAACACATTTGAGATTACAtt                                                      >  1:2344534/1‑70 (MQ=255)
                 tccatTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAAACACATTTGAGATTACATTGa                                                    <  1:760545/71‑1 (MQ=255)
                   catTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAAACACATTTGAGATTACATTGaa                                                   <  1:2360795/70‑1 (MQ=255)
                   catTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAAACACATTTGAGATTACATTGaa                                                   >  1:1839369/1‑70 (MQ=255)
                    atTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAAACACATTTGAGATTACATTGaaa                                                  >  1:1372078/1‑70 (MQ=255)
                      ttttCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTaaa                                                                        >  1:1746909/1‑46 (MQ=255)
                      ttttCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTaaa                                                                        >  1:1325237/1‑46 (MQ=255)
                      ttttCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTaaa                                                                        >  1:370848/1‑46 (MQ=255)
                      ttttCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAAACACATTTGAGATTACATTGAAAACt                                               >  1:3261298/1‑71 (MQ=255)
                      ttttCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAAACACATTTGAGATTACATTGAAAACt                                               >  1:2487386/1‑71 (MQ=255)
                      atttCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAAACACATTTGAGAtt                                                           >  1:1600374/2‑59 (MQ=255)
                       tttCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAAACACATTTGAGATTAc                                                         >  1:2132618/1‑60 (MQ=255)
                       tttCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAAACACATTTGAGATTACATTGAAAACTg                                              >  1:237225/1‑71 (MQ=255)
                        ttCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAAACACATTTGAGATTACATTGAAAACTg                                              <  1:1843632/70‑1 (MQ=255)
                          cATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAAACACATTTGAGATTACATTGAAAACTGGa                                            <  1:1914080/70‑1 (MQ=255)
                            tGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAAACACTTTTGAGATTACATTGAAAACTGGAggg                                         <  1:764872/71‑1 (MQ=255)
                              aGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAAACACATTTGAGATTACATTg                                                     <  1:198663/57‑1 (MQ=255)
                              aGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAAACACATTTGAGATTACATTg                                                     <  1:426085/57‑1 (MQ=255)
                               gTAAACTGCGACTCGTTTTCAAGGTAAAAACAGTAAACACATTTGAGATTACATTGAAAACt                                               <  1:3031792/62‑1 (MQ=255)
                                tAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAAACACATTTGAGATTACATTGAAAACTGGAGGGCTTg                                     >  1:374297/1‑71 (MQ=255)
                                  aaCTGCGACTCTTTTTCAAGGTAAAAACAGTAAACACATTTga                                                               >  1:1384820/1‑43 (MQ=255)
                                  aaCTGCGACTCTTTTTCAAGGTAAAAACAGTAAACACATTTGAGATTACATTGAAAACTGGAGGGCTTgg                                    <  1:1870268/70‑1 (MQ=255)
                                  aaCTGCAACTCTTTTTCAAGGTAAAAACAGTAAACACATTTGAGATTACATTGAAAACTGGAGGGCTTgg                                    <  1:3471616/70‑1 (MQ=255)
                                  aaCAGCGACTCTTTTTCAAGGTAAAAACAGTAAACACATTTGAGATTACATTGAAAACTGGAGGGCTTgg                                    <  1:588433/70‑1 (MQ=255)
                                    cTGCGACTCTTTTTCAAGGTAAAAACAGTAAACACATTTgaga                                                             <  1:2701781/43‑1 (MQ=255)
                                     tGCGACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTg                                 >  1:12140/1‑70 (MQ=255)
                                     tGCGACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTg                                 >  1:635983/1‑70 (MQ=255)
                                     tGCGACTCTTTTTCAAGGTAAAAACAGTAAACACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTg                                 >  1:2056943/1‑70 (MQ=255)
                                       cGACTCTTTTTCAATGTAAAAACAGTAAACACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGt                                <  1:1806328/69‑1 (MQ=255)
                                       cGACTCTTTTTCAAGGTAAAAACAGTAAACACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGt                                <  1:2402059/69‑1 (MQ=255)
                                       cGACTCTTTTTCAAGGTAAAAACAGTAAACACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTaa                              <  1:127945/71‑1 (MQ=255)
                                        gACTCTTTTTCAAGGTAAAAACAGTAAACACATTTGAGATTACATTGAAAACTGGAGGGCTTGTCTGTAAc                             <  1:2341941/71‑1 (MQ=255)
                                        gACTCTTTTTCAAGGTAAAAACAGTAAACACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAAc                             <  1:80184/71‑1 (MQ=255)
                                         aCTCTTTTTCAAGGTAAAAACAGTAAACACATTTGAGATTACATTg                                                     >  1:1253898/1‑46 (MQ=255)
                                          ctctTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGt                           <  1:1096625/71‑1 (MQ=255)
                                          ctctTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGt                           <  1:1005632/71‑1 (MQ=255)
                                               tttCAAGGTAAAAACAGTAAACACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCa                        <  1:1462752/69‑1 (MQ=255)
                                               tttCAAGGTAAAAACAGTAAACACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATc                      >  1:536779/1‑71 (MQ=255)
                                               tttCAAGGTAAAAACAGTAAACACATCTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCAt                       <  1:2963455/70‑1 (MQ=255)
                                                ttCAAGGTAAAAACAGTAAACACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCg                     >  1:1673503/1‑71 (MQ=255)
                                                 tCAAGGTAAAAACAGTAAACACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCg                     <  1:3532614/70‑1 (MQ=255)
                                                 tCAAGGTAAAAACAGTAAACACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCg                     <  1:3109076/70‑1 (MQ=255)
                                                   aaGGTAAAAACAGTAAACACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCt                   <  1:2208429/70‑1 (MQ=255)
                                                      gTAAAAACAGTAAACACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTAtct               >  1:96631/1‑71 (MQ=255)
                                                        aaaaaCAGTAAACACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCg             >  1:3219976/1‑71 (MQ=255)
                                                        aaaaaCAGTAAACACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCg             >  1:1910514/1‑71 (MQ=255)
                                                          aaaCAGTAAACACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCt           >  1:779310/1‑71 (MQ=255)
                                                           aaCAGTAAACACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCt           >  1:862898/1‑70 (MQ=255)
                                                           aaCAGTAAACACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTc          <  1:1047455/71‑1 (MQ=255)
                                                           aaCAGTAAACACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTc          <  1:931296/71‑1 (MQ=255)
                                                           aaCAGTAAACACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTc          <  1:682887/71‑1 (MQ=255)
                                                            aCAGTAAACACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCa         <  1:2487190/71‑1 (MQ=255)
                                                             cAGTAAACACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCaa        <  1:3088575/71‑1 (MQ=255)
                                                              aGTAAACACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCaa        >  1:3294792/1‑70 (MQ=255)
                                                              aGTAAACACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCaa        >  1:1130023/1‑70 (MQ=255)
                                                                  aaCACATTTGAGATTACATTGAAAACTGGAGGTCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAgg   <  1:2945443/71‑1 (MQ=255)
                                                                  aaCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAgg   <  1:61485/71‑1 (MQ=255)
                                                                   acacATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAggc  <  1:1421602/71‑1 (MQ=255)
                                                                   |                                                                      
CAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAAACACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAGGC  >  Exported/965777‑965914

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