Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,248,571 T→C 16.9% T182A (ACA→GCA)  dhaM ← putative dihydroxyacetone‑specific PTS enzymes: HPr, EI components

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,248,5710TC16.9% 288.8 / 49.5 135T182A (ACA→GCA) dhaMputative dihydroxyacetone‑specific PTS enzymes: HPr, EI components
Reads supporting (aligned to +/- strand):  ref base T (47/62);  new base C (9/13);  total (59/76)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.60e-01

Only 100 of 136 total aligned reads displayed.

CGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGG  >  NC_000913/1248473‑1248670
                                                                                                  |                                                                                                   
cGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTc                                                                                                    <  1:1284475/100‑1 (MQ=255)
cGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTc                                                                                                    <  2:3624642/100‑1 (MQ=255)
cGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTc                                                                                                    <  1:409790/100‑1 (MQ=255)
    tCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGACATTTAATGCCGATa                                                                                                <  2:1209799/100‑1 (MQ=255)
    tCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATAGCGGCATTAAATGTCGATa                                                                                                <  1:361015/100‑1 (MQ=255)
      gTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGACATTTAATGCCGATatg                                                                                              <  1:748413/100‑3 (MQ=255)
      gTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGACATTTAATGCCGATatg                                                                                              <  2:3380660/100‑3 (MQ=255)
       ttATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAgg                                                                                             <  1:1941339/100‑1 (MQ=255)
       ttATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATATCGATAAgg                                                                                             <  2:1044893/100‑1 (MQ=255)
        tataGCGAACTTGTAGTAACGTAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCATCAACATATCGGCATTAAATGTCGATAAg                                                                                              >  2:703964/1‑98 (MQ=255)
        tataGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCCACTGCCGTTTTTTTCCAGCAACATATCGCCATTAAATGTCGAGaa                                                                                               >  2:1490863/1‑97 (MQ=255)
        tataGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGt                                                                                            >  1:665259/1‑100 (MQ=255)
          taGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGtat                                                                                          >  1:577371/1‑100 (MQ=255)
          taGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGACATTTAATGCCGATatgttgc                                                                                          <  2:608932/100‑7 (MQ=38)
           aGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGACATTTAATGCCGATatgttgct                                                                                         <  1:377104/100‑8 (MQ=38)
            gcgcACTTGCAGTAACGCAATCTGGTTATTACTCTCTGGTGTGACGCATTGGCCGTTTTTTTCCAGCAACATATCGACATTTAATGCCGATatgttgct                                                                                         <  2:3772065/99‑8 (MQ=25)
            gCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACGTATCGACATTTAATGCCGATacgttgctg                                                                                        <  1:240281/100‑9 (MQ=255)
            gCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATaa                                                                                        >  2:3617843/1‑100 (MQ=255)
            gCGAACTTGCAGTAACGCAATCTGGTTAATACTCTCTGTTGTGACGCATTTTCCGTTTTTTTCCAGCATCATATCGACATTTAATGCCGATatgttgctg                                                                                        <  1:3579818/100‑9 (MQ=25)
              gAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAAc                                                                                      >  2:2892312/1‑100 (MQ=255)
              gAACTTGTAGTAACACAGTCTGGTTAATACTCTCTCGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGACATTTAATGCCGATatgttgctgg                                                                                       <  2:2668630/99‑10 (MQ=25)
               aaCTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAAcc                                                                                     >  1:808413/1‑100 (MQ=255)
                aCTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCa                                                                                    >  2:522857/1‑100 (MQ=255)
                   tGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGcc                                                                                 >  1:1802301/1‑100 (MQ=255)
                     tagtaACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCgg                                                                               <  1:733966/100‑1 (MQ=255)
                      agtaACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCggg                                                                              <  1:2375858/100‑1 (MQ=255)
                      agtaACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCggg                                                                              <  2:1949578/100‑1 (MQ=255)
                       gtaACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGa                                                                             >  2:66514/1‑100 (MQ=255)
                        taACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAg                                                                            <  1:1350560/100‑1 (MQ=255)
                         aaCGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCTGCAATATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAgg                                                                           >  2:3089417/1‑100 (MQ=255)
                         aaCGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAgg                                                                           >  1:271603/1‑100 (MQ=255)
                          aCGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAATTGGCGATAAGGTATAAGCCAGCCGGGAGGc                                                                          >  2:3622326/1‑100 (MQ=255)
                          aCGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGc                                                                          >  1:1302238/1‑100 (MQ=255)
                           ctcAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGcc                                                                         <  2:2313763/98‑1 (MQ=255)
                           cGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGcc                                                                         >  2:2057292/1‑100 (MQ=255)
                              aaTCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGa                                                                      >  2:3056015/1‑100 (MQ=255)
                              aaTCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGa                                                                      <  1:1160970/100‑1 (MQ=255)
                              aaTCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGa                                                                      <  1:329723/100‑1 (MQ=255)
                              aaTCTGGTTAATACTCTCTGGTGGGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCg                                                                        <  2:1314030/98‑1 (MQ=255)
                               aTCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGAc                                                                     <  2:1872499/100‑1 (MQ=255)
                               aTCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGAc                                                                     <  2:665252/100‑1 (MQ=255)
                                tCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACg                                                                    <  1:826306/100‑1 (MQ=255)
                                  tGGTTAATACTATGTGGAGTGACGCATTTCCCGTTTTTTTACAGTAACGTAACGGCAACGATTGACGACACGGTATAAACCAGCCGTGAGGCGGGACGt                                                                   >  2:3427070/1‑99 (MQ=255)
                                   ggTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTAc                                                                 <  1:590398/100‑1 (MQ=255)
                                   ggTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTAc                                                                 <  2:3786404/100‑1 (MQ=255)
                                   ggTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTAc                                                                 <  1:884240/100‑1 (MQ=255)
                                    gTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACa                                                                <  1:1309793/100‑1 (MQ=255)
                                    gTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACa                                                                <  2:3251001/100‑1 (MQ=255)
                                     ttAATACTCTCTGGTGTGGCGCATTTGCTGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAATCAGCCGGGAGGCCGGACGTACa                                                                <  2:1562406/99‑1 (MQ=255)
                                      tAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCATCATATCGGCATTAAATGTCGATAAGGTATAAACAAGCCGGGAGGCCGGACGTACATg                                                              >  2:857872/1‑100 (MQ=255)
                                      tAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATg                                                              >  2:310437/1‑100 (MQ=255)
                                      tAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATg                                                              >  1:1281787/1‑100 (MQ=255)
                                      tAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATg                                                              >  2:771990/1‑100 (MQ=255)
                                      tAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATg                                                              <  1:3806029/100‑1 (MQ=255)
                                       aaTACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGc                                                             <  2:2583124/100‑1 (MQ=255)
                                          aCTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAgg                                                          <  1:98631/100‑1 (MQ=255)
                                          aCTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAgg                                                          <  1:522863/100‑1 (MQ=255)
                                          aCTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAgg                                                          <  2:693584/100‑1 (MQ=255)
                                            tctctGGTGTGACGCATTTGCCGTTTTTTTTCTTTATCTTATCGGCAAAAAATGTCGATAAGGTATAAACCAGCCTGGAGGCCGGACGTACATg                                                              >  2:3304225/1‑94 (MQ=255)
                                             ctctGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCg                                                       <  1:1546841/100‑1 (MQ=255)
                                             ctctGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCg                                                       >  1:2567286/1‑100 (MQ=255)
                                              tctGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGt                                                      >  2:1650566/1‑100 (MQ=255)
                                                 ggTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTAc                                                   <  1:829579/100‑1 (MQ=255)
                                                      gACGCATTTGCCGTTTTTTTTTTTTTTTTTATAGGCATAAAATGACGCAAAGGTATAAACCAGCCGGGAGGCCGGACGTGCATTCCGGCCGTTACGGt                                                >  2:375023/1‑98 (MQ=255)
                                                      gACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGtt                                               <  2:3215321/99‑1 (MQ=255)
                                                        cGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGttttt                                            <  1:467051/100‑1 (MQ=255)
                                                         gCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATACAGGCCGTTACGGtttttt                                           <  2:50815/100‑1 (MQ=255)
                                                          cATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTa                                          <  1:3756752/100‑1 (MQ=255)
                                                          cATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTa                                          <  1:990652/100‑1 (MQ=255)
                                                          cATTTGCCGTTTTTTTCCAGCAACATACCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTa                                          <  1:2867114/100‑1 (MQ=255)
                                                           aTTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTAt                                         <  1:2596362/100‑1 (MQ=255)
                                                           aTTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTAt                                         >  1:1834919/1‑100 (MQ=255)
                                                            tttGCCGTTTTTTTCTTTTAATATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGtttttttt                                         >  2:3819388/1‑97 (MQ=255)
                                                             ttGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTAt                                         >  2:337650/1‑98 (MQ=255)
                                                                ccGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGacca                                    >  2:3127358/1‑100 (MQ=255)
                                                                     tttttCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGcc                               <  1:779617/100‑1 (MQ=255)
                                                                      ttttCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCa                              >  2:1277708/1‑100 (MQ=255)
                                                                      ttttCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCa                              <  1:423257/100‑1 (MQ=255)
                                                                        ttCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCaga                            >  2:3444604/1‑100 (MQ=255)
                                                                          ccAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAgaga                          >  1:374394/1‑100 (MQ=255)
                                                                            tataccttatCGACATTTAATGCCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAAc                        >  1:2514081/8‑100 (MQ=38)
                                                                             gCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACg                       <  1:1872277/100‑1 (MQ=255)
                                                                               tccttatCGACATTTAATGCCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACggg                     <  2:3821371/96‑1 (MQ=255)
                                                                               aaCATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACggg                     <  1:950652/100‑1 (MQ=255)
                                                                               aaCATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGAAGTACATGCAGGACGTTACGGTTTTTTATGACCACGGCCCGAGAACggg                     >  2:896065/1‑100 (MQ=255)
                                                                                ccttatCGACATTTAATGCCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACggga                    >  1:3821399/4‑99 (MQ=255)
                                                                                ccttatCGACATTTAATGCCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACggg                     >  2:2394303/4‑99 (MQ=255)
                                                                                ccttatCGACATTTAATGCCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGc                    >  1:75383/4‑100 (MQ=255)
                                                                                ccttatCGACATTTAATGCCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGc                    >  2:2864996/4‑100 (MQ=255)
                                                                                 cttatCGACATTTAATGCCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGct                   >  1:3085257/3‑100 (MQ=255)
                                                                                  atatCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGctt                  <  2:1305452/100‑1 (MQ=255)
                                                                                   tatCGGCATTAAATGTCGATAAGGTATAAACTAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGcttc                 >  2:456685/1‑100 (MQ=255)
                                                                                    atCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGcttct                <  2:1447745/100‑1 (MQ=255)
                                                                                     tCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGACGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGcttctt               <  1:2927662/100‑1 (MQ=255)
                                                                                       ggCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCa             >  1:3172229/1‑100 (MQ=255)
                                                                                       ggCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCa             <  2:3499000/100‑1 (MQ=255)
                                                                                          aTTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCg          <  2:1783056/100‑1 (MQ=255)
                                                                                            tAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTa        >  1:2073433/1‑100 (MQ=255)
                                                                                               aTGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGc     >  2:2668711/1‑100 (MQ=255)
                                                                                                tGTCGATAAGGTATAAACCAGCGGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCa    >  2:2211955/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
CGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGG  >  NC_000913/1248473‑1248670

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: