Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,547,976 T→C 10.2% V62A (GTT→GCT)  htrL → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,547,9760TC10.2% 152.7 / 9.3 59V62A (GTT→GCT) htrLhypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (23/30);  new base C (3/3);  total (26/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.23e-01

TACTTTGAGCGTCTGGCCGCACTTGAAAATGAGATGATTATTTTCACTTCTCCCGATTTGAAACCAAGGGTTGAAGCGATTCGCAACGGAAAACCAACAACGGTTATCGTTATCGATATAAAAAAGAAATTTAGATATAT  >  Exported/2547906‑2548045
                                                                      |                                                                     
tACTTTGAGCGTCTGGCCGCACTTGAAAATGAGATGATTATTTTCACTTCTCCCGATTTGAAACCAAGGGt                                                                       >  1:1521461/1‑71 (MQ=255)
  cTTTGAGCGTCTGGCCGCACTTGAAAATGAGATGATTATTTTCACTTCTCCCGATTTGAAACCAAGGGTTg                                                                     <  1:2299004/71‑1 (MQ=255)
  cTTTGAGCGTCTGGCCGCACTTGAAAATGAGATGATTATTTTCACTTCTCCCGATTTGAAACCAAGGGTTg                                                                     <  1:1199100/71‑1 (MQ=255)
    ttGAGCGTCTGGCCGCACTTGAAAATGAGATGATTATTTTCACTTCTCCCGATTTGAAACCAAGGGTTg                                                                     >  1:920740/1‑69 (MQ=255)
     tGAGCGTCTGGCCGCACTTGAAAATGAGATGATTATTTTCACTTCTCCCGATTTGAAACCAAGGGTTGa                                                                    <  1:1534327/69‑1 (MQ=255)
     tGAGCGTCTGGCCGCACTTGAAAATGAGATGATTATTTTCACTTCTCCCGATTTGAAACCAAGGGTTGAAg                                                                  >  1:1963337/1‑71 (MQ=255)
     tGAGCGTCTGGCCGCACTTGAAAATGAGATGATTATTTTCACTTCTCCCGATTTGAAACCAAGGGCTGAAg                                                                  >  1:1163719/1‑71 (MQ=255)
      gAGCGTCTGGCCGCACTTGAAAATGAGATGATTATTTTCACTTCTCCCGATTTGAAACCAAGGGCTGa                                                                    <  1:1093411/68‑1 (MQ=255)
       aGCGTCTGGCCGCACTTGAAAATGAGATGATTATTTTCACTTCTCCCGATTTGAAACCAAGGGTTGAAGCg                                                                <  1:1999098/71‑1 (MQ=255)
        gCGTCTGGCCGCACTTGAAAATGAGATGATTATTTTCACTTCTCCCGATTTGAAACCAAGGGTTGAAGCGa                                                               <  1:1792169/71‑1 (MQ=255)
          gTCTGGCCGCACTTGAAAATGAGATGATTATTTTCACTTCTCCCGATTTGAAACCAAGGGTTGAAGCGAtt                                                             >  1:3208507/1‑71 (MQ=255)
          gTCTGGCCGCACTTGAAAATGAGATGATTATTTTCACTTCTCCCGATTTGAAACCAAGGGTTGAAGCGAtt                                                             <  1:1127940/71‑1 (MQ=255)
               gCCGCACTTGAAAATGAGATGATTATTTTCACTTCTCCCGATTTGAAACCAAGGGTTGAAGCGATTCg                                                           <  1:3248024/68‑1 (MQ=255)
               gCCGCACTTGAAAATGAGATGATTATTTTCACTTCTCCCGATTTGAAACCAAGGGTTGAAGCGATTCg                                                           <  1:1004863/68‑1 (MQ=255)
               gCCGCACTTGAAAATGAGATGATTATTTTCACTTCTCCCGATTTGAAACCAAGGGCTGAAGCGATTCGCaa                                                        <  1:79609/71‑1 (MQ=255)
               gCCGCACTTGAAAATGAGATGATTATTTTCACTTCTCCCGATTTGAAACCAAGGGCTGAAGCGATTCGCaa                                                        <  1:506011/71‑1 (MQ=255)
                ccGCACTTGAAAATGAGATGATTATTTTCACTTCTCCCGATTTGAAACCAAGGGTTGaa                                                                   >  1:1134183/1‑59 (MQ=255)
                ccGCACTTGAAAATGAGATGATTATTTTCACTTCTCCCGATTTGAAACCAAGGGTTGAAGCGATTCGCa                                                         <  1:2833034/69‑1 (MQ=255)
                 cGCACTTGAAAATGAGATGATTATTTTCACTTCTCCCGATTTGAAACCAAGGGTTGAAGCGATTCGCAACg                                                      >  1:2923626/1‑71 (MQ=255)
                    aCTTGAAAATGAGATGATTATTTTCACTTCTCCCGATTTGAAACCAAGGGTTg                                                                     <  1:3624474/53‑1 (MQ=255)
                    aCTTGAAAATGAGATGATTATTTTCACTTCTCCCGATTTGAAACCAAGGGTTGAAGCGATTCGCAACGGa                                                    >  1:2639375/1‑70 (MQ=255)
                       tGAAAATGAGATGATTATTTTCACTTCTCCCGATTTGAAACCAAGGGTTGAAGCGATTCGCAACGGaaaa                                                 >  1:251755/1‑70 (MQ=255)
                         aaaaTGAGATGATTATTTTCACTTCTCCCGATTTGAAACCAAGGGTTGAAGCGATTCGCAACGGAAAAcc                                               <  1:2485127/70‑1 (MQ=255)
                          aaaTGAGATGATTATTTTCACTTCTCCCGATTTGAAACCAAGGGTTGAAGCGATTCGCAACGGAAAACcaa                                             >  1:688336/1‑71 (MQ=255)
                            aTGAGATGATTATTTTCACTTCTCCCGATTTGAAACCAAGGGTTGAAGCGATTCGCAACGGAAAACcaaca                                           <  1:2350586/71‑1 (MQ=255)
                             tGAGATGATTATTTTCACTTCTCCCGATTTGAAACCAAGGGTTGAAGCGATTCGCAACGGAAAACcaacaa                                          <  1:1317088/71‑1 (MQ=255)
                             tGAGATGATTATTTTCACTTCTCCCGATTTGAAACCAAGGGTTGAAGCGATTCGCAACGGAAAACcaacaa                                          <  1:1771488/71‑1 (MQ=255)
                                gatgatTATTTTCACTTCTCCCGATTTGAAACCAAGGGTTGAAGCGATTCGCAACGGaa                                                   <  1:1303391/59‑1 (MQ=255)
                                    attattTTCACTTCTCCCGATTTGAAACCAAGGGTTGAAGc                                                                 <  1:2051976/41‑1 (MQ=255)
                                     ttattTTCACTTCTCCCGATTTGAAACCAAGGGCTGAAGCGATTCGCAACGGAAAACCAACAACGgttatc                                  >  1:1087910/1‑71 (MQ=255)
                                     ttattTTCACTTCTCCCGATTTGAAACAAAGGGCTGAAGCGATTCGCAACGGAAAACCAACAACGgttat                                   >  1:3713881/1‑70 (MQ=255)
                                      tattTTCACTTCTCCCGATTTGAAACCAAGGGTTGAAGCGATTCGCAAc                                                       >  1:2304969/1‑49 (MQ=255)
                                      tattTTCACTTCTCCCGATTTGAAACCAAGGGTTGAAGCGATTCGCAAc                                                       >  1:399265/1‑49 (MQ=255)
                                        ttttCACTTCTCCCGATTTGAAACCAAGGGTTGAAGCGATTCGCAACGGAAAACCAACAACGgttatcgtt                               <  1:2274495/71‑1 (MQ=255)
                                             aCTTCTCCCGATTTGAAACCAAGGGTTGAAGCGATTCGCAACGGAAAACCAACAACGGttatcgttatcg                           <  1:2917412/70‑1 (MQ=255)
                                             aCTTCTCCCGATTTGAAACCAAGGGTTGAAGCGATTCGCAACGGAAAACCAACAACGGttatcgttatcg                           <  1:2976768/70‑1 (MQ=255)
                                              cTTCTCCCGATTTGAAACCAAGGGTTGAAGCGATTCGCAACGGaa                                                   <  1:2094215/45‑1 (MQ=255)
                                                 ctcCCGATTTGAAACCAAGGGTTGAAGCGATTCGCaa                                                        >  1:2286581/1‑37 (MQ=255)
                                                 ctcCCGATTTGAAACCAAGGGTTGAAGCGATTCGCaa                                                        >  1:2809559/1‑37 (MQ=255)
                                                   cccGATTTGAAACCAAGGGTTGAAGCGATTCGCAACGGAAAACCAACAACGGTTATCGTTATCGATATaa                     >  1:3235385/1‑70 (MQ=255)
                                                    ccGATTTGAAACCAAGGGTTGAAGCGATTCGCAACGGAAAACCAACAACGGTTATCGTTATCGATATaaa                    <  1:1104372/70‑1 (MQ=255)
                                                       aTTTGAAACCAAGGGTTGAAGCGATTCGCAACGGAAAACCAACAACGGTTATCGTTATCGATATaaaa                   >  1:1517563/1‑68 (MQ=255)
                                                              aCCAAGGGTTGAAGCGATTCGCAACGGAAAACCAACAACGGTTATCGTTATCGATATAAAAAAGaaa             >  1:3675307/1‑67 (MQ=255)
                                                              aCCAAGGGTTGAAGCGATTCGCAACGGAAAACCAACAACGGTTATCGTTATCGATATAAAAAAGaaa             >  1:1535055/1‑67 (MQ=255)
                                                              aCCAAGGGTTGAAGCGATTCGCAACGGAAAACCAACAACGGTTATCGTTATCGATATAAAAAAGAAATTTa         <  1:520090/71‑1 (MQ=255)
                                                              aCCAAGGGTTGAAGCGATTCGCAACGGAAAACCAACAACGGTTATCGTTATCGATATAAAAAAGAAATTTa         <  1:3706059/71‑1 (MQ=255)
                                                              aCCAAGGGTTGAAGCGATTCGCAACGGAAAACCAACAACGGTTATCGTTATCGATATAAAAAAGAAATTTa         <  1:174264/71‑1 (MQ=255)
                                                              aCCAAGGGTTGAAGCGATTCGCAACGGAAAACCAACAACGGTTATCGTTATCGATATAAAAAAGAAATTTa         <  1:1721465/71‑1 (MQ=255)
                                                               ccAAGGGTTGAAGCGATTCGCAACGGAAAACCAACAACGGTTATCGTTATCGATATAAAAAAGAAAtt           >  1:108095/1‑68 (MQ=255)
                                                               ccAAGGGTTGAAGCGATTCGCAACGGAAAACCAACAACGGTTATCGTTATCGATATAAAAAAGAAATTTAg        <  1:502940/71‑1 (MQ=255)
                                                                 aaGGGTTGAAGCGATTCGCAACGGAAAACCAACAACGgttatc                                  >  1:3081744/1‑43 (MQ=255)
                                                                 aaGGGTTGAAGCGATTCGCAACGGAAAACCAACAACGgttatc                                  >  1:1744794/1‑43 (MQ=255)
                                                                 aaGGGTTGAAGCGATTCGCAACGGAAAACCAACAACGGTTATCGTTATCGATATAAAAAAGaa              <  1:538062/63‑1 (MQ=255)
                                                                   gggTTGAAGCGATTCGCAACGGAAAACCAACAACGGTTATCGTTATCGATATAAAAAAGAAATTTAGata     >  1:1858329/1‑70 (MQ=255)
                                                                   gggTTGAAGCGATTCGCAACGGAAAACCAACAACGGTTATCGTTATCGATATAAAAAAGAAATTTAGata     >  1:1475894/1‑70 (MQ=255)
                                                                    ggTTGAAGCGATTCGCAACGGAAAACCAACAACGgttatcgttat                             <  1:3239907/45‑1 (MQ=255)
                                                                    ggTTGAAGCGATTCGCAACGGAAAACCAACAACGGttatcgttatcg                           >  1:2680820/1‑47 (MQ=255)
                                                                    ggTTGAAGCGATTCGCAACGGAAAACCAACAACGGttatcgttatcg                           >  1:2490042/1‑47 (MQ=255)
                                                                     gTTGAAGCGATTCGCAACGGAAAACCAACAACGGTTATCGTTATCGATATAAAAAAGAAATTTAGatatat  <  1:1924018/71‑1 (MQ=255)
                                                                     gTTGAAGCGATTCGCAACGGAAAACAACAACGGTTATCGTTATCGATATAAAAAAGAAATTTAGatatat  <  1:961313/70‑1 (MQ=255)
                                                                      ttGAAGCGATTCGCAACGGAAAACCAACAACGgttatcgtt                               <  1:2436147/41‑1 (MQ=255)
                                                                      ttGAAGCGATTCGCAACGGAAAACCAACAACGgttatc                                  <  1:2742607/38‑1 (MQ=255)
                                                                      |                                                                     
TACTTTGAGCGTCTGGCCGCACTTGAAAATGAGATGATTATTTTCACTTCTCCCGATTTGAAACCAAGGGTTGAAGCGATTCGCAACGGAAAACCAACAACGGTTATCGTTATCGATATAAAAAAGAAATTTAGATATAT  >  Exported/2547906‑2548045

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: