Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 2,622,878 T→C 6.5% G393G (GGT→GGC ppk → polyphosphate kinase, component of RNA degradosome

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090482,622,8780TC 6.5% 213.5 / 5.2 76G393G (GGT→GGCppkpolyphosphate kinase, component of RNA degradosome
Reads supporting (aligned to +/- strand):  ref base T (36/35);  new base C (3/2);  total (39/37)
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

AACATTTACCGCGTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGCCAACATTCACTGGGCGAAGC  >  AP009048/2622813‑2622948
                                                                 |                                                                      
aaCATTTACCGCGTGGCGAAAGATTCACGCATCATCGACTGGATGATCCACGCCGCACATAACGGTAAGa                                                                    <  1:981944/70‑1 (MQ=255)
      tACCGCGTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCa                                                              <  1:4145235/70‑1 (MQ=255)
      tACCGCGTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCa                                                              <  1:3419858/70‑1 (MQ=255)
      tACCGCGTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCa                                                              <  1:1126552/70‑1 (MQ=255)
         cgcgTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCgt                                                          >  1:434913/1‑71 (MQ=255)
          gcgTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCgtg                                                         >  1:2657853/1‑71 (MQ=255)
            gTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCgtggt                                                       <  1:1654729/71‑1 (MQ=255)
               gCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGaa                                                                   >  1:3769428/1‑56 (MQ=255)
               gCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGtggtgg                                                     >  1:985761/1‑70 (MQ=255)
               gCGAAAGATTCACGCATCATCGACTCGATGATCCACGC‑‑CACATAACGGTAAGAAAGTCACCGTGGTGGTTg                                                  <  1:812176/71‑1 (MQ=255)
                cGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGaaa                                                                  <  1:566251/56‑1 (MQ=255)
                cGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGtt                                                   <  1:1483402/71‑1 (MQ=255)
                cGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGtt                                                   <  1:519875/71‑1 (MQ=255)
                 gAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTg                                                  <  1:527807/71‑1 (MQ=255)
                 gAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTg                                                  <  1:3393536/71‑1 (MQ=255)
                  aaaGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGtt                                                   <  1:3082607/69‑1 (MQ=255)
                  aaaGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGtt                                                   >  1:672753/1‑69 (MQ=255)
                   aaGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAg                                                >  1:1319059/1‑71 (MQ=255)
                   aaGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAg                                                <  1:2550496/71‑1 (MQ=255)
                     gATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGtt                                              <  1:3944196/71‑1 (MQ=255)
                     gATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGtt                                              <  1:4126685/71‑1 (MQ=255)
                     gATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGt                                               >  1:3627821/1‑70 (MQ=255)
                      aTTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTa                                             >  1:186345/1‑71 (MQ=255)
                       ttCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTAc                                            <  1:1765473/71‑1 (MQ=255)
                          aCGCATCATCGACTCGATGATCCACGCCGCACATAACGGCAAGaaa                                                                  <  1:2531644/46‑1 (MQ=255)
                           cGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAgg                                         >  1:3083925/1‑70 (MQ=255)
                           cGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGgc                                        <  1:3224352/71‑1 (MQ=255)
                             catcatCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGtt                                                   <  1:3144375/58‑1 (MQ=255)
                              atcatcGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTc                                                               <  1:654207/45‑1 (MQ=255)
                              atcatcGACTCGATGATCCACGCCGCACATAACGGCAAGAAAGTCACCGTGGTGGTTGAGTTACAGGcgcg                                     >  1:3373837/1‑71 (MQ=255)
                               tcatcGACTCGATGATCCACGCCGCACATAACGGTAAGa                                                                    <  1:1482925/39‑1 (MQ=255)
                               tcatcGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTg                                                  <  1:3474608/57‑1 (MQ=255)
                                 atcGACTCGATGATCCACGCCTCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGttt                                  >  1:1964970/1‑71 (MQ=255)
                                 atcGACTCGATGATCCACGCCGCACATAACGGTAAGaaa                                                                  >  1:1405025/1‑39 (MQ=255)
                                 atcGACTCGATGATCCACGCCGCACATAACGGTAAGaaa                                                                  >  1:2571132/1‑39 (MQ=255)
                                 atcGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGtt                                   >  1:774663/1‑70 (MQ=255)
                                 atcGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGt                                    >  1:4099070/1‑69 (MQ=255)
                                  tcGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTc                                 <  1:717522/71‑1 (MQ=255)
                                  tcGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTc                                 <  1:3638312/71‑1 (MQ=255)
                                   cGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTc                                                               >  1:2484083/1‑40 (MQ=255)
                                   cGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGtt                                                   >  1:4009178/1‑52 (MQ=255)
                                    gACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCgtgg                                                        >  1:203430/1‑46 (MQ=255)
                                      cTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTcga                               >  1:4140114/1‑69 (MQ=255)
                                       tCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTcgacg                             >  1:87025/1‑70 (MQ=255)
                                        cGATGATCCACGCCGCACATAACGGTAAGAAAGTCa                                                              >  1:4191543/1‑36 (MQ=255)
                                        cGATGATCCACGCCGCACATAACGGCAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTcgacg                             >  1:1332325/1‑69 (MQ=255)
                                         gatgatCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGtt                                              <  1:4159446/51‑1 (MQ=255)
                                         gatgatCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTcgacg                             >  1:1599854/1‑68 (MQ=255)
                                         gatgatCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTcgacg                             >  1:2156042/1‑68 (MQ=255)
                                         gatgatCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTcgacg                             >  1:3997364/1‑68 (MQ=255)
                                           tgatCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACa                                           >  1:865860/1‑52 (MQ=255)
                                           tgatCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACa                                           >  1:610868/1‑52 (MQ=255)
                                           tgatCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAgg                                         >  1:4169358/1‑54 (MQ=255)
                                             atCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGc                      <  1:487384/71‑1 (MQ=255)
                                              tCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGcc                     <  1:2028376/71‑1 (MQ=255)
                                              tCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGcc                     <  1:1219279/71‑1 (MQ=255)
                                              tCCACGCCGCACATAACGGCAAGAAAGTCACCGTGGTGGTTGa                                                 <  1:567371/43‑1 (MQ=255)
                                               ccACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGcc                     >  1:2642849/1‑70 (MQ=255)
                                               ccACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGcc                     <  1:776144/70‑1 (MQ=255)
                                                 aCGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTa                                             <  1:4084732/44‑1 (MQ=255)
                                                  cgccgcACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTAc                                            >  1:3493584/1‑44 (MQ=255)
                                                   gccgcACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGCCAACAt                >  1:1718043/1‑71 (MQ=255)
                                                    ccgcACATAACGGTAAGAAAGTCACCGTGGTGGtt                                                   >  1:1325362/1‑35 (MQ=255)
                                                       cacaTAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACgaagaa                        <  1:970394/59‑1 (MQ=255)
                                                        acaTAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTAc                                            >  1:1517313/1‑38 (MQ=255)
                                                        acaTAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGCCAACATTCACt           >  1:3661409/1‑71 (MQ=255)
                                                        acaTAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGCCAACATTCACt           <  1:2001613/71‑1 (MQ=255)
                                                         caTAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTa                                             >  1:3810764/1‑36 (MQ=255)
                                                         caTAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGCCAGCATTCACTg          >  1:3514553/1‑71 (MQ=255)
                                                         caTAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGCCAACATTCACTg          <  1:3937120/71‑1 (MQ=255)
                                                          aTAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGCCAACATTCACTgg         <  1:398201/71‑1 (MQ=255)
                                                            aaCGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGCCAACATTCACTGGGc       >  1:1202821/1‑71 (MQ=255)
                                                             aCGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGCCAACATTCACTGGGCg      <  1:705039/71‑1 (MQ=255)
                                                              cGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTagacga                            >  1:2207292/1‑48 (MQ=255)
                                                              cGGCAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGCCAACATTCACTGGGCg      >  1:196462/1‑70 (MQ=255)
                                                               ggTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGCCAACATTCACTGGGCGaa    <  1:3438430/71‑1 (MQ=255)
                                                                 tAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGCCAACATTCACTGGGCGAAgc  <  1:2037655/71‑1 (MQ=255)
                                                                 tAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGCCAACATTCACTGGGCGAAgc  <  1:1228730/71‑1 (MQ=255)
                                                                 |                                                                      
AACATTTACCGCGTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGCCAACATTCACTGGGCGAAGC  >  AP009048/2622813‑2622948

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: