Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,195,102 G→A 40.4% intergenic (+277/‑39) ibsD → / → ibsE toxic membrane protein/toxic membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,195,1020GA40.4% 105.7 / 172.4 164intergenic (+277/‑39)ibsD/ibsEtoxic membrane protein/toxic membrane protein
Reads supporting (aligned to +/- strand):  ref base G (42/54);  new base A (33/32);  total (75/89)
Fisher's exact test for biased strand distribution p-value = 4.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

Only 100 of 171 total aligned reads displayed.

GTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGTTTCCCGACTTACTAAC  >  NC_000913/3195003‑3195201
                                                                                                    |                                                                                                   
gTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAg                                                                                                     <  1:823213/100‑1 (MQ=255)
gTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAg                                                                                                     >  1:2774940/1‑100 (MQ=255)
gTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAg                                                                                                     <  1:2705587/100‑1 (MQ=255)
  gACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGcc                                                                                                   <  2:2074754/100‑1 (MQ=255)
  gACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAAcc                                                                                                   <  2:2653186/100‑1 (MQ=255)
      gCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTc                                                                                               >  1:907304/1‑100 (MQ=255)
         cGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAACCCCTcc                                                                                              <  2:1389722/98‑1 (MQ=255)
          cccagaGGCAATTTTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCga                                                                                           >  2:3105461/6‑100 (MQ=21)
            ggAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCgaag                                                                                         <  2:3885185/100‑1 (MQ=255)
             gAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCgaaga                                                                                        >  2:3586365/1‑100 (MQ=255)
              aaGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCgaagaa                                                                                       <  1:319235/100‑1 (MQ=255)
              aaGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCgaagaa                                                                                       <  1:437509/100‑1 (MQ=255)
                ggCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAgg                                                                                     <  2:2760521/100‑1 (MQ=255)
                 tcAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAACCCCTCCCGAAGAAggg                                                                                    <  2:1209743/99‑1 (MQ=255)
                 tcAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAACCCCTCCCGAAGAAggg                                                                                    <  2:823513/99‑1 (MQ=255)
                  cAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAACCCCTCCCGAAGAAgggg                                                                                   <  2:4017862/100‑1 (MQ=255)
                   aaaaTTGCCTCTGGGAAAGCATTGCGCTGGGTCCGGCGCTCGTCAAGACTCGTGGGGCAGCAAGGGGCTGAAACGGGAAAACCCCTCCCGAAGAAGGGGc                                                                                  >  1:407998/1‑100 (MQ=21)
                    aaaCTGCCTCTGGGAAAGCATTGCGCGGGGTCTGACGCTCATCAACAATCGGGAAGTCAGCAAAGGGCTGAAACGGGAAAGCCCCTCCCGAGGAAGGGGc                                                                                  <  1:3715152/100‑1 (MQ=11)
                       ttGCCTCTGGGAAAGCATTGCGGGGGGTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAACCCCTCCCGAAGAAGGGGCCTTg                                                                              <  2:2430775/100‑1 (MQ=255)
                          ccTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGtat                                                                           >  1:279998/1‑100 (MQ=255)
                          ccTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGtat                                                                           <  2:516675/100‑1 (MQ=255)
                          ccTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGtat                                                                           >  2:637773/1‑100 (MQ=255)
                           gtctAGGAAAGCATTGCGCGGGGTCCGGCGCACATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAACCCCCTCCCGAAGAAGGGGCCTTGt                                                                             <  1:1128595/96‑1 (MQ=255)
                            nctGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAGACGGGAACGCCCCCCCCGAAGAAGGGGCCCTGtat                                                                           >  1:548344/2‑98 (MQ=18)
                             ctGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAg                                                                        <  1:3047197/100‑1 (MQ=255)
                               gggAATACATTGGGCGGGGGCCGGCGCTCATCAACTATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAACCCCTCCCGAAGAAGGGGCCTTGTATaan                                                                        <  2:1981714/98‑2 (MQ=18)
                                 gAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGaaa                                                                    <  1:3317476/100‑1 (MQ=255)
                                 gAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGAAAAACGCCTCCCGAAGAAGGGGCCTTGTATAAGGaa                                                                     >  1:2878646/1‑99 (MQ=255)
                                  aaggattggcggGGGGCCGGGCGCTCACCAAAAACGGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAg                                                                   <  1:2245695/92‑1 (MQ=18)
                                     gCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGt                                                                <  1:297775/100‑1 (MQ=255)
                                     gCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAACCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAggg                                                                 >  2:1461137/1‑99 (MQ=255)
                                     gCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAACCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGt                                                                >  2:1023997/1‑100 (MQ=255)
                                      cATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGtt                                                               >  2:3453893/1‑100 (MQ=255)
                                      cATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGtt                                                               >  2:2675189/1‑100 (MQ=255)
                                       aTTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAACCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTa                                                              >  2:3424497/1‑100 (MQ=255)
                                         tGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTatg                                                            <  2:1589803/100‑1 (MQ=255)
                                          gcgcGGGGTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAACCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTatga                                                           <  1:1386665/100‑1 (MQ=255)
                                          gcgcGGGGTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAACCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTatga                                                           <  1:2445205/100‑1 (MQ=255)
                                            gcgGGGTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAACCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTatgatg                                                         >  1:1493853/1‑100 (MQ=255)
                                                ggTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAACCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGc                                                     <  2:283333/100‑1 (MQ=255)
                                                ggTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAACCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGc                                                     >  1:1988930/1‑100 (MQ=255)
                                                 gTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAACCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCt                                                    >  1:3453635/1‑100 (MQ=255)
                                                   ccGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCGGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCg                                                  >  1:2738543/1‑100 (MQ=255)
                                                   ccGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCg                                                  <  2:672467/100‑1 (MQ=255)
                                                   ccGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAACCCCTCCCGAAGAAGGGGCCTTGGATAAGGAAAGGGTTATGATGAAGCTCg                                                  >  1:718883/1‑100 (MQ=255)
                                                     ggCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACTGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGtc                                                >  1:1302939/1‑100 (MQ=255)
                                                       cgcTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGtcat                                              <  1:106837/100‑1 (MQ=255)
                                                       cgcTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGtcat                                              >  1:3767302/1‑100 (MQ=255)
                                                        gcTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAACCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGtcatc                                             >  2:2421130/1‑100 (MQ=255)
                                                          tcatcaACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAACCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTcatcat                                           <  1:401053/100‑1 (MQ=255)
                                                           catcaACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATa                                          >  2:3449662/1‑100 (MQ=255)
                                                           catcaACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATa                                          >  2:3510170/1‑100 (MQ=255)
                                                           catcaACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAACCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATa                                          >  2:1031214/1‑100 (MQ=255)
                                                              caacaaTCGGGGGGCAGCAAGGGGCTGAAACGGGAAAACCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTg                                       >  2:4189265/1‑100 (MQ=255)
                                                              caacaaTCGGGGGGCAGCAAGGGGCTGAAACGGGAAAACCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTg                                       >  1:262928/1‑100 (MQ=255)
                                                               aacaaTCGGGGGGCAGCAAGGGGCTGAAACGGGAAAACCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGACGAGCTTCATCATaaccc                                      <  1:1055980/100‑5 (MQ=18)
                                                                acaaTCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGTGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgt                                     <  1:1472635/100‑1 (MQ=255)
                                                                acaaTCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTgg                                      <  2:1047565/99‑1 (MQ=255)
                                                                 caaTCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgtt                                    >  1:1773861/1‑100 (MQ=255)
                                                                 caaTCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgtt                                    >  2:734894/1‑100 (MQ=255)
                                                                  aaTCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgttg                                   <  2:544740/100‑1 (MQ=255)
                                                                  aaTCGGGGGGCAGCAAGGGGCTGAAACGGGAAAACCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgttg                                   <  2:3372305/100‑1 (MQ=255)
                                                                    tCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgttgtg                                 <  2:3065186/100‑1 (MQ=255)
                                                                     cGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgttgtgt                                <  1:2205188/100‑1 (MQ=255)
                                                                      ggggggCAGCAAGGGGCTGAAACGGGAAAACCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgttgtgtt                               <  2:3550608/100‑1 (MQ=255)
                                                                      ggggagAAGCAAGGGGTTGAAACGGGAATACTTGTCGCGAAGAAGGGGCTTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgttgtg                                 <  1:4237961/98‑1 (MQ=17)
                                                                       gggggCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgttgtgttg                              <  1:2195891/100‑1 (MQ=255)
                                                                        ggggCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgttgtgttgt                             <  1:848296/100‑1 (MQ=255)
                                                                           ggCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTa                           >  1:1790251/1‑100 (MQ=255)
                                                                            gcagcaAGGGGCTGAAACGGGAAAACCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTAc                          <  2:2234091/100‑1 (MQ=255)
                                                                            gcagcaAGGGGCTGAAACGGGAAAACCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTAc                          <  2:3157618/100‑1 (MQ=255)
                                                                             cagcaAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACt                         <  1:1930782/100‑1 (MQ=255)
                                                                             cagcaAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACt                         <  1:286638/100‑1 (MQ=255)
                                                                             cagcaAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACt                         >  2:2161130/1‑100 (MQ=255)
                                                                             cagcaAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACt                         <  1:2043864/100‑1 (MQ=255)
                                                                             cagcaAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACt                         >  1:801747/1‑100 (MQ=255)
                                                                             cagcaAGGGGCTGAAACGGGAAAACCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCAGACTGGTTGTGTTGTTAc                          >  2:3923419/1‑99 (MQ=255)
                                                                              agcaAGGGGCTGAAACGGGAAAACCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTg                        <  1:489295/100‑1 (MQ=255)
                                                                                caAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGtt                      <  2:2167373/100‑1 (MQ=255)
                                                                                 aaGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTa                     <  1:2418709/100‑1 (MQ=255)
                                                                                   ggggCTGAAACGGGAAAACCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTacg                   >  2:2486958/1‑98 (MQ=255)
                                                                                   ggggCTGAAACGGGAAAACCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAg                   <  2:3191692/100‑1 (MQ=255)
                                                                                     ggCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGtt                 >  2:1684925/1‑100 (MQ=255)
                                                                                      gCTGAAACGGGAAAACCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGttt                >  1:4041085/1‑100 (MQ=255)
                                                                                       cTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGttt                >  2:3306574/1‑99 (MQ=255)
                                                                                        tGAAACGGGAAAACCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTAATGTTAAGTTTcc              >  1:71163/1‑100 (MQ=255)
                                                                                         gAAACGGGAAAACCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGTTTccc             <  2:1318470/100‑1 (MQ=255)
                                                                                          aaaCGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGTTTCCCg            <  2:3795091/100‑1 (MQ=255)
                                                                                          aaaCGGGAAACCCCTCCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGTTTCCCg            <  1:2183311/100‑1 (MQ=255)
                                                                                           aaCGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGTTTCCCGa           >  1:3032692/1‑100 (MQ=255)
                                                                                             cGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGTTTCCCGACt         >  2:136105/1‑100 (MQ=255)
                                                                                             cGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGTTTCCCGACt         <  1:3581693/100‑1 (MQ=255)
                                                                                              gggAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGTTTCCCGACtt        <  2:642528/100‑1 (MQ=255)
                                                                                               ggAAATCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGTTTCCCGACTTa       <  2:1908645/100‑1 (MQ=255)
                                                                                                gAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGTTTCCCGAc          <  1:1526121/96‑1 (MQ=255)
                                                                                                gAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGTTTCCCGACTTAc      >  1:4077207/1‑100 (MQ=255)
                                                                                                 aaaaCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGTTTCCCGACTTACt     <  2:1770563/100‑1 (MQ=255)
                                                                                                  aaGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGTTTCCCGACTTACTa    <  1:3784493/100‑1 (MQ=255)
                                                                                                   aGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGTTTCCCGACTTACCaa   >  2:1864402/1‑100 (MQ=255)
                                                                                                    gCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGTTTCCCGACTTACTaac  >  2:3341257/1‑100 (MQ=255)
                                                                                                    |                                                                                                   
GTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGTTTCCCGACTTACTAAC  >  NC_000913/3195003‑3195201

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