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.G97.6% 300.9 / ‑0.2 82intergenic (‑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 . (2/0);  new base G (36/44);  total (38/44)
Fisher's exact test for biased strand distribution p-value = 2.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑G‑GC‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGA  >  NC_000913/3560357‑3560547
                                                                                                     |                                                                                            
tGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑G‑Gcg                                                                                              <  2:725117/100‑2 (MQ=255)
  cGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑G‑Gcgca                                                                                            >  1:1688579/1‑97 (MQ=255)
  cGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑G‑Gcgca                                                                                            <  1:1594861/100‑4 (MQ=255)
   gggCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑G‑Gcgcag                                                                                           <  2:309949/100‑5 (MQ=255)
     gCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑G‑GCGCAGGa                                                                                         >  2:1008094/1‑100 (MQ=255)
       gATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑G‑GCGCAGGATc                                                                                       <  2:1165532/100‑1 (MQ=255)
        aTGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑G‑GCGCAGGATCa                                                                                      <  2:80113/100‑1 (MQ=255)
        aTGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑G‑GCGCAGGATCa                                                                                      <  2:673066/100‑1 (MQ=255)
        aTGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑G‑GCGCAGGATCa                                                                                      <  1:1175714/100‑1 (MQ=255)
        aTGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑G‑GCGCAGGATCa                                                                                      <  1:2028777/100‑1 (MQ=255)
          gcgcTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑G‑GCGCAGGATCAgg                                                                                    >  2:321329/1‑100 (MQ=255)
           cgcTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑G‑GCGCAGGATCAgg                                                                                    >  2:1021503/1‑99 (MQ=255)
           cgcTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑G‑GCGCAGGATCAGGt                                                                                   <  1:1332013/100‑1 (MQ=255)
              tctTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑G‑GCGCAGGATCAGGtttt                                                                                <  1:1233448/100‑1 (MQ=255)
                tttttCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGc                                                                              >  1:627338/1‑100 (MQ=255)
                tttttCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGc                                                                              <  1:1619919/100‑1 (MQ=255)
                   ttcttcGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTcc                                                                           <  1:343333/100‑1 (MQ=255)
                     cttcGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGc                                                                         >  1:1163611/1‑100 (MQ=255)
                       tcGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCa                                                                       <  1:1418385/100‑1 (MQ=255)
                         ggTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑G‑GCGCCGGATAAGGCTTTGCCccgcccgc                                                                     >  2:1986780/1‑100 (MQ=255)
                         ggTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGc                                                                     <  1:1138303/100‑1 (MQ=255)
                             tGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCAt                                                                 >  2:320802/1‑100 (MQ=255)
                              gggTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCAtt                                                                >  1:368074/1‑100 (MQ=255)
                              gggTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCAtt                                                                >  1:1589552/1‑100 (MQ=255)
                              gggTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCAtt                                                                >  2:944574/1‑100 (MQ=255)
                               ggTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTg                                                               <  2:353734/100‑1 (MQ=255)
                               ggTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTg                                                               <  2:640379/100‑1 (MQ=255)
                               cttccgatctcGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTg                                                               <  2:1494344/90‑1 (MQ=255)
                                gTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGa                                                              >  1:286012/1‑100 (MQ=255)
                                gTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGa                                                              >  1:290375/1‑100 (MQ=255)
                                gTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGa                                                              <  1:116377/100‑1 (MQ=255)
                                    ccTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑GG‑CAGGATCAGGTTTAGCTCCGCCAGCTCATTGAGGTc                                                          >  2:1399264/1‑100 (MQ=255)
                                    ccTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTc                                                          >  1:1354658/1‑100 (MQ=255)
                                    ccTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTc                                                          >  2:1100279/1‑100 (MQ=255)
                                     cTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGATGGCAGCGCCGCACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTcg                                                         >  1:354462/1‑100 (MQ=255)
                                         cGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCgcggc                                                     <  2:348486/100‑1 (MQ=255)
                                         cGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGATcggaag                                                     >  1:1494354/1‑96 (MQ=255)
                                          gATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCgcggcg                                                    <  2:1290968/100‑1 (MQ=255)
                                          gATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCgcggc                                                     >  2:985916/1‑99 (MQ=255)
                                            tCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGaa                                                  <  1:1649943/100‑1 (MQ=255)
                                            tCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑G‑GCGCAGGAGCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGaa                                                  <  2:449876/100‑1 (MQ=255)
                                             cGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAAt                                                 <  1:477892/100‑1 (MQ=255)
                                               tGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAg                                               <  2:357719/100‑1 (MQ=255)
                                                gCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGt                                              >  1:36560/1‑100 (MQ=255)
                                                  cACGGCGTGTTAACCGAACTGGCAGGCAGCCCCGCACCGCCATGAT‑GCCC‑CAGGATCAGGTTTCGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCt                                            >  2:9415/1‑100 (MQ=255)
                                                    cGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGc                                          <  1:394762/100‑1 (MQ=255)
                                                     ggCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCg                                         <  2:125807/100‑1 (MQ=255)
                                                      gCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCgg                                        <  1:1906585/100‑1 (MQ=255)
                                                      gCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCgg                                        >  1:1065126/1‑100 (MQ=255)
                                                         tgtTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCt                                     <  1:619216/100‑1 (MQ=255)
                                                         tgtTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCt                                     <  2:183120/100‑1 (MQ=255)
                                                         tgtTAAACGAACTGGATGGTTGCGGCGCACCGCCATGAT‑G‑GCGCAGGAGCAATTATTTCCCTGCCAGCTCTTTGAGGTGGGGTCGAATTTTCTGCggg                                       >  2:690501/1‑98 (MQ=255)
                                                           ttAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGa                                   <  1:2012822/100‑1 (MQ=255)
                                                            tAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGAc                                  >  1:306321/1‑100 (MQ=255)
                                                            tAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGAc                                  >  2:846535/1‑100 (MQ=255)
                                                            tAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGAc                                  >  2:1669067/1‑100 (MQ=255)
                                                             aaCCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACg                                 >  2:793124/1‑100 (MQ=255)
                                                             aaCCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACg                                 <  1:1290371/100‑1 (MQ=255)
                                                              aCCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACgg                                >  1:714338/1‑100 (MQ=255)
                                                               ccGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGga                               >  1:16139/1‑100 (MQ=255)
                                                               ccGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGga                               >  1:413558/1‑100 (MQ=255)
                                                                   aCTGGAAGGCAGCGCCGCACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGaaa                           <  2:1207828/100‑1 (MQ=255)
                                                                    cTGGAAGGCAGCGCCGCACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAAt                          >  1:1972038/1‑100 (MQ=255)
                                                                      ggAAGGCAGCGCCGCACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGc                        <  1:1068292/100‑1 (MQ=255)
                                                                      ggAAGGCAGCGCCGCACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGc                        >  2:267905/1‑100 (MQ=255)
                                                                      ggAAGGCAGCGCCGCACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGc                        >  2:1197019/1‑100 (MQ=255)
                                                                      ggAAGGCAGCGCCGCACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGc                        <  1:1496526/100‑1 (MQ=255)
                                                                        aaGGCAGCGCCGCACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGctc                      <  1:1008097/100‑1 (MQ=255)
                                                                         aGGCAGCGCCGCACCGCCATGAT‑G‑GCGCAGGATCAGGGTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAGATGctct                     >  2:1116434/1‑100 (MQ=255)
                                                                          ggCAGCGCCGCACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTa                    <  1:671153/100‑1 (MQ=255)
                                                                          ggCAGCGCCGCACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTa                    >  2:705006/1‑100 (MQ=255)
                                                                          ggCAGCGCCGCACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTa                    <  1:321332/100‑1 (MQ=255)
                                                                            cAGCGCCGCACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTAcc                  <  1:1920095/100‑1 (MQ=255)
                                                                            cAGCGCCGCACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTAcc                  <  1:1272683/100‑1 (MQ=255)
                                                                               tgccgcACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGc               <  1:1050423/99‑1 (MQ=255)
                                                                               cgccgcACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGc               >  1:1360744/1‑100 (MQ=255)
                                                                                  cgcACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAATTGCTCTACCAGctct            >  1:2039529/1‑100 (MQ=255)
                                                                                  cgcACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGctct            >  2:1135249/1‑100 (MQ=255)
                                                                                    cACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTc          <  2:413555/100‑1 (MQ=255)
                                                                                    cACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTc          <  2:368071/100‑1 (MQ=255)
                                                                                    cACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTc          <  2:205529/100‑1 (MQ=255)
                                                                                    cACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTc          <  1:946673/100‑1 (MQ=255)
                                                                                       cGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGt       <  1:1215370/100‑1 (MQ=255)
                                                                                         ccATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTAc     >  1:1117358/1‑100 (MQ=255)
                                                                                           atgat‑g‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTg   <  1:999380/100‑1 (MQ=255)
                                                                                            tgat‑g‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGa  <  1:980861/100‑1 (MQ=255)
                                                                                                     |                                                                                            
TGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑G‑GC‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGA  >  NC_000913/3560357‑3560547

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑