Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,259,700 C→T 16.3% intergenic (‑404/+19) 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,7000CT16.3% 321.4 / 45.2 135intergenic (‑404/+19)psuK/fruApseudouridine kinase/fused fructose‑specific PTS enzymes: IIBcomponent/IIC components
Reads supporting (aligned to +/- strand):  ref base C (63/49);  new base T (12/10);  total (76/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.32e-01

Only 100 of 138 total aligned reads displayed.

TTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAAT  >  NC_000913/2259601‑2259799
                                                                                                    |                                                                                                   
 tCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGcc                                                                                                    <  1:2475818/100‑1 (MQ=255)
   tACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCt                                                                                                  >  2:2656578/1‑100 (MQ=255)
    aCTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTg                                                                                                 >  2:1104538/1‑100 (MQ=255)
     cTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTGTCTTTTCAGTCATCAAATCTTACAGGACAGGAAATTG‑CTTCCCTGt                                                                                                >  1:1833959/1‑100 (MQ=255)
     cTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGt                                                                                                >  2:392901/1‑100 (MQ=255)
     cTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCACTCATCACATGTTACAGGACAGGAACTTT‑CTGCCCTGt                                                                                                >  1:926114/1‑100 (MQ=255)
      tGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTa                                                                                               >  2:407311/1‑100 (MQ=255)
       gAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTaa                                                                                              >  1:1050010/1‑100 (MQ=255)
        aaaGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGGCAGAAATTTC‑CTGTCCTGTAac                                                                                             <  2:2499766/100‑1 (MQ=37)
        aaaGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGGCA‑GAAATTTCCTGTCCTGTAac                                                                                             <  2:1610969/100‑1 (MQ=37)
        aaaGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTa                                                                                               <  1:2664880/98‑1 (MQ=255)
        aaaGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAac                                                                                             <  2:2865657/100‑1 (MQ=255)
        aaaGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAac                                                                                             <  2:2103218/100‑1 (MQ=255)
        aaaGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAac                                                                                             <  2:1757392/100‑1 (MQ=255)
        aaaGATAAAACCATTTAAATTCATAAAATTACAGTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACTGGACAGGAAATTT‑CTGCCCTGTAac                                                                                             <  1:404596/100‑1 (MQ=255)
         aaGATAAACCCATTTCAATTCATACCATTACACTCACAAACCCGTTATTCTCTTTGCGCTCATCAAATGTTACAGCTCAGAAATTTC‑CTGTCCTGTAac                                                                                             <  1:268896/99‑1 (MQ=19)
            aTAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGGCA‑GAAATTTCCTGTCCTGTAacattt                                                                                         <  2:2646691/100‑4 (MQ=25)
            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)
               aaaCCATCTAATTTCATAATATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGGCAGAAATTTC‑CTGTCCTGTAacatttgat                                                                                      <  1:2473290/100‑7 (MQ=21)
                aaCCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGGCA‑GAAATTTCCTGTCCTGTAacatttgatg                                                                                     <  2:2454940/100‑8 (MQ=21)
                aaCCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCttt                                                                                     >  1:2704084/1‑100 (MQ=255)
                  ccATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTtta                                                                                   >  2:206901/1‑100 (MQ=255)
                   cATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTttat                                                                                  <  2:2898697/100‑1 (MQ=255)
                    aTTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTttatt                                                                                 <  1:936160/100‑1 (MQ=255)
                    aTTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTttatt                                                                                 >  1:606416/1‑100 (MQ=255)
                    aTTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTttatt                                                                                 >  2:2266637/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)
                                aaaaTTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCg                                                                     <  1:1284041/100‑1 (MQ=255)
                                 aaaTTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGc                                                                    <  1:166219/100‑1 (MQ=255)
                                    ttACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTAc                                                                 <  2:1992067/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                                                          <  1:829870/100‑1 (MQ=255)
                                           cacaAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTcc                                                          <  2:511170/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)
                                                      cATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGAc                                                <  1:220536/99‑1 (MQ=255)
                                                      cATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACg                                               >  2:273407/1‑100 (MQ=255)
                                                         tctcTTTGCCTTCAGCAAATGTAACAGGACAAGAAATTT‑CTGCTCTGGAACACACCTTTTATTACGCTGCTTTTTCTATTGCGTCCACTTCCGGAc                                                <  1:484152/96‑1 (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:1894926/1‑100 (MQ=255)
                                                            ctTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAg                                         >  1:565071/1‑100 (MQ=255)
                                                             tttGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAgg                                        <  1:155567/100‑1 (MQ=255)
                                                                gcgcTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAg                                     <  2:1593620/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)
                                                                       tcaAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTATTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAg                               >  1:547728/1‑99 (MQ=255)
                                                                       tcaAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAgg                              <  1:2245267/100‑1 (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)
                                                                            tGTTACAGGGCAGAAATTTC‑CTGTCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCaaa                         >  2:659947/1‑100 (MQ=37)
                                                                            tGTTACAGGGCA‑GAAATTTCCTGTCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCaaa                         >  2:2314805/1‑100 (MQ=37)
                                                                             gTTACAGGGCAGAAATTTC‑CTGTCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAAc                        >  1:2165507/1‑100 (MQ=37)
                                                                             gTTACAGGGCA‑GAAATTTCCTGTCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAAc                        >  1:1929168/1‑100 (MQ=37)
                                                                                gctcatgtggaAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCg                     <  2:2540639/92‑1 (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)
                                                                                  aggacaggaAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGcc                   <  2:944717/100‑1 (MQ=255)
                                                                                  aggacaggaAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGcc                   <  1:205820/100‑1 (MQ=255)
                                                                                   ggacaggaAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGcca                  >  2:2654059/1‑100 (MQ=255)
                                                                                    gacaggaAAGTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACGTCCGCACGTTTCAGGAAGGCATAGGCCAAACCCGccac                 <  1:1793620/100‑1 (MQ=255)
                                                                                      gcagaAATTTCCTGTCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCaccacc                <  1:1927355/97‑1 (MQ=255)
                                                                                      caggaAATTT‑CTGCCCTGTAACACACCTTTTATTACGGTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGccacca               >  1:1548131/1‑100 (MQ=255)
                                                                                       aggaAATTT‑CTTCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGccaccan              <  1:400589/100‑2 (MQ=255)
                                                                                       aggaAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGcc                   <  1:2266627/95‑1 (MQ=255)
                                                                                        ggaAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAATGCATAGGCCAAACCCGCCACCAGc             >  1:490611/1‑100 (MQ=255)
                                                                                        ggaAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGc             >  1:1649430/1‑100 (MQ=255)
                                                                                        ggaAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGc             >  1:490782/1‑100 (MQ=255)
                                                                                        ggaAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGc             >  2:1163864/1‑100 (MQ=255)
                                                                                        ggaAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGc             >  2:2187034/1‑100 (MQ=255)
                                                                                         gaAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCg            <  2:275744/100‑1 (MQ=255)
                                                                                          aaaTTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGt           <  1:805958/100‑1 (MQ=255)
                                                                                             ttt‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTAcc        <  2:565865/100‑1 (MQ=255)
                                                                                               t‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCgg      <  2:723356/100‑1 (MQ=255)
                                                                                                ccTGTCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCgg      <  1:2422458/99‑1 (MQ=255)
                                                                                                 cTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGc     >  1:2291883/1‑100 (MQ=255)
                                                                                                 cTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGc     >  2:2319084/1‑100 (MQ=255)
                                                                                                    cccTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCaat  <  1:2347351/100‑1 (MQ=255)
                                                                                                    |                                                                                                   
TTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAAT  >  NC_000913/2259601‑2259799

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