Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,811,2080TC50.0% 2.1 / 115.5 70V92A (GTT→GCT) tcyPcystine/sulfocysteine:cation symporter
Reads supporting (aligned to +/- strand):  ref base T (12/23);  new base C (16/19);  total (28/42)
Fisher's exact test for biased strand distribution p-value = 4.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CTGCAAATGATCGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGTTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCAGTTTTCTGA  >  NC_000913/1811159‑1811255
                                                 |                                               
ctgcAAATGATCGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGt                                                 >  1:1407601/1‑50 (MQ=255)
  gcAAATGATCGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGTTg                                               >  1:3806812/1‑50 (MQ=255)
   cAAATGATCGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGtt                                                <  1:2504149/48‑1 (MQ=255)
   cAAATGATCGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGtt                                                <  1:3206720/48‑1 (MQ=255)
   cAAATGATCGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGtt                                                <  1:213845/48‑1 (MQ=255)
   cAAATGATCGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGtt                                                <  1:651701/48‑1 (MQ=255)
   cAAATGATCGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGTTGc                                              <  1:1405993/50‑1 (MQ=255)
   cAAATGATCGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGTTGc                                              <  1:3560001/50‑1 (MQ=255)
   cAAATGATCGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGc                                              >  1:196065/1‑50 (MQ=255)
   cAAATGATCGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGc                                              <  1:1632619/50‑1 (MQ=255)
    aaaTGATCGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGTTGcc                                             >  1:2692277/1‑50 (MQ=255)
     aaTGATCGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGTTGccc                                            <  1:2178266/50‑1 (MQ=255)
     aaTGATCGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGTTGccc                                            >  1:2094785/1‑50 (MQ=255)
     aaTGATCGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGccc                                            <  1:3129325/50‑1 (MQ=255)
     aaTGATCGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGccc                                            <  1:1777657/50‑1 (MQ=255)
      aTGATCGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGTTGCCCg                                           <  1:3275510/50‑1 (MQ=255)
      aTGATCGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCg                                           >  1:2018717/1‑50 (MQ=255)
      aTGATCGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCg                                           <  1:2791349/50‑1 (MQ=255)
        gATCGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTc                                         >  1:4300922/1‑50 (MQ=255)
        gATCGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTc                                         >  1:2135987/1‑50 (MQ=255)
         aTCGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGTTGCCCGTCt                                        <  1:4055847/50‑1 (MQ=255)
         aTCGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTc                                         <  1:311304/49‑1 (MQ=255)
         aTCGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCt                                        <  1:2162094/50‑1 (MQ=255)
         aTCGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCt                                        >  1:1772069/1‑50 (MQ=255)
          tCGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTg                                       <  1:1701252/50‑1 (MQ=255)
           cGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTg                                       >  1:1818260/1‑49 (MQ=255)
            gTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGTTGCCCGTCTGc                                      <  1:769005/49‑1 (MQ=255)
            gTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGTTGCCCGTCTGCa                                     >  1:2740376/1‑50 (MQ=255)
            gTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGTTGCCCGTCTGCa                                     >  1:683997/1‑50 (MQ=255)
            gTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCa                                     >  1:4046512/1‑50 (MQ=255)
            gTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCa                                     >  1:1475655/1‑50 (MQ=255)
            gTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCa                                     >  1:3788916/1‑50 (MQ=255)
               aTGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATaa                                  >  1:2290907/1‑50 (MQ=255)
                tGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGTTGCCCGTCTGCATaa                                  <  1:772938/49‑1 (MQ=255)
                tGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATaa                                  <  1:1210097/49‑1 (MQ=255)
                tGCCGCTAGTCTTCGCCTCTATGCTGAGCGCGGTTGCCCGTCTGCATaa                                  <  1:3245501/49‑1 (MQ=255)
                  ccGTTAGTCTTCGCCTCTATTCTGAGCGCGGTTGCCCGTCTGCATAACg                                <  1:3727936/49‑1 (MQ=255)
                   cGTTAGTCTTCGCCTCTATTCTGAGCGCGGTTGCCCGTCTGCATAACGCa                              >  1:2608209/1‑50 (MQ=255)
                   cGTTAGTCTTCGCCTCTATTCTGAGCGCGGTTGCCCGTCTGCATAACGCa                              >  1:2218706/1‑50 (MQ=255)
                    gTTAGTCTTCGCCTCTATTCTGAGCGCGGTTGCCCGTCTGCATAACGCAt                             >  1:1504321/1‑50 (MQ=255)
                     ttAGTCTTCGCCTCTATTCTGAGCGCGGTTGCCCGTCTGCATAACGCAtc                            <  1:819562/50‑1 (MQ=255)
                     ttAGTCTTCGCCTCTATTCTGAGCGCGGTTGCCCGTCTGCATAACGCAtc                            <  1:123028/50‑1 (MQ=255)
                       aGTCTTCGCCTCTATTCTGAGCGCGGTTGCCCGTCTGCATAACGCAtct                           <  1:2948405/49‑1 (MQ=255)
                       aGTCTTCGCCTCTATTCTGAGCGCGGTTGCCCGTCTGCATAACGCAtct                           <  1:3838455/49‑1 (MQ=255)
                       aGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATAACGCAtct                           <  1:66782/49‑1 (MQ=255)
                        gTCTTCGCCTCTATTCTGAGCGCGGTTGCCCGTCTGCATAACGCAtctc                          <  1:325568/49‑1 (MQ=255)
                        gTCTTCGCCTCTATTCTGAGCGCGGTTGCCCGTCTGCATAACGCAtctc                          <  1:3963480/49‑1 (MQ=255)
                        gTCTTCGCCTCTATTCTGAGCGCGGTTGCCCGTCTGCATAACGCATCTCa                         <  1:1380642/50‑1 (MQ=255)
                        gTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATAACGCAtct                           >  1:2954712/1‑48 (MQ=255)
                        gTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATAACGCATCTCa                         <  1:1408255/50‑1 (MQ=255)
                          cTTCGCCTCTATTCTGAGCGCGGTTGCCCGTCTGCATAACGCATCTCAGt                       >  1:426251/1‑50 (MQ=255)
                             cGCCTCTATTCTGAGCGCGGTTGCCCGTCTGCATAACGCATCTCAGtt                      <  1:84094/48‑1 (MQ=255)
                             cGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATAACGCATCTCAGt                       <  1:1991328/47‑1 (MQ=255)
                             cGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATAACGCATCTCAGt                       <  1:1478996/47‑1 (MQ=255)
                             cGCCTCTATGCTGAGCGCGGTTGCCCGTCTGCATAACGCATCTCAGTTa                     <  1:2695662/49‑1 (MQ=255)
                              gCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATAACGCATCTCAGTTAgg                   >  1:20184/1‑50 (MQ=255)
                                  ctATTCTGAGCGCGGTTGCCCGTCTGCATAACGCATCTCAGTTAGGCaa                >  1:2379181/1‑49 (MQ=255)
                                   tATTCTGAGCGCGGCTGCCCGTCTGCATAACGCATCTCAGTTAGGCaaaa              >  1:1579132/1‑50 (MQ=255)
                                     ttCTGAGCGCGGTTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAAt             <  1:257078/49‑1 (MQ=255)
                                     ttCTGAGCGCGGCTGCCCGTCTGCATAACGCATCTCAGTTAGGCaaaa              <  1:2489503/48‑1 (MQ=255)
                                     ttCTGAGCGCGGCTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAAt             <  1:927787/49‑1 (MQ=255)
                                      tCTGAGCGCGGCTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAAt             >  1:671916/1‑48 (MQ=255)
                                       cTGAGCGCGGTTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCa           <  1:1967412/49‑1 (MQ=255)
                                        tGAGCGCGGTTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCAGt         >  1:346882/1‑50 (MQ=255)
                                          aGCGCTGTTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCAGt         <  1:3664451/48‑1 (MQ=255)
                                          aGCGCGGTTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCAGt         <  1:3278677/48‑1 (MQ=255)
                                           gcgcGGCTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCAGtttt      <  1:4307876/50‑1 (MQ=255)
                                           gcgcGGCTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCAGtttt      >  1:2842656/1‑50 (MQ=255)
                                            cgcgGTTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCAGTTTTc     <  1:461205/50‑1 (MQ=255)
                                            cgcgGCTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCAGTTTTc     <  1:4366009/50‑1 (MQ=255)
                                             gcgGCTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCAGTTTTCt    <  1:2289332/50‑1 (MQ=255)
                                             gcgGCTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCAGTTTTCt    >  1:102647/1‑50 (MQ=255)
                                             gcgGCTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCAGTTTTCt    <  1:781311/50‑1 (MQ=255)
                                              cgGCTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCAGTTTTCTg   <  1:3122735/50‑1 (MQ=255)
                                               ggTTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCAGTTTTCTGa  >  1:2840297/1‑50 (MQ=255)
                                                 |                                               
CTGCAAATGATCGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGTTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCAGTTTTCTGA  >  NC_000913/1811159‑1811255

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