Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 70,740 T→G 20.4% G118G (GGT→GGG araC → DNA‑binding transcriptional dual regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP00904870,7400TG20.4% 80.9 / 35.1 59G118G (GGT→GGGaraCDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base T (22/25);  new base G (6/6);  total (28/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.20e-01

TTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAACGCC  >  AP009048/70673‑70806
                                                                   |                                                                  
tttACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGtttc                                                                 <  1:1094871/71‑1 (MQ=255)
  tACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGtttctt                                                               >  1:943840/1‑71 (MQ=255)
  tACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGtttctt                                                               >  1:2666879/1‑71 (MQ=255)
  tACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGgttct                                                                <  1:2236746/70‑1 (MQ=255)
   aCTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGgttcttt                                                              >  1:4295748/1‑71 (MQ=255)
     tttCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGtttcttt                                                              <  1:1830874/69‑1 (MQ=255)
      ttCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGt                                                                    <  1:4484649/62‑1 (MQ=255)
      ttCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGt                                                                    <  1:1692918/62‑1 (MQ=255)
       tCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCg                                                            >  1:1358566/1‑69 (MQ=255)
       tCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGcc                                                          >  1:4320031/1‑71 (MQ=255)
       tCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGcc                                                          >  1:1793602/1‑71 (MQ=255)
       tCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGcc                                                          >  1:4757972/1‑71 (MQ=255)
       tCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGcc                                                          >  1:3083280/1‑71 (MQ=255)
        cGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGccc                                                         >  1:348627/1‑71 (MQ=255)
            cgcgcgCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGtttct                                                                <  1:2107922/60‑1 (MQ=255)
            cgcgcgCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGa                                                      <  1:1388628/70‑1 (MQ=255)
             gcgcgcCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATg                                                    <  1:3845053/71‑1 (MQ=255)
             gcgcgcCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATg                                                    <  1:4106264/71‑1 (MQ=255)
               gcgcCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATg                                                    >  1:4875297/1‑69 (MQ=255)
               gcgcCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGa                                                   <  1:2882129/70‑1 (MQ=255)
                 gcCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGtttct                                                                <  1:4545187/55‑1 (MQ=255)
                 gcCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAgc                                                >  1:405441/1‑71 (MQ=255)
                 gcCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAgc                                                <  1:4187779/71‑1 (MQ=255)
                 gcCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAg                                                 <  1:4733366/70‑1 (MQ=255)
                   cTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAgcgc                                              <  1:4351299/71‑1 (MQ=255)
                      cTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGt                                                                    >  1:2003815/1‑46 (MQ=255)
                      cTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGt                                                                    >  1:4461804/1‑46 (MQ=255)
                      cTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAgcgc                                              <  1:779543/68‑1 (MQ=255)
                       tGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGtttct                                                                <  1:4419985/49‑1 (MQ=255)
                         gCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGtttct                                                                >  1:9109/1‑47 (MQ=255)
                         gCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGaa                                                  <  1:3879029/61‑1 (MQ=255)
                         gCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCa                                          <  1:3665051/69‑1 (MQ=255)
                          cATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAg                                         >  1:3137125/1‑69 (MQ=255)
                           aTGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGtttct                                                                <  1:3848456/45‑1 (MQ=255)
                               aTGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAgcg                                               <  1:2915913/58‑1 (MQ=255)
                               aTGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCAtt                                  >  1:1336791/1‑71 (MQ=255)
                                tGGCTTATCTGGCCGTCAATATTTGCGAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCAttt                                 <  1:3677752/71‑1 (MQ=255)
                                tGGCTTAACTGGCCGTCAATATTTGCCAATACGGGgttctttct                                                            >  1:190117/1‑43 (MQ=255)
                                tGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCAttt                                 <  1:4101649/71‑1 (MQ=255)
                                  gCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCAc                                            >  1:5101282/1‑58 (MQ=255)
                                  gCTTAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGa                                                      >  1:4340819/1‑48 (MQ=255)
                                     tAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCa                               >  1:1005434/1‑68 (MQ=255)
                                      aaCTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCg                            >  1:3677726/1‑70 (MQ=255)
                                       aCTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGAc                          <  1:1546665/71‑1 (MQ=255)
                                       aCTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGAc                          >  1:1629556/1‑71 (MQ=255)
                                        cTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGc                             <  1:3557586/67‑1 (MQ=255)
                                        cTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGAcc                         >  1:3204386/1‑71 (MQ=255)
                                           gCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGt                      >  1:459224/1‑71 (MQ=255)
                                           gCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGt                      <  1:2475554/71‑1 (MQ=255)
                                           gCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGt                      >  1:4583343/1‑71 (MQ=255)
                                            ccGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGtt                     <  1:2087507/71‑1 (MQ=255)
                                             cGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGttt                    <  1:231834/71‑1 (MQ=255)
                                             cGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGttt                    <  1:3041848/71‑1 (MQ=255)
                                             cGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGttt                    <  1:4178955/71‑1 (MQ=255)
                                             cGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGttt                    <  1:4462148/71‑1 (MQ=255)
                                              gTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGAt                                                     <  1:4771188/37‑1 (MQ=255)
                                              gTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTg                   >  1:13040/1‑71 (MQ=255)
                                               tCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTgg                  <  1:5147075/71‑1 (MQ=255)
                                                 aaTATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTggg                 >  1:3210472/1‑70 (MQ=255)
                                                 aaTATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTgg                  >  1:3772194/1‑69 (MQ=255)
                                                  atatTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCAc                                            <  1:3990562/42‑1 (MQ=255)
                                                   tatTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGtt                     <  1:372991/64‑1 (MQ=255)
                                                   tatTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCa               >  1:2214464/1‑70 (MQ=255)
                                                       tGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATc           <  1:3772208/70‑1 (MQ=255)
                                                       tGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCa          >  1:5005015/1‑71 (MQ=255)
                                                       tGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCa          >  1:1063051/1‑71 (MQ=255)
                                                          cAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCAtt        <  1:1042933/70‑1 (MQ=255)
                                                          cAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTa       <  1:865112/71‑1 (MQ=255)
                                                            aTACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCAttt                                 <  1:242873/43‑1 (MQ=255)
                                                            aTACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGtt                     >  1:4664303/1‑55 (MQ=255)
                                                            aTACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAAc     >  1:2798433/1‑71 (MQ=255)
                                                            aTACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGtt                     <  1:414184/55‑1 (MQ=255)
                                                              aCGGGTTTCTTTCGCCCGGATGAAGCGCACCAGcc                                       >  1:2520072/1‑35 (MQ=255)
                                                              aCGGGTTTCTTTCGCCCGGATGAAGCGCACCAGcc                                       >  1:1206057/1‑35 (MQ=255)
                                                               cGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAACGcc  <  1:33778/71‑1 (MQ=255)
                                                                   |                                                                  
TTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAACGCC  >  AP009048/70673‑70806

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: