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.

TCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAG  >  NC_000913/2177233‑2177430
                                                                                                    |                                                                                                  
tCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAAAGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGcc                                                                                                    >  1:2262107/1‑100 (MQ=255)
cttccgatctaTGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGcc                                                                                                    <  2:280926/90‑1 (MQ=255)
 cAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGcc                                                                                                    >  1:2460594/1‑99 (MQ=255)
 cAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCt                                                                                                   >  1:2676029/1‑100 (MQ=255)
 cAAAATATAAGGAAATTATTTGTTTTTTAATATCGAGAAAAAGAAAAAAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGcc                                                                                                    <  2:3338992/99‑1 (MQ=255)
     atatAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAAATATTCGAGGAGTGAAAGGCAAAAAAACGGCCTcc                                                                                                 >  2:2769749/1‑98 (MQ=255)
        taATGAAATTATTTGTTTTTAAATATCGAGATAACGAGCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATa                                                                                            <  1:971899/100‑1 (MQ=255)
          aTGAAATTATTTGTTTTTAAATATCTAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGATTGTAAGTCAAAAAAACGGCCTCCCGATAgg                                                                                          <  2:397061/100‑1 (MQ=255)
          aTGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGccagatcggaag                                                                                          >  1:280931/1‑90 (MQ=255)
          aTGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAgg                                                                                          >  2:2889326/1‑100 (MQ=255)
          aTGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAgg                                                                                          >  1:2926923/1‑100 (MQ=255)
          aTGAAAGTCTTTTTTTTTAAATATTAGGATCATGATCACAATAACTACCATATGAAAATTATTCGTGTACTTATAGGCAAATAAACGGCCTCCCGATAg                                                                                           <  2:2275274/99‑1 (MQ=255)
             aaaTTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTAGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGaa                                                                                       >  1:402311/1‑100 (MQ=255)
              aaTTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAg                                                                                      <  1:1786925/100‑1 (MQ=255)
                ttattTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGcc                                                                                    <  1:1492597/100‑1 (MQ=255)
                 tattTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCTATCGGGAGGCCg                                                                                   <  1:2706461/100‑1 (MQ=255)
                  attTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGATGAGTGAAAGGCAAAAAAACGGCTTCCCTATCGGGAGGCCGt                                                                                  <  2:1225541/100‑1 (MQ=255)
                  attTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGt                                                                                  >  1:2848122/1‑100 (MQ=255)
                  attGGTGTTTGAATATCGAGATGACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGt                                                                                  <  2:3057250/100‑1 (MQ=255)
                   tttGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTa                                                                                 <  1:2856093/100‑1 (MQ=255)
                    ttGTTTTTAAATATCGAGATAACGGTCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAATCGGCCTCCCGATAGGGAAGCCGTa                                                                                 <  2:36974/99‑1 (MQ=255)
                    ttGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCTTCCCTATCGGGAGGCCGtt                                                                                 <  1:1903184/99‑2 (MQ=38)
                     tGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGc                                                                               <  1:1648155/100‑1 (MQ=255)
                      gTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCa                                                                              >  2:243998/1‑100 (MQ=255)
                       tttttAAATATGGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGCTAGGGAAGCCGTAGCaa                                                                             >  2:3160377/1‑100 (MQ=255)
                       tttttAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCTTCCCTATCGGGAGGCCGttttt                                                                              <  1:1021661/99‑5 (MQ=38)
                        ttttAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCTTCCCTATCGGGAGGCCGttttttt                                                                            <  1:2781132/100‑7 (MQ=37)
                        tagtAAATATAGAGATAACGATCATAAAAACGACAATATATAAATTATTCGAGAAGTGAAAGGTAAAAAAACTGTTTCCCGATAGGGAAGCCGTa                                                                                 <  2:1087160/92‑1 (MQ=255)
                         tttAAATATCGAGATAACGATCACAATAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCTTCCCTATCGGGAGGCCGtttttttg                                                                           <  1:1147001/100‑8 (MQ=37)
                         tttAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAAACGGCCTCCCGAAAGGAAACCCGAAGaaa                                                                             >  1:3681390/1‑99 (MQ=255)
                         tttAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCTTCCCTATCGGGAGGCCGtttttttg                                                                           <  2:2186636/100‑8 (MQ=37)
                         tttAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCTTCCCTATCGGGAGGCCGtttttttg                                                                           <  1:1477938/100‑8 (MQ=37)
                           tAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTg                                                                         <  2:3161493/100‑1 (MQ=255)
                           tAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTg                                                                         >  2:1731880/1‑100 (MQ=255)
                           tAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAAGGAAGCCGTAGCAAAGTg                                                                         >  1:1889190/1‑100 (MQ=255)
                            aaaTATCGAGATAACGATCACAAAAATGACAATATGAAAATTATTCGAGGAGAAAAAAGAAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAg                                                                           <  2:1061373/97‑1 (MQ=255)
                               tatCGAGATAACGATCACAAAAACGACATTATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGcgcg                                                                     <  1:2340463/100‑1 (MQ=255)
                               tatCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGcgcg                                                                     <  1:550710/100‑1 (MQ=255)
                                 tCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCgtg                                                                   <  2:2145141/100‑1 (MQ=255)
                                 tCGAGATAACAATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCgtg                                                                   <  2:3756175/100‑1 (MQ=255)
                                   gagaTAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGtt                                                                 >  1:1622358/1‑100 (MQ=255)
                                   gagaTAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGtt                                                                 >  1:808688/1‑100 (MQ=255)
                                      aTAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGttttt                                                              <  1:987833/100‑1 (MQ=255)
                                       tAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTa                                                             >  2:1116486/1‑100 (MQ=255)
                                         aCGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAATCGGCCTTTCGATAGGGAAGCCGTAGCAAAGTGCGCGTGtttt                                                               <  2:2689504/96‑1 (MQ=255)
                                          cGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGc                                                          >  2:3353323/1‑100 (MQ=255)
                                           gctCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGcc                                                         >  2:3427741/3‑100 (MQ=255)
                                               acaAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGAt                                                     >  2:3190246/1‑100 (MQ=255)
                                                  aaaaCGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAAAGGCCTCACGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGAt                                                     >  2:2298219/1‑97 (MQ=255)
                                                   aaaCGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCgg                                                 <  1:167877/100‑1 (MQ=255)
                                                   aaaCGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCgg                                                 <  1:244003/100‑1 (MQ=255)
                                                       gACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTa                                             <  2:2523261/100‑1 (MQ=255)
                                                        aCAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTaa                                            <  1:2266789/100‑1 (MQ=255)
                                                         cAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTaca                                           <  1:2650868/100‑3 (MQ=255)
                                                         cAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTaaa                                           <  1:3531408/100‑1 (MQ=255)
                                                         cAATAGGAAAATTATTTGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGAGAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTaaa                                           <  2:3701431/100‑1 (MQ=255)
                                                           atatGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACa                                         <  1:2289266/100‑1 (MQ=255)
                                                             atGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATa                                       <  1:1116497/100‑1 (MQ=255)
                                                             atGAAAATTAGTCGAGGAGTGAAAGTCATAATAACGGCCCCCCGATTGTGAAGCAGTATCGAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTa                                             <  2:1511764/94‑1 (MQ=25)
                                                              tGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATAc                                      >  1:903298/1‑100 (MQ=255)
                                                                   attattCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGc                                 <  1:2963608/100‑1 (MQ=255)
                                                                      attCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGt                              >  1:3331296/1‑100 (MQ=255)
                                                                       ttCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTg                             <  2:1143051/100‑1 (MQ=255)
                                                                       ttCGAGGAGTGAAAGGCAAAAAAACGGCCTACCGATATGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTg                             <  1:671986/100‑1 (MQ=255)
                                                                        tCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGc                            <  1:77702/100‑1 (MQ=255)
                                                                         cGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCt                           >  1:73269/1‑100 (MQ=255)
                                                                          gaggagTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTa                          <  1:3311682/100‑1 (MQ=255)
                                                                          gaggagTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTa                          <  1:3465474/100‑1 (MQ=255)
                                                                          gaggagTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTa                          <  1:3699052/100‑1 (MQ=255)
                                                                          caggaggGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTTAACATACGTGGCAGTGCTa                          <  2:2441651/99‑1 (MQ=255)
                                                                            ggagTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTAt                         <  1:1411927/99‑1 (MQ=255)
                                                                            ggagTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATg                        >  2:390882/1‑100 (MQ=255)
                                                                              agTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTa                      <  2:618073/100‑1 (MQ=255)
                                                                              agTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTa                      <  2:2845090/100‑1 (MQ=255)
                                                                              agTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTa                      <  2:2169626/100‑1 (MQ=255)
                                                                               gTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTa                      >  2:2948339/1‑99 (MQ=255)
                                                                                tGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAgg                    <  1:508451/100‑1 (MQ=255)
                                                                                tGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAgg                    <  1:778387/100‑1 (MQ=255)
                                                                                 gAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGc                   >  2:1702518/1‑100 (MQ=255)
                                                                                    aGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATg                >  2:532530/1‑100 (MQ=255)
                                                                                     agaaaaaaaaaCGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATg                <  2:584627/96‑1 (MQ=255)
                                                                                       cAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTTCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATa             <  1:890753/100‑1 (MQ=255)
                                                                                       cAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATa             >  1:3819311/1‑100 (MQ=255)
                                                                                       cAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATa             >  2:3836878/1‑100 (MQ=255)
                                                                                       cAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATa             >  1:1010373/1‑100 (MQ=255)
                                                                                       cAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATa             >  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)
                                                                                                    |                                                                                                  
TCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAG  >  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: