Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,370,614 T→A 22.3% intergenic (‑202/‑74) yjeH ← / → groS putative transporter/Cpn10 chaperonin GroES, small subunit of GroESL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,370,6140TA22.3% 363.8 / 92.7 173intergenic (‑202/‑74)yjeH/groSputative transporter/Cpn10 chaperonin GroES, small subunit of GroESL
Reads supporting (aligned to +/- strand):  ref base T (64/70);  new base A (13/26);  total (77/96)
Fisher's exact test for biased strand distribution p-value = 1.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.99e-01

Only 100 of 174 total aligned reads displayed.

TGTGAGGTGAATCAGGGTTTTCACCCGATTTTGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCGTCCATTGCATGATCG  >  NC_000913/4370516‑4370713
                                                                                                  |                                                                                                   
tgtgAGGTGAATCAGGGTTTTCACCCGATTTTGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTc                                                                                                    >  2:2846230/1‑100 (MQ=255)
gaaatgtgaggtgaatcagggttttcacccgattttgtgctgatcagaattttttttctttttcccccttgaaggggcgaagcctcatccccatttGGTc                                                                                                    <  1:1603773‑M2/4‑1 (MQ=255)
 gtgATGTGAATCAGGGTTTTCACCCGATTTTGTGCTGATCAGAACTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGACa                                                                                                   <  1:1161691/100‑1 (MQ=255)
    tgtgaggtgaatcagggttttcacccgattttgtgctgatcagaattttttttctttttcccccttgaaggggcgaagcctcatccccatttGGTCACCa                                                                                                >  2:3367949‑M2/93‑100 (MQ=255)
      tgaggtgaatcagggttttcacccgattttgtgctgatcagaattttttttttttttttcccttgaaggggcgaagcctcatccccatttGGTCACCAGc                                                                                              >  1:2655552‑M2/91‑100 (MQ=255)
      tgaggtgaatcagggttttcacccgattttgtgctgatcagaattttttttctttttcccccttgaaggggcgaagcctcatccccatttGGTCACCAGc                                                                                              <  1:667113‑M2/10‑1 (MQ=255)
       tGAATCAGGGTTTTCACCCGATTTTGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCACCAGcc                                                                                             >  2:1846429/1‑100 (MQ=255)
        gAATCAGGGTTTTCACCCGATTTTGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGACACCAGcc                                                                                             >  1:534895/1‑99 (MQ=255)
          gtgaatcagggttttcacccgattttgtgctgatcagaattttttttctttttcccccttgaaggggcgaagcctcatccccatttGGTCACCAGCCggg                                                                                          <  2:694542‑M2/14‑1 (MQ=255)
          aTCAGGGTTTTCACCCGATTTTGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGACACCAGCCggg                                                                                          >  1:3433400/1‑100 (MQ=255)
            gaatcagggttttcacccgattttgtgctgatcagaattttttttcttttttccccttgaaggggcgaagcctcatccccatttGGTCACCAGCCgg                                                                                           >  1:1545433‑M2/85‑97 (MQ=255)
            gaatcagggttttcacccgattttgtgctgatcagaattttttttctttttcccccttgcaggggcgaagcctcatccccatttCGTCACCAGccgggta                                                                                        >  2:3638715‑M2/85‑98 (MQ=255)
             aatcagggttttcacccgattttgtgctgatcagaattttttttctttttcccccttgaaggggcgaagcctcatccccatttGGTCACCAGCCGGGaaa                                                                                       <  1:1620308‑M2/17‑1 (MQ=255)
              gggTTTTCACCCGATTTTGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCACCAGCCGGGAAAc                                                                                      <  1:346967/100‑1 (MQ=255)
              atcagggttttcacccgattttgtgctgatcagaattttttttctttttcccccttgaaggggcgaagcctcatccccatttGGTCACCAGCCGGGAAAc                                                                                      <  2:2349445‑M2/18‑1 (MQ=255)
              atcagggttttcacccgattttgtgctgatcagaattttttttctttttcccccttgaaggggcgaagcctcatccccatttGGTCACCAGCCGGGAAAc                                                                                      >  1:945597‑M2/83‑100 (MQ=255)
               tcagggttttcacccgattttgggctgatcagaattttttttctttttcccccttgaaggggcgaagcctcatccccatttGGTCACCAGCCGGGAAAcc                                                                                     <  2:3116148‑M2/19‑1 (MQ=255)
               ggTTTTCACCGGGTTTTGTGCTGATCAGAATTTTTGTGCTGTTGCCCCCTTCAAGGGGCGAAGCCTCATCCCCAGT‑‑‑‑TGGTCACCAGCCGGGAAACCACGt                                                                                 <  1:3164418/100‑1 (MQ=255)
                cagggttttcacccgattttgtgctgatcagaattttttttctttttcccccttgaaggggcgaagcctcatccccatttGGTCACCAGCCGGGAAACCa                                                                                    <  2:1482832‑M2/20‑1 (MQ=255)
                cagggttttcacccgattttgtgctgatcagaattttttttctttttcccccttgaaggggcgaagcctcatccccatttGGTCACCAGCCGGGAAACCa                                                                                    >  1:1961401‑M2/81‑100 (MQ=255)
                cagggttttcacccgattttgtgctgatcagaattttttttctttttcccccttgaaggggcgaagcctcatccccatttGGTCACCAGCCGGGAAACCa                                                                                    >  2:1357437‑M2/81‑100 (MQ=255)
                  gggttttcacccgattttgtgctgatcagaattttttttttttttttctctttaaggggcgaagcctcatccccatttGTTCACCAGCCGGGAAACCACg                                                                                  >  1:2074614‑M2/79‑100 (MQ=255)
                    gttttcacccgattttgtgctgatcagaattttttttctttttcccccttgaaggggcgaagcctcatccccatttGGTCACCAGCCGGGAAACCACGTa                                                                                <  1:2030323‑M2/24‑1 (MQ=255)
                    gttttcacccgattttgtgctgatcagaattttttttctttttcccccttgaaggggcgaagcctcatccccatttGGTCACCAGCCGGGAAACCACGTa                                                                                >  2:376289‑M2/77‑100 (MQ=255)
                     ttttcacccgattttgtgctgatcagaattttttttctttttcccccttgaaggggcgaagcctcatccccatttGGTCACCAGCCGGGAAACCACGtan                                                                               <  1:3602951‑M2/25‑2 (MQ=255)
                     ttttcacccgattttgtgctgatcagaattttttttctttttcccccttgaaggggcgaagcctcatccccatttGGTCACCAGCCGGGAAACCACGTaa                                                                               <  2:501553‑M2/25‑1 (MQ=255)
                      tttcacccgattttgtgctgatcagaattttttttctttttcccccttgaaggggcgaagcctcatccccatttGGTCACCAGCCGGGAAACCACGTAAg                                                                              <  2:1414519‑M2/26‑1 (MQ=255)
                      tttcacccgattttgtgctgatcagaattttttttctttttcccccttgaaggggcgaagcctcatccccatttGGTCACCAGCCGGGAAACCACGTAAg                                                                              >  2:667051‑M2/75‑100 (MQ=255)
                      tatcacccgattttgtgctgatcagaattttttttctttttcccccttgaaggggcgaagcctcatccccatttGGTCACCAGCCGGGAAACCACGTAAg                                                                              >  2:2503876‑M2/75‑100 (MQ=255)
                       cccGATTTTGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGc                                                                             >  2:2793708/1‑100 (MQ=255)
                       cccGATTTTGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCCCATCCCCATTTCTCTGGTCACCAGCCGGGAAACCACGTAAGc                                                                             >  2:2272917/1‑100 (MQ=255)
                        tcacccgattttgtgctgatcagaattttttttctttttcccccttgaaggggcgaagcctcatccccatttGGTCACCAGCCGGGAAACCACGTAAGCt                                                                            >  1:1715343‑M2/73‑100 (MQ=255)
                         cacccgattttgtgctgatcagaattttttttctttttcccccttgaaggggcgaagcctcatccccatttGGTCACCAGCCGGGAAACCACGTAAGCTc                                                                           <  2:1200895‑M2/29‑1 (MQ=255)
                           aTTTTGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCACCAGCCGGGAAACCACGTAAGCTCCg                                                                         <  2:102524/100‑1 (MQ=255)
                            ttttGGGCTAATCAGAATTTTTTTTCTATTTCCCCCTTGAAGGCGCGAAGCCTCTTCCCCATTTCTCTGGTTACCAGCCGGGAAACCACGTAAg                                                                              <  1:3254211/94‑1 (MQ=255)
                            ccgattttgtgctgatcagaattttttttctttttcccccttgaaggggcgaagcctcatccccatttGGTCACCAGCCGGGAAACCACGTAAGCTCCgg                                                                        <  2:2882665‑M2/32‑1 (MQ=255)
                             tttGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGc                                                                       <  2:2482435/100‑1 (MQ=255)
                                ttttgtgctgatcagaattttttttctttttcccccttgaaggggcgaagcctcatccccatttGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTc                                                                    <  1:3103987‑M2/36‑1 (MQ=255)
                                 tgCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCa                                                                   <  2:3450551/100‑1 (MQ=255)
                                  gCTGATCAGAATTTTTTTTTTTTTTTTTTTTTTAAGGAGCGAAGCCTCATCTTCTTT‑‑‑TGTTTACCAGTCGTGAAACCACGTAACCTCCGGCTTCATCCa                                                               >  1:414822/1‑99 (MQ=255)
                                    gtgctgatcagaattttttttctttttcccccttgaaggggcgaagcctcatccccatttGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCAccc                                                                <  2:284287‑M2/40‑1 (MQ=255)
                                     gATCAGAATTTTTTTTCTTTTTTTCCTTTTAAGGTGCTAAGCCTCATCCCCATTTC‑‑‑‑GTCACCATCCGGGTAACCATCTAATCTCCGTCGTCACCCATAAc                                                           >  1:2680751/1‑100 (MQ=255)
                                      gctgatcagaattttttttctttttcccccttgaaggggcgaagcctcatccccatttGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCAt                                                              <  2:3270046‑M2/42‑1 (MQ=255)
                                      aTCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCAt                                                              >  2:568045/1‑100 (MQ=255)
                                        tgatcagaattttttttctttttcccccttgaaggggcgaagcctcatccccatttGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATaa                                                            >  2:1521471‑M2/57‑100 (MQ=255)
                                         gatcagaattttttttctttttcccccttgaaggggcgaagcctcatccccatttGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCAtaan                                                           <  1:2256397‑M2/45‑2 (MQ=255)
                                          gAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACa                                                          <  2:3054077/100‑1 (MQ=255)
                                          gAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACa                                                          <  1:2643777/100‑1 (MQ=255)
                                          atcagaattttttttctttttcccccttgaaggggcgaagcctcatccccatttGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACa                                                          >  1:2196379‑M2/55‑100 (MQ=255)
                                           tcagaatttttttttttttttccccttgaaggggcgaagcctcatccccatttGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAg                                                         >  1:454891‑M2/54‑100 (MQ=255)
                                              tttttttCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATa                                                      <  2:1384197/100‑1 (MQ=255)
                                              gaattttttttctttttcccccttgaaggggcgaagcctcatccccatttGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATa                                                      >  2:229939‑M2/51‑100 (MQ=255)
                                               aattttttttctttttcccccttgaaggggcgaagcctcatccccatttGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATAc                                                     <  1:1680869‑M2/51‑1 (MQ=255)
                                                 ttttctttttcCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACg                                                    <  2:173240/99‑1 (MQ=255)
                                                   ttctttttcCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGAc                                                 <  2:1529045/100‑1 (MQ=255)
                                                    ttttttttttttttccttgaaggggcgaagcctcatccccatttGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACTCATAACTGATACGGACt                                                >  1:507029‑M2/45‑100 (MQ=255)
                                                     ctttttcCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACtt                                               >  2:1697265/1‑100 (MQ=255)
                                                         ttcCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTtctc                                           <  1:2818027/100‑1 (MQ=255)
                                                         ctttttcccccttgaaggggcgaagcctcatccccatttGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTtctc                                           <  1:555906‑M2/61‑1 (MQ=255)
                                                          tttttcccccttgaaggggcgaagcctcatccccatttGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCa                                          >  2:1856132‑M2/39‑100 (MQ=255)
                                                           ttttcccccttgaaggggcgaagcctcatccccatttGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCaa                                         <  2:1023986‑M2/63‑1 (MQ=255)
                                                           ttttcccccttgaaggggcgaagcctcatccccatttGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCaa                                         >  1:2133528‑M2/38‑100 (MQ=255)
                                                            tttcccccttgaaggggcgaagcctcatccccatttGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCaaa                                        <  1:3317926‑M2/64‑1 (MQ=255)
                                                             ttcccccttgaaggggcgaagcctcatccccatttGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAg                                       <  2:1974411‑M2/65‑1 (MQ=255)
                                                             ttcccccttgaaggggcgaagcctcatccccatttGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAg                                       <  2:237938‑M2/65‑1 (MQ=255)
                                                                ttGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGgag                                    <  1:3232052/100‑1 (MQ=255)
                                                                ccccttgaaggggcgaagcctcatccccatttGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGgag                                    >  2:3090191‑M2/33‑100 (MQ=255)
                                                                 tGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGgaga                                   <  1:573462/100‑1 (MQ=255)
                                                                 cccttgaaggggcgaagcctcatccccatttGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGgaga                                   <  1:859080‑M2/69‑1 (MQ=255)
                                                                 cccttgaaggggcgaagcctcatccccatttGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGgaga                                   <  1:1498345‑M2/69‑1 (MQ=255)
                                                                  ccttgaaggggcgaagcctcatccccatttGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGagag                                  <  2:2074627‑M2/70‑1 (MQ=255)
                                                                  ccttgaaggggcgaagcctcatccccatatGGTCACCAGCCGGGAAACCACGTAAGCTACGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGagag                                  <  1:391123‑M2/70‑1 (MQ=255)
                                                                    ttgaaggggcgaagcctcatccccatttGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGtt                                <  1:1232176‑M2/72‑1 (MQ=255)
                                                                    aGGGGCGAAGCCTCATCCCCATT‑‑‑‑TGGTCACCTGCCTGGAATCCACGTATGCTCCGGCTGCACCCATTACAGATACTGACTTTTTCAATGGTGAGTTAt                              >  1:304865/1‑98 (MQ=255)
                                                                     gatctgggcgaagcctcatccccatttGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTa                               <  1:3084723‑M2/73‑1 (MQ=255)
                                                                      gggCGAAGCCTCATCCCCATTTCTCTGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTAt                              >  2:1073480/1‑100 (MQ=255)
                                                                      gggCGAAGCCTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTa                               <  1:1022719/99‑1 (MQ=255)
                                                                        gCGAAGCCTCATCCCCATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCa                            <  1:2793691/100‑1 (MQ=255)
                                                                        ccgatctcTCATCCCCATTTCTCTGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCa                            <  1:907103/93‑1 (MQ=255)
                                                                         cGAAGCCTCATCCCCATTTCTCTGGACACCAGCCGGGACACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCaa                           >  1:2286093/1‑100 (MQ=255)
                                                                          gggcgaagcctcatccccatttGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAgagttaagatc                          >  2:3084745‑M2/23‑95 (MQ=255)
                                                                          gggcaaagcctcatccccatttGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACACATAACAGAAACGGACTTTCTACAAGGAGCGCTATCaa                           >  1:3134493‑M2/23‑99 (MQ=255)
                                                                            gcgaagcctcatccccatttGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGa                        >  2:1304108‑M2/21‑100 (MQ=255)
                                                                             cgaagcctcatcccaatatGGTCACCAGCCGGGAAACCATGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTAGCa                            <  1:3566858‑M2/76‑1 (MQ=255)
                                                                               cTCATCCCCATTTCTCTGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCaagatcgg                     >  2:907108/1‑94 (MQ=255)
                                                                               aagcctcatccccatttGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAata                     <  2:186790‑M2/83‑1 (MQ=255)
                                                                                   tCCCCATTTCTCTGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCg                 <  1:3563690/100‑1 (MQ=255)
                                                                                   ctcatcccgatttGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCg                 <  1:2467625‑M2/87‑1 (MQ=255)
                                                                                   ctcatccccatttGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCg                 <  1:1527525‑M2/87‑1 (MQ=255)
                                                                                     cccATTTCTCTGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCGTc               <  1:2342517/100‑1 (MQ=255)
                                                                                     catccccatttGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCCAAGGAGAGTTATCAATGAATATTCGTc               >  1:109432‑M2/12‑100 (MQ=255)
                                                                                      ccATTTCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGAACTTTCTCAAAGAAGAGATATCAAAGAATATTCg                 >  1:1171690/1‑97 (MQ=255)
                                                                                          ttCTCTGGACACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCGTCCATTg          <  1:1020577/100‑1 (MQ=255)
                                                                                            ctctGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCGTcagatcgg        >  2:1960416/1‑93 (MQ=255)
                                                                                            atttGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCGTCCATTGCa        >  2:3152268‑M2/5‑100 (MQ=255)
                                                                                             tttGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCGTCCATTGCAt       <  2:2655567‑M2/97‑1 (MQ=255)
                                                                                               tGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCGTCCATTGCATGa     <  2:2024862‑M2/99‑1 (MQ=255)
                                                                                                ggTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCGTCCATTGCATGAt    >  2:790648/1‑100 (MQ=255)
                                                                                                 gTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCGTCCATTGCATGATc   >  2:401683/1‑100 (MQ=255)
                                                                                                  tCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCGTCCATTGCATGATcg  <  1:1787594/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
TGTGAGGTGAATCAGGGTTTTCACCCGATTTTGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCGTCCATTGCATGATCG  >  NC_000913/4370516‑4370713

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: