Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,177,332 C→T 16.6% intergenic (‑128/+180) 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,3320CT16.6% 251.3 / 33.7 107intergenic (‑128/+180)gatY/fbaBD‑tagatose 1,6‑bisphosphate aldolase 2, catalytic subunit/fructose‑bisphosphate aldolase class I
Reads supporting (aligned to +/- strand):  ref base C (42/47);  new base T (7/11);  total (49/58)
Fisher's exact test for biased strand distribution p-value = 6.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 109 total aligned reads displayed.

TCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGC  >  NC_000913/2177233‑2177431
                                                                                                   |                                                                                                   
tCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGTcc                                                                                                     >  1:98814/1‑100 (MQ=255)
    aaTATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTccc                                                                                                 <  1:681763/100‑1 (MQ=255)
        taATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATa                                                                                             <  2:97401/100‑1 (MQ=255)
        taATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATa                                                                                             >  1:2132379/1‑100 (MQ=255)
         aaTGAAATTATTTGTTTTTTTTTTTTTTGATAACGATAACAAAAACGACAATATGATATTTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAg                                                                                            >  1:616786/1‑100 (MQ=255)
          aTGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAgg                                                                                           <  1:408744/100‑1 (MQ=255)
          aTGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAgg                                                                                           >  2:1362777/1‑100 (MQ=255)
           tGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAggg                                                                                          >  2:1173952/1‑100 (MQ=255)
           tGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAggg                                                                                          >  1:154785/1‑100 (MQ=255)
            gAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGa                                                                                         <  2:2026348/100‑1 (MQ=255)
            gAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGA‑‑GAAAGGCAAAAAAACGGCCTCCCGATAGGGAAg                                                                                       <  1:2014921/100‑1 (MQ=255)
                  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                                                                                  <  1:428334/100‑2 (MQ=255)
                   tttGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCTTCCCTATCGGGAGGCCGtt                                                                                  <  2:198691/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)
                        ttttAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAGCGGCCTCCCGATAGGGAAGCCGTAGCaaa                                                                             >  1:927475/1‑100 (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)
                          ttAAATATCGAAATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAAAAAAAAAAAAAAAAACGGCTTCCCTATCGGGAGGCCGtt                                                                                  <  1:583954/93‑2 (MQ=19)
                            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                                                                   >  1:906325/1‑100 (MQ=255)
                                  cGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCgtgt                                                                   >  2:1185266/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)
                                       tAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTa                                                              >  2:2256311/1‑100 (MQ=255)
                                       tAACGATAACAAAAACGACAATATGAAAATTTTTCGAGGAGTGAAAGGCAAAAATACTGCTTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTa                                                              <  1:1932355/100‑1 (MQ=255)
                                         aCGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATg                                                            >  2:309005/1‑100 (MQ=255)
                                         aCGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATg                                                            <  2:1088780/100‑1 (MQ=255)
                                         aCGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATg                                                            >  2:79933/1‑100 (MQ=255)
                                         aCGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATg                                                            <  1:723364/100‑1 (MQ=255)
                                              cacaAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGa                                                       >  1:1581508/1‑100 (MQ=255)
                                               acaAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGAt                                                      >  2:112094/1‑100 (MQ=255)
                                               acaAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTACGCGTGTTTTTATGCCGGAt                                                      >  2:1921491/1‑100 (MQ=255)
                                                caAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGCTAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATg                                                     >  2:480607/1‑100 (MQ=255)
                                                caAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATg                                                     >  1:220947/1‑100 (MQ=255)
                                                caAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATg                                                     >  1:1149086/1‑100 (MQ=255)
                                                caAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATg                                                     >  2:862264/1‑100 (MQ=255)
                                                caAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATg                                                     >  1:1261663/1‑100 (MQ=255)
                                                   aaaCGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCgg                                                  <  2:1427624/100‑1 (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:2443369/100‑1 (MQ=255)
                                                          aaTATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAAc                                           <  1:1023763/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)
                                                              tGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGAGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATAc                                       >  2:1844254/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)
                                                                   attattCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGc                                  <  2:154786/100‑1 (MQ=255)
                                                                    ttattCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCa                                 <  2:1522744/100‑1 (MQ=255)
                                                                      attCGATGATTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGt                               <  1:1370572/100‑1 (MQ=255)
                                                                       ttCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTg                              >  1:450493/1‑100 (MQ=255)
                                                                         cGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCt                            >  1:2049738/1‑100 (MQ=255)
                                                                         cGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCt                            >  2:2024688/1‑100 (MQ=255)
                                                                         cGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCt                            >  2:1887767/1‑100 (MQ=255)
                                                                         cGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCt                            >  1:1912761/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)
                                                                                        aaaaaaaCGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATaa             >  1:436105/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:1252274/4‑100 (MQ=38)
                                                                                           gctaCGGCTTCCCTATCGGGAGGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGAc          >  1:1166173/4‑100 (MQ=38)
                                                                                           aaaaCGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGAc          <  2:2124816/100‑1 (MQ=255)
                                                                                           aaaaCGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGAc          <  2:2124873/100‑1 (MQ=255)
                                                                                           aaaaCGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTCTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGAc          <  1:1594704/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)
                                                                                                   |                                                                                                   
TCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGC  >  NC_000913/2177233‑2177431

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: