Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,259,686 A→G 17.1% intergenic (‑390/+33) psuK ← / ← fruA pseudouridine kinase/fused fructose‑specific PTS enzymes: IIBcomponent/IIC components

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,259,6860AG17.1% 255.8 / 41.6 114intergenic (‑390/+33)psuK/fruApseudouridine kinase/fused fructose‑specific PTS enzymes: IIBcomponent/IIC components
Reads supporting (aligned to +/- strand):  ref base A (50/43);  new base G (11/8);  total (61/53)
Fisher's exact test for biased strand distribution p-value = 8.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01

Only 100 of 132 total aligned reads displayed.

TGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCAC  >  NC_000913/2259587‑2259784
                                                                                                   |                                                                                                   
tgaACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTacagga                                                                                                     >  2:1562113/1‑100 (MQ=255)
tgaACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTacagga                                                                                                     <  1:1137409/100‑1 (MQ=255)
 gaACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTacaggac                                                                                                    >  2:622065/1‑100 (MQ=255)
 gaACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTacaggac                                                                                                    <  2:2638062/100‑1 (MQ=255)
 gaACGTTTGAAGTTTCTACTGAAAGATAAAACCATCTAAAATCATAAAATTACATTCACAAACCGGCCATTCTCTTTGAGCCAACCAAATGTTacaggac                                                                                                    >  1:2462773/1‑100 (MQ=255)
  aaCGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTacaggaca                                                                                                   >  2:2209818/1‑100 (MQ=255)
     gTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCGCAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTAcaggacagga                                                                                                >  2:1804378/1‑100 (MQ=255)
     gTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTAcaggacagga                                                                                                >  2:2664104/1‑100 (MQ=255)
     gTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTAcaggacagga                                                                                                >  2:1804652/1‑100 (MQ=255)
      tttGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACAGTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGaa                                                                                               <  1:2562084/100‑1 (MQ=255)
            gTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑Ct                                                                                        <  2:2591547/100‑1 (MQ=255)
             tttCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTg                                                                                       <  2:2465197/100‑1 (MQ=255)
             tttCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTg                                                                                       <  1:2678044/100‑1 (MQ=255)
               tCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGcc                                                                                     <  1:2475818/100‑1 (MQ=255)
                cTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGGCA‑GAAATTT‑cctgtcc                                                                                   <  2:1624511/100‑7 (MQ=38)
                 tACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCt                                                                                   >  2:2656578/1‑100 (MQ=255)
                  aCTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTg                                                                                  >  2:1104538/1‑100 (MQ=255)
                   cTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGt                                                                                 >  1:2579378/1‑100 (MQ=255)
                   cTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGt                                                                                 >  2:392901/1‑100 (MQ=255)
                   cTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCACTCATCACATGTTACAGGACAGGAACTTT‑CTGCCCTGt                                                                                 >  1:926114/1‑100 (MQ=255)
                   cTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTACGCTCATCAAATGTTACAGGACAGGAAATTT‑CTCCCCTGt                                                                                 >  1:1834681/1‑100 (MQ=255)
                    tGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTa                                                                                >  2:407311/1‑100 (MQ=255)
                     gAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTaa                                                                               >  1:1050010/1‑100 (MQ=255)
                     gAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTaa                                                                               >  1:481681/1‑100 (MQ=255)
                      aaaGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGGCA‑GAAATTTCCTGTCCTGTAac                                                                              <  1:2219552/100‑1 (MQ=37)
                      aaaGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGGCA‑GAAATTTCCTGTCCTGTAac                                                                              <  2:1610969/100‑1 (MQ=37)
                      aaaGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAac                                                                              <  1:2307145/100‑1 (MQ=255)
                      aaaGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAac                                                                              <  2:2103218/100‑1 (MQ=255)
                      aaaGATAAAACCATTTAAATTCATAAAATTACAGTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACTGGACAGGAAATTT‑CTGCCCTGTAac                                                                              <  1:404596/100‑1 (MQ=255)
                      aaaGATAAAACCATTCAAATTCATAAAATTACATTCACAAACCGGTCATTCCCTTTGCGCTCATCAAATGATACCGGACAGGAAATTT‑CTGCCCTGTAac                                                                              >  1:97205/1‑100 (MQ=255)
                       aaGATAAACCCATTTCAATTCATACCATTACACTCACAAACCCGTTATTCTCTTTGCGCTCATCAAATGTTACAGCTCAGAAATTTC‑CTGTCCTGTAac                                                                              <  1:268896/99‑1 (MQ=19)
                          aTAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAacacac                                                                          >  2:839269/1‑100 (MQ=255)
                           tAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACAcc                                                                         >  2:1469041/1‑100 (MQ=255)
                            aaaaCCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCt                                                                        >  2:2792576/1‑100 (MQ=255)
                             aaaCCATTTAAATTCCTAAAATAACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGGCAGAAATTTC‑CTGTCCTGTAacatttg                                                                         <  1:1133348/98‑5 (MQ=21)
                              aaCCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGGCA‑GAAATTTCCTGTCCTGTAacatttgatg                                                                      <  2:2454940/100‑8 (MQ=21)
                              aaCCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCttt                                                                      >  1:2704084/1‑100 (MQ=255)
                              aaCCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCttt                                                                      <  2:449633/100‑1 (MQ=255)
                                ccATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTtta                                                                    >  2:206901/1‑100 (MQ=255)
                                 cATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTttat                                                                   <  2:2898697/100‑1 (MQ=255)
                                  aTTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTttatt                                                                  >  1:606416/1‑100 (MQ=255)
                                  aTTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTttatt                                                                  <  1:812699/100‑1 (MQ=255)
                                   tttAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTttatta                                                                 <  1:1136522/100‑1 (MQ=255)
                                     tAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACg                                                               >  2:257405/1‑100 (MQ=255)
                                      aaaTTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACgc                                                              >  1:659694/1‑100 (MQ=255)
                                       aaTTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGATCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACgct                                                             >  1:687891/1‑100 (MQ=255)
                                        aTTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACgctg                                                            >  1:1439802/1‑100 (MQ=255)
                                        aTTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACgct                                                             >  1:878952/1‑99 (MQ=255)
                                              aaaaTTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCg                                                      <  1:1284041/100‑1 (MQ=255)
                                                    aCATTCGCAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTg                                                <  1:1804372/100‑1 (MQ=255)
                                                    aCATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTg                                                <  2:2454488/100‑1 (MQ=255)
                                                    aCATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTg                                                <  1:1804646/100‑1 (MQ=255)
                                                     cATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGc                                               >  2:163840/1‑100 (MQ=255)
                                                      aTTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCCGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCg                                              <  1:1562109/100‑1 (MQ=255)
                                                         cacaAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTcc                                           >  1:1974319/1‑100 (MQ=255)
                                                         cacaAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTcc                                           >  2:2778977/1‑100 (MQ=255)
                                                         cacaAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTcc                                           <  2:76110/100‑1 (MQ=255)
                                                             aaCCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACtt                                       >  2:321973/1‑100 (MQ=255)
                                                               ccGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTcc                                     >  1:1372472/1‑100 (MQ=255)
                                                               ccGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTcc                                     <  2:651791/100‑1 (MQ=255)
                                                                cGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCg                                    <  2:193950/100‑1 (MQ=255)
                                                                 ggTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCgg                                   <  1:769994/100‑1 (MQ=255)
                                                                 ggTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCgg                                   >  2:2316617/1‑100 (MQ=255)
                                                                 ggTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCg                                    >  1:540200/1‑99 (MQ=255)
                                                                    cATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGAc                                 <  1:220536/99‑1 (MQ=255)
                                                                    cATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACg                                >  2:273407/1‑100 (MQ=255)
                                                                        ctctTTGCGCTCATCAAATGTTACAGGACAGGACATTT‑CTGGCCGGTAACACACCTTTTATTACGCTGCTTTCCCTACTGCGTCCACTTCCGGCCGttt                             >  1:1205134/1‑99 (MQ=255)
                                                                        ctctTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTc                            >  1:2795586/1‑100 (MQ=255)
                                                                          ctTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGAc                                 >  1:1339260/1‑93 (MQ=255)
                                                                          ctTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAg                          >  2:1922057/1‑100 (MQ=255)
                                                                          ctTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAg                          >  1:565071/1‑100 (MQ=255)
                                                                          ctTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAg                          >  1:1894926/1‑100 (MQ=255)
                                                                           tttGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAgg                         <  1:155567/100‑1 (MQ=255)
                                                                               cgcTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAgg                     <  2:2271092/100‑1 (MQ=255)
                                                                                gcTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGc                    >  2:1135115/1‑100 (MQ=255)
                                                                                 cTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGCACGTTTCAGGAAGGCa                   >  1:436539/1‑100 (MQ=255)
                                                                                   catcaAATGTTATAGGTCAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATa                 <  1:2440475/100‑1 (MQ=255)
                                                                                    atcaAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAg                <  2:2907277/100‑1 (MQ=255)
                                                                                     tcaAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTATTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAg                >  1:547728/1‑99 (MQ=255)
                                                                                      ggtgtGTTACAGGGCA‑GAAATTTCCTGTCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCAACTTCCGGACGTTTCAGGAAGGCATAGGc              >  1:1226207/5‑100 (MQ=25)
                                                                                      caAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGc              >  2:1316516/1‑100 (MQ=255)
                                                                                      caAATGTTACAGGACAGGAAATCT‑CTGCCCTATAACAGACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGc              <  1:1957125/100‑1 (MQ=255)
                                                                                         gtGTTACAGGGCAGAAATTTC‑CTGTCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCaa           >  1:922888/2‑100 (MQ=37)
                                                                                         gtGTTACAGGGCA‑GAAATTTCCTGTCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCaa           >  2:1927363/2‑100 (MQ=37)
                                                                                          tGTTACAGGGCA‑GAAATTTCCTGTCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCaaa          >  2:2314805/1‑100 (MQ=37)
                                                                                           gTTACAGGGCAGAAATTTC‑CTGTCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAAc         >  2:648614/1‑100 (MQ=37)
                                                                                           gTTACAGGGCAGAAATTTC‑CTGTCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAAc         >  1:2165507/1‑100 (MQ=37)
                                                                                           gTTACAGGGCA‑GAAATTTCCTGTCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAAc         >  2:2479515/1‑100 (MQ=37)
                                                                                           gTTACAGGGCA‑GAAATTTCCTGTCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAAc         >  1:1929168/1‑100 (MQ=37)
                                                                                           gTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAAc         >  1:2790354/1‑100 (MQ=255)
                                                                                             tACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAAccc       >  1:809360/1‑100 (MQ=255)
                                                                                              acagggca‑gaAATTTCCTGTCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCg      >  2:180992/1‑100 (MQ=38)
                                                                                              acaggacaggAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCg      >  1:878004/1‑100 (MQ=255)
                                                                                              acaggacaggAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCg      <  1:2864650/100‑1 (MQ=255)
                                                                                               caggacaggaAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGc     >  1:2439577/1‑100 (MQ=255)
                                                                                                aggacaggaAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGcc    <  1:205820/100‑1 (MQ=255)
                                                                                                aggacaggaAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGcc    <  2:944717/100‑1 (MQ=255)
                                                                                                 ggacaggaAATTT‑CTGCCCTGTAACACACCTTTTATTACTCTTTTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGATGGCATAGGCCAAACCCGcca   >  1:1647336/1‑100 (MQ=255)
                                                                                                 ggacaggaAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGcca   >  2:2654059/1‑100 (MQ=255)
                                                                                                  gacaggaAAGTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACGTCCGCACGTTTCAGGAAGGCATAGGCCAAACCCGccac  <  1:1793620/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
TGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCAC  >  NC_000913/2259587‑2259784

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