Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,177,332 C→T 15.0% 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,3320CT15.0% 345.6 / 43.3 140intergenic (‑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 (53/65);  new base T (12/9);  total (66/74)
Fisher's exact test for biased strand distribution p-value = 3.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.92e-01

Only 100 of 146 total aligned reads displayed.

TCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAG  >  NC_000913/2177233‑2177430
                                                                                                   |                                                                                                  
tCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGcc                                                                                                    >  1:1797931/1‑100 (MQ=255)
  aaaaTATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTc                                                                                                  >  1:1824996/1‑100 (MQ=255)
  aaaaTATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGAAAAAAAACGGcc                                                                                                    <  2:2147302/98‑1 (MQ=255)
        taATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATa                                                                                            >  2:2275580/1‑100 (MQ=255)
        taATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATa                                                                                            <  2:1456205/100‑1 (MQ=255)
        taATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATa                                                                                            <  2:2766423/100‑1 (MQ=255)
        taATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGAAAAAAAAAAAAAAAAAACGGCCTCCCGATa                                                                                            <  2:2001531/100‑1 (MQ=255)
         aaTGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAg                                                                                           <  2:1122156/100‑1 (MQ=255)
          aTGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAgg                                                                                          >  2:2156190/1‑100 (MQ=255)
           tGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTGTTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAggg                                                                                         >  1:2324037/1‑100 (MQ=255)
           tGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAggg                                                                                         >  2:1337368/1‑100 (MQ=255)
            gAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGa                                                                                        <  2:1171176/100‑1 (MQ=255)
            gAAATTATGTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCTAGGAGTAAAAAGAAAAAAAACGGCCTCCCGATAggg                                                                                         <  2:2299927/99‑1 (MQ=255)
              aaTTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAggg                                                                                         <  1:770940/97‑1 (MQ=255)
               attattTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGc                                                                                     <  2:2996317/100‑1 (MQ=255)
                tgacgtgcttttAAATATCGAGATGACCATCACATACACGACAATATGAAAATTATGCTAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAgg                                                                                          <  2:3297601/86‑1 (MQ=255)
                 tattTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCg                                                                                   <  1:15978/100‑1 (MQ=255)
                  attGGTTTGTAAAGATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGTCTCCCGATAGGGAAGCCGt                                                                                  <  2:1622835/100‑1 (MQ=255)
                   tttGTTTTTAAATATCGAGATAACGATCACAAAAACGATAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCTTCCCTATCGGGAGGCCGtt                                                                                 <  2:1480907/100‑2 (MQ=255)
                   tttGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTa                                                                                 <  1:2619248/100‑1 (MQ=255)
                    ttGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAg                                                                                <  1:1337376/100‑1 (MQ=255)
                    ttGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAg                                                                                <  2:1496714/100‑1 (MQ=255)
                     tGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGc                                                                               >  1:1463462/1‑100 (MQ=255)
                      gttagtAAATATTGAAAAAACTAGCACAAAAACGACAATATGAAAGTTATTCGAGGAGTGAACGGCAAAAAACCGGCATCCCGCTCGGGCTGCCGtggac                                                                              >  2:670425/6‑96 (MQ=255)
                        ttttAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCTTCCCTATCGGGAGGCCGttttttt                                                                            <  2:1606513/100‑7 (MQ=37)
                        ttttAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGACCTCCCGATAGGGAAGCCGTAGCaaa                                                                            >  2:2550385/1‑100 (MQ=255)
                        gtttAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCTTCCCTATCGGGAGGCCGttttttt                                                                            <  2:1510885/99‑7 (MQ=37)
                         tttAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCTTCCCTATCGGGAGGCCGtttttttg                                                                           <  2:3586524/100‑8 (MQ=37)
                         tttAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCTTCCCTATCGGGAGGCCGtttttttg                                                                           <  1:3302656/100‑8 (MQ=37)
                         tttAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAg                                                                           <  2:1386786/100‑1 (MQ=255)
                          ttAAATATCTACATAACGAGTACAAAAACGACAATATGAAAATTATTCGATGAGTGAAAGGCAAATATACGGCTTCCCTATCGGGAGGCCGtttttttgc                                                                          <  2:2994218/100‑9 (MQ=16)
                          ttAAATATCGAGATAACGATCACCAAAACGACAATAGGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCTTCCCTATCGGGAGGCCGtttttttg                                                                           <  2:1436877/99‑8 (MQ=25)
                           tAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCTTCCCTATCGGGAGGCCGtttttttgcc                                                                         <  2:2638647/100‑10 (MQ=25)
                            aaaTATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTgc                                                                        >  1:900494/1‑100 (MQ=255)
                               tatCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGcgcg                                                                     >  1:2717550/1‑100 (MQ=255)
                               tatCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGcgcg                                                                     <  1:1160090/100‑1 (MQ=255)
                               tatCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGcgcg                                                                     >  2:1427225/1‑99 (MQ=255)
                                 tCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAAAAAAAAAAAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCgtg                                                                   <  2:497848/100‑1 (MQ=255)
                                  cGAGATAACGATCACAACAACGACACTATGAAAATTATTCGAGGAGTGAAAGGCAAAAAACCGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCgtgt                                                                  <  1:2545571/100‑1 (MQ=255)
                                       tAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTa                                                             >  1:2085214/1‑100 (MQ=255)
                                       tAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTa                                                             <  1:1759624/100‑1 (MQ=255)
                                       tAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTa                                                             >  1:2955684/1‑100 (MQ=255)
                                       tAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTa                                                             >  1:1061460/1‑100 (MQ=255)
                                        aaCGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTa                                                             <  2:2949587/99‑1 (MQ=255)
                                         tcGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATg                                                           <  2:1705800/99‑1 (MQ=255)
                                          cGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGc                                                          >  1:2929787/1‑100 (MQ=255)
                                          cGATCACAAAAACGACAAGATGAAAATTATTCGATGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATg                                                           <  2:635387/99‑1 (MQ=255)
                                             tCACAAAAACTACAATATGAAAATTATTCGAGGAGTAAAAAGCAAAAAAACTGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCg                                                        <  2:1966797/99‑1 (MQ=255)
                                             tCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCgg                                                       >  1:521841/1‑100 (MQ=255)
                                               ataaaaaCGACAAGATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGAt                                                     >  1:2474327/3‑100 (MQ=255)
                                                caAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATg                                                    >  1:1186273/1‑100 (MQ=255)
                                                caAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATg                                                    >  2:190318/1‑100 (MQ=255)
                                                    aaCGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGgt                                                <  2:1139465/100‑1 (MQ=255)
                                                     aCGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGgtg                                               <  1:1024207/100‑1 (MQ=255)
                                                      cGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGgtgt                                              >  2:3561543/1‑100 (MQ=255)
                                                       gACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTa                                             <  2:245470/100‑1 (MQ=255)
                                                            tatGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACAt                                        >  1:1476061/1‑100 (MQ=255)
                                                             atGAAGATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGGGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATa                                       <  2:3461542/100‑1 (MQ=255)
                                                              tGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATAc                                      >  2:2358076/1‑100 (MQ=255)
                                                               gAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACg                                     <  1:1095945/100‑1 (MQ=255)
                                                                 aaaTTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTg                                   <  2:718365/100‑1 (MQ=255)
                                                                 aaaTTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTg                                   >  1:2218616/1‑100 (MQ=255)
                                                                   attattCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGc                                 <  2:1202434/100‑1 (MQ=255)
                                                                   attattCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGc                                 <  1:1710389/100‑1 (MQ=255)
                                                                   attattCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGc                                 <  2:347307/100‑1 (MQ=255)
                                                                      attCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCa                                <  1:792488/98‑1 (MQ=255)
                                                                         cGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCt                           <  1:1459730/100‑1 (MQ=255)
                                                                         cGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCt                           <  2:594577/100‑1 (MQ=255)
                                                                         cGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCt                           >  1:3226097/1‑100 (MQ=255)
                                                                            ggagTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATg                        >  1:5997/1‑100 (MQ=255)
                                                                               gTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAg                     >  1:3261503/1‑100 (MQ=255)
                                                                                  aaaGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCa                  <  2:1894969/100‑1 (MQ=255)
                                                                                    aGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATg                <  1:1938283/100‑1 (MQ=255)
                                                                                      gCAAAAAAACGGCCTCCCGATAGGTAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCTGTGTAAACATACGTGGCAGTGCTATGTAGGCATGa               <  2:31857/99‑1 (MQ=255)
                                                                                      gCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCTGTGTAAACATACGTGGCAGTGCTATGTAGGCATGAt              <  1:2574551/100‑1 (MQ=255)
                                                                                      gCAAAAAAACGGCCTCCCGATAGGGAAGCACTAGCAAAGTGCGCGCGTTTTTATGCCGGATGCGGTGTAAACACACGCGCCAGTGCCTTGTAGGTAt                 >  2:2913012/1‑97 (MQ=37)
                                                                                       cAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTACACATACGTGGCAGTGCTATGTCGGCATGATa             >  1:3380888/1‑100 (MQ=255)
                                                                                       cAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATa             >  2:1108748/1‑100 (MQ=255)
                                                                                       cAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATa             >  1:1563955/1‑100 (MQ=255)
                                                                                       cAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATa             >  1:2191494/1‑100 (MQ=255)
                                                                                        aaaaaaaCGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATaa            <  1:1108757/100‑1 (MQ=255)
                                                                                        aaaaaaaCGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATaa            >  2:569879/1‑100 (MQ=255)
                                                                                        aaaaaaaCGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATaa            <  1:1303559/100‑1 (MQ=255)
                                                                                         ttgctaCGGCTTCCCTATCGGGAGGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAg           >  1:1505152/6‑100 (MQ=37)
                                                                                         aaaaaaCGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAg           >  2:687677/1‑100 (MQ=255)
                                                                                           aaaaCGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGAc         <  1:3083784/100‑1 (MQ=255)
                                                                                           aaaaCGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGAc         >  2:855120/1‑100 (MQ=255)
                                                                                            ctaCGGCTTCCCTATCGGGAGGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCCGTGCTATGTAGGCATGATAAGACg        >  2:2480974/3‑100 (MQ=37)
                                                                                            ctaCGGCTTCCCTATCGGGAGGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGAc         >  2:3608049/3‑99 (MQ=38)
                                                                                            ctaCGGCTTCCCTATCGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACg        >  2:870865/3‑100 (MQ=38)
                                                                                            ctaCGGCTTCCCTATCGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACg        >  2:2700065/3‑100 (MQ=38)
                                                                                             taCGGCTTCCCTATCGGGAGGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGc       >  1:885364/2‑100 (MQ=38)
                                                                                             aaCGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGc       >  1:3266329/1‑100 (MQ=255)
                                                                                              aCGGCTTCCCTATCGGGAGGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCt      >  1:1966138/1‑100 (MQ=38)
                                                                                                ggCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTc    >  1:1017047/1‑100 (MQ=255)
                                                                                                ggCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTc    <  1:3169564/100‑1 (MQ=255)
                                                                                                 gCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCa   <  1:2021780/100‑1 (MQ=255)
                                                                                                  ccTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAg  >  1:1623007/1‑100 (MQ=255)
                                                                                                  ccTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAg  >  1:2805181/1‑100 (MQ=255)
                                                                                                   cTCCCGATAGGGAAGCCGTAGCAAGGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAg  <  2:734621/99‑1 (MQ=255)
                                                                                                   |                                                                                                  
TCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAG  >  NC_000913/2177233‑2177430

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: