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.

GTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGG‑GGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGTTTCCCGACTTACTAAC  >  NC_000913/3195003‑3195201
                                                                                                    |                                                                                                   
gTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGG‑GGCAGCAAGGGGCTGAAACGGGAAAg                                                                                                     <  1:823213/100‑1 (MQ=255)
gTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGG‑GGCAGCAAGGGGCTGAAACGGGAAAg                                                                                                     >  1:2774940/1‑100 (MQ=255)
gTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGG‑GGCAGCAAGGGGCTGAAACGGGAAAg                                                                                                     <  1:2705587/100‑1 (MQ=255)
  gACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGG‑GGCAGCAAGGGGCTGAAACGGGAAAGcc                                                                                                   <  2:2074754/100‑1 (MQ=255)
  gACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGG‑GGCAGCAAGGGGCTGAAACGGGAAAAcc                                                                                                   <  2:2653186/100‑1 (MQ=255)
      gCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGG‑GGCAGCAAGGGGCTGAAACGGGAAAGCCCCTc                                                                                               >  1:907304/1‑100 (MQ=255)
         cGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGG‑GGCAGCAAGGGGCTGAAACGGGAAAACCCCTcc                                                                                              <  2:1389722/98‑1 (MQ=255)
          cccagaGGCAATTTTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGG‑GGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCga                                                                                           >  2:3105461/6‑100 (MQ=21)
            ggAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGG‑GGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCgaag                                                                                         <  2:3885185/100‑1 (MQ=255)
             gAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGG‑GGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCgaaga                                                                                        >  2:3586365/1‑100 (MQ=255)
              aaGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGG‑GGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCgaagaa                                                                                       <  1:319235/100‑1 (MQ=255)
              aaGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGG‑GGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCgaagaa                                                                                       <  1:437509/100‑1 (MQ=255)
                ggCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGG‑GGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAgg                                                                                     <  2:2760521/100‑1 (MQ=255)
                 tcAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGG‑GGCAGCAAGGGGCTGAAACGGGAAAACCCCTCCCGAAGAAggg                                                                                    <  2:1209743/99‑1 (MQ=255)
                 tcAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGG‑GGCAGCAAGGGGCTGAAACGGGAAAACCCCTCCCGAAGAAggg                                                                                    <  2:823513/99‑1 (MQ=255)
                  cAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGG‑GGCAGCAAGGGGCTGAAACGGGAAAACCCCTCCCGAAGAAgggg                                                                                   <  2:4017862/100‑1 (MQ=255)
                   aaaaTTGCCTCTGGGAAAGCATTGCGCTGGGTCCGGCGCTCGTCAAGACTCGTGG‑GGCAGCAAGGGGCTGAAACGGGAAAACCCCTCCCGAAGAAGGGGc                                                                                  >  1:407998/1‑100 (MQ=21)
                    aaaCTGCCTCTGGGAAAGCATTGCGCGGGGTCTGACGCTCATCAACAATCGGGAAGTCAGCAAAGGGCTGAAACGGGAAAGCCCCTCCCGAGGAAGGGGc                                                                                  <  1:3715152/100‑1 (MQ=11)
                       ttGCCTCTGGGAAAGCATTGCGGGGGGTCCGGCGCTCATCAACAATCGGGG‑GGCAGCAAGGGGCTGAAACGGGAAAACCCCTCCCGAAGAAGGGGCCTTg                                                                              <  2:2430775/100‑1 (MQ=255)
                          ccTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGG‑GGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGtat                                                                           >  1:279998/1‑100 (MQ=255)
                          ccTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGG‑GGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGtat                                                                           <  2:516675/100‑1 (MQ=255)
                          ccTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGG‑GGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGtat                                                                           >  2:637773/1‑100 (MQ=255)
                           gtctAGGAAAGCATTGCGCGGGGTCCGGCGCACATCAACAATCGGGG‑GGCAGCAAGGGGCTGAAACGGGAAACCCCCTCCCGAAGAAGGGGCCTTGt                                                                             <  1:1128595/96‑1 (MQ=255)
                            nctGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGG‑GGCAGCAAGGGGCTGAGACGGGAACGCCCCCCCCGAAGAAGGGGCCCTGtat                                                                           >  1:548344/2‑98 (MQ=18)
                             ctGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGG‑GGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAg                                                                        <  1:3047197/100‑1 (MQ=255)
                               gggAATACATTGGGCGGGGGCCGGCGCTCATCAACTATCGGGG‑GGCAGCAAGGGGCTGAAACGGGAAAACCCCTCCCGAAGAAGGGGCCTTGTATaan                                                                        <  2:1981714/98‑2 (MQ=18)
                                 gAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGG‑GGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGaaa                                                                    <  1:3317476/100‑1 (MQ=255)
                                 gAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGG‑GGCAGCAAGGGGCTGAAACGGAAAAACGCCTCCCGAAGAAGGGGCCTTGTATAAGGaa                                                                     >  1:2878646/1‑99 (MQ=255)
                                  aaggattggcggGGGGCCGGGCGCTCACCAAAAACGGGGG‑GGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAg                                                                   <  1:2245695/92‑1 (MQ=18)
                                     gCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGG‑GGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGt                                                                <  1:297775/100‑1 (MQ=255)
                                     gCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGG‑GGCAGCAAGGGGCTGAAACGGGAAAACCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAggg                                                                 >  2:1461137/1‑99 (MQ=255)
                                     gCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGG‑GGCAGCAAGGGGCTGAAACGGGAAAACCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGt                                                                >  2:1023997/1‑100 (MQ=255)
                                      cATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGG‑GGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGtt                                                               >  2:3453893/1‑100 (MQ=255)
                                      cATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGG‑GGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGtt                                                               >  2:2675189/1‑100 (MQ=255)
                                       aTTGCGCGGGGTCCGGCGCTCATCAACAATCGGGG‑GGCAGCAAGGGGCTGAAACGGGAAAACCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTa                                                              >  2:3424497/1‑100 (MQ=255)
                                         tGCGCGGGGTCCGGCGCTCATCAACAATCGGGG‑GGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTatg                                                            <  2:1589803/100‑1 (MQ=255)
                                          gcgcGGGGTCCGGCGCTCATCAACAATCGGGG‑GGCAGCAAGGGGCTGAAACGGGAAAACCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTatga                                                           <  1:1386665/100‑1 (MQ=255)
                                          gcgcGGGGTCCGGCGCTCATCAACAATCGGGG‑GGCAGCAAGGGGCTGAAACGGGAAAACCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTatga                                                           <  1:2445205/100‑1 (MQ=255)
                                            gcgGGGTCCGGCGCTCATCAACAATCGGGG‑GGCAGCAAGGGGCTGAAACGGGAAAACCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTatgatg                                                         >  1:1493853/1‑100 (MQ=255)
                                                ggTCCGGCGCTCATCAACAATCGGGG‑GGCAGCAAGGGGCTGAAACGGGAAAACCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGc                                                     <  2:283333/100‑1 (MQ=255)
                                                ggTCCGGCGCTCATCAACAATCGGGG‑GGCAGCAAGGGGCTGAAACGGGAAAACCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGc                                                     >  1:1988930/1‑100 (MQ=255)
                                                 gTCCGGCGCTCATCAACAATCGGGG‑GGCAGCAAGGGGCTGAAACGGGAAAACCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCt                                                    >  1:3453635/1‑100 (MQ=255)
                                                   ccGGCGCTCATCAACAATCGGGG‑GGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCGGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCg                                                  >  1:2738543/1‑100 (MQ=255)
                                                   ccGGCGCTCATCAACAATCGGGG‑GGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCg                                                  <  2:672467/100‑1 (MQ=255)
                                                   ccGGCGCTCATCAACAATCGGGG‑GGCAGCAAGGGGCTGAAACGGGAAAACCCCTCCCGAAGAAGGGGCCTTGGATAAGGAAAGGGTTATGATGAAGCTCg                                                  >  1:718883/1‑100 (MQ=255)
                                                     ggCGCTCATCAACAATCGGGG‑GGCAGCAAGGGGCTGAAACTGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGtc                                                >  1:1302939/1‑100 (MQ=255)
                                                       cgcTCATCAACAATCGGGG‑GGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGtcat                                              <  1:106837/100‑1 (MQ=255)
                                                       cgcTCATCAACAATCGGGG‑GGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGtcat                                              >  1:3767302/1‑100 (MQ=255)
                                                        gcTCATCAACAATCGGGG‑GGCAGCAAGGGGCTGAAACGGGAAAACCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGtcatc                                             >  2:2421130/1‑100 (MQ=255)
                                                          tcatcaACAATCGGGG‑GGCAGCAAGGGGCTGAAACGGGAAAACCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTcatcat                                           <  1:401053/100‑1 (MQ=255)
                                                           catcaACAATCGGGG‑GGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATa                                          >  2:3449662/1‑100 (MQ=255)
                                                           catcaACAATCGGGG‑GGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATa                                          >  2:3510170/1‑100 (MQ=255)
                                                           catcaACAATCGGGG‑GGCAGCAAGGGGCTGAAACGGGAAAACCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATa                                          >  2:1031214/1‑100 (MQ=255)
                                                              caacaaTCGGGG‑GGCAGCAAGGGGCTGAAACGGGAAAACCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTg                                       >  2:4189265/1‑100 (MQ=255)
                                                              caacaaTCGGGG‑GGCAGCAAGGGGCTGAAACGGGAAAACCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTg                                       >  1:262928/1‑100 (MQ=255)
                                                               aacaaTCGGGG‑GGCAGCAAGGGGCTGAAACGGGAAAACCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGACGAGCTTCATCATaaccc                                      <  1:1055980/100‑5 (MQ=18)
                                                                acaaTCGGGG‑GGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGTGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgt                                     <  1:1472635/100‑1 (MQ=255)
                                                                acaaTCGGGG‑GGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTgg                                      <  2:1047565/99‑1 (MQ=255)
                                                                 caaTCGGGG‑GGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgtt                                    >  1:1773861/1‑100 (MQ=255)
                                                                 caaTCGGGG‑GGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgtt                                    >  2:734894/1‑100 (MQ=255)
                                                                  aaTCGGGG‑GGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgttg                                   <  2:544740/100‑1 (MQ=255)
                                                                  aaTCGGGG‑GGCAGCAAGGGGCTGAAACGGGAAAACCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgttg                                   <  2:3372305/100‑1 (MQ=255)
                                                                    tCGGGG‑GGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgttgtg                                 <  2:3065186/100‑1 (MQ=255)
                                                                     cGGGG‑GGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgttgtgt                                <  1:2205188/100‑1 (MQ=255)
                                                                      gggg‑ggCAGCAAGGGGCTGAAACGGGAAAACCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgttgtgtt                               <  2:3550608/100‑1 (MQ=255)
                                                                      gggg‑agAAGCAAGGGGTTGAAACGGGAATACTTGTCGCGAAGAAGGGGCTTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgttgtg                                 <  1:4237961/98‑1 (MQ=17)
                                                                       ggg‑ggCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgttgtgttg                              <  1:2195891/100‑1 (MQ=255)
                                                                        gg‑ggCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgttgtgttgt                             <  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)
                                                                                                    |                                                                                                   
GTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGG‑GGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGTTTCCCGACTTACTAAC  >  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: ‑