Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,356,030	A→G	7.8%	intergenic (+297/‑95)	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_002947	1,356,030	0	A	G	7.8%	185.4 / 5.2	64	intergenic (+297/‑95)	PP_1180/PP_1181	hypothetical protein/two‑component system response regulator
Reads supporting (aligned to +/- strand):  ref base A (28/31);  new base G (2/3);  total (30/34)
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.36e-01

CCTTGTGCCAAGTGGTTTGCAGTATTTTTCTGCACCGTGAAGCCCTTGGGCAGTTTTTTTTGGGCGCGAAGCCTGTCACAAAAAACGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCA‑TCTACATTCA‑TACTCTCAGCTGCCCTGACACGATCCCAAGAACAACAAGGCCAAAGCCATGGACCACCCCATACCCCGCATCCTCATCGTCGAAGACGACCAGCGCCTC  >  NC_002947/1355884‑1356175                                                                                                                                                   |                                                                                                                                                    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                                                                                                                                >  4:63316/1‑125 (MQ=37)                                            ccTTGGGCAGTTTTTTTTGGGCGCGAAGCCTGTCACAAAAAACGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGCGAGACCGAGCGCCCCGCGGGCGGTGCTCGAtc                                                                                                                                <  3:63316/125‑1 (MQ=37)                                             cTTGGGCAGTTTTTTTTGGGCGCGAAGCCTGTCACAAAAAACGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCA‑TCTACATTCa                                                                                                      >  7:80361/1‑149 (MQ=255)                                             cTTGGGCAGTTTTTTTTGGGCGCGAAGCCTGTCACAAAAAACGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCA‑TCTACATTCa                                                                                                      >  3:251926/1‑149 (MQ=255)                                              ttGGGCAGTTTTTTTTGGGCGCGAAGCCTGTCACAAAAAACGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCA‑TCTACATTCa                                                                                                      >  3:239972/1‑148 (MQ=255)                                                   cAGTTTTTTTTGGGCGCGAAGCCTGTCACAAAAAACGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCt                                                                                                                                     >  2:347445/1‑113 (MQ=255)                                                   cAGTTTTTTTTGGGCGCGAAGCCTGTCACAAAAAACGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCt                                                                                                                                     <  1:347445/113‑1 (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:98378/149‑1 (MQ=255)                                                                                                gAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCA‑TCTACATTCA‑TACTCTCAGCTGCCCTGACACGATCCCAAGAACAACAAGGCCAAAGCCATg                                                  >  8:231204/1‑149 (MQ=255)                                                                                                gAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCA‑TCTACATTCA‑TACTCTCAGCTGCCCTGACACGATCCCAAGAACAACAAGGCCAAAGCCATg                                                  <  2:36958/149‑1 (MQ=255)                                                                                                        ggTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCA‑TCTACATTCA‑TACTCTCAGCTGCCCTGACACGATCCCAAGAACAACAAGGCCAAAGCCATGGACCAccc                                          >  2:358966/1‑149 (MQ=255)                                                                                                             ggCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCA‑TCTACATTCA‑TACTCTCAGCTGCCCTGACACGATCCCAAGAACAACAAGGCCAAAGCCATGGac                                               >  3:182290/1‑139 (MQ=255)                                                                                                             ggCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCA‑TCTACATTCA‑TACTCTCAGCTGCCCTGACACGATCCCAAGAACAACAAGGCCAAAGCCATGGac                                               <  4:182290/139‑1 (MQ=255)                                                                                                                   tCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCA‑TCTACATTCA‑TACTCTCAGCTGCCCTGACACGATCCCAAGAACAACAAGGCCAAAGCCATGGACCACCCCATACCCCGc                                <  4:239972/148‑1 (MQ=255)                                                                                                                     tGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCA‑TCTACATTCA‑TACTCTCAGCTGCCCTGACACGATCCCAAGAACAACAAGGCCAAAGCCATGGACCACCCCATACCCCGCATc                             >  7:27796/1‑149 (MQ=255)                                                                                                                     tGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCA‑TCTACATTCA‑TACTCTCAGCTGCCCTGACACGATCCCAAGAACAACAAGGCCAAAGCCATGGACCACCCCATACCCCGCATc                             >  6:57451/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)                                                                                                                                               gAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCA‑TCTACATTCA‑TACTCTCAGCTGCCCTGACACGATCCCAAGAACAACAAGGCCAAAGCCATGGACCACCCCATACCCCGCATCCTCATCGTCGAAGACGACCAGCGCCt   <  8:115443/149‑1 (MQ=255)                                                                                                                                                aGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCA‑TCTACATTCA‑TACTCTCAGCTGCCCTGACACGATCCCAAGAACAACAAGGCCAAAGCCATGGACCACCCCATACCCCGCATCCTCATCGTCGAAGACGACCAGCGCCTc  <  3:149791/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                    CCTTGTGCCAAGTGGTTTGCAGTATTTTTCTGCACCGTGAAGCCCTTGGGCAGTTTTTTTTGGGCGCGAAGCCTGTCACAAAAAACGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCA‑TCTACATTCA‑TACTCTCAGCTGCCCTGACACGATCCCAAGAACAACAAGGCCAAAGCCATGGACCACCCCATACCCCGCATCCTCATCGTCGAAGACGACCAGCGCCTC  >  NC_002947/1355884‑1356175

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