Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,356,024 T→C 7.8% intergenic (+291/‑101) PP_1180 → / → PP_1181 hypothetical protein/two‑component system response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,356,0240TC 7.8% 185.6 / 4.5 64intergenic (+291/‑101)PP_1180/PP_1181hypothetical protein/two‑component system response regulator
Reads supporting (aligned to +/- strand):  ref base T (29/30);  new base C (2/3);  total (31/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

CTGTAGGCCTTGTGCCAAGTGGTTTGCAGTATTTTTCTGCACCGTGAAGCCCTTGGGCAGTTTTTTTTGGGCGCGAAGCCTGTCACAAAAAACGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCA‑TCTACATTCA‑TACTCTCAGCTGCCCTGACACGATCCCAAGAACAACAAGGCCAAAGCCATGGACCACCCCATACCCCGCATCCTCATCGTCGAAGACGAC  >  NC_002947/1355877‑1356166
                                                                                                                                                   |                                                                                                                                                
cTGTAGGCCTTGTGCCAAGTGGTTTGCAGTATTTTTCTGCACCGTGAAGCCCTTGGGCAGTTTTTTTTGGGCGCGAAGCCTGTCACAAAAAACGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATc                                                                                                                                                 <  4:346836/149‑1 (MQ=255)
  gTAGGCCTTGTGCCAAGTGGTTTGCAGTATTTTTCTGCACCGTGAAGCCCTTGGGCAGTTTTTTTTGGGCGCGAAGCCTGTCACAAAAAACGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATc                                                                                                                                                 >  7:323862/1‑147 (MQ=255)
       ccTTGTGCCAAGTGGTTTGCAGTATTTTTCTGCACCGTGAAGCCCTTGGGCAGTTTTTTTTGGGCGCGAAGCCTGTCACAAAAAACGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCAcc                                                                                                                                          >  2:274341/1‑149 (MQ=255)
        cTTGTGCCAAGTGGTTTGCAGTATTTTTCTGCACCGTGAAGCCCTTGGGCAGTTTTTTTTGGGCGCGAAGCCTGTCACAAAAAACGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCAccg                                                                                                                                         >  7:98378/1‑149 (MQ=255)
                  gTGGTTTGCAGTATTTTTCTGCACCGTGAAGCCCTTGGGCAGTTTTTTTTGGGCGCGAAGCCTGTCACAAAAAACGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGgc                                                                                                                               <  8:323862/149‑1 (MQ=255)
                      tttGCAGTATTTTTCTGCACCGTGAAGCCCTTGGGCAGTTTTTTTTGGGCGCGAAGCCTGTCACAAAAAACGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGTGCACCGCCCGCGGGGCGCTc                                                                                                                           <  6:384858/149‑1 (MQ=255)
                        tGCAGTATTTTTCTGCACCGTGAAGCCCTTGGGCAGTTTTTTTTGGGCGCGAAGTCTGTCACAAAAAACGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCgg                                                                                                                         >  2:299183/1‑149 (MQ=255)
                                ttttCTGCACCGTGACGCCCTTGGGCAGTTTTTTTTGGGCGCGAAGCCTGTCACAAAAAACGGTGCAGTCGAGGGGGCGGTAGGGCGTGGCTGAAGATCTGGGGTGCCTGCGAGACCGAGCGCCCCGCGGGCGGTGCTCGATCTCCCAg                                                                                                                 <  5:377197/149‑1 (MQ=25)
                                ttttCAGCACCGTGATGCCCTTGGGCAGTTTTTTTTGGGCGCGAAGCCTGTCACAAAAAACGGTTCAGTCGTGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGCGAGACCGAGCGCCCCGCGGGCGGTGCTCGATCTCCCa                                                                                                                  <  8:231527/148‑1 (MQ=25)
                                   tCTGCACCGTGAAGCCCTTGGGCAGTTTTTTTTGGGCGCGAAGCCTGTCACAAAAAACGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCggg                                                                                                                                 <  5:345643/130‑1 (MQ=255)
                                   tCTGCACCGTGAAGCCCTTGGGCAGTTTTTTTTGGGCGCGAAGCCTGTCACAAAAAACGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCggg                                                                                                                                 >  6:345643/1‑130 (MQ=255)
                                   tCTGCACCGTGAAGCCCTTGGGCAGTTTTTTTTGGGCGCGAAGCCGGTCACAAAAAACGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAACATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAggg                                                                                                               <  4:344125/148‑1 (MQ=255)
                                       cACCGTGAAGCCCTTGGGCAGTTTTTTTTGGGCGCGAAGCCTGTCACAAAAAACGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGcccc                                                                                                           <  4:116068/148‑1 (MQ=255)
                                        aCCGTGAAGCCCTTGGGCAGTTTTTTTTGGGCGCGAAGCCTGTCACAAAAAACGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGccccc                                                                                                          >  3:177538/1‑148 (MQ=255)
                                          cGTGAAGCCCTTGGGCAGTTTTTTTTGGGCGCGAAGCCTGTCACAAAAAACGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCA‑t                                                                                                      <  5:362375/149‑1 (MQ=255)
                                           gTGAAGCCCTTGGGCAGTTTTTTTTGGGCGCGAAGCCTGTCACAAAAAACGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCA‑Tc                                                                                                     >  8:242742/1‑149 (MQ=255)
                                            tGAAGCCCTTGGGCAGTTTTTTTTGGGCGCGAAGCCTGTCACAAAAAACGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCA‑Tc                                                                                                     >  1:37475/1‑148 (MQ=255)
                                            tGAAGCCCTTGGGCAGTTTATTTTGGGCGCGAAGCCTGTCACAAAAAACGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCA‑TCt                                                                                                    <  1:124080/149‑1 (MQ=255)
                                              aaGCCCTTGGGCAGTTTTTTTTGGGCGCGAAGCCTGTCACAAAAAACGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCA‑TCTa                                                                                                   <  7:242742/148‑1 (MQ=255)
                                                  ccTTGGGCAGTTTTTTTTGGGCGCGAAGCCTGTCACAAAAAACGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGCGAGACCGAGCGCCCCGCGGGCGGTGCTCGAtc                                                                                                                       <  3:63316/125‑1 (MQ=37)
                                                  ccTTGGGCAGTTTTTTTTGGGCGCGAAGCCTGTCACAAAAAACGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGCGAGACCGAGCGCCCCGCGGGCGGTGCTCGAtc                                                                                                                       >  4:63316/1‑125 (MQ=37)
                                                   cTTGGGCAGTTTTTTTTGGGCGCGAAGCCTGTCACAAAAAACGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCA‑TCTACATTCa                                                                                             >  3:251926/1‑149 (MQ=255)
                                                   cTTGGGCAGTTTTTTTTGGGCGCGAAGCCTGTCACAAAAAACGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCA‑TCTACATTCa                                                                                             >  7:80361/1‑149 (MQ=255)
                                                    ttGGGCAGTTTTTTTTGGGCGCGAAGCCTGTCACAAAAAACGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCA‑TCTACATTCa                                                                                             >  3:239972/1‑148 (MQ=255)
                                                         cAGTTTTTTTTGGGCGCGAAGCCTGTCACAAAAAACGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCt                                                                                                                            <  1:347445/113‑1 (MQ=255)
                                                         cAGTTTTTTTTGGGCGCGAAGCCTGTCACAAAAAACGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCt                                                                                                                            >  2:347445/1‑113 (MQ=255)
                                                         cAGTTTTTTTTGGGCGCGAAGCCTGTCACAAAAAACGGTTCAGGCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCA‑TCTACATTCA‑TActct                                                                                      <  2:279773/149‑1 (MQ=255)
                                                          aGTTTTTTTTGGGCGCGAAGCCTGTCACAAAAAACGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCA‑TCTACATTCA‑TACtctc                                                                                     <  6:136018/149‑1 (MQ=255)
                                                           gTTTTTTTTTGGCGCGAAGCCTGTCACTAAAAACGGTTCAGTCTAGGGGGCGGTAGGGCTTGTCTGAAGATCTGGTGTGCCTGGGATATCGAGCACCGCCCTCGGGGCGCTCGGTCTAGCAGGGCCCCCCA‑TCTACATTAA‑TACTCTCa                                                                                    <  5:212591/149‑1 (MQ=255)
                                                            ttttttttGGGCGCGAAGCCTTTCACAAAAAACGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCt                                                                                                                            <  4:24997/110‑1 (MQ=255)
                                                            ttttttttGGGCGCGAAGCCTGTCACAAAAAACGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCt                                                                                                                            >  3:24997/1‑110 (MQ=255)
                                                                     ggCGCGAAGCCTGTCACAAAAAACGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCA‑TCTACATTCA‑TACTCTCAGCTGCCCTGa                                                                          <  8:87840/149‑1 (MQ=255)
                                                                             gCCTGTCACAAAAAACGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCA‑TCTACATTCA‑TACTCTCAGCTGCCCTGACACGATcc                                                                  <  2:37475/149‑1 (MQ=255)
                                                                              ccTGTCAGAAATAAGGGTTCTGTCGAGGGGGCGGTAGGGCGTGTCTGATGATCTGGGGTGCCTGTGATATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCA‑TCTACATTCT‑TACTCTCAGCTGCCCTGACACGATccc                                                                 <  2:250264/149‑1 (MQ=255)
                                                                               cTGTCACAAAAAACGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCA‑TCTACATTCA‑TACTCTCAGCTGCCCTGACACGATccc                                                                 >  4:245758/1‑148 (MQ=255)
                                                                                 gTCACAAAAAACGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCA‑TCTACATTCA‑TACTCTCAGCTGCCCTGACACGATCCCaa                                                               >  7:328565/1‑148 (MQ=255)
                                                                                    acaAAAAACGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGTTGCCTGTTATATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCA‑TCTACATTCA‑TACTCTCAGCTGCCCTGACACGATCCCAAGaa                                                            <  1:299183/148‑1 (MQ=255)
                                                                                        aaaaCGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCA‑TCTACATTCA‑TACtctc                                                                                     <  6:346115/119‑1 (MQ=255)
                                                                                        aaaaCGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCA‑TCTACATTCA‑TACtctc                                                                                     >  5:346115/1‑119 (MQ=255)
                                                                                          aaCGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGGCTCGCAGGGCCCCCCACTCTACATTCA‑TACTCTCAGCTGCCCTGACACGATCCCAAGAAcaacaa                                                      <  1:274341/149‑1 (MQ=255)
                                                                                           aCGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCA‑TCTACATTCA‑TACTCTCAGCTGCCCTGACACGATCCCAAGAACAACAAgg                                                    <  4:251926/149‑1 (MQ=255)
                                                                                                      gAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCA‑TCTACATTCA‑TACTCTCAGCTGCCCTGACACGATCCCAAGAACAACAAGGCCAAAGCCATg                                         >  8:231204/1‑149 (MQ=255)
                                                                                                      gAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCA‑TCTACATTCA‑TACTCTCAGCTGCCCTGACACGATCCCAAGAACAACAAGGCCAAAGCCATg                                         <  8:98378/149‑1 (MQ=255)
                                                                                                      gAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCA‑TCTACATTCA‑TACTCTCAGCTGCCCTGACACGATCCCAAGAACAACAAGGCCAAAGCCATg                                         <  2:36958/149‑1 (MQ=255)
                                                                                                              ggTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCA‑TCTACATTCA‑TACTCTCAGCTGCCCTGACACGATCCCAAGAACAACAAGGCCAAAGCCATGGACCAccc                                 >  2:358966/1‑149 (MQ=255)
                                                                                                                   ggCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCA‑TCTACATTCA‑TACTCTCAGCTGCCCTGACACGATCCCAAGAACAACAAGGCCAAAGCCATGGac                                      <  4:182290/139‑1 (MQ=255)
                                                                                                                   ggCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCA‑TCTACATTCA‑TACTCTCAGCTGCCCTGACACGATCCCAAGAACAACAAGGCCAAAGCCATGGac                                      >  3:182290/1‑139 (MQ=255)
                                                                                                                         tCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCA‑TCTACATTCA‑TACTCTCAGCTGCCCTGACACGATCCCAAGAACAACAAGGCCAAAGCCATGGACCACCCCATACCCCGc                       <  4:239972/148‑1 (MQ=255)
                                                                                                                           tGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCA‑TCTACATTCA‑TACTCTCAGCTGCCCTGACACGATCCCAAGAACAACAAGGCCAAAGCCATGGACCACCCCATACCCCGCATc                    >  6:57451/1‑149 (MQ=255)
                                                                                                                           tGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCA‑TCTACATTCA‑TACTCTCAGCTGCCCTGACACGATCCCAAGAACAACAAGGCCAAAGCCATGGACCACCCCATACCCCGCATc                    >  7:27796/1‑149 (MQ=255)
                                                                                                                            gAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCA‑TCTACATTCAGTACTCTCAGCTGCCCTGACACGATCCCAAGa                                                             <  1:39486/108‑1 (MQ=255)
                                                                                                                            gAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCA‑TCTACATTCA‑TACTCTCAGCTGCCCTGACACGATCCCAAGaacaac                                                        >  6:246860/1‑112 (MQ=255)
                                                                                                                            gAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCA‑TCTACATTCA‑TACTCTCAGCTGCCCTGACACGATCCCAAGaacaac                                                        <  5:246860/112‑1 (MQ=255)
                                                                                                                            gAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCA‑TCTACATTCA‑TACTCTCAGCTGCCCTGACACGATCCCAAGAACAACAAGGCCAAAGCCATg                                         <  6:339690/127‑1 (MQ=255)
                                                                                                                            gAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCA‑TCTACATTCA‑TACTCTCAGCTGCCCTGACACGATCCCAAGAACAACAAGGCCAAAGCCATGGACCACCCCATACCCCGCATcc                   >  5:302744/1‑149 (MQ=255)
                                                                                                                            gAAGATCTGGGGTGCCTCGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCA‑TCTACATTCA‑TACTCTCAGCTGCCCTGACACGATCCCAAGAACAACAAGGCCAAAGCCATGGACAACCCCATACCCCGCATc                    >  4:149791/1‑148 (MQ=255)
                                                                                                                             aaGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCA‑TCTACATTCA‑TACTCTCAGCTGCCCTGACACGATCCCAAGAACAACAAGGCCAAAGCCATGGACCACCCCATACCCCGCATCCt                  >  7:115443/1‑149 (MQ=255)
                                                                                                                              aGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCA‑TCTACATTCA‑TACTCTCAGCTGCCCTGACACGATCCCAAGAACAACAAGGCCAAAGCCATGGACCACCCCATACCCCGCATCCt                  >  6:19839/1‑148 (MQ=255)
                                                                                                                               gATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCA‑TCTACATTCA‑TACTCTCAGCTGCCCTGACACGATCCCAAGAACAACAAGGCCAAAGCCATGGACCACCCCATACCCCGCATCCTc                 >  8:31017/1‑148 (MQ=255)
                                                                                                                                 tCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCA‑TCTACATTCA‑TACTCTCAGCTGCCCTGACACGATCCCAAGAACAACAAGGCCAAAGCCATGGACCACCCCATACCCCGCCTCCTCAtc              >  4:175962/1‑149 (MQ=255)
                                                                                                                                      gggcccTGCGAGACCGAGCGCCCCGCGGGCGGTGCTCGATCTCCCA‑GGCACCCCA‑TCTACATTCA‑TACTCTCAGCTGCCCTGACACGATCCCAAGAACAACAAGGCCAAAGCCATGGACCACCCCATACCCCGCATCCTCAtcgtcg          >  6:312170/5‑147 (MQ=21)
                                                                                                                                      ggTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCA‑TCTACATTCA‑TACTCTCAGCTGCCCTGACACGATCCCAAGAACAACAAGGCCAAAGCCATGGACCACCCCATACCCCGCATCCTCATCGTCGa         <  1:358966/149‑1 (MQ=255)
                                                                                                                                            tGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCA‑TCTACATTCA‑TACTCTCAGCTGCCCTGACACGATCCCAAGAACAACAAGGCCAAAGCCATGGACCACCCCATACCCCGCATCCTCATCGTCGAAgacga   >  2:314589/1‑149 (MQ=255)
                                                                                                                                              ggAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCA‑TCTACATTCA‑TACTCTCAGCTGCCCTGACACGATCCCAAGAACAACAAGGCCAAAGCCATGGACCACCCCATACCCCGCATCCTCATCGTCGAAgacgac  <  7:31017/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
CTGTAGGCCTTGTGCCAAGTGGTTTGCAGTATTTTTCTGCACCGTGAAGCCCTTGGGCAGTTTTTTTTGGGCGCGAAGCCTGTCACAAAAAACGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCA‑TCTACATTCA‑TACTCTCAGCTGCCCTGACACGATCCCAAGAACAACAAGGCCAAAGCCATGGACCACCCCATACCCCGCATCCTCATCGTCGAAGACGAC  >  NC_002947/1355877‑1356166

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