Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,524,463 Δ1 bp 15.4% intergenic (+95/+18) yncE → / ← ansP ATP‑binding protein, periplasmic, function unknown/L‑asparagine transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,524,4630A.15.4% 386.8 / 69.3 136intergenic (+95/+18)yncE/ansPATP‑binding protein, periplasmic, function unknown/L‑asparagine transporter
Reads supporting (aligned to +/- strand):  ref base A (61/52);  new base . (13/8);  total (75/61)
Fisher's exact test for biased strand distribution p-value = 6.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 137 total aligned reads displayed.

TGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTC  >  NC_000913/1524364‑1524561
                                                                                                    |                                                                                                  
tGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTa                                                                                                    >  1:3290353/1‑100 (MQ=255)
 gTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAg                                                                                                   >  1:1524982/1‑100 (MQ=255)
    aTAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCg                                                                                                <  2:2954803/100‑1 (MQ=255)
    aTAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCg                                                                                                >  2:1390593/1‑100 (MQ=255)
      atttatGGCGAGGGTGAAAACGGTGAGACTTGCCTTCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCggg                                                                                              <  2:751938/97‑1 (MQ=255)
      aGATATGGCGAGGGTTAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCggg                                                                                              <  2:3765345/100‑1 (MQ=255)
       gATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCTACGT‑GCCGGGc                                                                                             <  2:617249/100‑1 (MQ=255)
       gATATGGCAAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGc                                                                                             >  2:2947019/1‑100 (MQ=255)
        atatGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCTACGT‑GCCGGGCa                                                                                            <  1:3036701/100‑1 (MQ=255)
        atatGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCTACGT‑GCCGGGCa                                                                                            <  1:2158613/100‑1 (MQ=255)
         tatGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCTACGT‑GCCGGGCaa                                                                                           <  1:3373691/100‑1 (MQ=255)
          atGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAAt                                                                                          >  1:3797034/1‑100 (MQ=255)
           tGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCGTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCTACGT‑GCCGGGCaaat                                                                                         <  2:707149/100‑3 (MQ=255)
               gAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCTGGGCAATTTTGc                                                                                     <  2:2514937/100‑1 (MQ=255)
               atctgTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGc                                                                                     <  1:1275887/96‑1 (MQ=255)
                 gggTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCTACGT‑GCCGGGCaaattcagat                                                                                   <  2:524893/100‑9 (MQ=37)
                 gggTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCaa                                                                                   <  2:2976134/100‑1 (MQ=255)
                  ggcgAAAACGGTGAGATTTCCCTGCCTTTTCTCACGGCGTATGCGGCACTTTCGGCAGTAATCTGAATTTGCTCGGCTACGT‑GCCGTGCAAattcagat                                                                                   <  1:2631222/96‑6 (MQ=255)
                  ggTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCaa                                                                                   >  1:980246/1‑99 (MQ=255)
                    tGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTa                                                                                <  2:256039/100‑1 (MQ=255)
                    tGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTa                                                                                >  2:2208755/1‑100 (MQ=255)
                    tGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATATGAATTAGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTa                                                                                >  1:2680065/1‑100 (MQ=255)
                    ngAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTa                                                                                >  1:2971625/2‑100 (MQ=255)
                      aaaaCGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAAc                                                                              >  2:2292635/1‑100 (MQ=255)
                       aaaCGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACt                                                                             >  2:1500349/1‑100 (MQ=255)
                         aCGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCg                                                                           <  2:901099/100‑1 (MQ=255)
                           ggTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGcc                                                                         <  1:3398628/100‑1 (MQ=255)
                            gTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCg                                                                           >  1:3422840/1‑97 (MQ=255)
                                gaCTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTc                                                                    >  2:472662/1‑100 (MQ=255)
                                 aCTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCa                                                                   <  1:521092/100‑1 (MQ=255)
                                    tgcctgccTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGc                                                                <  1:3015047/100‑1 (MQ=255)
                                    tgcctgccTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGc                                                                >  1:2998504/1‑100 (MQ=255)
                                     gcctgcctTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCt                                                               >  2:2488132/1‑100 (MQ=255)
                                      cctgcctTTTATCACGGCGTATACGGCACTTTCGTCAGTAATCTGAATTTGCCGGGC‑ACGTAGCCGTGCAATTTTGCAATTAACTCGCCGTTTCAGGCt                                                               <  1:1504810/99‑1 (MQ=255)
                                       ctgcctTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTa                                                             <  2:3352465/100‑1 (MQ=255)
                                       ctgcctTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTa                                                             <  2:2493502/100‑1 (MQ=255)
                                        tgcctTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTaa                                                            <  1:127135/100‑1 (MQ=255)
                                        tgcctTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTaa                                                            >  2:2545452/1‑100 (MQ=255)
                                          cctTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAac                                                          <  2:412798/100‑1 (MQ=255)
                                           ctTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAaca                                                         >  2:1688394/1‑100 (MQ=255)
                                            ttttATCACGGCGTATCCGGCACTTGCGTAAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAaca                                                         <  1:429656/99‑1 (MQ=255)
                                             tttATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAAcaca                                                       <  2:1349418/100‑1 (MQ=255)
                                              ttATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACaa                                                      >  2:1725673/1‑100 (MQ=255)
                                              ttATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACaa                                                      <  2:3520584/100‑1 (MQ=255)
                                               tATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAAt                                                     >  2:726034/1‑100 (MQ=255)
                                                aTCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAAtt                                                    >  1:1469433/1‑100 (MQ=255)
                                                 tCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTaaagatcggaa                                                   >  2:3708465/1‑92 (MQ=255)
                                                  cACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTc                                                  >  2:3800271/1‑100 (MQ=255)
                                                   aCGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTCGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTcc                                                 >  2:1561555/1‑100 (MQ=255)
                                                   aCGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTcc                                                 <  1:2871443/100‑1 (MQ=255)
                                                   aCGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTcc                                                 <  1:2845806/100‑1 (MQ=255)
                                                   aCGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTc                                                  >  1:2086416/1‑99 (MQ=255)
                                                   aCGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTACCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTcc                                                 >  2:1558138/1‑100 (MQ=255)
                                                    cGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCt                                                >  2:2620024/1‑100 (MQ=255)
                                                    cGGCGTATACGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCACTTAACTCGAAGTTTCAGGCATAAAAACAATTTcc                                                 >  1:2260276/1‑99 (MQ=255)
                                                      gCGTATCCTGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAAGTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTg                                               <  1:1014753/99‑1 (MQ=255)
                                                            ccGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTttct                                        >  2:927073/1‑100 (MQ=255)
                                                            ccGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTttct                                        <  1:3159163/100‑1 (MQ=255)
                                                             cGGCGCTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTttctt                                       <  1:1298873/100‑1 (MQ=255)
                                                             cGGCACTTTCGTCAGTAATCTGAATTTTCCCGGC‑ACGTAGCCGGTCAATTTTGCAATTAACTCTCCTTTTCAGGCTTAAACACAATTTCCTGTTTttctt                                       <  1:895243/100‑1 (MQ=255)
                                                              ggCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCATTTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTttctt                                       <  1:2972900/99‑1 (MQ=255)
                                                              ggCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTttcttc                                      >  2:3095060/1‑100 (MQ=255)
                                                               gCACTTTCGTTAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCATTTAACTCGCCGTTTCAGGCTTAAACACAGGTTCCTGTTTttct                                        <  1:404631/97‑1 (MQ=255)
                                                                  cTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCt                                  >  1:2874883/1‑100 (MQ=255)
                                                                     tCGTGAGTAATCTCAAGTTGCCGGGG‑ACTCTGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCAtcctcc                               <  1:2690831/100‑1 (MQ=255)
                                                                      cGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCt                              <  2:1316033/100‑1 (MQ=255)
                                                                       gTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTACCTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTc                             <  2:1905441/100‑1 (MQ=255)
                                                                       gTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTTAACACCATTTCCTGTTTTTCTTCATCTTCCt                              >  1:172496/1‑99 (MQ=255)
                                                                       gTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTc                             <  2:2990151/100‑1 (MQ=255)
                                                                          aGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTcgac                          >  2:2475887/1‑100 (MQ=255)
                                                                                 tGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGgc                   >  2:2489358/1‑100 (MQ=255)
                                                                                  gAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGgcg                  <  2:1188392/100‑1 (MQ=255)
                                                                                   aaattGCCCGGCTACGT‑GCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGgcgc                 >  2:2779406/4‑100 (MQ=39)
                                                                                   aaTTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGgcgc                 <  2:1699144/100‑1 (MQ=255)
                                                                                   aaTTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGgcgc                 <  1:1171193/100‑1 (MQ=255)
                                                                                   aaTTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGgcgc                 >  1:3637805/1‑100 (MQ=255)
                                                                                   aaTTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGgcgc                 >  2:2770305/1‑100 (MQ=255)
                                                                                     tttGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAg               >  2:3267698/1‑100 (MQ=255)
                                                                                     tttGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAg               <  1:2299604/100‑1 (MQ=255)
                                                                                     attGCCCGGCTACGT‑GCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCTTCGACGACTg                     >  1:3629804/2‑94 (MQ=38)
                                                                                     attGCCCGGCTACGT‑GCCGGGCAATTGTGCAATTAACTCGCCGTTTCAGGCTGAAACACAATTTACTGTTTTTCGTCATCTTCCTCGACGAATGGCGCAg               >  1:714199/2‑100 (MQ=255)
                                                                                     attGCCCGGCTACGT‑GCCGAGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTACCAGTATTTATTCATCTTCCTAGACGACt                      >  1:1091989/2‑93 (MQ=255)
                                                                                       tGCCCGGCTACGT‑GCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTg             >  2:3528544/1‑100 (MQ=255)
                                                                                       tGCCCGGCTACGT‑GCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTg             >  2:1812918/1‑100 (MQ=255)
                                                                                       tGCCCGGCTACGT‑GCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTATCTTCATCTTCCTCGACGACTGGCGCAGTg             >  1:738224/1‑100 (MQ=255)
                                                                                        gCCCGGCTACGT‑GCCGGGCACTTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTg             >  1:3026874/1‑99 (MQ=255)
                                                                                        gCCCGGCTACGT‑GCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGc            >  2:1745100/1‑100 (MQ=255)
                                                                                        gCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGc            <  2:2865636/100‑1 (MQ=255)
                                                                                         cccGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCt           >  1:2379090/1‑100 (MQ=255)
                                                                                          ccGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCtg          <  2:3552715/100‑1 (MQ=255)
                                                                                          ccGGC‑ACGTAGCCGGGCAATATTGCAATTAACTCGCCGTTTCAGCCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCtg          >  1:199778/1‑100 (MQ=255)
                                                                                           cGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCtgt         >  1:1677744/1‑100 (MQ=255)
                                                                                            ggC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCtgtg        >  1:2598722/1‑100 (MQ=255)
                                                                                            ggC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCtgtg        <  2:1034272/100‑1 (MQ=255)
                                                                                             gC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGa       >  2:1742555/1‑100 (MQ=255)
                                                                                                aCGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAAt     <  1:2644274/100‑1 (MQ=255)
                                                                                                aCGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAAt     <  2:2379102/100‑1 (MQ=255)
                                                                                                aCGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAAt     <  2:2034538/100‑1 (MQ=255)
                                                                                                   tAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTc  >  2:1043011/1‑100 (MQ=255)
                                                                                                   tAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTc  >  2:504625/1‑100 (MQ=255)
                                                                                                    |                                                                                                  
TGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTC  >  NC_000913/1524364‑1524561

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: ‑