Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 4,285,507 T→C 9.7% S387G (AGT→GGT) ‡ yjcF ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090484,285,5070TC 9.7% 338.7 / 23.6 123S387G (AGT→GGT) ‡yjcFconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (54/57);  new base C (7/5);  total (61/62)
Fisher's exact test for biased strand distribution p-value = 5.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

Only 100 of 124 total aligned reads displayed.

TTGTGACATGCCAGTATGAGTGGCACAATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCTGGACGCTAGCCTCATCGA  >  AP009048/4285439‑4285575
                                                                    |                                                                    
ttGTGACATGCCAGTATGAGTGGCACAATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCACTAt                                                                    >  1:2815005/1‑71 (MQ=255)
 tgtgACATGCCAGTATGAGTGGCACAATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCACTAt                                                                    <  1:7622797/70‑1 (MQ=255)
     aCATGCCAGTATGAGTGGCACAATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATg                                                               >  1:56096/1‑71 (MQ=255)
        tGCCAGTATGAGTGGCACAATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATg                                                               >  1:6509919/1‑68 (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:3179428/1‑71 (MQ=255)
           cAGTATGAGTGGCACAATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGata                                                         >  1:2329128/1‑71 (MQ=255)
           cAGTATGAGTGGCACAATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGata                                                         <  1:4351503/71‑1 (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)
              tATGAGTGGCACAATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCTCCATTAAtagt                                                             >  1:5699654/1‑61 (MQ=255)
              tATGAGTGGCACAATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCTCCATTAAtagt                                                             >  1:5013465/1‑61 (MQ=255)
              tATGAGTGGCACAATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCACTAtt                                                                   <  1:7141638/58‑1 (MQ=255)
              tATGAGTGGCACAATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCACTAtt                                                                   <  1:7114475/58‑1 (MQ=255)
              tATGAGTGGCACAATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCACTAtt                                                                   <  1:7574800/58‑1 (MQ=255)
              tATGAGTGGCACAATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATg                                                               <  1:1350562/62‑1 (MQ=255)
              tATGAGTGGCACAATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGATAttt                                                      <  1:1122661/71‑1 (MQ=255)
              tATGAGTGGCACAATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGATAttt                                                      <  1:6369810/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:6294341/71‑1 (MQ=255)
                tGAGTGGCACAATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGATATTTTc                                                    >  1:6953685/1‑71 (MQ=255)
                tGAGTGGCACAATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGATATTTTc                                                    >  1:1067910/1‑71 (MQ=255)
                  aGTGGCACAATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCCCTATTaa                                                                 <  1:3440635/56‑1 (MQ=255)
                  aGTGGCACAATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTaa                                                                 <  1:1057043/56‑1 (MQ=255)
                  aGTGGCACAATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTaa                                                                 <  1:4226715/56‑1 (MQ=255)
                    tGGCACAATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATg                                                               >  1:2659907/1‑56 (MQ=255)
                    tGGCACAATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATg                                                               >  1:277535/1‑56 (MQ=255)
                      gCACAATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTaa                                                                 >  1:3481194/1‑52 (MQ=255)
                      gCACAATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGATATTTTCCTcac                                               >  1:1262516/1‑70 (MQ=255)
                       cacaATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCTCCATTAAtag                                                              >  1:7150014/1‑52 (MQ=38)
                       cacaATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCTCCATTAAtag                                                              >  1:6413945/1‑52 (MQ=38)
                         caATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGATATTTTcc                                                   <  1:4564040/63‑1 (MQ=255)
                         caATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGATATTTTCCTCACCACg                                           >  1:2707330/1‑71 (MQ=255)
                         caATCTTCTTCAAAGAATTTCTGTAGTATTACCAGGATATCACTATTaa                                                                 >  1:5554179/1‑49 (MQ=255)
                           aTCTCCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGATATTTTCCTCACCACg                                           >  1:2061019/1‑69 (MQ=255)
                            tcttctTCAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGATATTTTCCTCACCACg                                           >  1:4234496/1‑68 (MQ=255)
                            tcttctTCAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGATATTTTCCTCACCACg                                           >  1:3860581/1‑68 (MQ=255)
                             cttcttCAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGat                                                          <  1:7301526/52‑1 (MQ=255)
                             cttcttCAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGat                                                          <  1:4494936/52‑1 (MQ=255)
                             cttcttCAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGATATTTTCCTCACCACGgag                                        >  1:5842440/1‑70 (MQ=255)
                              ttcttcAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGaa                                      >  1:1674843/1‑71 (MQ=255)
                                   cAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATg                                    <  1:6999719/68‑1 (MQ=255)
                                    aaaGAATTTCTGCAGTGTTACCAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTaa                                >  1:7070773/1‑71 (MQ=255)
                                     aaGAATTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATg                                    >  1:6624332/1‑66 (MQ=255)
                                      aGAATTTCTGTAGTGTTACCAGGATATCTCCATTAATAGTGatat                                                        <  1:7391900/45‑1 (MQ=18)
                                      aGAATTTCTGTAGTGTTACCAGGATATCTCCATTAATAGTGatat                                                        <  1:7921279/45‑1 (MQ=18)
                                        aaTTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATAcc                            <  1:668430/71‑1 (MQ=255)
                                         aTTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATAcc                            >  1:7142870/1‑70 (MQ=255)
                                               gTAGTGTTACCAGGATATCACTATTAATGGAGATATTTTCCTcac                                               <  1:194647/45‑1 (MQ=255)
                                               gTAGTGTTACCAGGATATCACTATTAATGGAGATATTTTCCTcac                                               <  1:5866781/45‑1 (MQ=255)
                                                tAGTGTTACCAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCt                    >  1:2195630/1‑71 (MQ=255)
                                                tAGTGTTACCAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCt                    >  1:4103437/1‑71 (MQ=255)
                                                  gtgtTACCAGGATATCACTATTAATGGAGGTATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCTg                   <  1:3564977/70‑1 (MQ=255)
                                                  gtgtTACCAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCTgg                  <  1:6290043/71‑1 (MQ=255)
                                                  gtgtTACCAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCTg                   >  1:4894910/1‑70 (MQ=255)
                                                  gtgtTACCAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCTg                   <  1:6889764/70‑1 (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)
                                                        ccAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTa                                 <  1:489698/50‑1 (MQ=255)
                                                        ccAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCTGGACGCTa            <  1:6857110/71‑1 (MQ=255)
                                                         cAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCTGGACGCTAg           <  1:3549779/71‑1 (MQ=255)
                                                         cAGGATATCACTATTAATGGAGATATTTTCCTAACCACGGAGAATGGTAATACCCATTCCCt                    >  1:3672195/1‑62 (MQ=255)
                                                          aGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCTGGAc                <  1:6919827/65‑1 (MQ=255)
                                                          aGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCTGGAc                <  1:3150919/65‑1 (MQ=255)
                                                          aGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACACATTCCCTGGAc                <  1:7492872/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)
                                                             atatCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCATGGACGCTAGCCTc       <  1:7168078/71‑1 (MQ=255)
                                                              tatCACTATTAATGGAGATATTTTCCTCACCACGGAGAATgg                                   <  1:5488762/42‑1 (MQ=255)
                                                              tatCACTATTAATGGAGATATTTTCCTCACCACGGAGAATgg                                   <  1:4317496/42‑1 (MQ=255)
                                                              tatCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCa                          >  1:5946545/1‑51 (MQ=255)
                                                               atCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGt                                  <  1:1426690/42‑1 (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)
                                                                  aCTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCTGGACGCTAGCCTCATCGa  <  1:3010553/71‑1 (MQ=255)
                                                                   cTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCTGGACGCTAGCCTCATc    >  1:3450363/1‑68 (MQ=255)
                                                                   cTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCTGGACGCTAGCCTCATc    >  1:401724/1‑68 (MQ=255)
                                                                   cTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCTGGACGCTAGCCTCATc    >  1:4287422/1‑68 (MQ=255)
                                                                   cTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCTGGACGCTAGCCTCATc    >  1:1128049/1‑68 (MQ=255)
                                                                   cTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCTGGACGCTAGCCTCATCg   >  1:3760525/1‑69 (MQ=255)
                                                                   cTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCTGGACGCTAGCCTCATCGa  >  1:495084/1‑70 (MQ=255)
                                                                   cTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCTGGACGCTAGCCTCATCGa  <  1:1122230/70‑1 (MQ=255)
                                                                    |                                                                    
TTGTGACATGCCAGTATGAGTGGCACAATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCTGGACGCTAGCCTCATCGA  >  AP009048/4285439‑4285575

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: