Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 2,632,686 A→G 7.4% V14A (GTT→GCT)  guaB ← IMP dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090482,632,6860AG 7.4% 264.4 / 11.3 95V14A (GTT→GCT) guaBIMP dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (43/45);  new base G (3/4);  total (46/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.86e-01

GTTTTCGTCAGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAATATCTCGACCAGAGTGGTTAATAA  >  AP009048/2632617‑2632756
                                                                     |                                                                      
gTTTTCGTCAGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAAc                                                                       <  1:936661/71‑1 (MQ=255)
gTTTTCGTCAGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAAc                                                                       >  1:3463362/1‑71 (MQ=255)
   ttCGTCAGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAAcgtc                                                                    >  1:843748/1‑71 (MQ=255)
   ttCGTCAGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAAcgtc                                                                    >  1:4198404/1‑71 (MQ=255)
   ttCGTCAGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAAcgtc                                                                    >  1:3293177/1‑71 (MQ=255)
   ttCGTCAGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAAcgtc                                                                    >  1:2387162/1‑71 (MQ=255)
   ttCGTCAGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAAcgt                                                                     <  1:3955218/70‑1 (MQ=255)
     cGTCAGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCATGAACGAGGAGAAcgtcg                                                                   <  1:971555/70‑1 (MQ=255)
      gTCAGCTGGGTTCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAAcgtcgt                                                                  >  1:4218421/1‑70 (MQ=255)
      gTCAGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAAcgtcgt                                                                  >  1:3240938/1‑70 (MQ=255)
        cAGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAGCGTCGTCa                                                                >  1:1870659/1‑70 (MQ=255)
         aGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCaaa                                                              >  1:120795/1‑71 (MQ=255)
         aGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCaaa                                                              >  1:2499518/1‑71 (MQ=255)
         aGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCaaa                                                              >  1:1326979/1‑71 (MQ=255)
         aGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCaa                                                               <  1:4543416/70‑1 (MQ=255)
          gCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAAc                                                             >  1:451970/1‑71 (MQ=255)
            tGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAAc                                                             >  1:2768371/1‑69 (MQ=255)
              ggTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGt                                                           <  1:620867/69‑1 (MQ=255)
              ggTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGt                                                           <  1:615437/69‑1 (MQ=255)
               gTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAGCGTCGTCAAACGTCa                                                         <  1:4282440/70‑1 (MQ=255)
                 gCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAGCGTCGTCAAACGTCa                                                         >  1:783124/1‑68 (MQ=255)
                 gCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAAcgtcg                                                                   <  1:2728138/58‑1 (MQ=255)
                  cTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAAcgtcgt                                                                  >  1:292388/1‑58 (MQ=255)
                     aGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGcttc                                                  >  1:3484445/1‑71 (MQ=255)
                       gTCATCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGcttctt                                                <  1:3285876/71‑1 (MQ=255)
                       gTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGcttctt                                                >  1:2418293/1‑71 (MQ=255)
                        tCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGcttc                                                  <  1:2068200/68‑1 (MQ=255)
                        tCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCttt                                               <  1:4140770/71‑1 (MQ=255)
                        tCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCttt                                               <  1:908402/71‑1 (MQ=255)
                         cagcagTATTCGGCAGAACTGTAGAGTGAGCAGGAACGAGGAGAGCGTCGTCAAACGTCAGAGcttctt                                                >  1:3319172/1‑69 (MQ=255)
                         cagcagTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGcttctt                                                >  1:1105944/1‑69 (MQ=255)
                         cagcagTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGcttctt                                                >  1:721170/1‑69 (MQ=255)
                          agcagTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTa                                              <  1:3373887/70‑1 (MQ=255)
                          agcagTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTa                                              <  1:4403360/70‑1 (MQ=255)
                          agcagTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAg                                             <  1:4268072/71‑1 (MQ=255)
                          agcagTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAg                                             >  1:3279101/1‑71 (MQ=255)
                           gcagTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAGCGTCGTCAAACGTCAGAGCTTCTTTAGc                                            <  1:571942/71‑1 (MQ=255)
                           gcagTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAGCGTCGTCAAACGTCAGAGCTTCTTTAGc                                            <  1:577346/71‑1 (MQ=255)
                           gcagTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAg                                             <  1:1393204/70‑1 (MQ=255)
                            cagTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGc                                            >  1:2565704/1‑70 (MQ=255)
                             agTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGa                                          <  1:3319515/71‑1 (MQ=255)
                              gTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGAt                                         <  1:2648459/71‑1 (MQ=255)
                               tATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAAcg                                                                      <  1:1972300/41‑1 (MQ=255)
                               tATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGct                                                    <  1:1377393/59‑1 (MQ=255)
                               tATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATa                                        <  1:3450484/71‑1 (MQ=255)
                                 ttCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAAcgtcg                                                                   >  1:2132854/1‑42 (MQ=255)
                                 ttCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAAc                                                             >  1:1938567/1‑48 (MQ=255)
                                 ttCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACg                                      >  1:2535161/1‑71 (MQ=255)
                                 ttCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACg                                      >  1:3869003/1‑71 (MQ=255)
                                  tCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGAt                                         >  1:2887336/1‑67 (MQ=255)
                                     gCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGt                                                           <  1:4396596/46‑1 (MQ=255)
                                     gCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGcttctt                                                <  1:4177489/57‑1 (MQ=255)
                                      cAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCagag                                                      >  1:119215/1‑50 (MQ=255)
                                       aGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCa                                                         >  1:4057737/1‑46 (MQ=255)
                                       aGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCAt                                <  1:3510066/71‑1 (MQ=255)
                                       aGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCAt                                >  1:552849/1‑71 (MQ=255)
                                        gAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATg                               <  1:3216381/71‑1 (MQ=255)
                                         aaCGGTAGAGTGAGCAGGAACAAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATgg                              >  1:2363340/1‑71 (MQ=255)
                                          aCGGTAGAGTGACCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATggg                             <  1:678377/71‑1 (MQ=255)
                                           cGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATggg                             >  1:4532563/1‑70 (MQ=255)
                                            ggTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCa                           <  1:2467714/71‑1 (MQ=255)
                                             gTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGc                                  <  1:1513876/63‑1 (MQ=255)
                                                gagTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAatat                       >  1:3170744/1‑71 (MQ=255)
                                                gagTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAata                        <  1:274602/70‑1 (MQ=255)
                                                 agTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCa                                 <  1:3226850/60‑1 (MQ=255)
                                                   tGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGa                                          <  1:1225564/49‑1 (MQ=255)
                                                   tGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAATAtctc                    <  1:658466/71‑1 (MQ=255)
                                                   tGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAATAtctc                    <  1:2649517/71‑1 (MQ=255)
                                                   tGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGAGATACGTAGCATGGGCAATAtctc                    <  1:4331747/71‑1 (MQ=255)
                                                     aGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGc                                            >  1:2172205/1‑45 (MQ=255)
                                                     aGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAATATCTCGa                  >  1:4035132/1‑71 (MQ=255)
                                                      gCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGAt                                         <  1:2148821/47‑1 (MQ=255)
                                                       cAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCaa                          >  1:691499/1‑61 (MQ=255)
                                                        aGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAATATCTCGa                  >  1:1167596/1‑68 (MQ=255)
                                                        aGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAATATCTCGAc                 >  1:1162234/1‑69 (MQ=255)
                                                        aGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAATATCTCGACCa               <  1:839762/71‑1 (MQ=255)
                                                           aaCGAGGAGAGCGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAATATCTCGACCagag            <  1:3182676/71‑1 (MQ=255)
                                                           aaCGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAat                         <  1:1807947/58‑1 (MQ=255)
                                                           aaCGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAATATCTCGACCagag            <  1:4485829/71‑1 (MQ=255)
                                                           aaCGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAATATCTCGACCagag            <  1:4105170/71‑1 (MQ=255)
                                                           aaCGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAATATCTCGACCagag            <  1:666743/71‑1 (MQ=255)
                                                            aCGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAg                                             >  1:3918792/1‑37 (MQ=255)
                                                            aCGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATa                                        >  1:1427476/1‑42 (MQ=255)
                                                            aCGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAatat                       >  1:2161411/1‑59 (MQ=255)
                                                              gaggagAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAATATCTCGACCAGAGTgg         <  1:2565444/71‑1 (MQ=255)
                                                               aggagAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAATATCTCGACCAGAGTGGt        <  1:4116599/71‑1 (MQ=255)
                                                                ggagAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAATATCTCGACCAGAGTGGtt       <  1:700083/71‑1 (MQ=255)
                                                                ggagAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAATATCTCGACCAGAGTGGt        >  1:1272937/1‑70 (MQ=255)
                                                                 gagaACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAATATCTCGACCAGAGTGGTta      <  1:3032602/71‑1 (MQ=255)
                                                                  agagCGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTa                                    >  1:2127958/1‑40 (MQ=255)
                                                                  agaACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGt                                     >  1:748184/1‑39 (MQ=255)
                                                                  agaACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAATATCTCg                   >  1:3761734/1‑57 (MQ=255)
                                                                   gaACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCAt                                >  1:514061/1‑43 (MQ=255)
                                                                   gaACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAATATCTCGACCAGAGTGGTtaat    <  1:1902516/71‑1 (MQ=255)
                                                                   gaACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAATATCTCGACCAGAGTGGTtaat    <  1:461411/71‑1 (MQ=255)
                                                                   gaACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAATATCTCGACCAGAGTGGTtaa     >  1:4496886/1‑70 (MQ=255)
                                                                    aaCGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAATATCTCGACCAGAGTGGTtaata   <  1:1421828/71‑1 (MQ=255)
                                                                     aCGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAATATCTCGACCAGAGTGGt        >  1:3019486/1‑65 (MQ=255)
                                                                     aCGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAATATCTCGACCAGAGTGGTtaataa  <  1:2911832/71‑1 (MQ=255)
                                                                     |                                                                      
GTTTTCGTCAGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAATATCTCGACCAGAGTGGTTAATAA  >  AP009048/2632617‑2632756

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑