Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 70,581 C→T 29.3% V65V (GTC→GTT araC → DNA‑binding transcriptional dual regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP00904870,5810CT29.3% 117.2 / 58.8 75V65V (GTC→GTTaraCDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base C (33/20);  new base T (16/6);  total (49/26)
Fisher's exact test for biased strand distribution p-value = 4.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

TGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCT  >  AP009048/70511‑70650
                                                                      |                                                                     
tGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTc                                                                       <  1:2563541/71‑1 (MQ=255)
  aaaGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTg                                                                     >  1:3880503/1‑71 (MQ=255)
   aaGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGc                                                                    >  1:366112/1‑71 (MQ=255)
    aGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGcc                                                                   >  1:4250894/1‑71 (MQ=255)
    aGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGc                                                                    >  1:4311839/1‑70 (MQ=255)
     ggTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGcc                                                                   >  1:1278682/1‑70 (MQ=255)
     ggTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGcc                                                                   >  1:2884538/1‑70 (MQ=255)
     ggTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCg                                                                  <  1:529983/71‑1 (MQ=255)
      gTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGa                                                                 <  1:953937/71‑1 (MQ=255)
       ttATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCt                                                                      >  1:4065915/1‑65 (MQ=255)
         atatTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCg                                                              >  1:590614/1‑71 (MQ=255)
         atatTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGAcc                                                               <  1:3117049/70‑1 (MQ=255)
         atatTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGAc                                                                >  1:502602/1‑69 (MQ=255)
          tatTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCg                                                              <  1:4476863/70‑1 (MQ=255)
              ctcAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGa                                                         >  1:3745008/1‑71 (MQ=255)
                cAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGata                                                       >  1:4875706/1‑71 (MQ=255)
                 aaTCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGatat                                                      <  1:1048507/71‑1 (MQ=255)
                  aTCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGata                                                       >  1:397074/1‑69 (MQ=255)
                   tctcACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCg                                                                  <  1:4562010/57‑1 (MQ=255)
                      cACCATTCGCGGTCATGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGc                                                 >  1:2191455/1‑71 (MQ=255)
                      cACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGAcc                                                               <  1:4427440/57‑1 (MQ=255)
                       aCCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATAttt                                                    <  1:4108131/67‑1 (MQ=255)
                         cATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTg                                                  >  1:2095143/1‑67 (MQ=255)
                          aTTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCgg                                                             >  1:4954919/1‑55 (MQ=255)
                           ttCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGa                                                                 >  1:2254462/1‑50 (MQ=255)
                              gcgGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCg                                         >  1:3489366/1‑71 (MQ=255)
                                ggTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGc                                        >  1:262648/1‑70 (MQ=255)
                                ggTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGc                                        >  1:1615428/1‑70 (MQ=255)
                                   cAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAg                                     >  1:4865216/1‑70 (MQ=255)
                                    aGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAg                                     >  1:424527/1‑69 (MQ=255)
                                    aGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAgg                                    <  1:3655747/70‑1 (MQ=255)
                                    aGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAg                                     >  1:1983318/1‑69 (MQ=255)
                                     gggggTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGt                                                           <  1:581177/46‑1 (MQ=255)
                                      ggggTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGa                                                         >  1:3630727/1‑47 (MQ=255)
                                      ggggTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGAc                                                                <  1:3990789/40‑1 (MQ=255)
                                       gggTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCg                                                              >  1:2897657/1‑41 (MQ=255)
                                       gggTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTg                                                  >  1:3231651/1‑53 (MQ=255)
                                           ggtgAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGAtt                             >  1:3744594/1‑70 (MQ=255)
                                            gtgAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATtca                           >  1:3186627/1‑71 (MQ=255)
                                            gtgAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATtc                            <  1:2573053/70‑1 (MQ=255)
                                             tgAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATtcat                          <  1:1857965/71‑1 (MQ=255)
                                              gAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATtcat                          >  1:1117994/1‑70 (MQ=255)
                                               aaaaaTCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGgagaga                               >  1:1818623/1‑64 (MQ=255)
                                               aaaaaTCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATtcatca                        >  1:5074428/1‑71 (MQ=255)
                                                aaaaTCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCAc                       >  1:2481030/1‑71 (MQ=255)
                                                aaaaTCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATtcatca                        <  1:878414/70‑1 (MQ=255)
                                                aaaaTCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATtcatca                        >  1:3002723/1‑70 (MQ=255)
                                                 aaaTCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGtt                                             >  1:1217733/1‑48 (MQ=255)
                                                 aaaTCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCAc                       >  1:2922700/1‑70 (MQ=255)
                                                   aTCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGgaga                                 >  1:627644/1‑58 (MQ=255)
                                                   aTCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTAc                    >  1:4153671/1‑71 (MQ=255)
                                                     cAGGGACGAGAATTTTTCTGCCGACCGGGTGATATTTTg                                                  >  1:1318595/1‑39 (MQ=255)
                                                     cAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCt                                                >  1:1972826/1‑41 (MQ=255)
                                                     cAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACgg                  >  1:3749215/1‑71 (MQ=255)
                                                      aGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTa                     >  1:5277844/1‑67 (MQ=255)
                                                       gggACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGgtc                <  1:3327281/71‑1 (MQ=255)
                                                           cGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTAATCACTACGGTCGTCa            >  1:1506489/1‑71 (MQ=255)
                                                            gagaATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGgtcgtc             >  1:3080912/1‑69 (MQ=255)
                                                            gagaATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGgagag                                <  1:544507/50‑1 (MQ=255)
                                                            gagaATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCAt           >  1:3646145/1‑71 (MQ=255)
                                                            gagaATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCAt           <  1:1591208/71‑1 (MQ=255)
                                                             agaATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACg                   <  1:1241039/62‑1 (MQ=255)
                                                             agaATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATc          <  1:2886835/71‑1 (MQ=255)
                                                              gaATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATtc                            >  1:2315395/1‑52 (MQ=255)
                                                              gaATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGgtcgt              >  1:4131214/1‑66 (MQ=255)
                                                               aaTTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCg                                         <  1:566476/38‑1 (MQ=255)
                                                               aaTTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGga                                   <  1:140469/44‑1 (MQ=255)
                                                               aaTTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGgtc                >  1:2795038/1‑63 (MQ=255)
                                                                 tttGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCAt           <  1:1697628/66‑1 (MQ=255)
                                                                 tttGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCg        >  1:5043854/1‑69 (MQ=255)
                                                                  ttGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGgagaga                               <  1:3589843/45‑1 (MQ=255)
                                                                  ttGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGgagaga                               <  1:2903514/45‑1 (MQ=255)
                                                                   tGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGgtcgt              >  1:2428617/1‑61 (MQ=255)
                                                                    gTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAgg    <  1:5012641/70‑1 (MQ=255)
                                                                    gTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGc   >  1:9150/1‑71 (MQ=255)
                                                                     tCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCt  <  1:2120143/71‑1 (MQ=255)
                                                                      |                                                                     
TGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCT  >  AP009048/70511‑70650

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: