Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 71,214 T→C 51.1% F276F (TTT→TTC araC → DNA‑binding transcriptional dual regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP00904871,2140TC51.1% 7.4 / 130.1 88F276F (TTT→TTCaraCDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base T (24/19);  new base C (19/26);  total (43/45)
Fisher's exact test for biased strand distribution p-value = 2.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

GGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAA  >  AP009048/71146‑71284
                                                                    |                                                                      
ggTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCg                                                                      >  1:2528592/1‑71 (MQ=255)
  ttttGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGt                                                                     >  1:1190881/1‑70 (MQ=255)
    ttGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGcc                                                                  >  1:1142817/1‑71 (MQ=255)
    ttGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGcc                                                                  >  1:4264195/1‑71 (MQ=255)
    ttGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGcc                                                                  >  1:3639008/1‑71 (MQ=255)
    ttGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGc                                                                   <  1:2561113/70‑1 (MQ=255)
     tGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGcc                                                                  <  1:1894802/70‑1 (MQ=255)
      gACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCgg                                                                >  1:3893690/1‑71 (MQ=255)
        cGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGtt                                                              <  1:2463305/71‑1 (MQ=255)
          aTCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTtgt                                                            <  1:515647/71‑1 (MQ=255)
          aTCAACTCTATTTCTCGCGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTtgtg                                                           <  1:1211173/71‑1 (MQ=255)
              aCTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCgg                                                                <  1:19709/63‑1 (MQ=255)
                tctATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTgaaga                                                       >  1:2458288/1‑70 (MQ=255)
                 ctATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTgaa                                                         >  1:438859/1‑67 (MQ=255)
                  tATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTgaagaa                                                      <  1:1419445/69‑1 (MQ=255)
                  tATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTgaagaa                                                      <  1:1360350/69‑1 (MQ=255)
                  tATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGaaa                                                     >  1:3086791/1‑70 (MQ=255)
                  tATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGaaa                                                     >  1:3030251/1‑70 (MQ=255)
                  tATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGaaa                                                     >  1:3149851/1‑70 (MQ=255)
                    tttCTCGCGGGTATTTAAATAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGaaaa                                                    <  1:563159/69‑1 (MQ=255)
                    tttCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGaaaa                                                    <  1:3958265/69‑1 (MQ=255)
                    tttCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAg                                                  >  1:4480953/1‑71 (MQ=255)
                     ttCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGt                                                 >  1:3502836/1‑71 (MQ=255)
                      tctcGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGt                                                 <  1:2687299/70‑1 (MQ=255)
                      tctcGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGt                                                 >  1:2275071/1‑70 (MQ=255)
                       ctcGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTg                                                >  1:3038363/1‑70 (MQ=255)
                          gcgTGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAat                                             <  1:3895845/70‑1 (MQ=255)
                          gcgAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAatg                                            <  1:2675538/71‑1 (MQ=255)
                          gcgAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAatg                                            <  1:3578753/71‑1 (MQ=255)
                                  ttAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGtt                                                              <  1:1930414/45‑1 (MQ=255)
                                   tAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGcc                                    >  1:1428562/1‑70 (MQ=255)
                                    aaaaaaTGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGtg                                                              >  1:4361205/1‑42 (MQ=255)
                                     aaaaaTGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTc                                 <  1:1754009/71‑1 (MQ=255)
                                          tGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCTGTCAAGtt                            <  1:2379224/71‑1 (MQ=255)
                                          tGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGtt                            <  1:1050342/71‑1 (MQ=255)
                                           gCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGa                                               >  1:68548/1‑51 (MQ=255)
                                            cACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTgaaga                                                       >  1:2554581/1‑42 (MQ=255)
                                            cACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTg                           <  1:3530426/70‑1 (MQ=255)
                                             aCCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGaaaa                                                    <  1:1784285/44‑1 (MQ=255)
                                              ccGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCa                        >  1:2817484/1‑71 (MQ=255)
                                               cGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAg                                                  <  1:1184051/44‑1 (MQ=255)
                                               cGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAg                                                  <  1:327579/44‑1 (MQ=255)
                                               cGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTTAATGATGTAGCCGTCAAGTTGTCAt                       <  1:3736590/71‑1 (MQ=255)
                                               cGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCa                        >  1:4188276/1‑70 (MQ=255)
                                               cGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGa                                               <  1:2527642/47‑1 (MQ=255)
                                               cGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCAt                       <  1:1086833/71‑1 (MQ=255)
                                                ggggCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCAt                       >  1:1862494/1‑70 (MQ=255)
                                                 gggCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTg                                                >  1:938993/1‑44 (MQ=255)
                                                 gggCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATaa                     >  1:1468122/1‑71 (MQ=255)
                                                 gggCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGaaa                                                     <  1:1455348/39‑1 (MQ=255)
                                                 gggCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGaaa                                                     >  1:1269222/1‑39 (MQ=255)
                                                 gggCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATaa                     <  1:3602746/71‑1 (MQ=255)
                                                  ggCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAAt                    <  1:3892948/71‑1 (MQ=255)
                                                  ggCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTc                         >  1:2722721/1‑66 (MQ=255)
                                                  ggCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAAa                    <  1:368998‑M1/71‑2 (MQ=255)
                                                    ccAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTg                  <  1:43006/71‑1 (MQ=255)
                                                    ccAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGaaaaa                                                   <  1:2087049/38‑1 (MQ=255)
                                                    ccAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAAagg                  <  1:4577005‑M1/71‑4 (MQ=255)
                                                     cAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTgg                 >  1:2700420/1‑71 (MQ=255)
                                                     cAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTgg                 >  1:1055326/1‑71 (MQ=255)
                                                      aGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTgg                 >  1:41173/1‑70 (MQ=255)
                                                      aGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGt                                                 <  1:963422/38‑1 (MQ=255)
                                                      aGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAAagggt                >  1:3601592‑M1/1‑66 (MQ=255)
                                                       gCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAatgat                                          <  1:1567908/44‑1 (MQ=255)
                                                       gCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAatg                                            <  1:815351/42‑1 (MQ=255)
                                                        cccGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATaa                     <  1:3972539/64‑1 (MQ=255)
                                                        cccGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAAagggtcc              >  1:40117‑M1/1‑64 (MQ=255)
                                                         ccGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAAagggtccc             >  1:4594555‑M1/1‑63 (MQ=255)
                                                          cgagcgagTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAAc             >  1:1189812/1‑70 (MQ=255)
                                                          cgagcgagTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAAc             >  1:88110/1‑70 (MQ=255)
                                                          cgagcgagTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACg            <  1:3387115/71‑1 (MQ=255)
                                                          cgagcgagTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAg                              >  1:3540285/1‑53 (MQ=255)
                                                          cgagcgagTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAAagggtccc             >  1:3476460‑M1/1‑62 (MQ=255)
                                                          cgagcgacTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAAc             >  1:409047/1‑70 (MQ=255)
                                                           gagcgagTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGt                <  1:1549478/66‑1 (MQ=255)
                                                            agcgagTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGaa          >  1:3617069/1‑71 (MQ=255)
                                                            agcgagTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGaa          <  1:1545164/71‑1 (MQ=255)
                                                             gcgagTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTc                         <  1:4345889/55‑1 (MQ=255)
                                                             gcgagTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAAagggtccctaag         >  1:1646799‑M1/1‑59 (MQ=255)
                                                              cgagTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTc                                 <  1:2447351/46‑1 (MQ=255)
                                                              cgagTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAAagggtccctaagt        <  1:1151462‑M1/71‑14 (MQ=255)
                                                              cgagTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAAag                   >  1:1621110‑M1/1‑58 (MQ=255)
                                                               gagTTCCGTGCCGGTTGTGAAGAAAAAGTGAatgat                                          <  1:1028980/36‑1 (MQ=255)
                                                                agTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGcc                                    >  1:1500478/1‑41 (MQ=255)
                                                                agTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGc                                     >  1:526623/1‑40 (MQ=255)
                                                                 gTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGa     >  1:2144731/1‑71 (MQ=255)
                                                                 gTTCCGTGCCGGTTGTGAAGAAAAAGTGAatgatg                                         <  1:3141213/35‑1 (MQ=255)
                                                                 gTTCCGTGCCGGTTGTGAAGAAAAGTGAATGATGTAGCCGTCAAGTTGTCATa                      <  1:672348/53‑1 (MQ=255)
                                                                  tttCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCaa                               <  1:2881378/44‑1 (MQ=255)
                                                                  tttCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTc                         <  1:1048387/50‑1 (MQ=255)
                                                                    tCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACaa  >  1:4095204/1‑71 (MQ=255)
                                                                    tCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACaa  >  1:1412089/1‑71 (MQ=255)
                                                                    |                                                                      
GGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAA  >  AP009048/71146‑71284

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: