Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,248,561 G→A 17.4% A185V (GCC→GTC)  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,5610GA17.4% 324.1 / 51.4 140A185V (GCC→GTC) dhaMputative dihydroxyacetone‑specific PTS enzymes: HPr, EI components
Reads supporting (aligned to +/- strand):  ref base G (48/66);  new base A (10/14);  total (60/80)
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.96e-01

Only 100 of 144 total aligned reads displayed.

ATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCA  >  NC_000913/1248462‑1248659
                                                                                                   |                                                                                                  
 tCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGc                                                                                                   <  2:3703433/100‑1 (MQ=255)
 tCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGc                                                                                                   <  1:909562/100‑1 (MQ=255)
  cAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCa                                                                                                  >  2:3015927/1‑100 (MQ=255)
  cAGACGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCa                                                                                                  <  1:737894/100‑1 (MQ=255)
    ggCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCAtt                                                                                                >  2:361013/1‑100 (MQ=255)
    ggCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCAtt                                                                                                >  2:897790/1‑100 (MQ=255)
    ggCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCAtt                                                                                                >  2:2590375/1‑100 (MQ=255)
     gcgcAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGACAttt                                                                                               >  1:2986030/1‑99 (MQ=255)
      cgcAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTaa                                                                                              <  2:159436/100‑1 (MQ=255)
         tgCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATg                                                                                           <  1:2910559/99‑1 (MQ=255)
         aGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATg                                                                                           <  2:2178149/100‑1 (MQ=255)
          gCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATga                                                                                          >  2:2910538/1‑99 (MQ=255)
           cGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTc                                                                                         <  1:808814/100‑1 (MQ=255)
           cGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTc                                                                                         <  2:3624642/100‑1 (MQ=255)
           cGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTc                                                                                         <  1:1284475/100‑1 (MQ=255)
           cGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGACATTTAATgcc                                                                                         <  2:2986008/100‑3 (MQ=255)
               tCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGACATTTAATGCCGATa                                                                                     <  2:1209799/100‑1 (MQ=255)
               tCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATAGCGGCATTAAATGTCGATa                                                                                     <  1:361015/100‑1 (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)
                   tataGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGt                                                                                 >  1:665259/1‑100 (MQ=255)
                     taGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAGGTtat                                                                               >  2:277308/1‑100 (MQ=255)
                     taGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGtat                                                                               >  1:577371/1‑100 (MQ=255)
                     taGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGtat                                                                               >  2:127112/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)
                       gCGAGCTGGTAGTAACGCAATCTTGTTATTACTCTCTGGTGTTACGCATTTGCCGTTTTTTTCCAGCAACATATCGACATTTAATGCCGATatgttgctg                                                                             <  2:3382176/100‑9 (MQ=25)
                       gCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATaa                                                                             >  2:3617843/1‑100 (MQ=255)
                       gCGAACTTGCAGTAACGCAATCTGGTTAATACTCTCTGTTGTGACGCATTTTCCGTTTTTTTCCAGCATCATATCGACATTTAATGCCGATatgttgctg                                                                             <  1:3579818/100‑9 (MQ=25)
                             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)
                                    aaCGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCTGCAATATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAgg                                                                >  2:3089417/1‑100 (MQ=255)
                                    aaCGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAgg                                                                >  1:271603/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:329723/100‑1 (MQ=255)
                                         aaTCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGa                                                           <  1:1160970/100‑1 (MQ=255)
                                         aaTCTGGTTAATACTCTCTGGTGGGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCg                                                             <  2:1314030/98‑1 (MQ=255)
                                          aTCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGAc                                                          <  2:1872499/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                                                      <  2:3786404/100‑1 (MQ=255)
                                              ggTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTAc                                                      >  1:2661292/1‑100 (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)
                                                  aaTACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGc                                                  <  2:2583124/100‑1 (MQ=255)
                                                     aCTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAgg                                               <  2:693584/100‑1 (MQ=255)
                                                     aCTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAgg                                               <  1:522863/100‑1 (MQ=255)
                                                        ctctGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCg                                            <  1:1546841/100‑1 (MQ=255)
                                                         tctGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGt                                           >  2:1650566/1‑100 (MQ=255)
                                                          ctGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGtt                                          <  1:1470294/100‑1 (MQ=255)
                                                           tGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTa                                         <  1:2506195/100‑1 (MQ=255)
                                                            ggTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTAc                                        <  1:829579/100‑1 (MQ=255)
                                                            ggTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTAc                                        >  1:575007/1‑100 (MQ=255)
                                                                    gCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATACAGGCCGTTACGGtttttt                                <  2:50815/100‑1 (MQ=255)
                                                                     cATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTa                               <  1:990652/100‑1 (MQ=255)
                                                                     cATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTa                               <  1:3756752/100‑1 (MQ=255)
                                                                     cATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTa                               <  2:577365/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)
                                                                         tGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGac                           <  1:413073/100‑1 (MQ=255)
                                                                           ccGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGacca                         >  2:3127358/1‑100 (MQ=255)
                                                                             gTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGaccacc                       <  2:808404/100‑1 (MQ=255)
                                                                             gTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGaccacc                       >  1:3006049/1‑100 (MQ=255)
                                                                                tttttCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGcc                    <  1:779617/100‑1 (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)
                                                                                          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         >  2:368341/4‑100 (MQ=255)
                                                                                           ccttatCGACATTTAATGCCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGc         >  1:75383/4‑100 (MQ=255)
                                                                                           ccttatCGACATTTAATGCCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGc         >  2:2864996/4‑100 (MQ=255)
                                                                                             atatCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGctt       <  2:1305452/100‑1 (MQ=255)
                                                                                              tatCGGCATTAAATGTCGATAAGGTATAAACTAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGcttc      >  2:456685/1‑100 (MQ=255)
                                                                                               atCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGcttct     <  2:1447745/100‑1 (MQ=255)
                                                                                                  ggCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCa  >  1:3172229/1‑100 (MQ=255)
                                                                                                  ggCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCa  <  2:3499000/100‑1 (MQ=255)
                                                                                                  ggCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCa  >  2:3084214/1‑100 (MQ=255)
                                                                                                   |                                                                                                  
ATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCA  >  NC_000913/1248462‑1248659

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: