Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,560,455:1 +G 100% intergenic (‑2/+1) glpR ← / ← glpR pseudogene, DNA‑binding transcriptional repressor;regulator; Energy metabolism, carbon: Anaerobic respiration; repressor of the glp operon/pseudogene, DNA‑binding transcriptional repressor;regulator; Energy metabolism, carbon: Anaerobic respiration; repressor of the glp operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,560,4551.G98.0% 354.6 / ‑0.0 100intergenic (‑2/+1)glpR/glpRpseudogene, DNA‑binding transcriptional repressor;regulator; Energy metabolism, carbon: Anaerobic respiration; repressor of the glp operon/pseudogene, DNA‑binding transcriptional repressor;regulator; Energy metabolism, carbon: Anaerobic respiration; repressor of the glp operon
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base G (36/62);  total (36/64)
Fisher's exact test for biased strand distribution p-value = 5.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 106 total aligned reads displayed.

TTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑C‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGA  >  NC_000913/3560356‑3560547
                                                                                                      |                                                                                            
ttGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑c                                                                                               <  2:1404229/100‑1 (MQ=255)
 tGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑cg                                                                                              <  2:1657102/100‑2 (MQ=255)
  gCGGGCGATGCGCGCTTTTTCTTCGGTCTGGGGGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑cgc                                                                                             <  1:1337490/100‑3 (MQ=255)
   cGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑cgca                                                                                            >  1:2167520/1‑97 (MQ=255)
     ggCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑cgcagg                                                                                          <  2:2011582/100‑6 (MQ=255)
     ggCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑cgcagg                                                                                          >  1:661679/1‑95 (MQ=255)
       cGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGAt                                                                                        >  2:768172/1‑100 (MQ=255)
        gatctgcTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATc                                                                                       <  2:2586654/95‑1 (MQ=255)
        gATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATc                                                                                       <  2:494742/100‑1 (MQ=255)
         aTGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCa                                                                                      <  2:2098331/100‑1 (MQ=255)
           gcgcTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAgg                                                                                    >  2:2456445/1‑100 (MQ=255)
            cgcTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGt                                                                                   <  1:1426755/100‑1 (MQ=255)
             gcTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAgatc                                                                                  >  1:2586639/1‑97 (MQ=255)
               tctTTGTCTTCGTTCTGGGTTGCCTTGCGATCGTGCCACTGCGTGTTAACCGAACTGGAAGGCAGCTCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGtttt                                                                                <  1:2111026/100‑1 (MQ=255)
                 tttttCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGc                                                                              >  2:1627204/1‑100 (MQ=255)
                 tttttCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGc                                                                              <  2:850184/100‑1 (MQ=255)
                    ttcttcGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTcc                                                                           <  2:1886710/100‑1 (MQ=255)
                     tcttcGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCg                                                                          <  2:254802/100‑1 (MQ=255)
                     tcttcGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCAGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCg                                                                          <  2:1199800/100‑1 (MQ=255)
                     tcgtcGGTCTGGCTGGCGTTACTATCGGGCCACGGCATGGTAAGGGAGCTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGAATCAGGt                                                                                   <  1:1518374/88‑1 (MQ=255)
                       ttcGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGcc                                                                        <  1:1119534/100‑1 (MQ=255)
                       ttcGGTCTGGGTGGCCTTGCGATCGTGCAACGGCGTGTTAACAGAACTGGAAGGCAGCGCCGAACCGCCATGAT‑GG‑CGCAGGATGAGGTTTTGCCCCGcc                                                                        <  1:944580/100‑1 (MQ=255)
                       tgcgctCTGTGTGGCCGTTCGATCGTGCCACGGCGTGTTAACCGAACTTCAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCgtc                                                                        <  1:122011/95‑3 (MQ=255)
                          ggTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGc                                                                     <  2:2825164/100‑1 (MQ=255)
                               gggTTGCCTTTCGCTCGTTCCACGGCGTGTTAATCGAACTGTAAGGCAGTTCCGTACCTTCATGAT‑GGTG‑CAGGATTATGTTTTGCTCCGCCAGCTCAtt                                                                <  1:582978/100‑1 (MQ=255)
                               gggTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCAtt                                                                >  1:1788651/1‑100 (MQ=255)
                               gggTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCAtt                                                                <  1:1515336/100‑1 (MQ=255)
                                ggTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTg                                                               <  1:2108148/100‑1 (MQ=255)
                                  tGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAg                                                             <  1:1201295/100‑1 (MQ=255)
                                   ggCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGatg                                                            <  2:586200/100‑3 (MQ=255)
                                   ggCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAgg                                                            <  1:2449156/100‑1 (MQ=255)
                                    gCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGt                                                           <  1:2696734/100‑1 (MQ=255)
                                    cccTTGCGATCGTGACACGCCGTGTTACCCGAACTAGGAGGCCGCGCCGCACCGCCATGAT‑GG‑CGCAGGACCAGGTTTTGCCCCGCCAGCACATTGAgg                                                            >  1:1020171/2‑99 (MQ=255)
                                     ccTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTc                                                          >  2:685438/1‑100 (MQ=255)
                                     ccTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAACTCATTGAgg                                                            >  1:1875535/1‑98 (MQ=255)
                                      cTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTcg                                                         <  2:2798171/100‑1 (MQ=255)
                                        tGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTcgcg                                                       <  2:1399909/100‑1 (MQ=255)
                                         gCGATCGTGCCCCGGCGTGTTAACCGCACTGGAAGGCCGCGCCCGACCGCCATGCT‑GG‑CGCCGGCTCGGGTTTTGCTCCGCCCGCTCACAGATGTcgc                                                        >  1:2689140/1‑98 (MQ=255)
                                          cGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCgcggc                                                     <  2:1771707/100‑1 (MQ=255)
                                          cGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCgcggc                                                     <  2:2579974/100‑1 (MQ=255)
                                          cGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCgcggc                                                     <  2:693178/100‑1 (MQ=255)
                                           gATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCgcggcg                                                    <  2:1713276/100‑1 (MQ=255)
                                             tCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGaa                                                  <  1:527715/100‑1 (MQ=255)
                                             tCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGaa                                                  >  2:1239057/1‑100 (MQ=255)
                                              cGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAAt                                                 >  2:2401235/1‑100 (MQ=255)
                                              cGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAAt                                                 <  2:1824274/100‑1 (MQ=255)
                                              cGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAAt                                                 >  2:930693/1‑100 (MQ=255)
                                               gTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATa                                                >  2:1873674/1‑100 (MQ=255)
                                                 gCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGt                                              >  1:2481147/1‑100 (MQ=255)
                                                 gCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGt                                              <  1:1017184/100‑1 (MQ=255)
                                                  ccACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTc                                             <  2:1278099/100‑1 (MQ=255)
                                                  ccACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTc                                             <  1:857246/100‑1 (MQ=255)
                                                   cACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCt                                            <  1:1263274/100‑1 (MQ=255)
                                                    aCGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTg                                           >  1:2514258/1‑100 (MQ=255)
                                                      ggCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCg                                         <  1:987413/100‑1 (MQ=255)
                                                       gCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCgg                                        <  1:2536163/100‑1 (MQ=255)
                                                        cGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCggg                                       <  1:2164639/100‑1 (MQ=255)
                                                         gtgtTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGc                                      <  1:450051/100‑1 (MQ=255)
                                                         gtgtTAACCGAACTATGAGTGAGCGCGTCACCGCCATGTT‑GG‑CGCAGTATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGc                                      <  1:2393209/100‑1 (MQ=255)
                                                          tgtTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCt                                     >  2:1145540/1‑100 (MQ=255)
                                                          tgtTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCACCTCATTGACGTCGCGGCGAATAGTCTGCGGGc                                      >  1:760695/1‑99 (MQ=255)
                                                             tAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGAc                                  >  2:135798/1‑100 (MQ=255)
                                                             tAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGAc                                  <  1:467740/100‑1 (MQ=255)
                                                              aaCCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACg                                 <  1:2255337/100‑1 (MQ=255)
                                                              aaCCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACg                                 >  2:629317/1‑100 (MQ=255)
                                                              aaCCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACg                                 >  2:773970/1‑100 (MQ=255)
                                                                ccGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACgg                                >  1:2300545/1‑99 (MQ=255)
                                                                ccGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGga                               >  2:2511751/1‑100 (MQ=255)
                                                                ccGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGga                               >  1:2704718/1‑100 (MQ=255)
                                                                ccGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGga                               >  2:2511571/1‑100 (MQ=255)
                                                                 cGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGgag                              >  2:284059/1‑100 (MQ=255)
                                                                 cGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGgag                              <  2:1982345/100‑1 (MQ=255)
                                                                  gAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGCCGgaga                             >  2:954892/1‑100 (MQ=255)
                                                                  gAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGgaga                             <  1:880424/100‑1 (MQ=255)
                                                                  gAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGgag                              <  1:1399023/99‑1 (MQ=255)
                                                                     cTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAAt                          >  2:1708168/1‑100 (MQ=255)
                                                                      tGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATg                         <  2:1858041/100‑1 (MQ=255)
                                                                      tGGAACGCAGCGCCGCAACGCCATGAT‑GG‑CGCAGGATCAGGTTTCGCTCCGCAAGCTAATTGACGTCGCGGCGAATAGTCTCCGGGCTGAAGGAGAAATg                         >  1:2038022/1‑100 (MQ=255)
                                                                       ggAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGc                        >  2:1590444/1‑100 (MQ=255)
                                                                       ggAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGc                        <  2:1908180/100‑1 (MQ=255)
                                                                       ggAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGc                        <  1:684389/100‑1 (MQ=255)
                                                                       ggAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGc                        <  2:1784327/100‑1 (MQ=255)
                                                                        gAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGct                       <  2:2595610/100‑1 (MQ=255)
                                                                          aGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGctct                     >  2:2893073/1‑100 (MQ=255)
                                                                           ggCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTa                    <  2:1909341/100‑1 (MQ=255)
                                                                               gcgcCGCACCGCCATGAT‑GG‑CGCAGGAGCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGACTAGTCTGCGGGCTGACGGAGAAATGctct                     <  1:2623069/95‑1 (MQ=255)
                                                                                cgccgcACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGc               >  2:1149650/1‑100 (MQ=255)
                                                                                cgccgcACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGc               <  1:51538/100‑1 (MQ=255)
                                                                                cgccgcACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGc               <  2:835960/100‑1 (MQ=255)
                                                                                 gccgcACCGCCATGATGGG‑G‑GAGGATCAGGTTTTGCTCCGCCAGCTCAGTGAGGTCGGGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGc               <  1:1819263/99‑1 (MQ=255)
                                                                                  ccgcACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGatc             >  1:115687/1‑100 (MQ=255)
                                                                                    cctgcGCCATCAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCtt           >  2:2023337/5‑100 (MQ=255)
                                                                                     cACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTc          >  1:1012118/1‑100 (MQ=255)
                                                                                       ccGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGTGCTGACGGAGAAATGCTCTACCAGCTCTTCgg        <  2:468886/100‑1 (MQ=255)
                                                                                       ccGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCgg        >  2:733156/1‑100 (MQ=255)
                                                                                        cGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCgg        <  1:1650750/99‑1 (MQ=255)
                                                                                        cGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGt       >  2:941440/1‑100 (MQ=255)
                                                                                           cATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTAc     <  2:873159/99‑1 (MQ=255)
                                                                                            atgat‑gG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTg   <  2:499836/100‑1 (MQ=255)
                                                                                             tgat‑gG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGa  <  2:2623088/100‑1 (MQ=255)
                                                                                                      |                                                                                            
TTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑C‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGA  >  NC_000913/3560356‑3560547

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: ‑