Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 184,000 G→T 5.4% I32I (ATC→ATA yaeH ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP009048184,0000GT 5.4% 311.8 / 2.2 92I32I (ATC→ATAyaeHconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (39/48);  new base T (3/2);  total (42/50)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

CTGACGCGGATATTTAAACTTAACGCTCTTGGCGAAAAACTCGCCTTTGTCTTTCTGGAAATAGATTTTCAGAATATCGTTGTTGGCTTCCTGCCGGAGGCTGTAACGATCAATTTCTTCAGGATTGGTAATA  >  AP009048/183937‑184069
                                                               |                                                                     
cTGACGCGGATATTTAAACTTAACGCTCTTGGCGAAAAACTCGCCTTTGTCTTTCTGGAAATAGATTTTCa                                                                <  1:5270403/71‑1 (MQ=255)
 tGACGCGGATATTTAAACTTAACGCTCTTGGCGAAAAACTCGCCTTTGTCTTTCTGAAAATCTAtttc                                                                  >  1:6774857/1‑67 (MQ=255)
   aCGCGGATATTTAAACTTAACGCTCTTGGCGAAAAACTCGCCTTTGTCTTTCTGGAAATAGATTTTCAg                                                               <  1:4574374/69‑1 (MQ=255)
   aCGCGGATATTTAAACTTAACGCTCTTGGCGAAAAACTCGCCTTTGTCTTTCTGGAAATAGATTTTCAGaa                                                             <  1:5717656/71‑1 (MQ=255)
   aCGCGGATATTTAAACTTAACGCTCTTGGCGAAAAACTCGCCTTTGTCTTTCTGGAAATAGATTTTCAGaa                                                             <  1:5140408/71‑1 (MQ=255)
         atatTTAAACTTAACGCTCTTGGCGAAAAACTCGCCTTTGTCTTTCTGGAAATAGATTTTCAg                                                               <  1:5482269/63‑1 (MQ=255)
         atatTTAAACTTAACGCTCTTGGCGAAAAACTCGCCTTTGTCTTTCTGGAAATAGATTTTCAg                                                               <  1:1828061/63‑1 (MQ=255)
          tatTTAAACTTAACGCTCTTGGCGAAAAACTCGCCTTTGTCTTTCTGGAAATAGATTTTCAGAATATCgtt                                                      <  1:7501758/71‑1 (MQ=255)
          tatTTAAACTTAACGCTCTTGGCGAAAAACTCGCCTTTGTCTTTCTGGAAATAGATTTTCAGAATATCgtt                                                      <  1:2813205/71‑1 (MQ=255)
           atTTAAACTTAACGCTCTTGGCGAAAAACTCGCCTTTGTCTTTCTGGAAATAGATTTTCAGAATATCgttg                                                     >  1:6345852/1‑71 (MQ=255)
           atTTAAACTTAACGCTCTTGGCGAAAAACTCGCCTTTGTCTTTCTGGAAATAGATTTTCAGAATATCgtt                                                      <  1:2302813/70‑1 (MQ=255)
            tttAAACTTAACGCTCTTGGCGAAAAACTCGCCTTTGTCTTTCTGGAAATAGATTTTCAGAata                                                           <  1:6219053/64‑1 (MQ=255)
            tttAAACTTAACGCTCTTGGCGAAAAACTCGCCTTTGTCTTTCTGGAAATAGATTTTCAGAATATCgttgt                                                    >  1:6766050/1‑71 (MQ=255)
              tAAACTTAACGCTCTTGGCGAAAAACTCGCCTTTGTCTTTCTGGAAATAGATTTTCAGAATATCGttgttg                                                  >  1:7615076/1‑71 (MQ=255)
               aaaCTTAACGCTCTTGGCGAAAAACTCGCCTTTGTCTTTCTGGAAATAGATTTTCAGAATATCgtt                                                      <  1:1994664/66‑1 (MQ=255)
               aaaCTTAACGCTCTTGGCGAAAAACTCGCCTTTGTCTTTCTGGAAATAGATTTTCAGAATATCGttgttg                                                  <  1:6655840/70‑1 (MQ=255)
               aaaCTTAACGCTCTTGGCGAAAAACTCGCCTTTGTCTTTCTGGAAATAGATTTTCAGAATATCGTTGTTgg                                                 >  1:458432/1‑71 (MQ=255)
                aaCTTAACGCTCTTGGCGAAAAACTCGCCTTTGTCTTTCTGGAAATAGATTTTCAGAATATCGttgttg                                                  <  1:6647019/69‑1 (MQ=255)
                    tAACGCTCTTGGCGAAAAACTCGCCTTTGTCTTTCTGGAAATAGATTTTCAGAATATCgttgt                                                    >  1:3175211/1‑63 (MQ=255)
                    tAACGCTCTTGGCGAAAAACTCGCCTTTGTCTTTCTGGAAATAGATTTTCAGAATATCgtt                                                      >  1:2041827/1‑61 (MQ=255)
                       cGCTCTTGGCGAAAAACTCGCCTTTGTCTTTCTGGAAATAGATTTTCAGAATATCGTTGTTGGCTTCCTg                                          >  1:2868358/1‑70 (MQ=255)
                        gCTCTTGGCGAAAAACTCGCCTTTGTCTTTCTGGAAATAGATTTTCAGAATATCGTTGTTGGCTTCCTGcc                                        >  1:3608911/1‑71 (MQ=255)
                         ctctTGGCGAAAAACTCGCCTTTGTCTTTCTGGAAATATATTTTCAGAATATCGTTGTTGGCTTCCTGCCg                                       <  1:4273174/71‑1 (MQ=255)
                         ctctTGGCGAAAAACTCGCCTTTGTCTTTCTGGAAATATATTTTCAGAATATCGTTGTTGGCTTCCTGCCg                                       <  1:4282350/71‑1 (MQ=255)
                         ctctTGGCGAAAAACTCGCCTTTGTCTTTCTGGAAATAGATTTTCAGAATATCgttgt                                                    >  1:788276/1‑58 (MQ=255)
                          tctTGGCGAAAAACTCGCCTTTGTGTTTCTGGAAATAGATTTTCAGAATATCGTTGTTGGCtt                                              <  1:6239270/63‑1 (MQ=255)
                          tctTGGCGAAAAACTCGCCTTTGTCTTTCTGGAAATAGATTTTCAGAATCTCgttgtt                                                   <  1:3569293/58‑1 (MQ=255)
                          tctTGGCGAAAAACTCGCCTTTGTCTTTCTGGAAATAGATTTTCAGAATATCgttgtt                                                   <  1:1812157/58‑1 (MQ=255)
                          tctTGGCGAAAAACTCGCCTTTGTCTTTCTGGAAATAGATTTTCAGAATATCGTTGTTGGCtt                                              <  1:2008621/63‑1 (MQ=255)
                            ttGGCGAAAAACTCGCCTTTGTCTTTCTGGAAATAGATTTTCAGAATATCGTTGTTGGCTTCCTGCCgg                                      <  1:7557905/69‑1 (MQ=255)
                             tGGCGAAAAACTCGCCTTTGTCTTTCTGGAAATAGATTTTCAGAATATCGTTGTTGGCTTCCTGCCGGAgg                                   >  1:7206580/1‑71 (MQ=255)
                              ggCGAAAAACTCGCCTTTGTCTTTCTGGAAATAGATTTTCAGAATATCgtt                                                      >  1:3586056/1‑51 (MQ=255)
                              ggCGAAAAACTCGCCTTTGTCTTTCTGGAAATAGATTTTCAGAATATCGTTGTTGGCTTCCTGCCGGAgg                                   <  1:246771/70‑1 (MQ=255)
                              ggCGAAAAACTCGCCTTTGTCTTTCTGGAAATAGATTTTCAGAATATCGTTGTTGGCTTCCTGCCGGAgg                                   <  1:5074434/70‑1 (MQ=255)
                              ggCGAAAAACTCGCCTTTGTCTTTCTGGAAATAGATTTTCAGAATATCGTTGTTGGCTTCCTGCCGGAgg                                   <  1:6511276/70‑1 (MQ=255)
                               gCGAAAAACTCGCCTTTGTCTTTCTGGAAATAGATTTTCAGAATATCg                                                        >  1:6653267/1‑48 (MQ=255)
                               gCGAAAAACTCGCCTTTGTCTTTCTGGAAATAGATTTTCAGAATATCGTTGTTGGCt                                               <  1:3308386/57‑1 (MQ=255)
                               gCGAAAAACTCGCCTTTGTCTTTCTGGAAATAGATTTTCAGAATATCGTTGTTGGCt                                               <  1:569425/57‑1 (MQ=255)
                               gCGAAAAACTCGCCTTTGTCTTTCTGGAAATAGATTTTCAGAATATCGTTGTTGGCt                                               <  1:5192047/57‑1 (MQ=255)
                               gCGAAAAACTCGCCTTTGTCTTTCTGGAAATAGATTTTCAGAATATCGTTGTTGGCt                                               <  1:2535639/57‑1 (MQ=255)
                               gCGAAAAACTCGCCTTTGTCTTTCTGGAAATAGATTTTAAGAatat                                                          >  1:5522937/1‑46 (MQ=255)
                                cGAAAAACTCGCCTTTGTCTTTCTTGAAATAGATTTTCAGAATATCGTTGTTGGCTTc                                             <  1:7844670/58‑1 (MQ=255)
                                   aaaaCTCGCCTTTGTCTTTCTGGAAATAGATTTTCAGAATATCGTTGTTGGCTTCCTGCCGGa                                     >  1:1297845/1‑63 (MQ=255)
                                   aaaaCTCGCCTTTGTCTTTCTGGAAATAGATTTTCAGAATATCGTTGTTGGCTTCCTGCCGGAGGCTGTaa                             >  1:7202448/1‑71 (MQ=255)
                                   aaaaCTCGCCTTTGTCTTTCTGGAAATAGATTTTCAGAATATCGTTGTTGGCTTCCTGCCGGAGGCTGTaa                             >  1:1725177/1‑71 (MQ=255)
                                    aaaCTCGCCTTTGTCTTTCTGGAAATAGATTTTCAGAat                                                            >  1:2605201/1‑39 (MQ=255)
                                    aaaCTCGCCTTTGTCTTTCTGGAAATAGATTTTCAGAATATCGTTGTTGGCTTCCTGCCGGAGGCTGTaa                             >  1:180560/1‑70 (MQ=255)
                                    aaaCTCGCCTTTGTCTTTCTGGAAATAGATTTTCAGAATATCGTTGTTGGCTTCCTGCCGGAGGCTGTAAc                            <  1:4191521/71‑1 (MQ=255)
                                     aaCTCGCCTTTGTCTTTCTGGAAATAGATTTTCAGAatat                                                          <  1:7419118/40‑1 (MQ=255)
                                     aaCTCGCCTTTGTCTTTCTGGAAATAGATTTTCAGAATATCGTTGTTGGCTTCCTGCCGGAGGCTGTAAc                            >  1:2772031/1‑70 (MQ=255)
                                         cGCCTTTGTCTTTCTGGAAATAGATTTTCAGAATATCGTTGTTGGCTTCCTGCCGGAGGCTGTa                              >  1:2526412/1‑64 (MQ=255)
                                          gCCTTTGTCTTTCTGGAAATAGATTTTCAGAATATCgttgtt                                                   >  1:6495482/1‑42 (MQ=255)
                                          gCCTTTGTCTTTCTGGAAATAGATTTTCAGAATATCGTTGTTGGCTTCCTGCCGGAGGCTGTAACGATc                        <  1:1980062/69‑1 (MQ=255)
                                           ccTTTGTCTTTCTGGAAATAGATTTTCAGAATATCg                                                        <  1:4965971/36‑1 (MQ=255)
                                           ccTTTGTCTTTCTGGAAATAGATTTTCAGAATATCGTTGTTGGCTTCCTGCCGGAGGCTGTAAc                            <  1:4224566/64‑1 (MQ=255)
                                            cTTTGTCTTTCTGGAAATAGATTTTCAGAATATCGTTGTTGGCTTCCTGCCGGAGGCTGTAACGATCAAtt                    <  1:7835443/71‑1 (MQ=255)
                                            cTTTGTCTTTCTGGAAATAGATTTTCAGAATATCGTTGTTGGCTTCCTGCCGGAGGCTGTAACGATCAAtt                    <  1:431469/71‑1 (MQ=255)
                                            cTTTGTCTTTCTGGAAATAGATTTTCAGAATATCGTTGTTGGCTTCCTGCCGGAGGCTGTAACGATCAAtt                    <  1:4029184/71‑1 (MQ=255)
                                            cTTTGTCTTTCTGGAAATAGATTTTCAGAATATCGTTGTTGGCTTCCTGCCGGAGGCTGTAACGATCAAtt                    <  1:1739079/71‑1 (MQ=255)
                                            cTTTGTCTTTCTGGAAATAGATTTTCAGAATATCGTTGTTGGCTTCCTGCCGGAGGCTGTAACGATCAAtt                    <  1:7580792/71‑1 (MQ=255)
                                            cTTTGTCTTTCTGGAAATAGATTTTCAGAATATCGTTGTTGGCTTCCTGCCGGAGGCTGTAACGATCAAt                     <  1:2083558/70‑1 (MQ=255)
                                             tttGTCTTTCTGGAAATAGATTTTCAGAATATCGTTGTTGGCTTcc                                            >  1:4489160/1‑46 (MQ=255)
                                             tttGTCTTTCTGGAAATAGATTTTCAGAATATCGTTGTTGGCTTcc                                            >  1:3396349/1‑46 (MQ=255)
                                              ttGTCTTTCTGGAAATAGATTTTCAGAATATCGTTGTTGGCTTCCTGCCGGAGGCTGTAACGATCAATttc                  >  1:5429720/1‑71 (MQ=255)
                                              ttGTCTTTCTGGAAATAGATTTTCAGAATATCGTTGTTGGCTTCCTGCCGGAGGCTGTAACGATCAATttc                  >  1:3172639/1‑71 (MQ=255)
                                               tGTCTTTCTGGAAATAGATTTTCAGAATATCgttg                                                     >  1:5646374/1‑35 (MQ=255)
                                               tGTCTTTCTGGAAATAGATTTTCAGAATATCGTTGTTGGCTTCCt                                           >  1:1694662/1‑45 (MQ=255)
                                               tGTCTTTCTGGAAATAGATTTTCAGAATATCGTTGTTGGCTTCCt                                           >  1:4276149/1‑45 (MQ=255)
                                               tGTCTTTCTGGAAATAGATTTTCAGAATATCGTTGTTGGCTTCCt                                           >  1:3355567/1‑45 (MQ=255)
                                               tGTCTTTCTGGAAATAGATTTTCAGAATATCGTTGTTGGCTTCCTGCCGGAgg                                   >  1:3732304/1‑53 (MQ=255)
                                                  tattCTGAAAATCTATTTCCAGAATATCGTTGTTGGCTTCCTGCCGGAGGCTGTAACGATCAATttctt                >  1:648375/3‑69 (MQ=37)
                                                   tttCTGGAAATAGATTTTCAGAATATCGTTGTTGGCtt                                              <  1:6636238/38‑1 (MQ=255)
                                                     tCTGGAAATAGATTTTCAGAATATCGTTGTTGGCTTCCTGCCGGAGGCTGTAACGATCAATTTCTTCAgg            <  1:7879826/70‑1 (MQ=255)
                                                     tCTGGAAATAGATTTTCAGAATATCGTTGTTGGCTTCCTGCCGGAGGCTGTAACGATCAATTTCTTCAGGa           >  1:3950032/1‑71 (MQ=255)
                                                     tCTGAAAATCTATTTCCAGAATATCGTTGTTGGCTTCCTGCCGGAGGCTGTAACGATCAATTTCTTCAGGa           >  1:647666/1‑71 (MQ=38)
                                                        ggAAATAGATTTTCAGAATATCGTTGTTGGCTTCCTGCCGGAGGCTGTAACGATCAATTTCTTCAGGATTg        >  1:3833548/1‑71 (MQ=255)
                                                         gAAATAGATTTTCAGAATATCGTTGTTGGCTTCCTGCCGGAGGCTGTAACGATCAATTTCTTCAGGAt          <  1:5735807/68‑1 (MQ=255)
                                                         gAAATAGATTTTCAGAATATCGTTGTTGGCTTCCTGCCGGAGGCTGTAACGATCAATTTCTTCAGGAt          <  1:3300254/68‑1 (MQ=255)
                                                         gAAATAGATTTTCAGAATATCGTTGTTGGCTTCCTGCCGGAGGCTGTAACGATCAATTTCTTCAGGATTgg       <  1:3799796/71‑1 (MQ=255)
                                                         gAAATAGATTTTCAGAATATCGTTGTTGGCTTCCTGCCGGAGGCTGTAACGATCAATTTCTTCAGGATTgg       <  1:1733249/71‑1 (MQ=255)
                                                          aaaTAGATTTTCAGAATATCGTTGTTGGCTTCCTGCCGGAGGCTGTAACGATCAATTTCTTCAGGATTGGt      <  1:3795247/71‑1 (MQ=255)
                                                            aTAGATTTTCAGAATATCGTTGTTGGCTTCCTGCCGGAGGCTGTAACGATCAATTTCTTCAGGATTGGTaa    <  1:1951533/71‑1 (MQ=255)
                                                            aTAGATTTTCAGAATATCGTTGTTGGCTTCCTGCCGGAGGCTGTAACGATCAATTTCTTCAGGATTGGTaa    <  1:3114838/71‑1 (MQ=255)
                                                            aTAGATTTTCAGAATATCGTTGTTGGCTTCCTGCCGGAGGCTGTAACGATCAATTTCTTCAGGATTGGTaa    <  1:90862/71‑1 (MQ=255)
                                                             tAGATTTTCAGAATATCGTTGTTGGCTTCCTGCCGGAg                                    >  1:6241821/1‑38 (MQ=255)
                                                             tAGATTTTCAGAATATCGTTGTTGGCTTCCTGCCGGAGGCTGTAACGATCAATttctt                >  1:545033/1‑58 (MQ=255)
                                                             tAGATTTTCAGAATATCGTTGTTGGCTTCCTGCCGGAGGCTGTAACGATCAATTTCTTCa              >  1:2321477/1‑60 (MQ=255)
                                                             tAGATTTTCAGAATATCGTTGTTGGCTTCCTGCCGGAGGCTGTAACGATCAATTTCTTCa              >  1:6803469/1‑60 (MQ=255)
                                                             tAGATTTTCAGAATATCGTTGTTGGCTTCCTGCCGGAGGCTGTAACGATCAATTTCTTCAGGATTGGTAAt   <  1:4370124/71‑1 (MQ=255)
                                                              aGATTTTCAGAATATCGTTGTTGGCTTCCTGCCgg                                      >  1:2299264/1‑35 (MQ=255)
                                                              aGATTTTCAGAATATCGTTGTTGGCTTCCTGCCgg                                      >  1:4548747/1‑35 (MQ=255)
                                                              aGATTTTCAGAATATCGTTGTTGGCTTCCTGCCGGAGGCTGTaa                             >  1:3446807/1‑44 (MQ=255)
                                                               gATTTTCAGAATATCGTTGTTGGCTTCCTGCCGTAGGCTGTAACGATCAATTTCTTCAGGATTGGTAATa  <  1:7479612/70‑1 (MQ=255)
                                                               gATTTTCAGAATATCGTTGTTGGCTTCCTGCCGGAGGCTGTAACGATCAATTTCTTCAGGAtt         >  1:3587770/1‑63 (MQ=255)
                                                               gATTTTCAGAATATCGTTGTTGGCTTCCTGCCGGAGGCTGTAACGATCAATTTCTTCAGGATTGGTAATa  <  1:4463700/70‑1 (MQ=255)
                                                               |                                                                     
CTGACGCGGATATTTAAACTTAACGCTCTTGGCGAAAAACTCGCCTTTGTCTTTCTGGAAATAGATTTTCAGAATATCGTTGTTGGCTTCCTGCCGGAGGCTGTAACGATCAATTTCTTCAGGATTGGTAATA  >  AP009048/183937‑184069

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