Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,524,458:1 +T 15.0% intergenic (+90/+23) 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,4581.T15.0% 421.1 / 83.4 133intergenic (+90/+23)yncE/ansPATP‑binding protein, periplasmic, function unknown/L‑asparagine transporter
Reads supporting (aligned to +/- strand):  ref base . (72/41);  new base T (12/8);  total (84/49)
Fisher's exact test for biased strand distribution p-value = 8.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.94e-01

Only 100 of 136 total aligned reads displayed.

TGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTG  >  NC_000913/1524359‑1524555
                                                                                                    |                                                                                                  
 gcgcTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑a                                                                                                   >  1:2454388/1‑100 (MQ=255)
  cgcTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑Ac                                                                                                  >  2:435244/1‑100 (MQ=255)
    cTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTCATCTGAATTTGCCCGGC‑CCGt                                                                                                >  2:2173147/1‑100 (MQ=255)
      gTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAg                                                                                              >  2:1639346/1‑100 (MQ=255)
      gTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAg                                                                                              >  1:3235363/1‑100 (MQ=255)
       tAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGc                                                                                             <  2:3741069/100‑1 (MQ=255)
         aTAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCg                                                                                           >  1:3719343/1‑100 (MQ=255)
          tAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCgg                                                                                          >  1:1703311/1‑100 (MQ=255)
          tAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCgg                                                                                          >  1:1943764/1‑100 (MQ=255)
           aGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCggg                                                                                         <  2:1368276/100‑1 (MQ=255)
           aGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCggg                                                                                         <  1:1674400/100‑1 (MQ=255)
             atatGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCTACGT‑GCCGGGC‑a                                                                                      <  2:1460090/100‑1 (MQ=255)
              tatGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCTACGT‑GCCGGGC‑aa                                                                                     <  2:1914584/100‑1 (MQ=255)
                tGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AAtt                                                                                   >  1:679050/1‑100 (MQ=255)
                tGGCGAGGGGGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCTACGT‑GCCGGGC‑aaat                                                                                   <  2:1752321/100‑3 (MQ=255)
                 ggCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AAttt                                                                                  >  1:1767264/1‑100 (MQ=255)
                  gCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCTACGT‑GCCGGGC‑AAattc                                                                                 <  1:1753191/100‑2 (MQ=38)
                   cGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTg                                                                                >  1:520998/1‑100 (MQ=255)
                   cGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTg                                                                                <  2:520992/100‑1 (MQ=255)
                   cGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTg                                                                                >  2:2734671/1‑100 (MQ=255)
                     aGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCTACGT‑GCCGGGC‑AAatt                                                                                  <  2:1692037/96‑1 (MQ=38)
                      gggTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCaa                                                                             >  2:1035589/1‑100 (MQ=255)
                       ggTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAAt                                                                            >  1:552324/1‑100 (MQ=255)
                        gTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAAtt                                                                           <  1:2674476/100‑1 (MQ=255)
                         tGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTa                                                                          >  1:3354506/1‑100 (MQ=255)
                          gAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTaa                                                                         >  2:121541/1‑100 (MQ=255)
                          gAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTAGCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTaa                                                                         >  2:2350232/1‑100 (MQ=255)
                           acaaCGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAAc                                                                        >  1:3534578/3‑100 (MQ=255)
                           aaaaCGGTGAGACTTGCCTGCCTTTTATCGCGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTCCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAAc                                                                        <  2:2298252/100‑1 (MQ=255)
                           aaaaCGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAAc                                                                        >  2:3435652/1‑100 (MQ=255)
                              aCGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCg                                                                     >  2:2624596/1‑100 (MQ=255)
                               cGGTGAGACTTGCCTGCCTTTTAAAAAAGCGTATAAGGCACTTTCGTCAGTAATCTGAATTTACGCAGC‑ACGTACCCGGGC‑AAATTTGCAATTAAc                                                                        >  2:1313504/1‑96 (MQ=255)
                                ggTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGcc                                                                   >  2:1100922/1‑100 (MQ=255)
                                     gaCTTGTCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTc                                                              >  1:3447163/1‑100 (MQ=255)
                                     gaCTTGCCTGCCTATTATCACGGCGTATCCGGCACTTTAGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTc                                                              >  2:2363269/1‑100 (MQ=255)
                                        ttGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAgg                                                           <  1:292295/100‑1 (MQ=255)
                                        ttGCCTGACTTTTATCACGCCGCATCCGGCACCTTCGTCAGTAGTCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAgg                                                           >  1:3206682/1‑100 (MQ=255)
                                        ttGCCTGACTTTTATCACGCCGCATCCGGCACCTTCGTCAGTAGTCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAgg                                                           >  1:1291285/1‑100 (MQ=255)
                                         tgcctgccTTTTATCACGGCGTATCCGGGACTTGCGTCAGTAATCTGAATTTGCTCGGC‑ACGTAGCCGGGC‑AATGTTGCAATTAACTCGCCGTTTCAGGc                                                          <  2:3233790/100‑1 (MQ=255)
                                          gcctgcctTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCt                                                         >  1:77137/1‑100 (MQ=255)
                                           cctgcctTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCtt                                                        >  1:897798/1‑100 (MQ=255)
                                               cctTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAac                                                    <  1:1944324/100‑1 (MQ=255)
                                               cctTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAac                                                    >  1:1133319/1‑100 (MQ=255)
                                                ctTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAaca                                                   <  2:1973911/100‑1 (MQ=255)
                                                ctTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAaca                                                   >  2:2702328/1‑100 (MQ=255)
                                                 ttttATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAacac                                                  >  1:1076284/1‑100 (MQ=255)
                                                    tATCACGGCGTATCCGGCACTTTCGTCAGCAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAAt                                               <  2:2298874/100‑1 (MQ=255)
                                                     aTCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAAtt                                              >  2:3094419/1‑100 (MQ=255)
                                                        aCGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCACTTAACTCGCCGTTTCAGGCTTAAACACCATTTcc                                           >  1:2945295/1‑100 (MQ=255)
                                                        aCGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTcc                                           >  2:2105545/1‑100 (MQ=255)
                                                          ggCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATCTCCTg                                         >  2:2753387/1‑100 (MQ=255)
                                                           gCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGt                                        <  2:834315/100‑1 (MQ=255)
                                                             gTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGtt                                       <  2:1676815/99‑1 (MQ=255)
                                                              tATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGtttt                                     >  1:743302/1‑100 (MQ=255)
                                                                 ccGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTttct                                  >  1:3091279/1‑100 (MQ=255)
                                                                  cGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTttctt                                 >  1:1438859/1‑100 (MQ=255)
                                                                  cGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTttctt                                 <  2:2207680/100‑1 (MQ=255)
                                                                  cGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTttctt                                 >  1:869272/1‑100 (MQ=255)
                                                                   ggCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTttcttc                                >  2:2625915/1‑100 (MQ=255)
                                                                    gCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCa                               <  1:304090/100‑1 (MQ=255)
                                                                    gCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCa                               <  2:3458410/100‑1 (MQ=255)
                                                                    gCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCa                               <  1:3054724/100‑1 (MQ=255)
                                                                     cACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCAt                              <  2:18881/100‑1 (MQ=255)
                                                                     cACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCAt                              >  2:3658230/1‑100 (MQ=255)
                                                                       cTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCt                            <  1:3031146/100‑1 (MQ=255)
                                                                        tttCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCtt                           <  1:2122988/100‑1 (MQ=255)
                                                                        tttCGTCAGTAAGCTGAATTTGCCTGGA‑ACGTAGCCGGGC‑AATTTTTCTATTAACTCTCCGTTTCATGCTTAAACACAATTTCCTGTTTTTCTTCATCt                            <  2:3362696/99‑1 (MQ=255)
                                                                         ttCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGTGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTc                          >  2:1755709/1‑100 (MQ=255)
                                                                           cGTCAGTAATCTGAATTTGCCCGGC‑AAGTAGCCGGGA‑AATTTTGCTATTAACTCGCCGTTTCGGGCGTAAACAAAATATCCTGTTGTTATTCTTCTTACt                        >  2:3648514/1‑100 (MQ=255)
                                                                            gTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTc                       >  2:604545/1‑100 (MQ=255)
                                                                            gTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTc                       <  1:231560/100‑1 (MQ=255)
                                                                            gTCAGTAATCTGAATAAGCCCGGC‑ACGTAGCCGTAC‑AGTGTTGCAATTAACTCGCCCTTTCACGCTTAAACACCAGTTCCTTTTATTCTTCATATTCCt                        >  2:1590793/1‑99 (MQ=255)
                                                                               aGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTcgac                    <  1:2390859/100‑1 (MQ=255)
                                                                                gTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTcgacg                   >  2:2323112/1‑100 (MQ=255)
                                                                                gTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTcgacg                   <  2:1816554/100‑1 (MQ=255)
                                                                                   aTCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACt                >  1:231846/1‑100 (MQ=255)
                                                                                    tCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTg               >  1:876344/1‑100 (MQ=255)
                                                                                     gcaaaattGCCCGGCTACGT‑GCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTgg              >  1:152639/7‑100 (MQ=38)
                                                                                     gcaaaattGCCCGGCTACGT‑GCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTgg              >  1:3619189/7‑100 (MQ=38)
                                                                                      tGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGATGACTGgc             <  2:2969369/100‑1 (MQ=255)
                                                                                      tGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGgc             >  2:1371376/1‑100 (MQ=255)
                                                                                      tGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGgc             >  2:450037/1‑100 (MQ=255)
                                                                                       gAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGgcg            <  1:3049423/100‑1 (MQ=255)
                                                                                        aaattGCCCGGCTACGT‑GCCGGGCAAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGgcg            >  2:3732738/4‑100 (MQ=37)
                                                                                        aaTTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAAATCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGgcgc           <  1:1032857/100‑1 (MQ=255)
                                                                                         aattGCCCGGCTACGT‑GCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCa          >  1:3774744/3‑100 (MQ=255)
                                                                                         aTTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCa          >  1:2000095/1‑100 (MQ=255)
                                                                                          tttGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAg         <  2:420406/100‑1 (MQ=255)
                                                                                            tGCCCGGCTACGT‑GCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCCGTg       >  2:744487/1‑100 (MQ=255)
                                                                                            tGCCCGGCTACGT‑GCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAg         >  2:870612/1‑98 (MQ=255)
                                                                                            tGCCCGGCTACGT‑GCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTg       >  2:3662378/1‑100 (MQ=255)
                                                                                            tGCCCGGCTACGT‑GCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTg       >  1:1820508/1‑100 (MQ=255)
                                                                                            tGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTg       >  2:1010774/1‑100 (MQ=255)
                                                                                             gCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGc      <  2:679043/100‑1 (MQ=255)
                                                                                              cccGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGc      >  2:2432337/1‑99 (MQ=255)
                                                                                              cccGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCt     >  2:517853/1‑100 (MQ=255)
                                                                                              cccGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCt     >  1:2228664/1‑100 (MQ=255)
                                                                                              cccGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCt     >  2:1532790/1‑100 (MQ=255)
                                                                                                cGGC‑ACGTAGCCGTGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGTCTGGCGCAGTGCtgt   <  2:681497/100‑1 (MQ=255)
                                                                                                 ggC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCtgtg  >  1:1464733/1‑100 (MQ=255)
                                                                                                    |                                                                                                  
TGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTG  >  NC_000913/1524359‑1524555

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