Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,232,364 C→A 10.2% G131V (GGA→GTA)  ptrB ← protease II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,232,3640CA10.2% 220.4 / 7.9 68G131V (GGA→GTA) ptrBprotease II
Reads supporting (aligned to +/- strand):  ref base C (31/30);  new base A (3/4);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.48e-01

CCGTACTGGCGTCGGGAAAGAAAATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTGGCATCGAGCAATGTTTCCCACTCATCCCA  >  Exported/1232295‑1232434
                                                                     |                                                                      
ccGTACTGGCGTCGGGAAAGAAAATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTcc                                                                       <  1:3652393/71‑1 (MQ=255)
ccGTACTGGCGTCGGGAAAGAAAATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTcc                                                                       <  1:1742057/71‑1 (MQ=255)
  gTACTGGCGTCGGGAAAGAAAATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGc                                                                     >  1:242790/1‑71 (MQ=255)
  gTACTGGCGTCGGGAAAGAAAATCTTCTGCCAGCGCCATAAAGGTGTTATCGGGCGTAATCGCCATTCCGc                                                                     >  1:2515542/1‑71 (MQ=255)
             gggAAAGAAAATCTTCTGGCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATaa                                                          <  1:2173542/71‑1 (MQ=255)
             gggAAAGAAAATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATaa                                                          <  1:1189974/71‑1 (MQ=255)
               gaaagaaaATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATaaaa                                                        <  1:1435691/71‑1 (MQ=255)
                  agaaaATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGc                                                                     >  1:1766811/1‑55 (MQ=255)
                   gaaaATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCa                                                    >  1:546977/1‑71 (MQ=255)
                   gaaaATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCa                                                    <  1:1310724/71‑1 (MQ=255)
                   gaaaATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCa                                                    <  1:2323041/71‑1 (MQ=255)
                   gaaaATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCa                                                    >  1:1637556/1‑71 (MQ=255)
                    aaaaTCTTCTGCCAGCGACATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCAc                                                   >  1:2319817/1‑71 (MQ=255)
                     aaaTCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCAc                                                   >  1:2356153/1‑70 (MQ=255)
                     aaaTCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCAc                                                   <  1:1746360/70‑1 (MQ=255)
                     aaaTCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACt                                                  <  1:1673387/71‑1 (MQ=255)
                     aaaTCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACt                                                  <  1:14212/71‑1 (MQ=255)
                        tcttctGCCAGCGCCATAATGGTGTTATCGGGCGGAATGGCGATTACGccc                                                                   >  1:2025436/1‑51 (MQ=21)
                        tcttctGCCAGCGCCATAATGGTGTTATCGGGCGGAATGGCGATTACGccc                                                                   >  1:3135235/1‑51 (MQ=21)
                        tcttctGCCAGCGCCATAATGGTGTTATCGGGCGGAATGGCGATTACGccc                                                                   >  1:1842334/1‑51 (MQ=21)
                         cttctGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTAt                                                >  1:107123/1‑69 (MQ=255)
                           tctGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGc                                            >  1:1840121/1‑71 (MQ=255)
                            ctGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATg                                               <  1:2257400/67‑1 (MQ=255)
                            ctGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATg                                               <  1:1772999/67‑1 (MQ=255)
                            ctGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATg                                               <  1:43657/67‑1 (MQ=255)
                                cAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTgcgc                                       <  1:1553584/71‑1 (MQ=255)
                                  gcgcCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCa                                                    <  1:3221466/56‑1 (MQ=255)
                                   cgcCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCt                                    >  1:329111/1‑71 (MQ=255)
                                   cgcCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCt                                    >  1:355919/1‑71 (MQ=255)
                                    gcCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCtt                                   <  1:2452894/71‑1 (MQ=255)
                                    gcCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCtt                                   <  1:3243234/71‑1 (MQ=255)
                                     ccATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTgcgcgc                                     <  1:2305181/68‑1 (MQ=255)
                                     ccATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTa                                  >  1:2349938/1‑71 (MQ=255)
                                      cATAATGGTGTTATCGGGCGTAATCGCCATTCCGCACATCGAATAAAACTCACTATGAGCTGCGCGCTTAt                                 <  1:3102388/71‑1 (MQ=255)
                                       aTAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTAtt                                >  1:664892/1‑71 (MQ=255)
                                        tAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTa                                                 <  1:957888/53‑1 (MQ=255)
                                        tAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTg                               >  1:3282933/1‑71 (MQ=255)
                                        tAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTg                               >  1:3287211/1‑71 (MQ=255)
                                        tAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTg                               >  1:1004755/1‑71 (MQ=255)
                                        tAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTg                               >  1:1857135/1‑71 (MQ=255)
                                        tAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTg                               >  1:1891808/1‑71 (MQ=255)
                                         aaTGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGGGCGCTTATTgg                              >  1:2246104/1‑71 (MQ=255)
                                         aaTGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTgg                              >  1:472771/1‑71 (MQ=255)
                                         aaTGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTgg                              >  1:3656216/1‑71 (MQ=255)
                                         aaTGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTgg                              >  1:3619948/1‑71 (MQ=255)
                                         aaTGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTgg                              >  1:2117673/1‑71 (MQ=255)
                                         aaTGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTgg                              >  1:3205421/1‑71 (MQ=255)
                                                  aTCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTGGCATCGa                        >  1:1101784/1‑68 (MQ=255)
                                                   tCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTGGCATCGAGCaa                    <  1:2859813/71‑1 (MQ=255)
                                                   tCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTGGCATCGAGCaa                    <  1:1338074/71‑1 (MQ=255)
                                                    cGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTGGCATCGAGCaa                    <  1:470352/70‑1 (MQ=255)
                                                     gggCGGAATGGCGATTACGCCCATCGAATAAAACTCACTATGAGc                                            <  1:1042718/45‑1 (MQ=18)
                                                     gggCGGAATGGCGATTACGCCCATCGAATAAAACTCACTATGAGc                                            <  1:113402/45‑1 (MQ=18)
                                                     gggCGGAATGGCGATTACGCCCATCGAATAAAACTCACTATGAGc                                            <  1:3324596/45‑1 (MQ=18)
                                                     gggCGGAATGGCGATTACGCCCATCGAATAAAACTCACTATGAGc                                            <  1:1601209/45‑1 (MQ=18)
                                                        cGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGa                                              <  1:1999930/40‑1 (MQ=255)
                                                         gTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTGGCATCGAGCAATGTTTc              <  1:2189975/71‑1 (MQ=255)
                                                         gTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTGGCATCGAGCAATGTTTc              >  1:1735470/1‑71 (MQ=255)
                                                          tAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGc                                            >  1:1805924/1‑40 (MQ=255)
                                                          tAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTGGCATCGAGCAATGTTTcc             <  1:2043933/71‑1 (MQ=255)
                                                          tAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTGGCATCGAGCAATGTTTcc             >  1:238452/1‑71 (MQ=255)
                                                             tCGCCATTCCGCCCATCGAATAAAACTCACTATGa                                              >  1:2978945/1‑35 (MQ=255)
                                                               gCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTGGCATCGAGc                      <  1:3155017/57‑1 (MQ=255)
                                                               gCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTGGCATCGAGc                      <  1:1708764/57‑1 (MQ=255)
                                                               gCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTGGCATCGAGCAATGTTTCCCACTc        <  1:2489568/71‑1 (MQ=255)
                                                                ccATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTGGCATCGAGCAATGTTTCCCACTc        >  1:1117393/1‑70 (MQ=255)
                                                                ccATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTGGCATCGAGCAATGTTTCCCACTCa       >  1:3390240/1‑71 (MQ=255)
                                                                   ttCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTa                                  <  1:1467919/41‑1 (MQ=255)
                                                                   ttCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTa                                  <  1:952355/41‑1 (MQ=255)
                                                                   ttCCGCCCATCGAATAAAACTCACTATGAGCTGCCCGCTTa                                  <  1:2792/41‑1 (MQ=255)
                                                                     ccGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTAtt                                >  1:601354/1‑41 (MQ=255)
                                                                     ccGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTGGCATCGAGCAAt                   >  1:2418162/1‑54 (MQ=255)
                                                                     ccGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTGGCATCGAGCAAt                   >  1:2815644/1‑54 (MQ=255)
                                                                     ccGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTGGCATCGAGCAAt                   >  1:3458154/1‑54 (MQ=255)
                                                                     ccGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTGGCATCGAGCAATGTTTCCCACTCATCCCa  >  1:2107557/1‑71 (MQ=255)
                                                                     |                                                                      
CCGTACTGGCGTCGGGAAAGAAAATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTGGCATCGAGCAATGTTTCCCACTCATCCCA  >  Exported/1232295‑1232434

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: