Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,177,332 C→T 15.3% 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.3% 313.0 / 40.5 136intergenic (‑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 (49/66);  new base T (11/10);  total (60/76)
Fisher's exact test for biased strand distribution p-value = 4.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.73e-01

Only 100 of 144 total aligned reads displayed.

TCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAG  >  NC_000913/2177233‑2177430
                                                                                                    |                                                                                                  
tCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAAAGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGcc                                                                                                    >  1:2262107/1‑100 (MQ=255)
cttccgatctaTGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGcc                                                                                                    <  2:280926/90‑1 (MQ=255)
 cAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGcc                                                                                                    >  1:2460594/1‑99 (MQ=255)
 cAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCt                                                                                                   >  1:2676029/1‑100 (MQ=255)
 cAAAATATAAGGAAATTATTTGTTTTTTAATATCGAGAAAAAGAAAAAAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGcc                                                                                                    <  2:3338992/99‑1 (MQ=255)
     atatAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAAATATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTcc                                                                                                 >  2:2769749/1‑98 (MQ=255)
        taATGAAATTATTTGTTTTTAAATATCGAGATAACGAGCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATa                                                                                            <  1:971899/100‑1 (MQ=255)
          aTGAAATTATTTGTTTTTAAATATCTAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGATTGTAAGTC‑AAAAAAACGGCCTCCCGATAgg                                                                                          <  2:397061/100‑1 (MQ=255)
          aTGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGccagatcggaag                                                                                          >  1:280931/1‑90 (MQ=255)
          aTGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAgg                                                                                          >  2:2889326/1‑100 (MQ=255)
          aTGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAgg                                                                                          >  1:2926923/1‑100 (MQ=255)
          aTGAAAGTCTTTTTTTTTAAATATTAGGATCATGATCACAATAACTACCATATGAAAATTATTCGTGTACTTATAGGC‑AAATAAACGGCCTCCCGATAg                                                                                           <  2:2275274/99‑1 (MQ=255)
             aaaTTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTAGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGaa                                                                                       >  1:402311/1‑100 (MQ=255)
              aaTTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAg                                                                                      <  1:1786925/100‑1 (MQ=255)
                ttattTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGcc                                                                                    <  1:1492597/100‑1 (MQ=255)
                 tattTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCTATCGGGAGGCCg                                                                                   <  1:2706461/100‑1 (MQ=255)
                  attTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGATGAGTGAAAGGC‑AAAAAAACGGCTTCCCTATCGGGAGGCCGt                                                                                  <  2:1225541/100‑1 (MQ=255)
                  attTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGt                                                                                  >  1:2848122/1‑100 (MQ=255)
                  attGGTGTTTGAATATCGAGATGACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGt                                                                                  <  2:3057250/100‑1 (MQ=255)
                   tttGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTa                                                                                 <  1:2856093/100‑1 (MQ=255)
                    ttGTTTTTAAATATCGAGATAACGGTCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAATCGGCCTCCCGATAGGGAAGCCGTa                                                                                 <  2:36974/99‑1 (MQ=255)
                    ttGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCTTCCCTATCGGGAGGCCGtt                                                                                 <  1:1903184/99‑2 (MQ=38)
                     tGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGc                                                                               <  1:1648155/100‑1 (MQ=255)
                      gTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCa                                                                              >  2:243998/1‑100 (MQ=255)
                       tttttAAATATGGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGCTAGGGAAGCCGTAGCaa                                                                             >  2:3160377/1‑100 (MQ=255)
                       tttttAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCTTCCCTATCGGGAGGCCGttttt                                                                              <  1:1021661/99‑5 (MQ=38)
                        ttttAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCTTCCCTATCGGGAGGCCGttttttt                                                                            <  1:2781132/100‑7 (MQ=37)
                        tagtAAATATAGAGATAACGATCATAAAAACGACAATATATAAATTATTCGAGAAGTGAAAGGT‑AAAAAAACTGTTTCCCGATAGGGAAGCCGTa                                                                                 <  2:1087160/92‑1 (MQ=255)
                         tttAAATATCGAGATAACGATCACAATAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCTTCCCTATCGGGAGGCCGtttttttg                                                                           <  1:1147001/100‑8 (MQ=37)
                         tttAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAAACGGCCTCCCGAAAGGAAACCCGAAGaaa                                                                             >  1:3681390/1‑99 (MQ=255)
                         tttAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCTTCCCTATCGGGAGGCCGtttttttg                                                                           <  2:2186636/100‑8 (MQ=37)
                         tttAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCTTCCCTATCGGGAGGCCGtttttttg                                                                           <  1:1477938/100‑8 (MQ=37)
                           tAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTg                                                                         <  2:3161493/100‑1 (MQ=255)
                           tAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTg                                                                         >  2:1731880/1‑100 (MQ=255)
                           tAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAAGGAAGCCGTAGCAAAGTg                                                                         >  1:1889190/1‑100 (MQ=255)
                            aaaTATCGAGATAACGATCACAAAAATGACAATATGAAAATTATTCGAGGAGAAAAAAGA‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAg                                                                           <  2:1061373/97‑1 (MQ=255)
                               tatCGAGATAACGATCACAAAAACGACATTATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGcgcg                                                                     <  1:2340463/100‑1 (MQ=255)
                               tatCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGcgcg                                                                     <  1:550710/100‑1 (MQ=255)
                                 tCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCgtg                                                                   <  2:2145141/100‑1 (MQ=255)
                                 tCGAGATAACAATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCgtg                                                                   <  2:3756175/100‑1 (MQ=255)
                                   gagaTAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGtt                                                                 >  1:1622358/1‑100 (MQ=255)
                                   gagaTAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGtt                                                                 >  1:808688/1‑100 (MQ=255)
                                      aTAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGttttt                                                              <  1:987833/100‑1 (MQ=255)
                                       tAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTa                                                             >  2:1116486/1‑100 (MQ=255)
                                         aCGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAATCGGCCTTTCGATAGGGAAGCCGTAGCAAAGTGCGCGTGtttt                                                               <  2:2689504/96‑1 (MQ=255)
                                          cGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGc                                                          >  2:3353323/1‑100 (MQ=255)
                                           gctCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGcc                                                         >  2:3427741/3‑100 (MQ=255)
                                               acaAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGAt                                                     >  2:3190246/1‑100 (MQ=255)
                                                  aaaaCGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAAAGGCCTCACGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGAt                                                     >  2:2298219/1‑97 (MQ=255)
                                                   aaaCGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCgg                                                 <  1:167877/100‑1 (MQ=255)
                                                   aaaCGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCgg                                                 <  1:244003/100‑1 (MQ=255)
                                                       gACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTa                                             <  2:2523261/100‑1 (MQ=255)
                                                        aCAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTaa                                            <  1:2266789/100‑1 (MQ=255)
                                                         cAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTaca                                           <  1:2650868/100‑3 (MQ=255)
                                                         cAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTaaa                                           <  1:3531408/100‑1 (MQ=255)
                                                         cAATAGGAAAATTATTTGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGAGAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTaaa                                           <  2:3701431/100‑1 (MQ=255)
                                                           atatGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACa                                         <  1:2289266/100‑1 (MQ=255)
                                                             atGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATa                                       <  1:1116497/100‑1 (MQ=255)
                                                             atGAAAATTAGTCGAGGAGTGAAAGTC‑ATAATAACGGCCCCCCGATTGTGAAGCAGTATCGAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTa                                             <  2:1511764/94‑1 (MQ=25)
                                                              tGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATAc                                      >  1:903298/1‑100 (MQ=255)
                                                                   attattCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGc                                 <  1:2963608/100‑1 (MQ=255)
                                                                      attCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGt                              >  1:3331296/1‑100 (MQ=255)
                                                                       ttCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTg                             <  2:1143051/100‑1 (MQ=255)
                                                                       ttCGAGGAGTGAAAGGC‑AAAAAAACGGCCTACCGATATGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTg                             <  1:671986/100‑1 (MQ=255)
                                                                        tCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGc                            <  1:77702/100‑1 (MQ=255)
                                                                         cGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCt                           >  1:73269/1‑100 (MQ=255)
                                                                          gaggagTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTa                          <  1:3311682/100‑1 (MQ=255)
                                                                          gaggagTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTa                          <  1:3465474/100‑1 (MQ=255)
                                                                          gaggagTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTa                          <  1:3699052/100‑1 (MQ=255)
                                                                          caggaggGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTTAACATACGTGGCAGTGCTa                          <  2:2441651/99‑1 (MQ=255)
                                                                            ggagTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTAt                         <  1:1411927/99‑1 (MQ=255)
                                                                            ggagTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATg                        >  2:390882/1‑100 (MQ=255)
                                                                              agTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTa                      <  2:618073/100‑1 (MQ=255)
                                                                              agTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTa                      <  2:2845090/100‑1 (MQ=255)
                                                                              agTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTa                      <  2:2169626/100‑1 (MQ=255)
                                                                               gTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTa                      >  2:2948339/1‑99 (MQ=255)
                                                                                tGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAgg                    <  1:508451/100‑1 (MQ=255)
                                                                                tGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAgg                    <  1:778387/100‑1 (MQ=255)
                                                                                 gAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGc                   >  2:1702518/1‑100 (MQ=255)
                                                                                    aGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATg                >  2:532530/1‑100 (MQ=255)
                                                                                     agaaaaaaaaaCGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATg                <  2:584627/96‑1 (MQ=255)
                                                                                       c‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTTCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATa             <  1:890753/100‑1 (MQ=255)
                                                                                       c‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATa             >  1:3819311/1‑100 (MQ=255)
                                                                                       c‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATa             >  2:3836878/1‑100 (MQ=255)
                                                                                       c‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATa             >  1:1010373/1‑100 (MQ=255)
                                                                                       c‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATa             >  2:538922/1‑100 (MQ=255)
                                                                                         aaaaaaaCGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATaa            >  2:2557970/1‑100 (MQ=255)
                                                                                          ttgctaCGGCTTCCCTATCGGGAGGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAg           >  1:2907154/6‑100 (MQ=37)
                                                                                           aaaaaCGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGa          <  2:3492443/100‑1 (MQ=255)
                                                                                            gctaCGGCTTCCCTATCGGGAGGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGAc         >  1:1152226/4‑100 (MQ=38)
                                                                                            gctaCGGCTTCCCTATCGGGAGGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACATAGCAGTGCTATGTAGGCATGATAAGAc         >  1:380591/4‑100 (MQ=37)
                                                                                            aaaaCGTCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGAGGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGAc         <  2:386225/100‑1 (MQ=255)
                                                                                             ctaCGGCTTCCCTATCGGGAGGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACg        >  1:2171210/3‑100 (MQ=38)
                                                                                             aaaCGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCTTGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGAc         <  1:696099/99‑1 (MQ=255)
                                                                                                cGGCTTCCCTATCGGGAGGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCtt     >  1:2250034/1‑100 (MQ=38)
                                                                                                cGGCTTCCCTATCGGGAGGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTATAAACATACATGGCAGTGCTATGTAGACATGATAAGACGCtt     >  1:67612/1‑100 (MQ=37)
                                                                                                cGGCTTCCCTATCGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATCAGACGCtt     >  2:1938239/1‑100 (MQ=39)
                                                                                                cGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCtt     <  1:3157370/100‑1 (MQ=255)
                                                                                                   ccTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTTTCAGTGCTATATAGGCATGATAAGAAGCTTCa   >  1:1301612/1‑99 (MQ=38)
                                                                                                   ccTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAg  >  1:2989228/1‑100 (MQ=255)
                                                                                                    |                                                                                                  
TCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAG  >  NC_000913/2177233‑2177430

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑