Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 197,160 G→A 34.9% F147F (TTC→TTT)  proS ← prolyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE197,1600GA34.9% 93.4 / 70.3 123F147F (TTC→TTT) proSprolyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base G (45/35);  new base A (23/20);  total (68/55)
Fisher's exact test for biased strand distribution p-value = 8.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.87e-01

Only 100 of 123 total aligned reads displayed.

GGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCGAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTGGATCTGATAGAAG  >  minE/197093‑197208
                                                                   |                                                
gggATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCGa                                                 <  1:1291818/69‑1 (MQ=255)
gggATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCGa                                                 <  1:1858119/69‑1 (MQ=255)
gggATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCGa                                                 <  1:415035/69‑1 (MQ=255)
gggATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCGa                                                 <  1:3250010/69‑1 (MQ=255)
gggATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCGa                                                 <  1:552317/69‑1 (MQ=255)
gggATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCGa                                                 <  1:2095374/69‑1 (MQ=255)
gggATTCCTGAGAAGTATGGAAAGAGTAAGCATCCTTCATCAGGAATTCGCGGGAACGCATGACGCCGa                                                 <  1:411849/69‑1 (MQ=255)
      ccTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCGAAAcgcg                                           <  1:2074442/69‑1 (MQ=255)
      ccTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCGAAAcgcg                                           <  1:1108103/69‑1 (MQ=255)
         gagaAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCGAAACGCGGa                                         <  1:2713494/68‑1 (MQ=255)
         gagaAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCGAAACGCGGa                                         <  1:508210/68‑1 (MQ=255)
         gagaAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCAAAAcgcg                                           >  1:2203045/1‑66 (MQ=255)
         gagaAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCAAAAcgcg                                           >  1:2862033/1‑66 (MQ=255)
         gagaAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCAAAAcgcg                                           >  1:2267739/1‑66 (MQ=255)
         gagaAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCAAAAcgcg                                           >  1:243879/1‑66 (MQ=255)
         gagaAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCAAAAcgcg                                           >  1:2979702/1‑66 (MQ=255)
         gagaAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCAAAAcgcg                                           >  1:1923964/1‑66 (MQ=255)
         gagaAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCAAAAcgcg                                           >  1:3211799/1‑66 (MQ=255)
         gagaAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCAAAAcgcg                                           >  1:2702674/1‑66 (MQ=255)
         gagaAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCAAAAcgcg                                           >  1:2623699/1‑66 (MQ=255)
         gagaAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCAAAAcgcg                                           >  1:2565762/1‑66 (MQ=255)
         gagaAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCAAAAcgcg                                           >  1:3383090/1‑66 (MQ=255)
         gagaAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCAAAAcgcg                                           >  1:956075/1‑66 (MQ=255)
         gagaAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCAAAAcgcg                                           >  1:830081/1‑66 (MQ=255)
         gagaAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCAAAAcgcg                                           >  1:791239/1‑66 (MQ=255)
         gagaAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCAAAAcgcg                                           >  1:1224352/1‑66 (MQ=255)
         gagaAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCAAAAcgcg                                           >  1:651657/1‑66 (MQ=255)
         gagaAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCAAAAcgcg                                           >  1:1377931/1‑66 (MQ=255)
         gagaAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCAAAAcgcg                                           >  1:140099/1‑66 (MQ=255)
         gagaAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCAAAAcgcg                                           >  1:1477446/1‑66 (MQ=255)
         gagaAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGAATGACGCCAAAAcgag                                           >  1:1305609/1‑64 (MQ=255)
            aaGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCGAAACGCGGACGc                                      <  1:2711037/68‑1 (MQ=255)
                        agTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCGAAACGCGGACGCACTTCGTCGCGGa                         <  1:1308501/69‑1 (MQ=255)
                        agTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCGAAACGCGGACGCACTTCGTCGCGGa                         <  1:316312/69‑1 (MQ=255)
                        agTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCGAAACGCGGACGCACTTCGTCGCGGa                         <  1:27128/69‑1 (MQ=255)
                            aGCATCTTTCATCACGAATTCGCGGGAACGCATGACGCCGAAACGCGGa                                         >  1:336122/1‑49 (MQ=255)
                              cATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCGAAACGCGGACGCACTTCGTCGCGGAACTTg                    <  1:2249905/68‑1 (MQ=255)
                              cATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCGAAACGCGGACGCACTTCGTCGCGGAACTTg                    <  1:839522/68‑1 (MQ=255)
                               aTCTTTCATCAGGAATTCGCGGGAACGCATGACGCCGAAACGCGGACGCACTTCGTCGCGGAACTTg                    <  1:5648/67‑1 (MQ=255)
                                tCTTTCATCAGGAATTCGCGGGAACGCATGACGCCGa                                                 >  1:2908930/1‑37 (MQ=255)
                                tCTTTCATCAGGAATTCGCGGGAACGCATGACGCCGa                                                 >  1:1224967/1‑37 (MQ=255)
                                  tttCATCAGGAATTCGCGGGAACGCATGACGCCGAAAcg                                             >  1:204512/1‑39 (MQ=255)
                                  tttCATCAGGAATTCGCGGGAACGCATGACGCCGAAACGCGGACGCACTTCGTCGCGGa                         >  1:2804359/1‑59 (MQ=255)
                                  tttCATCAGGAATTCGCGGGAACGCATGACGCCGAAACGCGGACGCACTTCGTCGCGGAACtt                     >  1:2618143/1‑63 (MQ=255)
                                  tttCATCAGGAATTCGCGGGAACGCATGACGCCGAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTg               >  1:1289277/1‑69 (MQ=255)
                                  tttCATCAGGAATTCGCGGGAACGCATGACGCCGAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTg               >  1:1138544/1‑69 (MQ=255)
                                  tttCATCAGGAATTCGCGGGAACGCATGACGCCGAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTg               >  1:1737913/1‑69 (MQ=255)
                                  tttCATCAGGAATTCGCGGGAACGCATGACGCCGAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTg               >  1:381483/1‑69 (MQ=255)
                                  tttCATCAGGAATTCGCGGGAACGCATGACGCCGAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTg               >  1:3070181/1‑69 (MQ=255)
                                  tttCATCAGGAATTCGCGGGAACGCATGACGCCGAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTg               >  1:1469980/1‑69 (MQ=255)
                                  tttCATCAGGAATTCGCGGGAACGCATGACGCCGAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTg               >  1:110119/1‑69 (MQ=255)
                                  tttCATCAGGAATTCGCGGGAACGCATGACGCCGAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTg               >  1:296192/1‑69 (MQ=255)
                                  tttCATCAGGAATTCGCGGGAACGCATGACGCCGAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTg               >  1:2134842/1‑69 (MQ=255)
                                  tttCATCAGGAATTCGCGGGAACGCATGACGCCGAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTg               >  1:2154390/1‑69 (MQ=255)
                                  tttCATCAGGAATTCGCGGGAACGCATGACGCCGAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTg               >  1:566041/1‑69 (MQ=255)
                                  tttCATCAGGAATTCGCGGGAACGCATGACGCCGAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTg               >  1:258045/1‑69 (MQ=255)
                                  tttCATCAGGAATTCGCGGGAACGCATGACGCCGAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTg               >  1:90640/1‑69 (MQ=255)
                                  tttCATCAGGAATTCGCGGGAACGCATGACGCCGAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTg               >  1:1017645/1‑69 (MQ=255)
                                      atcaGGAATTCGCGGGAACGCATGACGCCGAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTGGa             <  1:298373/67‑1 (MQ=255)
                                      atcaGGAATTCGCGGGAACGCATGACGCCGAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTGGa             <  1:716309/67‑1 (MQ=255)
                                      atcaGGAATTCGCGGGAACGCATGACGCCGAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTGGa             <  1:947296/67‑1 (MQ=255)
                                      atcaGGAATTCGCGGGAACGCATGACGCCGAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTGGa             <  1:3149962/67‑1 (MQ=255)
                                      atcaGGAATTCGCGGGAACGCATGACGCCGAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTGGa             <  1:666165/67‑1 (MQ=255)
                                      atcaGGAATTCGCGGGAACGCATGACGCCGAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTGGa             <  1:1229640/67‑1 (MQ=255)
                                      atcaGGAATTCGCGGGAACGCATGACGCCGAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTGGa             <  1:2367598/67‑1 (MQ=255)
                                      atcaGGAATTCGCGGGAACGCATGACGCCGAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTGGa             <  1:2174712/67‑1 (MQ=255)
                                           gAATTCGCGGGAACGCATGACGCCAAAACGCGGACGCACTTCGTCGCGTAACTTGGTCTg               <  1:913340/60‑1 (MQ=255)
                                           gAATTCGCGGGAACGCATGACGCCAAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTg               <  1:1449446/60‑1 (MQ=255)
                                           gAATTCGCGGGAACGCATGACGCCAAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTg               <  1:751353/60‑1 (MQ=255)
                                           gAATTCGCGGGAACGCATGACGCCAAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTg               <  1:687549/60‑1 (MQ=255)
                                           gAATTCGCGGGAACGCATGACGCCAAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTg               <  1:3071813/60‑1 (MQ=255)
                                           gAATTCGCGGGAACGCATGACGCCAAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTg               <  1:1639008/60‑1 (MQ=255)
                                           gAATTCGCGGGAACGCATGACGCCAAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTg               <  1:623063/60‑1 (MQ=255)
                                           gAATTCGCGGGAACGCATGACGCCAAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTg               <  1:2424315/60‑1 (MQ=255)
                                           gAATTCGCGGGAACGCATGACGCCAAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTg               <  1:3243094/60‑1 (MQ=255)
                                           gAATTCGCGGGAACGCATGACGCCAAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTg               <  1:2497229/60‑1 (MQ=255)
                                           gAATTCGCGGGAACGCATGACGCCAAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTg               <  1:2373402/60‑1 (MQ=255)
                                           gAATTCGCGGGAACGCATGACGCCAAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTg               <  1:1934265/60‑1 (MQ=255)
                                           gAATTCGCGGGAACGCATGACGCCAAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTg               <  1:2287673/60‑1 (MQ=255)
                                           gAATTCGCGGGAACGCATGACGCCAAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTg               <  1:2234087/60‑1 (MQ=255)
                                            acttCGCGAGAACGCATGACGCCAAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTg               <  1:2639955/57‑1 (MQ=255)
                                            aaTTCGCGGGAACGCATGACGCCAAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTg               <  1:3032149/59‑1 (MQ=255)
                                             aTTCGCGGGAACGCATGACGCCGAAACGCGGACGCACTTCGTCGCGGAACTTGGTc                 >  1:1042781/1‑56 (MQ=255)
                                             aTTCGCGGGAACGCATGACGCCGAAACGCGGACGCACTTCGTCGCGGAACTTGGTc                 >  1:570373/1‑56 (MQ=255)
                                             aTTCGCGGGAACGCATGACGCCGAAACGCGGACGCACTTCGTCGCGGAACTTGGTc                 >  1:2835140/1‑56 (MQ=255)
                                             aTTCGCGGGAACGCATGACGCCGAAACGCGGACGCACTTCGTCGCGGAACTTGGTc                 >  1:2981521/1‑56 (MQ=255)
                                             aTTCGCGGGAACGCATGACGCCGAAACGCGGACGCACTTCGTCGCGGAACTTGGTc                 >  1:3014822/1‑56 (MQ=255)
                                             aTTCGCGGGAACGCATGACGCCGAAACGCGGACGCACTTCGTCGCGGAACTTGGTc                 >  1:1193216/1‑56 (MQ=255)
                                             aTTCGCGGGAACGCATGACGCCGAAACGCGGACGCACTTCGTCGCGGAACTTGGTc                 >  1:3355583/1‑56 (MQ=255)
                                             aTTCGCGGGAACGCATGACGCCGAAACGCGGACGCACTTCGTCGCGGAACTTGGTc                 >  1:3285194/1‑56 (MQ=255)
                                             aTTCGCGGGAACGCATGACGCCGAAACGCGGACGCACGTCGTCGCGGAACTTGGTc                 >  1:3192664/1‑56 (MQ=255)
                                              ttCGCGGGAACGCATGACGCCGAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTGGa             <  1:1900703/59‑1 (MQ=255)
                                                cgcgGGAACGCATGACGCCGAAACGCGGACGCACTTCGTCGCGGaa                        >  1:2378213/1‑46 (MQ=255)
                                                cgcgGGAACGCATGACGCCGAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTGGATCTGATAGAAg  >  1:1528676/1‑68 (MQ=255)
                                                cgcgGGAACGCATGACGCCGAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTGGATCTGATAGAAg  >  1:1549310/1‑68 (MQ=255)
                                                cgcgGGAACGCATGACGCCGAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTGGATCTGATAGAAg  >  1:188490/1‑68 (MQ=255)
                                                cgcgGGAACGCATGACGCCGAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTGGATCTGATAGAAg  >  1:1945234/1‑68 (MQ=255)
                                                cgcgGGAACGCATGACGCCGAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTGGATCTGATAGAAg  >  1:2259677/1‑68 (MQ=255)
                                                cgcgGGAACGCATGACGCCGAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTGGATCTGATAGAAg  >  1:238985/1‑68 (MQ=255)
                                                cgcgGGAACGCATGACGCCGAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTGGATCTGATAGAAg  >  1:2436098/1‑68 (MQ=255)
                                                                   |                                                
GGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCGAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTGGATCTGATAGAAG  >  minE/197093‑197208

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: ‑