Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,707,335	Δ1 bp	10.0%	intergenic (+111/+31)	PP_1501 → / ← PP_1502	hypothetical protein/OmpA family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,707,335	0	T	.	10.0%	174.7 / 13.7	60	intergenic (+111/+31)	PP_1501/PP_1502	hypothetical protein/OmpA family protein
Reads supporting (aligned to +/- strand):  ref base T (23/31);  new base . (3/3);  total (26/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.85e-01

GGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGC‑GGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCGACGGTTCTGCGCACGTCCCGAGTTGCTCGAGTTGTCCGCTACCGGG  >  NC_002947/1707195‑1707483                                                                                                                                              |                                                                                                                                                      ggCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGC‑GGGGCt                                                                                                                                             >  1:404419/1‑150 (MQ=255)  gCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGC‑GGGGCtt                                                                                                                                            <  1:12223/150‑1 (MQ=255)   cGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGGAAAGCCCCGCGATCGA‑CGC‑GGGGCttt                                                                                                                                           <  2:1560841/150‑2 (MQ=38)      ggAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGC‑GGGGCTTGTcc                                                                                                                                        <  2:1881016/150‑1 (MQ=255)         aaCACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGC‑GGGGCTTGTCCttt                                                                                                                                     <  1:1110675/150‑1 (MQ=255)         aaCACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGC‑GGGGCTTGTCCttt                                                                                                                                     <  1:1183859/150‑1 (MQ=255)          acacCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGC‑GGGGCTTGTCCTTTc                                                                                                                                    <  2:404419/150‑1 (MQ=255)           cacCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGC‑GGGGCTTGTCCTTTCa                                                                                                                                   >  1:1302791/1‑150 (MQ=255)            acCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGC‑gggg                                                                                                                                               <  2:1855173/137‑1 (MQ=255)            acCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGC‑gggg                                                                                                                                               >  1:1855173/1‑137 (MQ=255)            acCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGC‑GGGGCTTGTCCTTTCAg                                                                                                                                  >  2:1517035/1‑150 (MQ=255)            acCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGC‑GGGGCTTGTCCTTTCAg                                                                                                                                  >  1:23407/1‑150 (MQ=255)            acCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGC‑GGGGCTTGTCCTTTCAg                                                                                                                                  >  1:997932/1‑150 (MQ=255)              cgcgGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGACGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAGGCCCCGCG‑TCGATCGCGGGGGCTTGTCCTTTCAGc                                                                                                                                 >  1:167399/1‑150 (MQ=255)               gcgGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGC‑GGGGCTTGTCCTTTCAGCTa                                                                                                                               <  1:1221938/150‑1 (MQ=255)                           gTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGC‑GGGGCTTGTCCTTTCATCTAGCCAGCCTTAgc                                                                                                                   >  2:1844312/1‑150 (MQ=255)                             gAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGC‑GGGGCTTGTCCTTTCAGCTAGCCAGCCTTAgcgc                                                                                                                 <  2:1788106/150‑1 (MQ=255)                                            caTTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGC‑GGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCAcc                                                                                                  <  2:997932/150‑1 (MQ=255)                                              tttttGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGC‑GGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGa                                                                                                <  1:1813087/150‑1 (MQ=255)                                                  tGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCTATCTCACAGGCGCCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑GCGATCGC‑GGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCgt                                                                                            <  1:580437/150‑1 (MQ=255)                                                   ggCGCCCGCCATATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGC‑GGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCgtg                                                                                           <  1:1516010/150‑1 (MQ=255)                                                             agaTCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGC‑GGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCgg                                                                                 <  2:485108/150‑1 (MQ=255)                                                             agaTCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGC‑GGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCgg                                                                                 <  2:485099/150‑1 (MQ=255)                                                                tCGAGCGCCGCCCGCGCGGCGCTCGATCTCGCGGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGC‑GGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTg                                                                              >  1:388667/1‑150 (MQ=39)                                                                tCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATcgc‑g                                                                                                                                                  <  1:764594/82‑1 (MQ=255)                                                                tCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATcgc‑g                                                                                                                                                  >  2:764594/1‑82 (MQ=255)                                                                 cGAGCGCCGCCCGCGCGGCGCTCGATCTCGCAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGC‑GGGGCTTGTCCttt                                                                                                                                     <  2:1240473/94‑1 (MQ=38)                                                                 cGAGCGCCGCCCGCGCGGCGCTCGATCTCGCAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGC‑GGGGCTTGTCCttt                                                                                                                                     >  1:1240473/1‑94 (MQ=38)                                                                   aGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGC‑GGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGtt                                                                           <  1:345175/150‑1 (MQ=255)                                                                       ccgcGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGC‑GGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCGGGCTCACCGGGCGTGGTGCCACCGGCTGGTTGTCg                                                                       >  2:1535347/1‑150 (MQ=255)                                                                         gcgcgGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGC‑GGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGtt                                                                     >  2:583766/1‑150 (MQ=255)                                                                         gcgcgGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGC‑GGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGtt                                                                     <  1:1494178/150‑1 (MQ=255)                                                                            cgggcggCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGC‑GGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGga                                                                  <  2:499156/150‑1 (MQ=255)                                                                                cggcgCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGC‑GGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATg                                                              <  1:1319337/150‑1 (MQ=255)                                                                                     cTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGC‑GGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCAc                                                         <  1:327166/150‑1 (MQ=255)                                                                                      tCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGC‑GGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGAt                                                               >  1:808721/1‑143 (MQ=255)                                                                                      tCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGC‑GGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGAt                                                               <  2:808721/143‑1 (MQ=255)                                                                                          tctcACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGC‑GGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCa                                                    >  2:771529/1‑150 (MQ=255)                                                                                             cacaGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGC‑GGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAACTGGTCACCTCCAccc                                                 >  1:282124/1‑150 (MQ=255)                                                                                              acaGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGC‑GGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCg                                                >  1:381821/1‑150 (MQ=255)                                                                                                aGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCGa                                               <  2:167399/150‑1 (MQ=255)                                                                                                  gcggcACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGC‑GGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCGACgg                                            >  1:1144289/1‑150 (MQ=255)                                                                                                  gcggcACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGC‑GGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCGACgg                                            >  1:560054/1‑150 (MQ=255)                                                                                                               aTCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGC‑GGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCGAAGGTTCTGCGCACGTc                               <  1:1844312/150‑1 (MQ=255)                                                                                                                         aaaGAAAAGCCCCGCG‑TCGATCGC‑GGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCGACGGTTCTGCGCACGTCCCGAGTTGCt                     >  1:345339/1‑150 (MQ=255)                                                                                                                            gAAAAGCCCCGCG‑TCGATCGC‑GGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGTTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCGACGGTTCTGCGCACGTCCCGAGTTGCTCGa                  <  2:42468/150‑1 (MQ=255)                                                                                                                              aaaGCCCCGCG‑TCGATCGC‑GGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCGACGGTTCTGCGCACGTCCCGAGTTGCTCGAGt                <  2:1560443/150‑1 (MQ=255)                                                                                                                               aaGCCCCGCGATCGA‑CGC‑GGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCTAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATg                                                              >  1:1148038/1‑103 (MQ=255)                                                                                                                               aaGCCCCGCGATCGA‑CGC‑GGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATg                                                              <  2:1148038/103‑1 (MQ=255)                                                                                                                               aaGCCCCGCGATCGA‑CGC‑GGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCGACGGTTCTGCGCACGTCCCGAGTTGCTCGAGtt               >  1:1620452/1‑150 (MQ=255)                                                                                                                               aaGCCCCGCG‑TCGATCGC‑GGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCGACGGTTCTGCGCACGTCCCGAGTTGCTCGAGtt               <  2:381821/150‑1 (MQ=255)                                                                                                                                 gCCCCGCGATCGA‑CGC‑GGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGAt                                                               <  2:1908221/100‑1 (MQ=255)                                                                                                                                 gCCCCGCGATCGA‑CGC‑GGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGAt                                                               >  1:1908221/1‑100 (MQ=255)                                                                                                                                   cccGCG‑TCGATCGC‑GGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCGACGGTTCTGCGCACGTCCCGAGTTGCTCGAGTTGTcc           <  1:771529/150‑1 (MQ=255)                                                                                                                                    ccGCG‑TCGATCGC‑GGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCGACGGTTCTGCGCACGTCCCGAGTTGCTCGAGTTGTCCg          >  1:324412/1‑150 (MQ=255)                                                                                                                                    ccGCG‑TCGATCGC‑GGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCGACGGTTCTGCGCACGTCCCGAGTTGCTCGAGTTGTCCg          >  2:1557438/1‑150 (MQ=255)                                                                                                                                       cg‑TCGATCGC‑GGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTgg                                                                                                        <  1:69994/53‑1 (MQ=255)                                                                                                                                       cg‑TCGATCGC‑GGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTgg                                                                                                        >  2:69994/1‑53 (MQ=255)                                                                                                                                           cGATCGC‑GGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCATAATCTGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCGACGGTTCTGCGCACGTCCCGAGTTGCTCGAGTTGTCCGCTACCg    <  1:497676/150‑1 (MQ=255)                                                                                                                                             aTCGC‑GGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCGACGGTTCTGCGCACGTCCCGAGTTGCTCGAGTTGTCCGCTACCggg  <  1:583766/150‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                      GGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGC‑GGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCGACGGTTCTGCGCACGTCCCGAGTTGCTCGAGTTGTCCGCTACCGGG  >  NC_002947/1707195‑1707483

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