Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,524,463 Δ1 bp 15.3% 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.3% 370.4 / 64.2 131intergenic (+95/+18)yncE/ansPATP‑binding protein, periplasmic, function unknown/L‑asparagine transporter
Reads supporting (aligned to +/- strand):  ref base A (70/41);  new base . (12/8);  total (82/49)
Fisher's exact test for biased strand distribution p-value = 8.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

Only 100 of 138 total aligned reads displayed.

GTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCA  >  NC_000913/1524365‑1524562
                                                                                                   |                                                                                                    
gTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAg                                                                                                     >  1:3235363/1‑100 (MQ=255)
gTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAg                                                                                                     >  2:1639346/1‑100 (MQ=255)
 tAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGc                                                                                                    <  2:3741069/100‑1 (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)
      gATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGc                                                                                               >  2:2664741/1‑100 (MQ=255)
       atatGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCTACGT‑GCCGGGC‑a                                                                                             <  2:1460090/100‑1 (MQ=255)
        tatGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCTACGT‑GCCGGGC‑aa                                                                                            <  2:1914584/100‑1 (MQ=255)
        tatGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑a                                                                                             >  1:922664/1‑99 (MQ=255)
          tGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AAtt                                                                                          >  1:679050/1‑100 (MQ=255)
          tGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTAGCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCTACGT‑GCCGGGC‑aaat                                                                                          <  2:504706/100‑3 (MQ=255)
          tGGCGAGGGGGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCTACGT‑GCCGGGC‑aaat                                                                                          <  2:1752321/100‑3 (MQ=255)
            gCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCTACGT‑GCCGGGC‑AAattc                                                                                        <  1:1486943/100‑2 (MQ=38)
            gCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCTACGT‑GCCGGGC‑AAattc                                                                                        <  1:1753191/100‑2 (MQ=38)
             cGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTg                                                                                       >  1:520998/1‑100 (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)
                  gTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAAtt                                                                                  <  1:2674476/100‑1 (MQ=255)
                   tGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTa                                                                                 >  1:3068788/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)
                                  ttGCCTGCCTTTTATCACGGTGTATCCGGCAGTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAg                                                                   <  2:782622/99‑1 (MQ=255)
                                  ttGCCTGACTTTTATCACGCCGCATCCGGCACCTTCGTCAGTAGTCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAgg                                                                  >  1:1291285/1‑100 (MQ=255)
                                   tgcctgccTTTTATCACGGCGTATCCGGGACTTGCGTCAGTAATCTGAATTTGCTCGGC‑ACGTAGCCGGGC‑AATGTTGCAATTAACTCGCCGTTTCAGGc                                                                 <  2:3233790/100‑1 (MQ=255)
                                     cctgcctTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCtt                                                               >  1:897798/1‑100 (MQ=255)
                                         cctTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAac                                                           <  1:2261926/100‑1 (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)
                                           ttttATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAacac                                                         >  1:1076284/1‑100 (MQ=255)
                                               aTCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAAtt                                                     >  2:3094419/1‑100 (MQ=255)
                                               aTCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAAtt                                                     >  1:3253836/1‑100 (MQ=255)
                                                  aCGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCACTTAACTCGCCGTTTCAGGCTTAAACACCATTTcc                                                  >  1:2945295/1‑100 (MQ=255)
                                                  aCGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTcc                                                  >  2:2105545/1‑100 (MQ=255)
                                                  aCGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTcc                                                  >  1:1279323/1‑100 (MQ=255)
                                                  aCGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTcc                                                  >  1:663287/1‑100 (MQ=255)
                                                    ggCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATCTCCTg                                                >  2:2753387/1‑100 (MQ=255)
                                                      cGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACATTTTCCTg                                                <  1:83451/98‑1 (MQ=255)
                                                       gTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGtt                                              <  2:1676815/99‑1 (MQ=255)
                                                        tATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGtttt                                            >  1:743302/1‑100 (MQ=255)
                                                         aTCCGGCACGTTCGGCAGTCATCTGCATTTGCCCGGC‑ACGTAGCCGTGC‑ACTTTGGCAATTAACTCGCCGCTACAGGCCTAAACACAATTTCCTg                                                >  2:3749662/1‑95 (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)
                                                             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)
                                                               cACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCAt                                     <  2:18881/100‑1 (MQ=255)
                                                               cACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCAt                                     >  2:3658230/1‑100 (MQ=255)
                                                                  tttCGTCAGTAAGCTGAATTTGCCTGGA‑ACGTAGCCGGGC‑AATTTTTCTATTAACTCTCCGTTTCATGCTTAAACACAATTTCCTGTTTTTCTTCATCt                                   <  2:3362696/99‑1 (MQ=255)
                                                                   ttCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGTGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTc                                 >  2:1755709/1‑100 (MQ=255)
                                                                      gTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTc                              >  2:604545/1‑100 (MQ=255)
                                                                      gTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTc                              >  2:604667/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:2932053/100‑1 (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:1744379/1‑100 (MQ=255)
                                                                              tCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTg                      >  1:876344/1‑100 (MQ=255)
                                                                               gcaaaattGCCCGGCTACGT‑GCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTgg                     >  1:152639/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:303339/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)
                                                                                      tGCCCGGCTACGT‑GCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCCGTg              >  2:744487/1‑100 (MQ=255)
                                                                                      tGCCCGGCTACGT‑GCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAg                >  2:870612/1‑98 (MQ=255)
                                                                                      tGCCCGGCTACGT‑GCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTg              >  1:1820508/1‑100 (MQ=255)
                                                                                      tGCCCGGCTACGT‑GCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTg              >  2:233331/1‑100 (MQ=255)
                                                                                      tGCCCGGCTACGT‑GCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTg              >  1:3586916/1‑100 (MQ=255)
                                                                                      tGCCCGGCTACGT‑GCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTg              >  2:3662378/1‑100 (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‑AAGTTGGCAATTAACGCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGc             <  1:3822601/99‑1 (MQ=255)
                                                                                          cGGC‑ACGTAGCCGTGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGTCTGGCGCAGTGCtgt          <  2:681497/100‑1 (MQ=255)
                                                                                                cGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAAtt     >  1:3673092/1‑100 (MQ=255)
                                                                                                cGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTCTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAAtt     >  1:3673182/1‑100 (MQ=255)
                                                                                                   aGCCGGGC‑AATTTTGTAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTCTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCa  <  1:2037339/100‑1 (MQ=255)
                                                                                                   aGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGTGCAGTGCGTTGAATTTCa  >  2:2497649/1‑100 (MQ=255)
                                                                                                   aGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCa  <  1:3298662/100‑1 (MQ=255)
                                                                                                   aGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCa  <  1:2664760/100‑1 (MQ=255)
                                                                                                   aGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCa  <  1:848356/100‑1 (MQ=255)
                                                                                                   aGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCa  <  2:1252802/100‑1 (MQ=255)
                                                                                                   aGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCa  <  2:876335/100‑1 (MQ=255)
                                                                                                   |                                                                                                    
GTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGC‑AATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCA  >  NC_000913/1524365‑1524562

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