Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 1,218,216 C→A 6.1% intergenic (+4/‑109) ymgB → / → ymgC hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090481,218,2160CA 6.1% 318.4 / 3.7 97intergenic (+4/‑109)ymgB/ymgChypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (40/51);  new base A (3/3);  total (43/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 103 total aligned reads displayed.

TCTCAAAGCTGATGTGATTAGAAAAACACTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTATTTTTGTCTGTCGCTTTAGACTATACA  >  AP009048/1218146‑1218284
                                                                      |                                                                    
tctcAAAGCTGATGTGATTAGAAAAACACTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTc                                                                      >  1:644299/1‑71 (MQ=255)
tctcAAAGCTGATGTGATTAGAAAAACACTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTc                                                                      >  1:3105498/1‑71 (MQ=255)
tctcAAAGCTGATGTGATTAGAAAAACACTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTc                                                                      >  1:4939475/1‑71 (MQ=255)
 ctcAAAGCTGATGTGATTAGAAAAACACTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTc                                                                      <  1:2756468/70‑1 (MQ=255)
  tcAAAGCTGATGTGATTAGAAAAACACTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCt                                                                     >  1:364161/1‑70 (MQ=255)
   cAAAGCTGATGTGATTAGAAAAACACTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTa                                                                    <  1:6835672/70‑1 (MQ=255)
    aaaGCTGATGTGATTAGAAAAACACTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCt                                                                     <  1:6577390/68‑1 (MQ=255)
    aaaGCTGATGTGATTAGAAAAACACTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCt                                                                     <  1:7280764/68‑1 (MQ=255)
    aaaGCTGATGTGATTAGAAAAACACTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAat                                                                  >  1:5834819/1‑71 (MQ=255)
     aaGCTGATGTGATTAGAAAAACACTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTa                                                                    >  1:1232565/1‑68 (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:6312619/1‑71 (MQ=255)
                   aGAAAAACACTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGtt                                                   >  1:4260427/1‑71 (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:5623980/1‑69 (MQ=255)
                       aaaCACTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTta                                                 >  1:7762676/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)
                         acacTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTa                                                                    <  1:483245/48‑1 (MQ=255)
                         acacTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTa                                                                    <  1:2789865/48‑1 (MQ=255)
                         acacTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTa                                                                    <  1:2101808/48‑1 (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)
                           acTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTa                                                                    >  1:387391/1‑46 (MQ=255)
                             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:4396429/71‑1 (MQ=255)
                                 aaTCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACtcat                                     <  1:7699811/71‑1 (MQ=255)
                                  aTCGTGTTGCGATACACAGCTGATGATATGTAACTTCt                                                                     <  1:7105363/38‑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)
                                             aTACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTAtttt                         <  1:4666071/71‑1 (MQ=255)
                                             aTACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTAtttt                         <  1:1369947/71‑1 (MQ=255)
                                              tACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACtca                                      <  1:4278364/57‑1 (MQ=255)
                                              tACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTAttttt                        <  1:2331540/71‑1 (MQ=255)
                                              tACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTAttttt                        <  1:6856924/71‑1 (MQ=255)
                                              tACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTAttttt                        <  1:7722735/71‑1 (MQ=255)
                                               acacAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGt                                                    <  1:3529732/42‑1 (MQ=255)
                                               acacAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGt                                                    <  1:5071444/42‑1 (MQ=255)
                                               acacAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTATTTTtg                       >  1:2030201/1‑71 (MQ=255)
                                                cacaGCTGATGATATGTAACTTCTAATTATTAAGTATAAg                                                     >  1:3927517/1‑40 (MQ=255)
                                                cacaGCTGATGATATGTAACTTCTAATTATTAAGTATAAg                                                     >  1:6351234/1‑40 (MQ=255)
                                                cacaGCTGATGATATGTAACTTCTAATTATTAAGTATAAg                                                     >  1:1589774/1‑40 (MQ=255)
                                                cacaGCTGATGATATGTAACTTCTAATTATTAAGTATAAg                                                     >  1:2177517/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:3037038/1‑36 (MQ=255)
                                                               gTAACTTCTAATTATTAAGTATAAGTTTATAGaaaa                                          >  1:3027862/1‑36 (MQ=255)
                                                               gTAACTTCTAATTATTAAGTATAAGTTTATAGaaaa                                          >  1:1478641/1‑36 (MQ=255)
                                                               gTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCtt                              <  1:4420718/48‑1 (MQ=255)
                                                                  aCTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTa                             >  1:5557294/1‑46 (MQ=255)
                                                                  aCTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTa                             >  1:5295045/1‑46 (MQ=255)
                                                                  aCTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTa                             >  1:1259737/1‑46 (MQ=255)
                                                                  aCTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTATTTTTGTCTGTCGCTTTAGACtata    <  1:5555411/71‑1 (MQ=255)
                                                                  aCTTATACTTAATAATTAGAAGTTTATAGAAAACTCATTCATCTTATTTTTGTCTGTCGCTTTAGACtata    <  1:3614554/71‑1 (MQ=37)
                                                                  aCTTATACTTAATAATTAGAAGTTTATAGAAAACTCATTCATCTTATTTTTGTCTGTCGCTTTAGACtata    <  1:3605096/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)
                                                                      |                                                                    
TCTCAAAGCTGATGTGATTAGAAAAACACTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTATTTTTGTCTGTCGCTTTAGACTATACA  >  AP009048/1218146‑1218284

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: