Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 1,218,219 A→C 7.3% intergenic (+7/‑106) ymgB → / → ymgC hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090481,218,2190AC 7.3% 209.0 / 10.0 82intergenic (+7/‑106)ymgB/ymgChypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (38/38);  new base C (3/3);  total (41/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAGCTGATGTGATTAGAAAAACACTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTATTTTTGTCTGTCGCTTTAGACTATACA  >  AP009048/1218150‑1218284
                                                                     |                                                                 
aaaGCTGATGTGATTAGAAAAACACTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAat                                                                  >  1:5834819/1‑71 (MQ=255)
    cTGATGTGATTAGAAAAACACTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAattatt                                                              <  1:4829954/71‑1 (MQ=255)
    cTGATGTGATTAGAAAAACACTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAattatt                                                              <  1:1637612/71‑1 (MQ=255)
       aTGTGATTAGAAAAACACTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAg                                                           <  1:3309601/71‑1 (MQ=255)
         gtgATTAGAAAAACACTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTaa                                                                   <  1:7076652/61‑1 (MQ=255)
         gtgATTAGAAAAACACTGGAAATCGTGTTGAGATACACAGCTGATGATATGTAACTTCTAATTATTAAGta                                                         >  1:6170359/1‑71 (MQ=255)
           gATTAGAAAAACACTTGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAattat                                                               <  1:5499315/63‑1 (MQ=255)
           gATTAGAAAAACACTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAattat                                                               <  1:6259390/63‑1 (MQ=255)
           gATTAGAAAAACACTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAattat                                                               <  1:2605406/63‑1 (MQ=255)
           gATTAGAAAAACACTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAatt                                                                 >  1:5367309/1‑61 (MQ=255)
            aTTAGAAAAACACTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAattatt                                                              <  1:6936826/63‑1 (MQ=255)
               aGAAAAACACTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGtt                                                   >  1:4260427/1‑71 (MQ=255)
               aGAAAAACACTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGtt                                                   >  1:6312619/1‑71 (MQ=255)
                gAAAAACACTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGtt                                                   <  1:6869511/70‑1 (MQ=255)
                gAAAAACACTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGtt                                                   <  1:3456177/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:5623980/1‑69 (MQ=255)
                   aaaCACTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTta                                                 >  1:4043939/1‑69 (MQ=255)
                    aaCACTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAg                                              <  1:1436108/71‑1 (MQ=255)
                    aaCACTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAg                                              <  1:72822/71‑1 (MQ=255)
                     acacTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGta                                                         >  1:975691/1‑59 (MQ=255)
                     acacTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGta                                                         >  1:699095/1‑59 (MQ=255)
                     acacTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGta                                                         >  1:7841438/1‑59 (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:339713/1‑71 (MQ=255)
                         tGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAAc                                         >  1:3684896/1‑71 (MQ=255)
                           gAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGtt                                                   <  1:6151840/59‑1 (MQ=255)
                           gAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACtc                                       <  1:2550131/71‑1 (MQ=255)
                             aaTCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTta                                                 <  1:2918116/59‑1 (MQ=255)
                             aaTCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTta                                                 <  1:5647634/59‑1 (MQ=255)
                             aaTCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTta                                                 <  1:2153894/59‑1 (MQ=255)
                             aaTCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACtcat                                     <  1:7699811/71‑1 (MQ=255)
                             aaTCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACtcat                                     <  1:4396429/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)
                                      gCGATACACAGCTGATGATATGTAACTTCTAAttatta                                                             >  1:2598859/1‑38 (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:7722735/71‑1 (MQ=255)
                                          tACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTAttttt                        <  1:2331540/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:6351234/1‑40 (MQ=255)
                                            cacaGCTGATGATATGTAACTTCTAATTATTAAGTATAAg                                                     >  1:2177517/1‑40 (MQ=255)
                                            cacaGCTGATGATATGTAACTTCTAATTATTAAGTATAAg                                                     >  1:3927517/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:5702830/1‑71 (MQ=255)
                                                    atgatATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTATTTTTGTCTGTCGCt              <  1:4561354/71‑1 (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:3027862/1‑36 (MQ=255)
                                                           gTAACTTCTAATTATTAAGTATAAGTTTATAGaaaa                                          >  1:3037038/1‑36 (MQ=255)
                                                           gTAACTTCTAATTATTAAGTATAAGTTTATAGaaaa                                          >  1:1478641/1‑36 (MQ=255)
                                                           gTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCtt                              <  1:4420718/48‑1 (MQ=255)
                                                              aCTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTa                             >  1:1259737/1‑46 (MQ=255)
                                                              aCTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTa                             >  1:5557294/1‑46 (MQ=255)
                                                              aCTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTa                             >  1:5295045/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:1770393/1‑69 (MQ=255)
                                                                  cTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTATTTTTGTCTGTCGCTTTAGACTATaca  <  1:1429108/69‑1 (MQ=255)
                                                                   tAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTATTTTTGTCTGTCGc               >  1:570792/1‑55 (MQ=255)
                                                                   tAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTATTTTTGTCTGTCGc               >  1:5041284/1‑55 (MQ=255)
                                                                     attattAACTATAAGTTTATAGAAAACTCATTCATCtt                              <  1:2362052/38‑1 (MQ=255)
                                                                     |                                                                 
AAAGCTGATGTGATTAGAAAAACACTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTATTTTTGTCTGTCGCTTTAGACTATACA  >  AP009048/1218150‑1218284

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: