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.

TGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCAC  >  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)
            gTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCt                                                                                        <  2:2591547/100‑1 (MQ=255)
             tttCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTg                                                                                       <  2:2465197/100‑1 (MQ=255)
             tttCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTg                                                                                       <  1:2678044/100‑1 (MQ=255)
               tCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGcc                                                                                     <  1:2475818/100‑1 (MQ=255)
                cTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGGCAGAAATTTcctgtcc                                                                                   <  2:1624511/100‑7 (MQ=38)
                 tACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCt                                                                                   >  2:2656578/1‑100 (MQ=255)
                  aCTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTg                                                                                  >  2:1104538/1‑100 (MQ=255)
                   cTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGt                                                                                 >  1:2579378/1‑100 (MQ=255)
                   cTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGt                                                                                 >  2:392901/1‑100 (MQ=255)
                   cTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCACTCATCACATGTTACAGGACAGGAACTTTCTGCCCTGt                                                                                 >  1:926114/1‑100 (MQ=255)
                   cTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTACGCTCATCAAATGTTACAGGACAGGAAATTTCTCCCCTGt                                                                                 >  1:1834681/1‑100 (MQ=255)
                    tGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTa                                                                                >  2:407311/1‑100 (MQ=255)
                     gAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTaa                                                                               >  1:1050010/1‑100 (MQ=255)
                     gAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTaa                                                                               >  1:481681/1‑100 (MQ=255)
                      aaaGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGGCAGAAATTTCCTGTCCTGTAac                                                                              <  1:2219552/100‑1 (MQ=37)
                      aaaGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGGCAGAAATTTCCTGTCCTGTAac                                                                              <  2:1610969/100‑1 (MQ=37)
                      aaaGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAac                                                                              <  1:2307145/100‑1 (MQ=255)
                      aaaGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAac                                                                              <  2:2103218/100‑1 (MQ=255)
                      aaaGATAAAACCATTTAAATTCATAAAATTACAGTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACTGGACAGGAAATTTCTGCCCTGTAac                                                                              <  1:404596/100‑1 (MQ=255)
                      aaaGATAAAACCATTCAAATTCATAAAATTACATTCACAAACCGGTCATTCCCTTTGCGCTCATCAAATGATACCGGACAGGAAATTTCTGCCCTGTAac                                                                              >  1:97205/1‑100 (MQ=255)
                       aaGATAAACCCATTTCAATTCATACCATTACACTCACAAACCCGTTATTCTCTTTGCGCTCATCAAATGTTACAGCTCAGAAATTTCCTGTCCTGTAac                                                                              <  1:268896/99‑1 (MQ=19)
                          aTAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAacacac                                                                          >  2:839269/1‑100 (MQ=255)
                           tAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACAcc                                                                         >  2:1469041/1‑100 (MQ=255)
                            aaaaCCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCt                                                                        >  2:2792576/1‑100 (MQ=255)
                             aaaCCATTTAAATTCCTAAAATAACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGGCAGAAATTTCCTGTCCTGTAacatttg                                                                         <  1:1133348/98‑5 (MQ=21)
                              aaCCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGGCAGAAATTTCCTGTCCTGTAacatttgatg                                                                      <  2:2454940/100‑8 (MQ=21)
                              aaCCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCttt                                                                      >  1:2704084/1‑100 (MQ=255)
                              aaCCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCttt                                                                      <  2:449633/100‑1 (MQ=255)
                                ccATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTtta                                                                    >  2:206901/1‑100 (MQ=255)
                                 cATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTttat                                                                   <  2:2898697/100‑1 (MQ=255)
                                  aTTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTttatt                                                                  >  1:606416/1‑100 (MQ=255)
                                  aTTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTttatt                                                                  <  1:812699/100‑1 (MQ=255)
                                   tttAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTttatta                                                                 <  1:1136522/100‑1 (MQ=255)
                                     tAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACg                                                               >  2:257405/1‑100 (MQ=255)
                                      aaaTTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACgc                                                              >  1:659694/1‑100 (MQ=255)
                                       aaTTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGATCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACgct                                                             >  1:687891/1‑100 (MQ=255)
                                        aTTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACgctg                                                            >  1:1439802/1‑100 (MQ=255)
                                        aTTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACgct                                                             >  1:878952/1‑99 (MQ=255)
                                              aaaaTTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCg                                                      <  1:1284041/100‑1 (MQ=255)
                                                    aCATTCGCAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTg                                                <  1:1804372/100‑1 (MQ=255)
                                                    aCATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTg                                                <  2:2454488/100‑1 (MQ=255)
                                                    aCATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTg                                                <  1:1804646/100‑1 (MQ=255)
                                                     cATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGc                                               >  2:163840/1‑100 (MQ=255)
                                                      aTTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCCGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCg                                              <  1:1562109/100‑1 (MQ=255)
                                                         cacaAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTcc                                           >  1:1974319/1‑100 (MQ=255)
                                                         cacaAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTcc                                           >  2:2778977/1‑100 (MQ=255)
                                                         cacaAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTcc                                           <  2:76110/100‑1 (MQ=255)
                                                             aaCCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACtt                                       >  2:321973/1‑100 (MQ=255)
                                                               ccGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTcc                                     >  1:1372472/1‑100 (MQ=255)
                                                               ccGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTcc                                     <  2:651791/100‑1 (MQ=255)
                                                                cGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCg                                    <  2:193950/100‑1 (MQ=255)
                                                                 ggTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCgg                                   <  1:769994/100‑1 (MQ=255)
                                                                 ggTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCgg                                   >  2:2316617/1‑100 (MQ=255)
                                                                 ggTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCg                                    >  1:540200/1‑99 (MQ=255)
                                                                    cATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGAc                                 <  1:220536/99‑1 (MQ=255)
                                                                    cATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACg                                >  2:273407/1‑100 (MQ=255)
                                                                        ctctTTGCGCTCATCAAATGTTACAGGACAGGACATTTCTGGCCGGTAACACACCTTTTATTACGCTGCTTTCCCTACTGCGTCCACTTCCGGCCGttt                             >  1:1205134/1‑99 (MQ=255)
                                                                        ctctTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTc                            >  1:2795586/1‑100 (MQ=255)
                                                                          ctTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGAc                                 >  1:1339260/1‑93 (MQ=255)
                                                                          ctTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAg                          >  2:1922057/1‑100 (MQ=255)
                                                                          ctTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAg                          >  1:565071/1‑100 (MQ=255)
                                                                          ctTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAg                          >  1:1894926/1‑100 (MQ=255)
                                                                           tttGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAgg                         <  1:155567/100‑1 (MQ=255)
                                                                               cgcTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAgg                     <  2:2271092/100‑1 (MQ=255)
                                                                                gcTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGc                    >  2:1135115/1‑100 (MQ=255)
                                                                                 cTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGCACGTTTCAGGAAGGCa                   >  1:436539/1‑100 (MQ=255)
                                                                                   catcaAATGTTATAGGTCAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATa                 <  1:2440475/100‑1 (MQ=255)
                                                                                    atcaAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAg                <  2:2907277/100‑1 (MQ=255)
                                                                                     tcaAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTATTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAg                >  1:547728/1‑99 (MQ=255)
                                                                                      ggtgtGTTACAGGGCAGAAATTTCCTGTCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCAACTTCCGGACGTTTCAGGAAGGCATAGGc              >  1:1226207/5‑100 (MQ=25)
                                                                                      caAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGc              >  2:1316516/1‑100 (MQ=255)
                                                                                      caAATGTTACAGGACAGGAAATCTCTGCCCTATAACAGACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGc              <  1:1957125/100‑1 (MQ=255)
                                                                                         gtGTTACAGGGCAGAAATTTCCTGTCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCaa           >  1:922888/2‑100 (MQ=37)
                                                                                         gtGTTACAGGGCAGAAATTTCCTGTCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCaa           >  2:1927363/2‑100 (MQ=37)
                                                                                          tGTTACAGGGCAGAAATTTCCTGTCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCaaa          >  2:2314805/1‑100 (MQ=37)
                                                                                           gTTACAGGGCAGAAATTTCCTGTCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAAc         >  2:648614/1‑100 (MQ=37)
                                                                                           gTTACAGGGCAGAAATTTCCTGTCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAAc         >  1:2165507/1‑100 (MQ=37)
                                                                                           gTTACAGGGCAGAAATTTCCTGTCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAAc         >  2:2479515/1‑100 (MQ=37)
                                                                                           gTTACAGGGCAGAAATTTCCTGTCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAAc         >  1:1929168/1‑100 (MQ=37)
                                                                                           gTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAAc         >  1:2790354/1‑100 (MQ=255)
                                                                                             tACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAAccc       >  1:809360/1‑100 (MQ=255)
                                                                                              acagggcagaAATTTCCTGTCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCg      >  2:180992/1‑100 (MQ=38)
                                                                                              acaggacaggAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCg      >  1:878004/1‑100 (MQ=255)
                                                                                              acaggacaggAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCg      <  1:2864650/100‑1 (MQ=255)
                                                                                               caggacaggaAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGc     >  1:2439577/1‑100 (MQ=255)
                                                                                                aggacaggaAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGcc    <  1:205820/100‑1 (MQ=255)
                                                                                                aggacaggaAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGcc    <  2:944717/100‑1 (MQ=255)
                                                                                                 ggacaggaAATTTCTGCCCTGTAACACACCTTTTATTACTCTTTTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGATGGCATAGGCCAAACCCGcca   >  1:1647336/1‑100 (MQ=255)
                                                                                                 ggacaggaAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGcca   >  2:2654059/1‑100 (MQ=255)
                                                                                                  gacaggaAAGTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACGTCCGCACGTTTCAGGAAGGCATAGGCCAAACCCGccac  <  1:1793620/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
TGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCAC  >  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: