Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,177,345 A→G 15.6% intergenic (‑141/+167) gatY ← / ← fbaB D‑tagatose 1,6‑bisphosphate aldolase 2, catalytic subunit/fructose‑bisphosphate aldolase class I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,177,3450AG15.6% 248.5 / 36.9 109intergenic (‑141/+167)gatY/fbaBD‑tagatose 1,6‑bisphosphate aldolase 2, catalytic subunit/fructose‑bisphosphate aldolase class I
Reads supporting (aligned to +/- strand):  ref base A (43/49);  new base G (7/10);  total (50/59)
Fisher's exact test for biased strand distribution p-value = 7.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.86e-01

Only 100 of 112 total aligned reads displayed.

GAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCA  >  NC_000913/2177245‑2177444
                                                                                                    |                                                                                                   
      attTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCTTCCCTATCGGGAGGCCGt                                                                                                <  1:57384/100‑1 (MQ=255)
      attTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCTTCCCTATCGGGAGGCCGt                                                                                                <  2:1374549/100‑1 (MQ=255)
      attTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGt                                                                                                >  2:423844/1‑100 (MQ=255)
      atatGTCTTTAAAAATCGAGAAAACAAATACAAAAACGACAATATGAAAATTATTCGAGGAGAGAAAGGAAAAAAAACGGCCTCCCGATAGGGAAg                                                                                                    <  1:489836/93‑1 (MQ=255)
      agttGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCTTCCCTATCGGGAGGCCGt                                                                                                <  1:2234218/98‑1 (MQ=255)
       tttGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCTTCCCTATCGGGAGGCCGtt                                                                                               <  2:198691/100‑2 (MQ=255)
       tttGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCTTCCCTATCGGGAGGCCGtt                                                                                               <  1:428334/100‑2 (MQ=255)
       tttGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTa                                                                                               >  2:182524/1‑100 (MQ=255)
       tttGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGAAAAAAGAAAAAAAACGGCTTCCCTATCGGGAGGCCGt                                                                                                <  1:2564527/99‑1 (MQ=255)
        ttGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAg                                                                                              <  2:2355757/100‑1 (MQ=255)
         tGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGc                                                                                             <  1:1937634/100‑1 (MQ=255)
          gTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCa                                                                                            >  2:1391988/1‑100 (MQ=255)
            ttttAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAGCGGCCTCCCGATAGGGAAGCCGTAGCaaa                                                                                          >  1:927475/1‑100 (MQ=255)
             tttAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAg                                                                                         <  1:536028/100‑1 (MQ=255)
             tttAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAg                                                                                         <  2:298047/100‑1 (MQ=255)
             tttAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAg                                                                                         <  2:1908367/100‑1 (MQ=255)
              ttAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCTTCCCTATCGGGAGGCCGtttttttgc                                                                                        <  1:1469206/100‑9 (MQ=25)
              ttAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCTTCCCTATCGGGAGGCCGtttttttgc                                                                                        <  2:2525524/100‑9 (MQ=25)
              ttAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCTATCGGGAGGCCGtttttttgc                                                                                        <  1:1565835/100‑9 (MQ=38)
                aaaGATTGAGATAAAGAGCAGAAAAAAGACAATATGAAAATTATTCGAGGAAAAAAAAAAAAAAAAACGGCCTCCCGATAGGGAAGCCGt                                                                                                <  1:1445559/90‑1 (MQ=255)
                   tatCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGcgcg                                                                                   <  1:2264541/100‑1 (MQ=255)
                    atCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCgt                                                                                  <  1:511687/100‑1 (MQ=255)
                     tCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCgtg                                                                                 >  1:2005988/1‑100 (MQ=255)
                     tCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCgtg                                                                                 >  2:1291439/1‑100 (MQ=255)
                      cGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCgtgt                                                                                >  2:1185266/1‑100 (MQ=255)
                      cGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCgtgt                                                                                >  1:906325/1‑100 (MQ=255)
                       gagaTAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGtt                                                                               >  2:1470199/1‑100 (MQ=255)
                         gaTAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGtttt                                                                             <  2:1269226/100‑1 (MQ=255)
                           tAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTa                                                                           <  2:1997580/100‑1 (MQ=255)
                           tAACGATAACAAAAACGACAATATGAAAATTTTTCGAGGAGTGAAAGGCAAAAATACTGCTTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTa                                                                           <  1:1932355/100‑1 (MQ=255)
                             aCGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATg                                                                         <  1:723364/100‑1 (MQ=255)
                             aCGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATg                                                                         >  2:79933/1‑100 (MQ=255)
                             aCGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATg                                                                         >  2:309005/1‑100 (MQ=255)
                             aCGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATg                                                                         <  2:1088780/100‑1 (MQ=255)
                                aTCACAAAAACGACAATATGAAAATTATTCGAGGAGTAAAAGGAAAAAAAACGACCTCCCGATAGGGAAGACGTAGCAAAATGACCGTGTTTTTATGCCg                                                                      >  1:2004201/1‑100 (MQ=255)
                                  cacaAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGa                                                                    >  1:1581508/1‑100 (MQ=255)
                                   acaAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGAt                                                                   >  2:112094/1‑100 (MQ=255)
                                   acaAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTACGCGTGTTTTTATGCCGGAt                                                                   >  2:1921491/1‑100 (MQ=255)
                                    caAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATg                                                                  >  2:862264/1‑100 (MQ=255)
                                    caAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATg                                                                  >  1:1149086/1‑100 (MQ=255)
                                    caAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATg                                                                  >  1:1515903/1‑100 (MQ=255)
                                    caAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATg                                                                  >  1:220947/1‑100 (MQ=255)
                                    caAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATg                                                                  >  1:1261663/1‑100 (MQ=255)
                                       aaaCGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCgg                                                               <  1:1673719/100‑1 (MQ=255)
                                          cGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGgtgt                                                            >  1:1321096/1‑100 (MQ=255)
                                            aCAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGgtgt                                                            <  2:675138/98‑1 (MQ=255)
                                            aCAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTaa                                                          <  2:1515335/100‑1 (MQ=255)
                                              aaTATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAAc                                                        <  2:1034459/100‑1 (MQ=255)
                                              aaTATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAAc                                                        <  1:1023763/100‑1 (MQ=255)
                                              aaTATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAAc                                                        <  2:2443369/100‑1 (MQ=255)
                                               atatGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACa                                                       <  1:1751607/100‑1 (MQ=255)
                                               atatGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACa                                                       <  1:2140783/100‑1 (MQ=255)
                                               atatGAAAATTATGCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACa                                                       <  1:1316697/100‑1 (MQ=255)
                                                tatGAACATTATTCTAGGAAAAAAAAAAAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTaa                                                          <  1:1138250/96‑1 (MQ=37)
                                                tatGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACAt                                                      >  2:1036554/1‑100 (MQ=255)
                                                   gAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACg                                                   >  2:584598/1‑100 (MQ=255)
                                                   gAAAATTATTCGAGGAGTGAAAGGCAAAAAAAAGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACg                                                   >  2:1267343/1‑100 (MQ=255)
                                                     aaaTTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTg                                                 <  2:803509/100‑1 (MQ=255)
                                                     aaaTTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTg                                                 >  2:2571560/1‑100 (MQ=255)
                                                       attattCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGc                                               <  2:154786/100‑1 (MQ=255)
                                                        ttattCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCa                                              <  2:1522744/100‑1 (MQ=255)
                                                           ttCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTg                                           >  1:450493/1‑100 (MQ=255)
                                                             cGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCt                                         >  2:2024688/1‑100 (MQ=255)
                                                             cGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCt                                         >  1:2049738/1‑100 (MQ=255)
                                                             cGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCt                                         >  1:1912761/1‑100 (MQ=255)
                                                             cGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCt                                         >  2:1887767/1‑100 (MQ=255)
                                                              gaggagTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTa                                        <  2:2263966/100‑1 (MQ=255)
                                                                 gagTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGt                                     <  1:564814/100‑1 (MQ=255)
                                                                    tGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAgg                                  <  1:500495/100‑1 (MQ=255)
                                                                      aaaGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCa                                <  1:2448939/100‑1 (MQ=255)
                                                                         ggCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGa                             >  1:594765/1‑100 (MQ=255)
                                                                           naaaaaaaCGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATa                           >  1:624601/2‑100 (MQ=255)
                                                                           ctttgctaCGGCTTCCCTATCGGGAGGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATa                           >  1:2000763/8‑100 (MQ=37)
                                                                           cAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATa                           >  2:1575942/1‑100 (MQ=255)
                                                                             atacgccGGCTTCCCGATAGGGAAGCCGTAGCACAGGGCGCGTGTTGTAATGCCGAATGAGGGGTGCACATATGTGGGAGTGCAATGTAGGCATGa                             <  1:2566288/90‑1 (MQ=255)
                                                                             aaaaaaCGGCCTCCCGATAGTGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGTTGTAAACATACGTGGCACTGCTATGTAGGCATGATAAg                         <  1:1583134/100‑1 (MQ=255)
                                                                              aaaaaCGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGa                        <  1:2116135/100‑1 (MQ=255)
                                                                               gctaCGGCTTCCCTATCGGGAGGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGAc                       >  1:1166173/4‑100 (MQ=38)
                                                                               gctaCGGCTTCCCTATCGGGAGGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGAc                       >  1:1252274/4‑100 (MQ=38)
                                                                               aaaaCGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGAc                       <  2:2124816/100‑1 (MQ=255)
                                                                               aaaaCGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGAc                       <  2:2124873/100‑1 (MQ=255)
                                                                               aaaaCGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGAc                       <  2:1975944/100‑1 (MQ=255)
                                                                                ctgcGGCTTCCCTATCGGGAGGCCGCAGCAGAGTGCGCGTGTGTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACg                      >  1:2163614/4‑100 (MQ=37)
                                                                                ctaCGGCTTCCCTATCGGGAGGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACg                      >  2:645399/3‑100 (MQ=38)
                                                                                aaaCGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACg                      <  1:1459805/100‑1 (MQ=255)
                                                                                 taCGGCTTCCCTATCGGGAGGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGc                     >  1:2454212/2‑100 (MQ=38)
                                                                                  aCGGCTTCCCTATCGGGAGGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCt                    >  2:1055970/1‑100 (MQ=38)
                                                                                       cTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGc               >  1:2134883/1‑100 (MQ=255)
                                                                                        tCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTCTGTAGGCATGATAAGACGCTTCAGCg              >  1:2312369/1‑100 (MQ=255)
                                                                                         cccGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGt             >  1:944769/1‑100 (MQ=255)
                                                                                          ccGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCg              <  1:510774/98‑1 (MQ=255)
                                                                                          ccGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGt             >  1:2588896/1‑99 (MQ=255)
                                                                                           cGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGTCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCg           <  2:1829667/100‑1 (MQ=255)
                                                                                           cGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCg           >  2:878825/1‑100 (MQ=255)
                                                                                            gATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGc          >  1:1453719/1‑100 (MQ=255)
                                                                                               aGGGAAGCCGTAGCGAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATc       <  2:567260/100‑1 (MQ=255)
                                                                                                  gAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGTTATGTAGGCATGATAAGACGCTGAAGCGTcttatcatg    <  1:1335640/100‑9 (MQ=38)
                                                                                                   aaGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGc   <  1:302669/100‑1 (MQ=255)
                                                                                                    aGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCa  <  1:423842/100‑1 (MQ=255)
                                                                                                    aGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCa  >  2:1759325/1‑100 (MQ=255)
                                                                                                    |                                                                                                   
GAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCA  >  NC_000913/2177245‑2177444

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: