Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 4,285,505 A→T 10.5% S387R (AGT→AGA) ‡ yjcF ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090484,285,5050AT10.5% 309.2 / 25.2 116S387R (AGT→AGA) ‡yjcFconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (47/56);  new base T (7/5);  total (54/62)
Fisher's exact test for biased strand distribution p-value = 5.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.17e-01

Only 100 of 119 total aligned reads displayed.

AGTTTGTGACATGCCAGTATGAGTGGCACAATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCTGGACGCTAGCCTCATCGA  >  AP009048/4285436‑4285575
                                                                     |                                                                      
aGTTTGTGACATGCCAGTATGAGTGGCACAATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCAc                                                                       >  1:14519/1‑71 (MQ=255)
 gTTTGTGACATGCCAGTATGAGTGGCACAATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCAc                                                                       <  1:4349142/70‑1 (MQ=255)
   ttGTGACATGCCAGTATGAGTGGCACAATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCACTAt                                                                    >  1:2815005/1‑71 (MQ=255)
        aCATGCCAGTATGAGTGGCACAATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATg                                                               >  1:56096/1‑71 (MQ=255)
           tGCCAGTATGAGTGGCACAATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATGgag                                                            <  1:2746764/71‑1 (MQ=255)
           tGCCAGTATGAGTGGCACAATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATGgag                                                            <  1:6681249/71‑1 (MQ=255)
             ccAGTATGAGTGGCACAATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCTCCATTaa                                                                 <  1:6169822/64‑1 (MQ=255)
              cAGTATGAGTGGCACAATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGata                                                         >  1:6558834/1‑71 (MQ=255)
              cAGTATGAGTGGCACAATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGata                                                         >  1:2329128/1‑71 (MQ=255)
              cAGTATGAGTGGCACAATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGata                                                         <  1:4351503/71‑1 (MQ=255)
              cAGTATGAGTGGCACAATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGata                                                         >  1:3179428/1‑71 (MQ=255)
              cAGTATGAGTGGCACAATCTTCTTCAAAGAATTTCTCTAGTGTTACCAGGATATCACTATTAATGGAGata                                                         >  1:7610370/1‑71 (MQ=255)
               aGTATGAGTTGCACAATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCTCCATTAAta                                                               <  1:5661808/64‑2 (MQ=255)
               aGTATGAGTGGCACAATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCTCCATTAAta                                                               <  1:4331607/64‑2 (MQ=255)
               aGTATGAGTGGCACAATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATg                                                               >  1:606329/1‑64 (MQ=255)
                gTATGAGTGGCACAATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCAc                                                                       <  1:7609264/55‑1 (MQ=255)
                 tATGAGTGGCACAATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCTCCATTAAtagt                                                             >  1:5699654/1‑61 (MQ=255)
                 tATGAGTGGCACAATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCTCCATTAAtagt                                                             >  1:5013465/1‑61 (MQ=255)
                 tATGAGTGGCACAATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCTCCATTAAtagt                                                             >  1:5827131/1‑61 (MQ=255)
                 tATGAGTGGCACAATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCACTAtt                                                                   <  1:7114475/58‑1 (MQ=255)
                 tATGAGTGGCACAATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCACTAtt                                                                   <  1:7574800/58‑1 (MQ=255)
                 tATGAGTGGCACAATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCACTAtt                                                                   <  1:7141638/58‑1 (MQ=255)
                 tATGAGTGGCACAATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATg                                                               <  1:1350562/62‑1 (MQ=255)
                 tATGAGTGGCACAATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATGGATATAttt                                                      <  1:16043/71‑1 (MQ=255)
                 tATGAGTGGCACAATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGATAttt                                                      <  1:6369810/71‑1 (MQ=255)
                 tATGAGTGGCACAATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGATAttt                                                      <  1:4921141/71‑1 (MQ=255)
                 tATGAGTGGCACAATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGATAttt                                                      <  1:1122661/71‑1 (MQ=255)
                  aTGAGTGGCACAATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGATAtttt                                                     >  1:6561846/1‑71 (MQ=255)
                   tGAGTGGCACAATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCTCCATTAATAGTGatatcct                                                     >  1:1599006/1‑67 (MQ=25)
                   tGAGTGGCACAATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGATATTTTc                                                    >  1:7109179/1‑71 (MQ=255)
                   tGAGTGGCACAATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGATATTTTc                                                    >  1:1067910/1‑71 (MQ=255)
                   tGAGTGGCACAATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGATATTTTc                                                    >  1:6953685/1‑71 (MQ=255)
                   tGAGTGGCACAATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGATATTTTc                                                    <  1:6294341/71‑1 (MQ=255)
                     aGTGGCACAATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCCCTATTaa                                                                 <  1:3440635/56‑1 (MQ=255)
                     aGTGGCACAATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTaa                                                                 <  1:4226715/56‑1 (MQ=255)
                     aGTGGCACAATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTaa                                                                 <  1:1057043/56‑1 (MQ=255)
                       tGGCACAATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCACTAtt                                                                   <  1:363986/52‑1 (MQ=255)
                       tGGCACAATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATg                                                               >  1:277535/1‑56 (MQ=255)
                       tGGCACAATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGATATTTTcc                                                   <  1:2837621/68‑1 (MQ=255)
                         gCACAATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGATATTTTCCTcac                                               >  1:1262516/1‑70 (MQ=255)
                          cacaATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCTCCATTAAtag                                                              >  1:6413945/1‑52 (MQ=38)
                           acaATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGATATTTTCCTcaccac                                            <  1:4236146/71‑1 (MQ=255)
                           acaATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGATATTTTCCTcaccac                                            <  1:7693754/71‑1 (MQ=255)
                            caATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGATATTTTCCTCACCACg                                           >  1:2707330/1‑71 (MQ=255)
                            caATCTTCTTCAAAGAATTTCTGTAGTATTACCAGGATATCACTATTaa                                                                 >  1:5554179/1‑49 (MQ=255)
                             aaTCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGATATTTTc                                                    >  1:718016/1‑61 (MQ=255)
                              aTCTCCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGATATTTTCCTCACCACg                                           >  1:2061019/1‑69 (MQ=255)
                               tcttctTCAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGATATTTTCCTCACCACg                                           >  1:4234496/1‑68 (MQ=255)
                                cttcttCAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGat                                                          <  1:4494936/52‑1 (MQ=255)
                                cttcttCAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGat                                                          <  1:7301526/52‑1 (MQ=255)
                                 ttcttcAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGaa                                      >  1:4829980/1‑71 (MQ=255)
                                 ttcttcAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGaa                                      >  1:1674843/1‑71 (MQ=255)
                                     tcAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGATAtt                                                       <  1:3916003/50‑1 (MQ=255)
                                      cAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATg                                    <  1:1214963/68‑1 (MQ=255)
                                       aaaGAATTTCTGCAGTGTTACCAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTaa                                >  1:7070773/1‑71 (MQ=255)
                                        aaGAATTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATg                                    >  1:6624332/1‑66 (MQ=255)
                                         aGAATTTCTGTAGTGTTACCAGGATATCTCCATTAATAGTGatat                                                        <  1:7921279/45‑1 (MQ=18)
                                         aGAATTTCTGTAGTGTTACCAGGATATCTCCATTAATAGTGatat                                                        <  1:7391900/45‑1 (MQ=18)
                                           aaTTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATAcc                            <  1:668430/71‑1 (MQ=255)
                                            aTTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATAcc                            >  1:7142870/1‑70 (MQ=255)
                                                cTGTAGTGTTACCAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTc                       >  1:4352828/1‑71 (MQ=255)
                                                  gTAGTGTTACCAGGATATCACTATTAATGGAGATATTTTCCTcac                                               <  1:5866781/45‑1 (MQ=255)
                                                   tAGTGTTACCAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCt                    >  1:4103437/1‑71 (MQ=255)
                                                   tAGTGTTACCAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCt                    >  1:2195630/1‑71 (MQ=255)
                                                    aGTGTTACCAGGATATCACTATTAATGGAGATATTTTCCTcacca                                             >  1:2681847/1‑45 (MQ=255)
                                                    aGTGTTACCAGGATATCACTATTAATGGAGATATTTTCCTcacca                                             >  1:7135051/1‑45 (MQ=255)
                                                     gtgtTACCAGGATATCACTATTAATGGAGGTATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCTg                   <  1:3564977/70‑1 (MQ=255)
                                                     gtgtTACCAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCTg                   <  1:6889764/70‑1 (MQ=255)
                                                     gtgtTACCAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCTg                   >  1:4894910/1‑70 (MQ=255)
                                                      ttttACCAGGATATCACTATTAATGGAGATATTTTCCTcacca                                             <  1:6379262/41‑1 (MQ=255)
                                                       gtTACCAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCTGGAc                >  1:1867847/1‑71 (MQ=255)
                                                        ttACCAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCa                          <  1:868079/60‑1 (MQ=255)
                                                        ttACCAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCa                          <  1:4532046/60‑1 (MQ=255)
                                                        ttACCAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCa                          <  1:4785981/60‑1 (MQ=255)
                                                        ttACCAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCTGGACg               <  1:6757487/71‑1 (MQ=255)
                                                         tACCAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCTGGACGc              >  1:5844083/1‑71 (MQ=255)
                                                          aCCAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATAccc                           <  1:3851622/57‑1 (MQ=255)
                                                          aCCAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCTg                   >  1:7911982/1‑65 (MQ=255)
                                                           ccAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCTGGACGCTa            <  1:6361016/71‑1 (MQ=255)
                                                           ccAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCTGGACGCTa            <  1:6857110/71‑1 (MQ=255)
                                                            cAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCTGGACGCTAg           <  1:7517571/71‑1 (MQ=255)
                                                            cAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCTGGACGCTAg           <  1:3549779/71‑1 (MQ=255)
                                                            cAGGATATCACTATTAATGGAGATATTTTCCTAACCACGGAGAATGGTAATACCCATTCCCt                    >  1:3672195/1‑62 (MQ=255)
                                                             aGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGaa                                      >  1:3384170/1‑43 (MQ=255)
                                                             aGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCTGGAc                <  1:3150919/65‑1 (MQ=255)
                                                             aGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCTGGAc                <  1:6919827/65‑1 (MQ=255)
                                                             aGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACACATTCCCTGGACGCt             >  1:485478/1‑68 (MQ=255)
                                                               gATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCTGGACGCTAGCCt        <  1:7715919/71‑1 (MQ=255)
                                                                atatCTCCATTAATAGTGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCTGGACGCTAGCCt        >  1:7620196/1‑70 (MQ=38)
                                                                atatCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCTGGACGCTAGCCt        >  1:5434930/1‑70 (MQ=255)
                                                                atatCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCTGGACGCTAGCCTc       <  1:6000118/71‑1 (MQ=255)
                                                                atatCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCATGGACGCTAGCCTc       <  1:7168078/71‑1 (MQ=255)
                                                                 tatCACTATTAATGGAGATATTTTCCTCACCACGGAGAATgg                                   <  1:5488762/42‑1 (MQ=255)
                                                                 tatCACTATTAATGGAGATATTTTCCTCACCACGGAGAATgg                                   <  1:2183154/42‑1 (MQ=255)
                                                                 tatCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCa                          >  1:5946545/1‑51 (MQ=255)
                                                                  atCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTaa                                >  1:5476717/1‑44 (MQ=255)
                                                                   tCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTa                                 >  1:3342691/1‑42 (MQ=255)
                                                                    cACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATa                              >  1:4875473/1‑44 (MQ=255)
                                                                     actactAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCTGGACGCTAGCCTCATCGa  <  1:747283/71‑1 (MQ=255)
                                                                     aCTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCTGGACGCTAGCCTCATCg   >  1:2592368/1‑70 (MQ=255)
                                                                     |                                                                      
AGTTTGTGACATGCCAGTATGAGTGGCACAATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCTGGACGCTAGCCTCATCGA  >  AP009048/4285436‑4285575

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: