Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,811,208 T→C V92A (GTT→GCT)  tcyP → cystine/sulfocysteine:cation symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,811,2080TC98.5% 759.2 / 0.8 197V92A (GTT→GCT) tcyPcystine/sulfocysteine:cation symporter
Reads supporting (aligned to +/- strand):  ref base T (0/3);  new base C (103/91);  total (103/94)
Fisher's exact test for biased strand distribution p-value = 1.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.94e-01

Only 100 of 199 total aligned reads displayed.

CTGCAAATGATCGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGTTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCAGTTTTCTGAC  >  NC_000913/1811159‑1811256
                                                 |                                                
  gcAAATGATCGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTg                                                >  2:4233750/1‑50 (MQ=255)
  gcAAATGATCGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTg                                                >  1:1670373/1‑50 (MQ=255)
   cAAATGATCGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGc                                               <  2:1743250/50‑1 (MQ=255)
   cAAATGATCGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGc                                               >  1:4237810/1‑50 (MQ=255)
   cAAATGATCGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGc                                               <  1:2666412/50‑1 (MQ=255)
   cAAATGATCGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGc                                               >  2:5815084/1‑50 (MQ=255)
    aaaTGATCGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGcc                                              >  1:97274/1‑50 (MQ=255)
     aaTGATCGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGccc                                             <  2:1514059/50‑1 (MQ=255)
      aTGATCGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCg                                            <  2:4928355/50‑1 (MQ=255)
       tGATCGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGt                                           >  2:714827/1‑50 (MQ=255)
        gATCGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTc                                          >  1:2929076/1‑50 (MQ=255)
        gATCGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTc                                          >  2:182896/1‑50 (MQ=255)
         aTCGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCt                                         >  1:5265545/1‑50 (MQ=255)
         aTCGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCt                                         <  2:108452/50‑1 (MQ=255)
         aTCGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCt                                         <  2:4044515/50‑1 (MQ=255)
          tCGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTg                                        <  2:3715903/50‑1 (MQ=255)
             ttATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCa                                      >  2:171956/1‑49 (MQ=255)
             ttATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCa                                      >  1:3426566/1‑49 (MQ=255)
             ttATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCAt                                     <  2:3200717/50‑1 (MQ=255)
             ttATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCAt                                     <  1:5078961/50‑1 (MQ=255)
             ttATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCAt                                     <  2:402503/50‑1 (MQ=255)
             ttATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCAt                                     <  2:1421838/50‑1 (MQ=255)
               aTGCCGTTTGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATaa                                   >  2:3293324/1‑50 (MQ=39)
               aTGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATaa                                   >  1:3133867/1‑50 (MQ=255)
               aTGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATaa                                   >  1:374601/1‑50 (MQ=255)
               aTGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATaa                                   >  1:5192139/1‑50 (MQ=255)
               aTGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATaa                                   >  1:2316003/1‑50 (MQ=255)
               aTGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATaa                                   >  2:5099646/1‑50 (MQ=255)
               aTGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATa                                    <  1:1246059/49‑1 (MQ=255)
               aTGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATa                                    <  1:4851177/49‑1 (MQ=255)
               aTGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATa                                    <  1:3285304/49‑1 (MQ=255)
                tGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATaa                                   <  1:3193780/49‑1 (MQ=255)
                 gCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATAACg                                 >  2:3813024/1‑50 (MQ=255)
                  ccGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATAACg                                 >  2:5357268/1‑49 (MQ=255)
                  ccGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATAACGc                                <  1:2697251/50‑1 (MQ=255)
                   cGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATAACGCa                               >  2:4190085/1‑50 (MQ=255)
                   cGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATAACGCa                               <  2:2231486/50‑1 (MQ=255)
                   cGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATAACGCa                               >  1:2232949/1‑50 (MQ=255)
                   cGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATAACGCa                               >  2:2775182/1‑50 (MQ=255)
                    gTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATAACGCAt                              >  2:5240137/1‑50 (MQ=255)
                    gTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATAACGCAt                              <  2:2318144/50‑1 (MQ=255)
                    gTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATAACGCAt                              >  2:2481027/1‑50 (MQ=255)
                     ttAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATAACGCAtc                             >  2:5273469/1‑50 (MQ=255)
                      tAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATAACGCAtct                            >  1:2205724/1‑50 (MQ=255)
                       aGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATAACGCAtct                            <  2:3257988/49‑1 (MQ=255)
                       aGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATAACGCAtct                            >  2:2541694/1‑49 (MQ=255)
                       aGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATAACGCAtct                            <  1:3753815/49‑1 (MQ=255)
                       aGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATAACGCAtct                            <  2:814410/49‑1 (MQ=255)
                        gTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATAACGCAtct                            >  1:3347493/1‑48 (MQ=255)
                        gTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATAACGCATCTCa                          <  1:5512918/50‑1 (MQ=255)
                        gTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATAACGCATCTCa                          <  1:1205496/50‑1 (MQ=255)
                         tCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATAACGCATCTCAg                         <  1:5845614/50‑1 (MQ=255)
                         tCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATAACGCATCTCAg                         >  1:4476229/1‑50 (MQ=255)
                          cTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATAACGCATCTCAGt                        >  1:394897/1‑50 (MQ=255)
                           ttCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATAACGCATCTCAGt                        >  2:4740245/1‑49 (MQ=255)
                           ttCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATAACGCATCTCAGt                        >  2:594881/1‑49 (MQ=255)
                             cGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATAACGCATCTCAGtt                       <  1:4583910/48‑1 (MQ=255)
                             cGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATAACGCATCTCAGtt                       <  1:726010/48‑1 (MQ=255)
                             cGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATAACGCATCTCAGt                        <  1:4943102/47‑1 (MQ=255)
                             cGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATAACGCATCTCAGt                        <  2:3021745/47‑1 (MQ=255)
                             cGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATAACGCATCTCAGTTa                      <  1:4376490/49‑1 (MQ=255)
                             cGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATAACGCATCTCAGTTAg                     >  2:5452345/1‑50 (MQ=255)
                             cGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATAACGCATCTCAGTTAg                     >  1:1992470/1‑50 (MQ=255)
                             cGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATAACGCATCTCAGTTAg                     >  1:4882340/1‑50 (MQ=255)
                               ccTCTATTCTGAGCGCGGCTGCCCGTCTGCATAACGCATCTCAGTTAgg                    >  1:3779352/1‑49 (MQ=255)
                                ctctATTCTGAGCGCGGCTGCCCGTCTGCATAACGCATCTCAGTTAGGc                   >  1:604993/1‑49 (MQ=255)
                                ctctATTCTGAGCGCGGCTGCCCGTCTGCATAACGCATCTCAGTTAGGc                   >  1:3457265/1‑49 (MQ=255)
                                ctctATTCTGAGCGCGGCTGCCCGTCTGCATAACGCATCTCAGTTAGGc                   >  2:5323362/1‑49 (MQ=255)
                                ctctATTCTGAGCGCGGCTGCCCGTCTGCATAACGCATCTCAGTTAGGCa                  <  2:4821635/50‑1 (MQ=255)
                                  ctATTCTGAGCGCGGCTGCCCGTCTGCATAACGCATCTCAGTTAGGCaa                 >  2:2267107/1‑49 (MQ=255)
                                   tATTCTGAGCGCGGCTGCCCGTCTGCATAACGCATCTCAGTTAGGCaaaa               >  1:67531/1‑50 (MQ=255)
                                   tATTCTGAGCGCGGCTGCCCGTCTGCATAACGCATCTCAGTTAGGCaaaa               >  1:149013/1‑50 (MQ=255)
                                   tATTCTGAGCGCGGCTGCCCGTCTGCATAACGCATCTCAGTTAGGCaa                 <  2:5638529/48‑1 (MQ=255)
                                   tATTCTGAGCGCGGCTGCCCGTCTGCATAACGCATCTCAGTTAGGCaa                 <  2:1295927/48‑1 (MQ=255)
                                    aTTCTGAGCGCGGCTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAAt              >  1:1931656/1‑50 (MQ=255)
                                     ttCTGAGCGCGGCTGCCCGTCTGCATAACGCATCTCAGTTAGGCaaaa               <  1:649305/48‑1 (MQ=255)
                                     ttCTGAGCGCGGCTGCCCGTCTGCATAACGCATCTCAGTTAGGCaaaa               <  1:3100101/48‑1 (MQ=255)
                                     ttCTGAGCGCGGCTGCCCGTCTGCATAACGCATCTCAGTTAGGCaaaa               <  2:1768417/48‑1 (MQ=255)
                                     ttCTGAGCGCGGCTGCCCGTCTGCATAACGCATCTCAGTTAGGCaaa                <  1:3826186/47‑1 (MQ=255)
                                     ttCTGAGCGCGGCTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAAt              <  1:4773710/49‑1 (MQ=255)
                                     ttCTGAGCGCGGCTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAAt              >  1:2773858/1‑49 (MQ=255)
                                     ttCTGAGCGCGGCTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAAt              <  1:1302846/49‑1 (MQ=255)
                                     ttCTGAGCGCGGCTGCCCGTCTGCATAACGCAGCTCAGTTAGGCAAAATc             <  2:1535262/50‑1 (MQ=255)
                                       cTGAGCGCGGTTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCa            <  2:5826713/49‑1 (MQ=255)
                                       cTGAGCGCGGCTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCa            <  2:4637206/49‑1 (MQ=255)
                                       cTGAGCGCGGCTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCa            <  1:2852020/49‑1 (MQ=255)
                                       cTGAGCGCGGCTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCa            <  1:26027/49‑1 (MQ=255)
                                        tGAGCGCGGCTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCAGt          >  2:351659/1‑50 (MQ=255)
                                         gAGCGCGGCTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCAGtt         >  1:2102959/1‑50 (MQ=255)
                                         gAGCGCGGCTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCAGtt         >  2:833379/1‑50 (MQ=255)
                                          aGCGCGGCTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCAGttt        >  1:2447963/1‑50 (MQ=255)
                                          aGCGCGGCTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCAGttt        <  2:1928885/50‑1 (MQ=255)
                                          aGCGCGGCTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCAGttt        <  2:653432/50‑1 (MQ=255)
                                          aGCGCGGCTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCAGttt        <  1:5201792/50‑1 (MQ=255)
                                          aGCGCGGCTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCAGtt         <  2:20507/49‑1 (MQ=255)
                                           gcgcGGCTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCAGtttt       <  1:5655286/50‑1 (MQ=255)
                                             gcgGCTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCAGTTTTCt     <  2:3610467/50‑1 (MQ=255)
                                             gcgGCTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCAGTTTTCt     >  1:4099601/1‑50 (MQ=255)
                                              cgGCTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCAGTTTTCt     >  2:5489401/1‑49 (MQ=255)
                                              cgGCTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCAGTTTTCTg    <  2:2579532/50‑1 (MQ=255)
                                                 |                                                
CTGCAAATGATCGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGTTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCAGTTTTCTGAC  >  NC_000913/1811159‑1811256

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: