Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 3,822,707 T→A 6.1% intergenic (‑52/+169) dinD ← / ← yicC DNA‑damage‑inducible protein/conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090483,822,7070TA 6.1% 228.1 / 6.1 82intergenic (‑52/+169)dinD/yicCDNA‑damage‑inducible protein/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (31/46);  new base A (3/2);  total (34/48)
Fisher's exact test for biased strand distribution p-value = 6.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

TTGATGATGTTCGTTCATTGCCCACTCCATGTGATATCAATTGCAAAGAAAGTAAATCTGTAACTGTATTTATATACAGTTGTGAGTATCTCACGTTCTCAGTGTGCTGGGCAAAAAAACATCGCAAAATCTGGAGCA  >  AP009048/3822638‑3822775
                                                                     |                                                                    
ttGATGATGTTCGTTCATTGCCCACTCCATGTGATATCAATTGCAAAGAAAGTAAATCTGTAACTGTAtt                                                                      <  1:1399450/70‑1 (MQ=255)
    tgatgTTCGTTCATTGCCCACTCCATGTGATATCAATTGCAAAGAAAGTAAATCTGTAACTGTATTtata                                                                  >  1:2972687/1‑70 (MQ=255)
        gttcgttcATTGCCCACTCCATGTGATATCAATTGCAAAGAAAGTAAATCTGTAACTGTATTtatat                                                                 >  1:318321/1‑67 (MQ=255)
         ttcgttcATTGCCCACTCCATGTGATATCAATTGCAAAGAAAGTAAATCTGTAACTGTATTTATATACAg                                                             <  1:3939998/70‑1 (MQ=255)
           cgttcATTGCCCACTCCATGTGATATCAATTGCAAAGAAAGTAAATCTGTAACTGTATTTATATACAGtt                                                           <  1:3927985/70‑1 (MQ=255)
           cgttcATTGCCCACTCCATGTGATATCAATTGCAAAGAAAGTAAATCTGTAACTGTATTTATATACAGTtg                                                          >  1:4127320/1‑71 (MQ=255)
           cgttcATTGCCCACTCCATGTGATATCAATTGCAAAGAAAGTAAATCTGTAACTGTATTTATATACAGTtg                                                          >  1:4207877/1‑71 (MQ=255)
            gttcATTGCCCACTCCATGTGATATCAATTGCAAAGAAAGTAAATCTGTAACTGTATTTATATACAGTtgt                                                         <  1:2014833/71‑1 (MQ=255)
            gttcATTGCCCACTCCATGTGATATCAATTGCAAAGAAAGTAAATCTGTAACTGTATTTATATACAGTtgt                                                         <  1:4282876/71‑1 (MQ=255)
            gttcATTGCCCACTCCATGTGATATCAATTGCAAAGAAAGTAAATCTGTAACTGTATTTATATACAGTtgt                                                         <  1:1854831/71‑1 (MQ=255)
             ttcATTGCCCACTCCATGTGATATCAATTGCAAAGAAAGTAAATCTGTAACTGTAttt                                                                     <  1:379813/58‑1 (MQ=255)
             ttcATTGCCCACTCCATGTGATATCAATTGCAAAGAAAGTAAATCTGTAACTGTATTTATATACAGTtgt                                                         >  1:3196697/1‑70 (MQ=255)
                   gCCCACTCCATGTGATATCAATTGCAAAGAAAGTAAATCTGTAACTGTATTTATATACAGTTGTGAGTAtc                                                  >  1:769959/1‑71 (MQ=255)
                   gCCCACTCCATGTGATATCAATTGCAAAGAAAGTAAATCTGTAACTGTATTTATATACAGTTGTGAGTAtc                                                  >  1:3759083/1‑71 (MQ=255)
                   gCCCACTCCATGTGATATCAATTGCAAAGAAAGTAAATCTGTAACTGTATTTATATACAGTTGTGAGTAtc                                                  >  1:4298974/1‑71 (MQ=255)
                   gCCCACTCCATGTGATATCAATTGCAAAGAAAGTAAATCTGTAACTGTATATAAATACAGttaca                                                        >  1:4931307/1‑62 (MQ=38)
                    cccACTCCATGTGATATCAATTGCAAAGAAAGTAAATCTGTAACTGTATTTATATACAGTTGTGAGTAtct                                                 <  1:2420059/71‑1 (MQ=255)
                    cccACTCCATGTGATATCAATTGCAAAGAAAGTAAATCTGTAACTGTATTTATATACAGTTGTGAGTAtct                                                 <  1:2176335/71‑1 (MQ=255)
                    cccACTCCATGTGATATCAATTGCAAAGAAAGTAAATCTGTAACTGTATTTATATACAGTTGTGAGTAtc                                                  <  1:2815607/70‑1 (MQ=255)
                      cACTCCATGTGATATCAATTGCAAAGAAAGTAAATCTGTAACTGTATTTATATACAGTtgt                                                         <  1:3330602/61‑1 (MQ=255)
                      cACTCCATGTGATATCAATTGCAAAGAAAGTAAATCTGTAACTGTATTTATATACAGTTGTGAGTATCTCa                                               <  1:380158/71‑1 (MQ=255)
                       aCTCCATGTGATATCAATTGCAAAGAAAGTAAATCTGTAACTGTATTTATATACAGTTGTGAGTATCTCAc                                              <  1:1359706/71‑1 (MQ=255)
                        cTCCATGTGATATCAATTGCAAAGAAAGTAAATCTGTAACTGTAttt                                                                     <  1:2553804/47‑1 (MQ=255)
                        cTCCATGTGATATCAATTGCAAAGAAAGTAAATCTGTAACTGTATTTATATACAGTTGTGAGTAt                                                   <  1:4941574/65‑1 (MQ=255)
                        cTCCATGTGATATCAATTGCAAAGAAAGTAAATCTGTAACTGTATATAAATACAGtta                                                          >  1:3421051/1‑57 (MQ=38)
                         tCCATGTGATATCAATTGCAAAGAAAGTAAATCTGTAACTGTATTTATATACa                                                              >  1:3953569/1‑53 (MQ=255)
                           cATGTGATATCAATTGCAAAGAAAGTAAATCTGTAACTGTATTTATATACAGTTGTGAg                                                      >  1:1869591/1‑59 (MQ=255)
                           cATGTGATATCAATTGCAAAGAAAGTAAATCTGTAACTGTATTTATATACAGTTGTGAGTATCTCACGtt                                           <  1:1653698/70‑1 (MQ=255)
                            aTGTGATATCAATTGCAAAGAAAGTAAATCTGTAACTGTAttt                                                                     <  1:1064641/43‑1 (MQ=255)
                            aTGTGATATCAATTGCAAAGAAAGTAAATCTGTAACTGTAttt                                                                     <  1:343315/43‑1 (MQ=255)
                            aTGTGATATCAATTGCAAAGAAAGTAAATCTGTAACTGTATTTATATACAGTTGTGAGTATCTCACGTtct                                         <  1:1420895/71‑1 (MQ=255)
                             tgtgATATCAATTGCAAAGAAAGTAAATCTGTAACTGTATTtatata                                                                >  1:1268833/1‑47 (MQ=255)
                               tgATATCAATTGCAAAGAAAGTAAATCTGTAACTGTATTTATATACAGTTGTGAGTATCTCACGTTCTCa                                       <  1:3237464/70‑1 (MQ=255)
                               tgATATCAATTGCAAAGAAAGTAAATCTGTAACTGTATTTATATACAGTTGTGAGTATCTCACGTTCTCa                                       <  1:3121878/70‑1 (MQ=255)
                               tgATATCAATTGCAAAGAAAGTAAATCTGTAACTGTATTTATATACAGTTGTGAGTATCTCACGTTCTCa                                       <  1:1598299/70‑1 (MQ=255)
                                gATATCAATTGCAAAGAAAGTAAATCTGTAACTGTATTTATATACAGTTGTGAGTATCTCACGTTCTCAgt                                     >  1:1789966/1‑71 (MQ=255)
                                  tatCAATTGCAAAGAAAGTAAATCTGTAACTGTATTTATATACAGTTGTGAGTATCTCACGTTCTCAgtgt                                   <  1:2206555/71‑1 (MQ=255)
                                  tatCAATTGCAAAGAAAGTAAATCTGTAACTGTATTTATATACAGTTGTGAGTATCTCACGTTCTCAgtgt                                   >  1:3202871/1‑71 (MQ=255)
                                  tatCAATTGCAAAGAAAGTAAATCTGTAACTGTATTTATATACAGTTGTGAGTATCTCACGTTCTCAgtgt                                   <  1:2321082/71‑1 (MQ=255)
                                   atCAATTGCAAAGAAAGTAAATCTGTAACTGTATTTATATACAGTtgt                                                         <  1:1114584/48‑1 (MQ=255)
                                    tCAATTGCAAAGAAAGTAAATCTGTAACTGTATTTATATACAGTTGTGAGTATCTCACGTTCTCa                                       <  1:4406789/65‑1 (MQ=255)
                                     cAATTGCAAAGAAAGTAAATCTGTAACTGTATTTATATAc                                                               >  1:1884330/1‑40 (MQ=255)
                                     cAATTGCAAAGAAAGTAAATCTGTAACTGTATTTATATACAGTTGTGAGt                                                     >  1:5139111/1‑50 (MQ=255)
                                     cAATTGCAAAGAAAGTAAATCTGTAACTGTATTTATATACAGTTGTGAGTATCTCACGTTCTCAGTGTGCt                                >  1:2197137/1‑71 (MQ=255)
                                     cAATTGCAAAGAAAGTAAATCTGTAACTGTATTTATATACAGTTGTGAGTATCTCACGTTCTCAGTGTGCt                                >  1:2029350/1‑71 (MQ=255)
                                     cAATTGCAAAGAAAGTAAATCTGTAACTGTATAtaaat                                                                 <  1:161974/38‑1 (MQ=37)
                                      aaTTGCAAAGAAAGTAAATCTGTAACTGTATTTATATACAGTTGTGAGTATCTCACGTTCTCAGTGTGCt                                <  1:2643253/70‑1 (MQ=255)
                                      aaTTGCAAAGAAAGTAAATCTGTAACTGTATTTATATACAGTTGTGAGTATCTCACGTTCTCAGTGTGCt                                <  1:147040/70‑1 (MQ=255)
                                      aaTTGCAAAGAAAGTAAATCTGTAACTGTATTTATATACAGTTGTGAGTATCTCACGTTCTCAGTGTGCt                                <  1:4750845/70‑1 (MQ=255)
                                      aaTTGCAAAGAAAGTAAATCTGTAACTGTATTTATATACAGTTGTGAGTATCTCACGTTCTCACTGTGCt                                <  1:1905/70‑1 (MQ=255)
                                       aTTGCAAAGAAAGTAAATCTGTAACTGTATTTATATACAGTTGTGAGTATCTCACGTTCTCAGTGTGCTgg                              <  1:2693136/71‑1 (MQ=255)
                                        ttGCAAAGAAAGTAAATCTGTAACTGTATTTATATACAGTtgtg                                                        <  1:4052312/44‑1 (MQ=255)
                                        ttGCAAAGAAAGTAAATCTGTAACTGTATTTATATACAGTtgtg                                                        <  1:138747/44‑1 (MQ=255)
                                        ttGCAAAGAAAGTAAATCTGTAACTGTATTTATATACAGTTGTGAGTATCTCACGTTCTCAGTGTGCTgg                              >  1:2790653/1‑70 (MQ=255)
                                         tGCAAAGAAAGTAAATCTGTAACTGTATTTATATACAGTTGTGAGTATCTCACGTTCTCAGTGTGCTggg                             >  1:433961/1‑70 (MQ=255)
                                          gCAAAGAAAGTAAATCTGTAACTGTATTTATATACAGTTGTGAGTATCTCACGTTCTCAGTGTGCTGGGc                            >  1:4853247/1‑70 (MQ=255)
                                          gCAAAGAAAGTAAATCTGTAACTGTATATAAATAc                                                               <  1:4141512/35‑1 (MQ=25)
                                           cAAAGAAAGTAAATCTGTAACTGTATTTATATACAGTTGTGAGTATCTCACGTTCTCAGTGTGCTGGGCaa                          <  1:2756221/71‑1 (MQ=255)
                                            aaagaaagTAAATCTGTAACTGTATTTATATACAGTTGTGAGTATCTCACGTTCTCAGTGTGCTGGGCaa                          >  1:2935288/1‑70 (MQ=255)
                                             aagaaagTAAATCTGTAACTGTATTTATATACAGTTGTGAGTATCTCACGTtct                                         >  1:1708355/1‑54 (MQ=255)
                                              agaaagTAAATCTGTAACTGTATTTATATACAGTTGTGAGTATCTCACGTTCTCAGTGTGCTGGGCaaaaa                       <  1:624766/71‑1 (MQ=255)
                                              agaaagTAAATCTGTAACTGTATTTATATACAGTTGTGAGTATCTCACGTTCTCAGTGTGCTGGGCaaaaa                       <  1:123869/71‑1 (MQ=255)
                                              agaaagTAAATCTGTAACTGTATTTATATACAGTTGTGAGTATCTCACGTTCTCAGTGTGCTGGGCaaaa                        >  1:4409622/1‑70 (MQ=255)
                                              agaaagTAAATCTGTAACTGTATTTATATACAGTTGTGAGTATCTCACGTTCTCAGTGTGCTGGGCaaaa                        >  1:553242/1‑70 (MQ=255)
                                               gaaagTAAATCTGTAACTGTATTTATATACAGTTGTGAGTATCTCACGTTCTCAGTGTGCTGGGCaaaaaa                      <  1:2876694/71‑1 (MQ=255)
                                               gaaagTAAATCTGTAACTGTATTTATATACAGTTGTGAGTATCTCACGTTCTCAGTGTGCTGGGCaaaaaa                      <  1:4318295/71‑1 (MQ=255)
                                               gaaagTAAATCTGTAACTGTATTTATATACAGTTGTGAGTATCTCACGTTCTCAGTGTGCTGGGCaaaaa                       >  1:259040/1‑70 (MQ=255)
                                                aaagTAAATCTGTAACTGTATTTATATACAGTTGTGAGTATCTCACGTTCTCAGTGTGCTGGGCaaaaaaa                     <  1:2476593/71‑1 (MQ=255)
                                                aaagTAAATCTGTAACTGTACTTATATACAGTTGTGAGTATCTCACGTTCTCAGTGTGCTGGGCaaaaaaa                     >  1:1082854/1‑71 (MQ=255)
                                                 aagTAAATCTGTAACTGTATTTATATACAGTTGTGAGTATCTCACGTTCTCAgt                                     >  1:5113015/1‑54 (MQ=255)
                                                 aagTAAATCTGTAACTGTATTTATATACAGTTGTGAGTATCTCACGTTCTCAGTGTGCTGGGCAAAAAAAc                    <  1:1448745/71‑1 (MQ=255)
                                                  agTAAATCTGTAACTGTATTTATATACAGTTGTGAGTATCTCACGTTCTCAGTGTGCTGGGCAAAAAAACa                   <  1:4361878/71‑1 (MQ=255)
                                                  agTAAATCTGTAACTGTATTTATATACAGTTGTGAGTATCTCACGTTCTCAGTGTGCTGGGCAAAAAAACa                   <  1:4735178/71‑1 (MQ=255)
                                                   gTAAATCTGTAACTGTATTTATATACAGTTGTGAGTATCTCACGTTCTCAGTGTGCTGGGCAAAAAAACAt                  <  1:3917765/71‑1 (MQ=255)
                                                    tAAATCTGTAACTGTATTTATATACAGTTGTGAGTAtct                                                 <  1:3532729/39‑1 (MQ=255)
                                                    tAAATCTGTAACTGTATTTATATACAGTTGTGAGTATCTCACGTTCTCAGTGTGCTGGGCa                           <  1:2701052/61‑1 (MQ=255)
                                                    tAAATCTGTAACTGTATTTATATACAGTTGTGAGTATCTCACGTTCTCAGTGTGCTGGGCa                           <  1:4719837/61‑1 (MQ=255)
                                                    tAAATCTGTAACTGTATTTATATACAGTTGTGAGTATCTCACGTTCTCAGTGTGCTGGGCAAAAAAACATc                 <  1:1238805/71‑1 (MQ=255)
                                                         cTGTAACTGTATTTATATACAGTTGTGAGTATCTCACGTTCTCAGTGTGCTGGGCAAAAAAACATCGCaa             >  1:303173/1‑70 (MQ=255)
                                                         cTGTAACTGTATTTATATACAGTTGTGAGTATCTCACGTTCTCAGTGTGCTGGGCAAAAAAACATCGCaa             >  1:4651628/1‑70 (MQ=255)
                                                           acaaCTGTATATAAATACAGTTGTGAGTATCTCTCGTTCTCAGTGTGCTGGGCAAAAAAACATCGCAAAAt          >  1:4928947/3‑71 (MQ=38)
                                                            tAACTGTATTTATATACAGTTGTGAGTATCTCACGTTCTCAGTGTGCTGGGCAAAAAAACATCGCaaa            <  1:1192302/68‑1 (MQ=255)
                                                            tAACTGTATTTATATACAGTTGTGAGTATCTCACGTTCTCAGTGTGCTGGGCAAAAAAACATCGCAAAATc         <  1:3610597/71‑1 (MQ=255)
                                                              aCTGTATTTATATACAGTTGTGAGTATCTCACGTTCTCAg                                      >  1:2409159/1‑40 (MQ=255)
                                                               cTGTATTTATATACAGTTGTGAGTATCTCACGTTCTCAGTGTGCTg                               <  1:971958/46‑1 (MQ=255)
                                                                 gTATTTATATACAGTTGTGAGTATCTCACGTTCTCAGTGTGCTggg                             >  1:5093385/1‑46 (MQ=255)
                                                                  tATTTATATACAGTTGTGAGTATCTCACGTTCTCAGtgtg                                  <  1:2249243/40‑1 (MQ=255)
                                                                    tttATATACAGTTGTGAGTATCTCACGTTCTCAGTGTGCTGGGCaaaaaa                      <  1:3770157/50‑1 (MQ=255)
                                                                    tttATATACAGTTGTGAGTATCTCACGTTCTCAGTGTGCTGGGCaaaaaa                      <  1:1302493/50‑1 (MQ=255)
                                                                    tttATATACAGTTGTGAGTATCTCACGTTCTCAGTGTGCTGGGCAAAAAAACATCGCAAAATCTGGAGCa  <  1:4848962/70‑1 (MQ=255)
                                                                     ttATATACAGTTGTGAGTATCTCACGTTCTCAGTGTGCt                                <  1:1123590/39‑1 (MQ=255)
                                                                     |                                                                    
TTGATGATGTTCGTTCATTGCCCACTCCATGTGATATCAATTGCAAAGAAAGTAAATCTGTAACTGTATTTATATACAGTTGTGAGTATCTCACGTTCTCAGTGTGCTGGGCAAAAAAACATCGCAAAATCTGGAGCA  >  AP009048/3822638‑3822775

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: