Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 4,027,857 A→G 20.0% L232L (TTG→CTG)  yhhT ← predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090484,027,8570AG20.0% 59.6 / 11.5 30L232L (TTG→CTG) yhhTpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (10/14);  new base G (4/2);  total (14/16)
Fisher's exact test for biased strand distribution p-value = 3.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

ACAGCAGCACCTGAATCATTGGCGGTACGGCGGAAATTACCGCGCCGATATTGGGCACGTAGTTGAGCAAAAACGCCAGTACTGCCCACATCAGCGCAAACTGCACCCCCATCAGCTCCAGCCCCAGCCAGACGATGACA  >  AP009048/4027788‑4027927
                                                                     |                                                                      
aCAGCAGCACCTGAATCATTGGCGGTACGGCGGAAATTACCGCGCCGATATTGGGCACGTAGTTGAGCaaa                                                                       <  1:697009/71‑1 (MQ=255)
aCAGCAGCACCTGAATCATTGGCGGTACGGCGGAAATTACCGCGCCGATATTGGGCACGTAGTTGAGCaaa                                                                       <  1:4236110/71‑1 (MQ=255)
aCAGCAGCACCTGAATCATTGGCGGTACGGCGGAAATTACCGCGCCGATATTGGGCACGTAGTTGAGCaaa                                                                       <  1:4200481/71‑1 (MQ=255)
aCAGCAGCACCTGAATCATTGGCGGTACGGCGGAAATTACCGCGCCGATATTGGGCACGTAGTTGAGCaa                                                                        >  1:1872572/1‑70 (MQ=255)
aCAGCAGCACCTGAATCATTGGCGGTACGGCGGAAATTACCGCGCCGATATTGGGCACGTAGTTGAGCaa                                                                        >  1:3259429/1‑70 (MQ=255)
   gcagcaCCTGAATCATTGGCGGTACGGCGGAAATTACCGCGCCGATATTGGGCACGTAGTTGAGCAAAAAc                                                                    <  1:305931/71‑1 (MQ=255)
     agcaCCTGAATCATTGGCGGTACGGCGGAAATTACCGCGCCGATATTGGGCACGTAGTTGAGCAAAAACGc                                                                  <  1:2683627/71‑1 (MQ=255)
     agcaCCTGAATCATTGGCGGTACGGCGGAAATTACCGCGCCGATATTGGGCACGTAGTTGAGCAAAAACGc                                                                  <  1:3217222/71‑1 (MQ=255)
         ccTGAATCATTGGCGGTACGGCGGAAATTACCGCGCCGATATTGGGCACGTAGTTGAGCAAAAACGCCAg                                                               >  1:1383434/1‑70 (MQ=255)
            gAATCATTGGCGGTACGGCGGAAATTACCGCGCCGATATTGGGCACGTAGTTGAGCAAAAACGCCAGTACt                                                           >  1:1724691/1‑71 (MQ=255)
              aTCATTGGCGGTACGGCGGAAATTACCGCGCCGATATTGGGCACGTAGTTGAGCAAAAACGCCAGt                                                              <  1:1951011/66‑1 (MQ=255)
                     gCGGTACGGCGGAAATTACCGCGCCGATATTGGGCACGTAGTTGAGCAAAAACGCCAGTACTGCCCACATc                                                  <  1:2155699/71‑1 (MQ=255)
                     gCGGTACGGCGGAAATTACCGCGCCGATATTGGGCACGTAGTTGAGCAAAAACGCCAGTACTGCCCACATc                                                  <  1:243602/71‑1 (MQ=255)
                           cggcggAAATTACCGCGCCGATATTGGGCACGTAGTTGAGCAAAAACGCCAg                                                               >  1:1413769/1‑52 (MQ=255)
                            ggcggAAATTACCGCGCCGATATTGGGCACGTAGTTGAGCAAAAACGCCAGTACTGCCCACATCAGCGCaa                                           <  1:2667752/71‑1 (MQ=255)
                             gcggAAATTACCGCGCCGATATTGGGCACGTAGTTGAGCaa                                                                        <  1:3809793/41‑1 (MQ=255)
                              cggAAATTACCGCGCCGATATTGGGCACGTAGTTGAGCAAAAACGCCAGTACTGCCCACATCAGCGCAAAc                                         >  1:2550605/1‑71 (MQ=255)
                                gAAATTACCGCGCCGATATTGGGCACGTAGTTGAGCAAAAACGCCa                                                                >  1:2003386/1‑46 (MQ=255)
                                 aaaTTACCGCGCCGATATTGGGCACGTAGTTGAGCAGAAACGCCAGTACTGCCCACATCAGCGCAAACTg                                       >  1:1017173/1‑70 (MQ=255)
                                   aTTACCGCGCCGATATTGGGCTCGTAGTTGAGCAAAAACGCCAGTACTGCCCACATCAGCGCAAACTGCAc                                    >  1:2332097/1‑71 (MQ=255)
                                   aTTACCGCGCCGATATTGGGCACGTAGTTGAGCAAAAACGCCAGTACTg                                                          >  1:3733914/1‑49 (MQ=255)
                                     tACCGCGCCGATATTGGGCACGTAGTTGAGCAAAAACGCCAGTACTGCCCACATCAGCGCAAACTGCAccc                                  >  1:979584/1‑71 (MQ=255)
                                      aCCGCGCCGATATTGGGCACGTAGTTGAGCAAAAACGCCAGTACTGCCCACATCAGCGCAAACTGCAccc                                  <  1:4542924/70‑1 (MQ=255)
                                       ccGCGCCGATATTGGGCACGTAGTTGAGCAGAAACGCCAGTACTGCCCACATCAGCGCAAACTGCAcccc                                 >  1:352459/1‑70 (MQ=255)
                                         gcgcCGATATTGGGCACGTAGTTGAGCAGAAACGCCAGTACTGCCCACATCAGCGCAAACTGCACCCCCAt                              <  1:1986081/71‑1 (MQ=255)
                                          cgcCGATATTGGGCACGTAGTTGAGCAAAAACGCCAGTACTGCCCACATCAGCGCAAACTGCACCCCCATc                             <  1:3942581/71‑1 (MQ=255)
                                            ccGATATTGGGCACGTAGTTGAGCAGAAACGCCAGTACTGCCCACATCAGCGCAAACTGCACCCCCATCAg                           >  1:3813558/1‑71 (MQ=255)
                                               atatTGGGCACGTAGTTGAGCAAAAACGCCAGTACTGCCCACATCAGTGCAAACTGCACCCCCATCAGCTc                        <  1:1175079/71‑1 (MQ=255)
                                                        aCGTAGTTGAGCAGAAACGCCAGTACTGCCCACATCAGCGCAAACTGCACCCCCATCAGCTccagccccag               >  1:2123163/1‑71 (MQ=255)
                                                        aCGTAGTTGAGCAAAAACGCCAGTACTGCCCACATCAGCGCAAACTGCACCCCCATCAGCTccagccccag               <  1:4431995/71‑1 (MQ=255)
                                                        aCGTAGTTGAGCAAAAACGCCAGTACTGCCCACATCAGCGCAAACTGCACCCCAATCAGCt                         >  1:501837/1‑61 (MQ=255)
                                                            aGTTGAGCAAAAACGCCAGTACTGCCCACATCAGCGCAAACTGCACCCCCATCa                            <  1:2823403/54‑1 (MQ=255)
                                                              atGAGCAGAAACGCCAGTACTGCCCACATCAGCGCAAACTGCACCCCCATCAGCTCCAGCCCCAGCCAGAc         <  1:1864637/70‑1 (MQ=255)
                                                                gAGCAAAAACGCCAGTACTGCCCACATCAGCGCAAACTGCACCCCCATCAGCTCCAGCCCCAGCCAGACg        <  1:1264898/70‑1 (MQ=255)
                                                                gAGCAAAAACGCCAGTACTGCCCACATCAGCGCAAACTGCACCCCCATCAGCTCCAGCCCCAGCCAGACg        <  1:2843928/70‑1 (MQ=255)
                                                                  gCAAAAACGCCAGTACTGCCCACATCAGCGCAAACTGCACCCCCATCAGCTCCAGCCCCAGCCAGACGAt      >  1:4365293/1‑70 (MQ=255)
                                                                    aaaaaCGCCAGTACTGCCCACATCAGCGCAAACTGCACCCCCATCAGCTCCAGCCCCAGCCAGACGATGac   <  1:205071/71‑1 (MQ=255)
                                                                    aaaaaCGCCAGTACTGCCCACATCAGCGCAAACTGCACCCCCATCAGCTCCAGCCCCAGCCAGACGATGa    >  1:2871073/1‑70 (MQ=255)
                                                                     aaaaCGCCAGTACTGCCCACATCAGCGCAAACTGCa                                     <  1:2458417/36‑1 (MQ=255)
                                                                     aaaaCGCCAGTACTGCCCACATCAGCGCAAACTGCACCCCCa                               >  1:1196597/1‑42 (MQ=255)
                                                                     aaaaCGCCAGTACTGCCCACATCAGCGCAAACTGCACCCCCATCa                            >  1:3056793/1‑45 (MQ=255)
                                                                     aaaaCGCCAGTACTGCCCACATCAGCGCAAACTGCACCCCCATCAGCTCCAGCCCCAGCCAGACGATGaca  >  1:2094150/1‑71 (MQ=255)
                                                                     |                                                                      
ACAGCAGCACCTGAATCATTGGCGGTACGGCGGAAATTACCGCGCCGATATTGGGCACGTAGTTGAGCAAAAACGCCAGTACTGCCCACATCAGCGCAAACTGCACCCCCATCAGCTCCAGCCCCAGCCAGACGATGACA  >  AP009048/4027788‑4027927

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: