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.

TTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAAT  >  NC_000913/2259601‑2259799
                                                                                                    |                                                                                                   
 tCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGcc                                                                                                    <  1:2475818/100‑1 (MQ=255)
   tACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCt                                                                                                  >  2:2656578/1‑100 (MQ=255)
    aCTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTg                                                                                                 >  2:1104538/1‑100 (MQ=255)
     cTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTGTCTTTTCAGTCATCAAATCTTACAGGACAGGAAATTGCTTCCCTGt                                                                                                >  1:1833959/1‑100 (MQ=255)
     cTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGt                                                                                                >  2:392901/1‑100 (MQ=255)
     cTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCACTCATCACATGTTACAGGACAGGAACTTTCTGCCCTGt                                                                                                >  1:926114/1‑100 (MQ=255)
      tGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTa                                                                                               >  2:407311/1‑100 (MQ=255)
       gAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTaa                                                                                              >  1:1050010/1‑100 (MQ=255)
        aaaGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGGCAGAAATTTCCTGTCCTGTAac                                                                                             <  2:2499766/100‑1 (MQ=37)
        aaaGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGGCAGAAATTTCCTGTCCTGTAac                                                                                             <  2:1610969/100‑1 (MQ=37)
        aaaGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTa                                                                                               <  1:2664880/98‑1 (MQ=255)
        aaaGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAac                                                                                             <  2:2865657/100‑1 (MQ=255)
        aaaGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAac                                                                                             <  2:2103218/100‑1 (MQ=255)
        aaaGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAac                                                                                             <  2:1757392/100‑1 (MQ=255)
        aaaGATAAAACCATTTAAATTCATAAAATTACAGTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACTGGACAGGAAATTTCTGCCCTGTAac                                                                                             <  1:404596/100‑1 (MQ=255)
         aaGATAAACCCATTTCAATTCATACCATTACACTCACAAACCCGTTATTCTCTTTGCGCTCATCAAATGTTACAGCTCAGAAATTTCCTGTCCTGTAac                                                                                             <  1:268896/99‑1 (MQ=19)
            aTAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGGCAGAAATTTCCTGTCCTGTAacattt                                                                                         <  2:2646691/100‑4 (MQ=25)
            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)
               aaaCCATCTAATTTCATAATATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGGCAGAAATTTCCTGTCCTGTAacatttgat                                                                                      <  1:2473290/100‑7 (MQ=21)
                aaCCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGGCAGAAATTTCCTGTCCTGTAacatttgatg                                                                                     <  2:2454940/100‑8 (MQ=21)
                aaCCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCttt                                                                                     >  1:2704084/1‑100 (MQ=255)
                  ccATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTtta                                                                                   >  2:206901/1‑100 (MQ=255)
                   cATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTttat                                                                                  <  2:2898697/100‑1 (MQ=255)
                    aTTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTttatt                                                                                 <  1:936160/100‑1 (MQ=255)
                    aTTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTttatt                                                                                 >  1:606416/1‑100 (MQ=255)
                    aTTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTttatt                                                                                 >  2:2266637/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)
                                aaaaTTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCg                                                                     <  1:1284041/100‑1 (MQ=255)
                                 aaaTTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGc                                                                    <  1:166219/100‑1 (MQ=255)
                                    ttACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTAc                                                                 <  2:1992067/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                                                          <  1:829870/100‑1 (MQ=255)
                                           cacaAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTcc                                                          <  2:511170/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)
                                                      cATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGAc                                                <  1:220536/99‑1 (MQ=255)
                                                      cATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACg                                               >  2:273407/1‑100 (MQ=255)
                                                         tctcTTTGCCTTCAGCAAATGTAACAGGACAAGAAATTTCTGCTCTGGAACACACCTTTTATTACGCTGCTTTTTCTATTGCGTCCACTTCCGGAc                                                <  1:484152/96‑1 (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:1894926/1‑100 (MQ=255)
                                                            ctTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAg                                         >  1:565071/1‑100 (MQ=255)
                                                             tttGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAgg                                        <  1:155567/100‑1 (MQ=255)
                                                                gcgcTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAg                                     <  2:1593620/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)
                                                                       tcaAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTATTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAg                               >  1:547728/1‑99 (MQ=255)
                                                                       tcaAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAgg                              <  1:2245267/100‑1 (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:659947/1‑100 (MQ=37)
                                                                            tGTTACAGGGCAGAAATTTCCTGTCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCaaa                         >  2:2314805/1‑100 (MQ=37)
                                                                             gTTACAGGGCAGAAATTTCCTGTCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAAc                        >  1:2165507/1‑100 (MQ=37)
                                                                             gTTACAGGGCAGAAATTTCCTGTCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAAc                        >  1:1929168/1‑100 (MQ=37)
                                                                                gctcatgtggaAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCg                     <  2:2540639/92‑1 (MQ=255)
                                                                                acagggcagaAATTTCCTGTCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCg                     >  2:180992/1‑100 (MQ=38)
                                                                                acaggacaggAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCg                     >  1:878004/1‑100 (MQ=255)
                                                                                acaggacaggAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCg                     <  1:2864650/100‑1 (MQ=255)
                                                                                  aggacaggaAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGcc                   <  2:944717/100‑1 (MQ=255)
                                                                                  aggacaggaAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGcc                   <  1:205820/100‑1 (MQ=255)
                                                                                   ggacaggaAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGcca                  >  2:2654059/1‑100 (MQ=255)
                                                                                    gacaggaAAGTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACGTCCGCACGTTTCAGGAAGGCATAGGCCAAACCCGccac                 <  1:1793620/100‑1 (MQ=255)
                                                                                      gcagaAATTTCCTGTCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCaccacc                <  1:1927355/97‑1 (MQ=255)
                                                                                      caggaAATTTCTGCCCTGTAACACACCTTTTATTACGGTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGccacca               >  1:1548131/1‑100 (MQ=255)
                                                                                       aggaAATTTCTTCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGccaccan              <  1:400589/100‑2 (MQ=255)
                                                                                       aggaAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGcc                   <  1:2266627/95‑1 (MQ=255)
                                                                                        ggaAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAATGCATAGGCCAAACCCGCCACCAGc             >  1:490611/1‑100 (MQ=255)
                                                                                        ggaAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGc             >  1:1649430/1‑100 (MQ=255)
                                                                                        ggaAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGc             >  1:490782/1‑100 (MQ=255)
                                                                                        ggaAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGc             >  2:1163864/1‑100 (MQ=255)
                                                                                        ggaAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGc             >  2:2187034/1‑100 (MQ=255)
                                                                                         gaAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCg            <  2:275744/100‑1 (MQ=255)
                                                                                          aaaTTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGt           <  1:805958/100‑1 (MQ=255)
                                                                                             tttCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTAcc        <  2:565865/100‑1 (MQ=255)
                                                                                               tCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCgg      <  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)
                                                                                                    |                                                                                                   
TTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAAT  >  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: