Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 1,218,230 T→G 6.8% intergenic (+18/‑95) ymgB → / → ymgC hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090481,218,2300TG 6.8% 229.2 / 10.5 88intergenic (+18/‑95)ymgB/ymgChypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (34/48);  new base G (3/3);  total (37/51)
Fisher's exact test for biased strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

AGAAAAACACTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTATTTTTGTCTGTCGCTTTAGACTATACACAGCATAATTTTATT  >  AP009048/1218165‑1218299
                                                                 |                                                                     
aGAAAAACACTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGtt                                                                  >  1:6312619/1‑71 (MQ=255)
aGAAAAACACTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGtt                                                                  >  1:4260427/1‑71 (MQ=255)
 gAAAAACACTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGtt                                                                  <  1:3456177/70‑1 (MQ=255)
 gAAAAACACTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGtt                                                                  <  1:6869511/70‑1 (MQ=255)
 gAAAAACACTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTATACTTAATAATTAGAAGt                                                                   <  1:5904878/69‑1 (MQ=25)
  aaaaaCACTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAagtttaagttta                                                                <  1:5905985/71‑1 (MQ=255)
    aaaCACTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTta                                                                >  1:7762676/1‑69 (MQ=255)
    aaaCACTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTta                                                                >  1:4043939/1‑69 (MQ=255)
    aaaCACTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTta                                                                >  1:5623980/1‑69 (MQ=255)
     aaCACTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAg                                                             <  1:1436108/71‑1 (MQ=255)
     aaCACTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAg                                                             <  1:72822/71‑1 (MQ=255)
      acacTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGa                                                            <  1:5537798/71‑1 (MQ=255)
       cacTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGttt                                                                 >  1:7819930/1‑65 (MQ=255)
       cacTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGaa                                                           <  1:227899/71‑1 (MQ=255)
       cacTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTATACTTAATAATTAGAAGttac                                                                >  1:5316384/1‑64 (MQ=21)
       cacTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTATACTTAATAATTAGAAGttac                                                                >  1:4536046/1‑64 (MQ=21)
          tGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGaaaa                                                         >  1:3762897/1‑70 (MQ=255)
          tGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAAc                                                        >  1:3684896/1‑71 (MQ=255)
          tGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAAc                                                        >  1:339713/1‑71 (MQ=255)
            gAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGtt                                                                  <  1:6151840/59‑1 (MQ=255)
            gAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACtc                                                      <  1:2550131/71‑1 (MQ=255)
              aaTCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTta                                                                <  1:2153894/59‑1 (MQ=255)
              aaTCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTta                                                                <  1:5647634/59‑1 (MQ=255)
              aaTCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTta                                                                <  1:2918116/59‑1 (MQ=255)
              aaTCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACtcat                                                    <  1:4396429/71‑1 (MQ=255)
              aaTCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACtcat                                                    <  1:7699811/71‑1 (MQ=255)
                  gtgtTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACtcattca                                                 >  1:2820014/1‑70 (MQ=255)
                   tgtTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACtcattcat                                                >  1:4402613/1‑70 (MQ=255)
                    gtTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGttt                                                                 <  1:4651752/52‑1 (MQ=255)
                    gtTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGttt                                                                 <  1:6688921/52‑1 (MQ=255)
                     ttGCGATACACAGCTGATTATATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCt                                              <  1:7543406/70‑1 (MQ=255)
                     ttGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGt                                                                   <  1:2728834/49‑1 (MQ=255)
                     ttGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCt                                              >  1:615799/1‑70 (MQ=255)
                         gATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTAttt                                         <  1:5623838/71‑1 (MQ=255)
                          aTACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTAtttt                                        <  1:1369947/71‑1 (MQ=255)
                          aTACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTAtttt                                        <  1:4666071/71‑1 (MQ=255)
                           tACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACtca                                                     <  1:4278364/57‑1 (MQ=255)
                           tACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTAttttt                                       <  1:6856924/71‑1 (MQ=255)
                           tACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTAttttt                                       <  1:2331540/71‑1 (MQ=255)
                           tACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTAttttt                                       <  1:7722735/71‑1 (MQ=255)
                            acacAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGt                                                                   <  1:5071444/42‑1 (MQ=255)
                            acacAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGt                                                                   <  1:3529732/42‑1 (MQ=255)
                            acacAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTATTTTtg                                      >  1:2030201/1‑71 (MQ=255)
                             cacaGCTGATGATATGTAACTTCTAATTATTAAGTATAAg                                                                    >  1:1589774/1‑40 (MQ=255)
                             cacaGCTGATGATATGTAACTTCTAATTATTAAGTATAAg                                                                    >  1:3927517/1‑40 (MQ=255)
                             cacaGCTGATGATATGTAACTTCTAATTATTAAGTATAAg                                                                    >  1:2177517/1‑40 (MQ=255)
                             cacaGCTGATGATATGTAACTTCTAATTATTAAGTATAAg                                                                    >  1:6351234/1‑40 (MQ=255)
                             cacaGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTATTTTtgt                                     >  1:5354050/1‑71 (MQ=255)
                              acaGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTtata                                                              >  1:1951783/1‑45 (MQ=255)
                              acaGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTATTTTtgtc                                    <  1:2697953/71‑1 (MQ=255)
                                aGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACtcattc                                                  >  1:1633118/1‑55 (MQ=255)
                                 gCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACt                                                       <  1:2991533/49‑1 (MQ=255)
                                 gCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTATTTTtgtctg                                  >  1:6986766/1‑70 (MQ=255)
                                     atgatATGTAACTTCTAATTATTAAGTATAAGTTta                                                                <  1:6489890/36‑1 (MQ=255)
                                     atgatATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTATTTTTGTCTGTCGCt                             <  1:4561354/71‑1 (MQ=255)
                                     atgatATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTATTTTTGTCTGTCGCt                             >  1:5702830/1‑71 (MQ=255)
                                      tgatATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTATTTTTGTCTGTCGCt                             <  1:3016612/70‑1 (MQ=255)
                                         tatGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTATTTTtgtctgt                                 >  1:4541322/1‑63 (MQ=255)
                                            taaaCTTATACTTAATAATTAGAAGTTTATAGAAAACTCATTCATCTTATTTTTGTCTGTCGCTTTAGAc                       >  1:2404236/3‑70 (MQ=25)
                                            gTAACTTCTAATTATTAAGTATAAGTTTATAGaaaa                                                         >  1:1478641/1‑36 (MQ=255)
                                            gTAACTTCTAATTATTAAGTATAAGTTTATAGaaaa                                                         >  1:3037038/1‑36 (MQ=255)
                                            gTAACTTCTAATTATTAAGTATAAGTTTATAGaaaa                                                         >  1:3027862/1‑36 (MQ=255)
                                            gTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCtt                                             <  1:4420718/48‑1 (MQ=255)
                                               aCTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTa                                            >  1:1259737/1‑46 (MQ=255)
                                               aCTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTa                                            >  1:5295045/1‑46 (MQ=255)
                                               aCTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTa                                            >  1:5557294/1‑46 (MQ=255)
                                               aCTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTATTTTTGTCTGTCGCTTTAGACtata                   <  1:5555411/71‑1 (MQ=255)
                                               aCTTATACTTAATAATTAGAAGTTTATAGAAAACTCATTCATCTTATTTTTGTCTGTCGCTTTAGACtata                   <  1:3605096/71‑1 (MQ=37)
                                               aCTTATACTTAATAATTAGAAGTTTATAGAAAACTCATTCATCTTATTTTTGTCTGTCGCTTTAGACtata                   <  1:3614554/71‑1 (MQ=37)
                                                cTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTATTTTTGTCTGTCGCttt                           >  1:4967469/1‑62 (MQ=255)
                                                cTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTATTTTTGTCTGTCGCTTTAGACtata                   <  1:6298013/70‑1 (MQ=255)
                                                cTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTATTTTTGTCTGTCGCTTTAGACTATac                  <  1:434514/71‑1 (MQ=255)
                                                   cTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTATTTTTGTCTGTCGCTTTAGACTATaca                 <  1:1429108/69‑1 (MQ=255)
                                                   cTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTATTTTTGTCTGTCGCTTTAGACTATaca                 >  1:1770393/1‑69 (MQ=255)
                                                    tAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTATTTTTGTCTGTCGc                              >  1:570792/1‑55 (MQ=255)
                                                    tAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTATTTTTGTCTGTCGc                              >  1:5041284/1‑55 (MQ=255)
                                                      attattAACTATAAGTTTATAGAAAACTCATTCATCtt                                             <  1:2362052/38‑1 (MQ=255)
                                                        tattaAGTATAAGTTTATAGAAAACTCATTCATCTTAtt                                          <  1:3629268/39‑1 (MQ=255)
                                                        tattaAGTATAAGTTTATAGAAAACTCATTCATCTTAtt                                          <  1:4180937/39‑1 (MQ=255)
                                                        tattaAGTATAAGTTTATAGAAAACTCATTCATCTTAtt                                          <  1:7076874/39‑1 (MQ=255)
                                                        tattaAGTATAAGTTTATAGAAAACTCATTCATCTTAtt                                          <  1:3408556/39‑1 (MQ=255)
                                                        tattaAGTATAAGCTTATAGAAAACTCATTCATCTTAtt                                          <  1:705084/39‑1 (MQ=255)
                                                          ttaAGTATAAGTTTATAGAAAACTCATTCATCTTATTTTtgt                                     >  1:6396115/1‑42 (MQ=255)
                                                          ttaAGTATAAGTTTATAGAAAACTCATTCATCTTATTTTTGTCTGTCGCTTTAGACTATACACAGCATAAt        <  1:1220540/71‑1 (MQ=255)
                                                             aGTATAAGTTTATAGAAAACTCATTCATCTTATTTTTGTCTGTCGCTTTAGACTATACACAGCATAAtttt     <  1:4161463/71‑1 (MQ=255)
                                                             aGTATAAGTTTATAGAAAACTCATTCATCTTATTTTTGTCTGTCGCTTTAGACTATACACAGCATAAtttt     <  1:6100145/71‑1 (MQ=255)
                                                              gTATAAGTTTATAGAAAACTCATTCATCTTATTTTTGTCTGTCGCTTTAGACTATACACAGCATAATTTTa    <  1:4528763/71‑1 (MQ=255)
                                                              gTATAAGTTTATAGAAAACTCATTCATCTTATTTTTGTCTGTCGCTTTAGACTATACACAGCATAATTTTa    <  1:5760332/71‑1 (MQ=255)
                                                              gTATAAGTTTATAGAAAACTCATTCATCTTATTTTTGTCTGTCGCTTTAGACTATACACAGCATAATTTTa    <  1:4595333/71‑1 (MQ=255)
                                                               tataAGTTTATAGAAAACTCATTCATCTTATTTTTGTCTGTCGCTTTAGACTATACACAGCATAATTTTa    <  1:1967806/70‑1 (MQ=255)
                                                                 taAGTTTATAGAAAACTCATTCATCTTATTTTtgt                                     >  1:173776/1‑35 (MQ=255)
                                                                 taAGTTTATAGAAAACTCATTCATCTTATTTTTGTCTGTCGCTTTAGACTATACACAGCATAATTTTAtt  <  1:587127/70‑1 (MQ=255)
                                                                 taAGTTTATAGAAAACTCATTCATCTTATTTTTGTCTGTCGCTTTAGACTATACACAGCATAATTTTAtt  <  1:3283064/70‑1 (MQ=255)
                                                                 |                                                                     
AGAAAAACACTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTATTTTTGTCTGTCGCTTTAGACTATACACAGCATAATTTTATT  >  AP009048/1218165‑1218299

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: