Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,945,378 A→G 8.9% intergenic (+693/+105) PP_5673 → / ← flhB D‑alanine‑‑D‑alanine ligase/flagellin export apparatus substrate specificity protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,945,3780AG 8.9% 191.9 / 4.0 64intergenic (+693/+105)PP_5673/flhBD‑alanine‑‑D‑alanine ligase/flagellin export apparatus substrate specificity protein
Reads supporting (aligned to +/- strand):  ref base A (31/27);  new base G (4/2);  total (35/29)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGC  >  NC_002947/4945229‑4945516
                                                                                                                                                     |                                                                                                                                          
tcgcgatgcgcCGCGCGGGCGGCGCTCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGcac                                                                                                                                           >  2:1365720/8‑150 (MQ=38)
      tgcgcCGCGCGGGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACaggag                                                                                                                                      >  2:581032/2‑150 (MQ=255)
      tgcgcCGCGCGGGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCACGCGCGCCGGGCACaggag                                                                                                                                      >  1:1047043/2‑150 (MQ=255)
         gccgccCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGa                                                                                                                                   <  1:244724/150‑1 (MQ=255)
          ccgcccgcGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCNTACAAACAAGGCNGTCGCGCGCGCCGGGCACAGGAGGNAc                                                                                                                                  <  1:1424442/150‑1 (MQ=255)
          ccgcccgcGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCCCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGAc                                                                                                                                  <  1:686495/150‑1 (MQ=255)
            gcgcgggcggcgCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCACGCGCGCCGGGCACAGGAGGGACTg                                                                                                                                <  2:1047043/144‑1 (MQ=255)
            gcccgcGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTg                                                                                                                                >  1:1651332/1‑150 (MQ=255)
                gcgcgGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGccc                                                                                                                            <  2:2271984/150‑1 (MQ=255)
                   cggcgCTCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGaaa                                                                                                                        <  1:1365720/150‑1 (MQ=255)
                          tCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATg                                                                                                                  <  2:1651332/150‑1 (MQ=255)
                           cgcgAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGt                                                                                                                 <  1:581032/150‑1 (MQ=255)
                            tgctgCAAGCCCGCTCCTACGGTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTa                                                                                                                <  2:2190657/146‑1 (MQ=255)
                             cgAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAg                                                                                                               <  2:1529362/150‑1 (MQ=255)
                                 cAAGCTCGCTCCTTCATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGgagc                                                                                                           <  2:2243135/150‑1 (MQ=255)
                                 cAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCCTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGCGCACAGGAGGGACTGGCCCGAAACAAATGTAgagat                                                                                                           >  1:2190657/1‑146 (MQ=255)
                                 cAAGCTCGCGCCTACCTTTGTTTCGGGCCACTTATCCCTGCCGGTTTTACGCGCGACCGCCTTCGCGTACCCCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGgagc                                                                                                           <  1:9658/150‑1 (MQ=25)
                                    gctcgctcCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGAACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACACATGTAGgagcgag                                                                                                        >  2:613144/1‑150 (MQ=255)
                                      ccgctcCTACGTTTGTTTCGGGCCAATTATTCCTGCGGGATTTGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACaaac                                                                                                                   >  2:648728/2‑135 (MQ=2)
                                      ccgctcCTACGTTTGTTTCGGGCCAATTATTCCTGCGGGATTTGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACaaac                                                                                                                   <  1:648728/135‑2 (MQ=2)
                                      ccgctcCTACGTTTGTTTCGGGCCAATTATTCCTGCGGGATTTGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAAAGTAGgagcgggatt                                                                                                     >  1:218211/2‑145 (MQ=0)
                                            ntACATTTGTTTCGGGCCAGTTATTCCTNCGGGATTTACGCGCGACCNCCTTCGCGTACACCTCGATATTGCGCCGTGCAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTCGATcg                                                                                                >  1:159392/2‑150 (MQ=255)
                                                aTTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACagga                                                                                                                                       <  1:2201028/107‑1 (MQ=255)
                                                aTTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACagga                                                                                                                                       >  2:2201028/1‑107 (MQ=255)
                                                aTTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGAt                                                                                            >  1:1679458/1‑150 (MQ=255)
                                                     tttCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGga                                                                                                             <  2:1960409/128‑1 (MQ=255)
                                                     tttCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGga                                                                                                             >  1:1960409/1‑128 (MQ=255)
                                                        cGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGcgccgcc                                                                                    >  2:1540155/1‑149 (MQ=255)
                                                         gggCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCGCGCGATGCGCcgcgc                                                                                   >  2:1683659/1‑150 (MQ=255)
                                                         gggCCAGTTATTACTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAAAAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCACGAAACAAATGTAGGAGCGAGCGTGCTCGCGATGCGCcgcgc                                                                                   >  1:179409/1‑150 (MQ=255)
                                                            ccAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCgcggg                                                                                >  2:1566998/1‑150 (MQ=255)
                                                            ccAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCgcggg                                                                                >  2:1518523/1‑150 (MQ=255)
                                                                ttATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGggcggc                                                                            >  2:937464/1‑150 (MQ=255)
                                                                    tCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCCCGCGCGgcgcat                                                                        <  1:1540155/150‑3 (MQ=255)
                                                                      cTGCGGGGTTTACGCGCGGCCGGCTGCGCGGACCGCTCGATATTGCGCCGTACAAGCAAGGCGGTCGGGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCCCGa                                                                      <  2:159392/150‑1 (MQ=21)
                                                                      cTGCGGGATTTGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACaaac                                                                                                                   >  2:2258007/1‑103 (MQ=2)
                                                                      cTGCGGGATTTGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACaaac                                                                                                                   <  1:2258007/104‑2 (MQ=2)
                                                                       tGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTg                                                                                                                                >  2:1954239/1‑91 (MQ=255)
                                                                       tGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTg                                                                                                                                <  1:1954239/91‑1 (MQ=255)
                                                                       tGCGGGATTTACGCGCGAACGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGgagcgag                                                                                                        >  2:919614/1‑115 (MQ=255)
                                                                       tGCGGGATTTACGCGCGAACGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGgagcgag                                                                                                        <  1:919614/115‑1 (MQ=255)
                                                                       tGCGGGATTTACGCCCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGTCCGTAAACAAATGTATTAGCGAGGTTTGTCGCGATGCGCCGCGCGGGCGGCGCTCGAt                                                                     >  1:398125/1‑150 (MQ=21)
                                                                               ttACGCGCGACCGCCATCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAgg                                                             <  2:1679458/150‑1 (MQ=255)
                                                                                          cGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGc                                                                                                                              <  1:452883/74‑1 (MQ=255)
                                                                                          cGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGc                                                                                                                              >  2:452883/1‑74 (MQ=255)
                                                                                            ccTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGggcggc                                                                            <  2:185049/122‑1 (MQ=255)
                                                                                            ccTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGggcggc                                                                            >  1:185049/1‑122 (MQ=255)
                                                                                                 gcgTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTca                                                                 <  1:2255963/128‑1 (MQ=255)
                                                                                                 gcgTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTca                                                                 >  2:2255963/1‑128 (MQ=255)
                                                                                                     acacCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAAc                                       >  2:393680/1‑150 (MQ=255)
                                                                                                         cTCGATATTGCGCCGTACAAACAAGGCGGGCGCGCGCGCCGGGCGCAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGAt                                   >  1:810690/1‑150 (MQ=255)
                                                                                                            gATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCaa                                >  1:497380/1‑150 (MQ=255)
                                                                                                             atatTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAAc                               <  2:497380/150‑1 (MQ=255)
                                                                                                             atatTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCAAGCTcgctc                                                                                             >  1:1042208/1‑86 (MQ=18)
                                                                                                             atatTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCAAGCTcgctc                                                                                             <  2:1042208/88‑3 (MQ=255)
                                                                                                                         acaaacaaGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCa                   >  1:1565601/1‑150 (MQ=255)
                                                                                                                              gaaGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGgcgc                    >  2:223502/2‑144 (MQ=255)
                                                                                                                              gaaGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGgcgc                    <  1:223502/143‑1 (MQ=255)
                                                                                                                               aaGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCcg             >  2:1762706/1‑150 (MQ=255)
                                                                                                                               aaGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCCCGCGCGgcgcatc                                                                       >  1:1506116/1‑89 (MQ=21)
                                                                                                                               aaGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCCCGCGCGgcgcatc                                                                       <  2:1506116/92‑4 (MQ=21)
                                                                                                                                    ggTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAgg                                                             <  2:2224615/97‑1 (MQ=255)
                                                                                                                                    ggTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAgg                                                             >  1:2224615/1‑97 (MQ=255)
                                                                                                                                     ttCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGa       >  1:1385852/2‑150 (MQ=255)
                                                                                                                                          gcgcgcCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGc  >  1:52535/1‑150 (MQ=255)
                                                                                                                                                     acaGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCgg   <  1:1538468/138‑1 (MQ=255)
                                                                                                                                                     acaGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCgg   >  2:1538468/1‑138 (MQ=255)
                                                                                                                                                     |                                                                                                                                          
AGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGC  >  NC_002947/4945229‑4945516

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: