Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 70,581 C→T 44.6% 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,5810CT44.6% 28.5 / 96.4 74V65V (GTC→GTTaraCDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base C (22/19);  new base T (17/16);  total (39/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.17e-01

TGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCT  >  AP009048/70511‑70650
                                                                      |                                                                     
tGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTc                                                                       >  1:3469461/1‑71 (MQ=255)
 gAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGTGGTGAAAAATCAGGGACGAGAAtttgtttg                                                                     >  1:3084593/1‑71 (MQ=255)
    aGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGcc                                                                   >  1:2799646/1‑71 (MQ=255)
    aGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGcc                                                                   >  1:4337396/1‑71 (MQ=255)
         atatTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGAcc                                                               <  1:2937488/70‑1 (MQ=255)
            ttCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGAcc                                                               <  1:242369/67‑1 (MQ=255)
             tctcAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGt                                                           >  1:2593217/1‑70 (MQ=255)
               tcAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTc                                                                       >  1:3156742/1‑56 (MQ=255)
               tcAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGc                                                                    >  1:4112899/1‑59 (MQ=255)
               tcAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGc                                                                    >  1:1184915/1‑59 (MQ=255)
               tcAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGa                                                         >  1:397998/1‑70 (MQ=255)
                cAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGata                                                       >  1:868576/1‑71 (MQ=255)
                 aaTCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGata                                                       <  1:3834933/70‑1 (MQ=255)
                 aaTCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGata                                                       <  1:4471382/70‑1 (MQ=255)
                  aTCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATAtt                                                     <  1:4512820/71‑1 (MQ=255)
                  aTCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATAtt                                                     <  1:1991548/71‑1 (MQ=255)
                  aTCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATAtt                                                     <  1:755675/71‑1 (MQ=255)
                      cACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGc                                                 >  1:1550968/1‑71 (MQ=255)
                       aCCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATAttt                                                    <  1:117282/67‑1 (MQ=255)
                       aCCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATAtttt                                                   <  1:977001/68‑1 (MQ=255)
                       aCCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATAttt                                                    <  1:3573905/67‑1 (MQ=255)
                       aCCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTg                                                  <  1:1137138/69‑1 (MQ=255)
                        ccATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTg                                               >  1:770170/1‑71 (MQ=255)
                           ttCGCGTTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTAGCCGACCGGGTGATAttt                                                    <  1:3186638/63‑1 (MQ=255)
                           ttCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGtt                                             <  1:2001120/70‑1 (MQ=255)
                              gcgGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCg                                         >  1:926372/1‑71 (MQ=255)
                                  tCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGcc                                       <  1:4519573/69‑1 (MQ=255)
                                  tCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAg                                     <  1:4386767/71‑1 (MQ=255)
                                   cAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAg                                     >  1:1597527/1‑70 (MQ=255)
                                    aGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTg                                                          <  1:2504753/48‑1 (MQ=255)
                                    aGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGga                                   <  1:516455/71‑1 (MQ=255)
                                    aGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGAcc                                                               >  1:3988227/1‑43 (MQ=255)
                                    aGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGa                                                         >  1:3109409/1‑49 (MQ=255)
                                    aGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATAttt                                                    <  1:3344158/54‑1 (MQ=255)
                                     gggggTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGtt                                             <  1:2978233/60‑1 (MQ=255)
                                     gggggTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAg                                     >  1:4384859/1‑68 (MQ=255)
                                      ggggTGGTGAAAAATCAGGGACGAGATTTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAgg                                    >  1:3679070/1‑68 (MQ=255)
                                      ggggTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCa                                      <  1:1611938/66‑1 (MQ=255)
                                        ggtggtGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGata                                                       <  1:2759234/47‑1 (MQ=255)
                                          tggtgAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGc                                        <  1:57839/60‑1 (MQ=255)
                                          tggtgAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGgaga                                 <  1:2905905/67‑1 (MQ=255)
                                          tggtgAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGtt                                             <  1:2877828/55‑1 (MQ=255)
                                            gtgAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATAtt                                                     >  1:260340/1‑45 (MQ=255)
                                             ttaaaaaTCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGgagag                                <  1:3867463/63‑1 (MQ=255)
                                             tgAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATtcat                          <  1:3478492/71‑1 (MQ=255)
                                              gAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGt                                                           >  1:2399970/1‑37 (MQ=255)
                                              gAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGga                                   >  1:2505309/1‑61 (MQ=255)
                                               aaaaaTCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATtcatca                        >  1:1375464/1‑71 (MQ=255)
                                               aaaaaTCAGGGACGAGAATTTGTCTGCCGACCGGGTGATAtttt                                                   >  1:770999/1‑44 (MQ=255)
                                               aaaaaTCAGGGACGAGAATTTGTCTGCCGACCGGGTGATAtt                                                     >  1:3616409/1‑42 (MQ=255)
                                               aaaaaTCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTTCCGCCAGGAGAGATtcatca                        >  1:491962/1‑71 (MQ=255)
                                               aaaaaTCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATtcatca                        >  1:3571898/1‑71 (MQ=255)
                                               aaaaaTCAGGGACGAGAATTTCTCTGCCGACCGGGTGATATTTTg                                                  <  1:4293974/45‑1 (MQ=255)
                                                 aaaTCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATtcat                          >  1:4182739/1‑67 (MQ=255)
                                                   aTCAGGGACGAGAATTTGTCTGCCGACCGGGTGATAtt                                                     >  1:911031/1‑38 (MQ=255)
                                                    tCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTa                     >  1:1180546/1‑69 (MQ=255)
                                                    tCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACg                   >  1:961243/1‑71 (MQ=255)
                                                    tCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGACATTCATCACTACg                   <  1:1784185/71‑1 (MQ=255)
                                                    tCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACg                   <  1:98175/71‑1 (MQ=255)
                                                      aGGTACGAGAATTTGTTTGCCGACCGGGTGATAttt                                                    >  1:4008769/1‑36 (MQ=37)
                                                      aGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGgt                 >  1:1163611/1‑71 (MQ=255)
                                                      aGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGAtt                             >  1:1060234/1‑59 (MQ=255)
                                                      aGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGAtt                             >  1:1810075/1‑59 (MQ=255)
                                                        ggACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGgtcg               >  1:4101785/1‑71 (MQ=255)
                                                         gACGAGAATTTGTTTGCCGACCGGGTGATATTTTg                                                  <  1:3948593/35‑1 (MQ=255)
                                                         gACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTccc                                          >  1:2385223/1‑43 (MQ=255)
                                                         gACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGgtcgt              <  1:1819016/71‑1 (MQ=255)
                                                           cGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATtc                            <  1:3148290/55‑1 (MQ=255)
                                                           cGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCa            >  1:3750529/1‑71 (MQ=255)
                                                            gagaATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCAt           >  1:3930979/1‑71 (MQ=255)
                                                            gagaATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCAt           <  1:1409514/71‑1 (MQ=255)
                                                               aaTTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATtc                            <  1:4410197/51‑1 (MQ=255)
                                                               aaTTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACgg                  >  1:366747/1‑61 (MQ=255)
                                                                aTTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCg        >  1:1959867/1‑70 (MQ=255)
                                                                 tttgtttgCCGACCGGGTGATATTTTGCTGTTCCCGCCAGgag                                  >  1:4196500/1‑43 (MQ=255)
                                                                 tttgtttgCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCg        >  1:2366369/1‑69 (MQ=255)
                                                                   tGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTa                     >  1:741201/1‑54 (MQ=255)
                                                                    gTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGga                                   >  1:2013148/1‑39 (MQ=255)
                                                                     tCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCt  <  1:4098528/71‑1 (MQ=255)
                                                                      |                                                                     
TGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCT  >  AP009048/70511‑70650

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: