Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,248,566 A→T 16.9% F183L (TTT→TTA 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,5660AT16.9% 340.4 / 55.1 137F183L (TTT→TTAdhaMputative dihydroxyacetone‑specific PTS enzymes: HPr, EI components
Reads supporting (aligned to +/- strand):  ref base A (48/65);  new base T (9/14);  total (58/79)
Fisher's exact test for biased strand distribution p-value = 8.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 138 total aligned reads displayed.

CGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTA  >  NC_000913/1248468‑1248664
                                                                                                  |                                                                                                  
cgcAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTaa                                                                                                   <  2:159436/100‑1 (MQ=255)
   tgCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATg                                                                                                <  1:2910559/99‑1 (MQ=255)
    gCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATga                                                                                               >  2:2910538/1‑99 (MQ=255)
     cGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTc                                                                                              <  1:409790/100‑1 (MQ=255)
     cGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTc                                                                                              <  1:808814/100‑1 (MQ=255)
     cGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTc                                                                                              <  2:3624642/100‑1 (MQ=255)
     cGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTc                                                                                              <  2:355505/100‑1 (MQ=255)
     cGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTc                                                                                              <  1:1284475/100‑1 (MQ=255)
         tCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGACATTTAATGCCGATa                                                                                          <  2:1209799/100‑1 (MQ=255)
          cGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGACATTTAATGCCGATat                                                                                         <  1:753311/100‑2 (MQ=255)
           gTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGACATTTAATGCCGATatg                                                                                        <  2:3202479/100‑3 (MQ=255)
           gTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGACATTTAATGCCGATatg                                                                                        <  1:748413/100‑3 (MQ=255)
           gTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGACATTTAATGCCGATatg                                                                                        <  2:3380660/100‑3 (MQ=255)
            ttATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAgg                                                                                       <  1:1941339/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                                                                                    >  2:127112/1‑100 (MQ=255)
               taGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGtat                                                                                    >  1:577371/1‑100 (MQ=255)
                 gcgcACTTGCAGTAACGCAATCTGGTTATTACTCTCTGGTGTGACGCATTGGCCGTTTTTTTCCAGCAACATATCGACATTTAATGCCGATatgttgct                                                                                   <  2:3772065/99‑8 (MQ=25)
                 gCGAGCTGGTAGTAACGCAATCTTGTTATTACTCTCTGGTGTTACGCATTTGCCGTTTTTTTCCAGCAACATATCGACATTTAATGCCGATatgttgctg                                                                                  <  2:3382176/100‑9 (MQ=25)
                 gCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACGTATCGACATTTAATGCCGATacgttgctg                                                                                  <  1:240281/100‑9 (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)
                       ttGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGc                                                                            >  2:1008595/1‑100 (MQ=255)
                        tGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGcc                                                                           >  1:1802301/1‑100 (MQ=255)
                           agtaACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCggg                                                                        <  2:1949578/100‑1 (MQ=255)
                           agtaACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCggg                                                                        <  1:2375858/100‑1 (MQ=255)
                            gtaACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGa                                                                       >  2:66514/1‑100 (MQ=255)
                             taACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGa                                                                       >  2:1874477/1‑99 (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                                                                <  1:1160970/100‑1 (MQ=255)
                                   aaTCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGa                                                                >  2:3056015/1‑100 (MQ=255)
                                   aaTCTGGTTAATACTCTCTGGTGGGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCg                                                                  <  2:1314030/98‑1 (MQ=255)
                                    aTCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGAc                                                               <  2:665252/100‑1 (MQ=255)
                                     tCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACg                                                              <  1:1838103/100‑1 (MQ=255)
                                     tCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACg                                                              <  1:826306/100‑1 (MQ=255)
                                       tgtttAATACTCTCGGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTa                                                            <  1:703971/97‑1 (MQ=255)
                                        ggTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTTTTTTTTTTTATCGTCATTAAATGTCGATAAGGTATTAACCAGCCGGGAGTTCGGACGt                                                             >  2:3302878/1‑98 (MQ=255)
                                        ggTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTAc                                                           <  1:884240/100‑1 (MQ=255)
                                        ggTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTAc                                                           <  2:2239090/100‑1 (MQ=255)
                                        ggTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTAc                                                           <  2:3786404/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                                                        >  1:1281787/1‑100 (MQ=255)
                                           tAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATg                                                        <  1:3806029/100‑1 (MQ=255)
                                           tAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATg                                                        >  2:771990/1‑100 (MQ=255)
                                            aaTACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGc                                                       <  2:2583124/100‑1 (MQ=255)
                                               aCTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAgg                                                    <  2:693584/100‑1 (MQ=255)
                                               aCTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAgg                                                    <  1:98631/100‑1 (MQ=255)
                                               aCTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAgg                                                    <  1:522863/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)
                                                     tGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTa                                              <  1:2506195/100‑1 (MQ=255)
                                                      ggTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTAc                                             >  1:575007/1‑100 (MQ=255)
                                                       gtgtGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTAc                                             <  1:3490712/99‑1 (MQ=255)
                                                              gCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATACAGGCCGTTACGGtttttt                                     <  2:50815/100‑1 (MQ=255)
                                                               cATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTa                                    <  2:577365/100‑1 (MQ=255)
                                                               cATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTa                                    <  1:990652/100‑1 (MQ=255)
                                                               cATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTa                                    <  1:3756752/100‑1 (MQ=255)
                                                               cATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTa                                    <  1:3617870/100‑1 (MQ=255)
                                                               cATTTGCCGTTTTTTTCCAGCAACATACCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTa                                    <  1:2867114/100‑1 (MQ=255)
                                                                aTTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTAt                                   <  1:2596362/100‑1 (MQ=255)
                                                                 tttGCCGTTTTTTTCTTTTAATATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGtttttttt                                   >  2:3819388/1‑97 (MQ=255)
                                                                   tGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGac                                <  1:413073/100‑1 (MQ=255)
                                                                     ccGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGacca                              >  2:3127358/1‑100 (MQ=255)
                                                                       gTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGaccacc                            >  1:3006049/1‑100 (MQ=255)
                                                                           ttttCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCa                        <  1:423257/100‑1 (MQ=255)
                                                                            tttCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCag                       <  1:3780660/100‑1 (MQ=255)
                                                                             ttCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCaga                      >  2:3444604/1‑100 (MQ=255)
                                                                               ccAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAgaga                    >  1:374394/1‑100 (MQ=255)
                                                                                  gCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACg                 <  1:1872277/100‑1 (MQ=255)
                                                                                  ataccttatCGACATTTAATGCCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGaa                   >  2:392955/7‑98 (MQ=38)
                                                                                    tccttatCGACATTTAATGCCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACggg               <  2:3821371/96‑1 (MQ=255)
                                                                                    aaCATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGAAGTACATGCAGGACGTTACGGTTTTTTATGACCACGGCCCGAGAACggg               >  2:896065/1‑100 (MQ=255)
                                                                                     ccttatCGACATTTAATGCCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACggg               >  2:2394303/4‑99 (MQ=255)
                                                                                     ccttatCGACATTTAATGCCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGc              >  2:368341/4‑100 (MQ=255)
                                                                                     ccttatCGACATTTAATGCCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGc              >  2:2864996/4‑100 (MQ=255)
                                                                                     ccttatCGACATTTAATGCCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGc              >  1:75383/4‑100 (MQ=255)
                                                                                      cttatCGACATTTAATGCCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGct             >  1:3085257/3‑100 (MQ=255)
                                                                                       atatCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGctt            <  2:1305452/100‑1 (MQ=255)
                                                                                            ggCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCa       <  2:3499000/100‑1 (MQ=255)
                                                                                            ggCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCa       >  2:3084214/1‑100 (MQ=255)
                                                                                            ggCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCa       >  1:3172229/1‑100 (MQ=255)
                                                                                              cATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATc     <  1:891372/100‑1 (MQ=255)
                                                                                               aTTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCg    <  2:1783056/100‑1 (MQ=255)
                                                                                               aTTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCg    >  2:1337458/1‑100 (MQ=255)
                                                                                                  |                                                                                                  
CGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTA  >  NC_000913/1248468‑1248664

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: