Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,050,088	Δ1 bp	8.7%	intergenic (+17/‑124)	PP_2662 → / → PP_2663	hypothetical protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,050,088	0	G	.	8.7%	201.4 / 15.1	69	intergenic (+17/‑124)	PP_2662/PP_2663	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (29/33);  new base . (3/3);  total (33/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.18e-01

GCTGTTCCATGACATCAACGACAACCTGAAACTGGTGGCCGAGTACAACCAGTTTTCCATCGACGGGCATGACACCAGTGCGCAGAACGAAGACACCGATACCTTTGCGGTGGGCGCGGTGTTGACCTGGTAAGCGGCCGGGGCCG‑CTAGGCGGCCCCGAATGACACGGTCTTTTACGCCCATGGAACCGCCCACCCGCTACCCAAGTAGCATTCGTCACCCACCACCACCTCCTTAAACTGATCACTCAGCCACCCCGTGCCGGAGACCACGATGCAGCAGGCTCATAACGACGGT  >  NC_002947/3049939‑3050235                                                                                                                                                       |                                                                                                                                                    gcTGTTCCATGACATCAACGACAACCTGAAACTGGTGGCCGAGTACAACCAGTTTTCCATCGACGGGCATGACACCAGTGCGCAGAACGAAGACACCGATACCTTTGCGGTGTGCGCGGTGTTGACCTGGTAAGCGGCCGGGGCCG‑CTAg                                                                                                                                                     >  1:865220/1‑150 (MQ=255) gcTGTTCCATGACATCAACGACAACCTGAAACTGGTGGCCGAGTACAACCAGTTTTCCATCGACGGGCATGACACCAGTGCGCAGAACGAAGACACCGATACCTTTGCGGTGGGCGCGGTGTTGACCTGGTAAGCGGCCGGGGCCG‑CTAg                                                                                                                                                     >  2:888761/1‑150 (MQ=255)    gTTCCATGACATCAACGACAACCTGAAACTGGTGGCCGAGTACAACCAGTTTTCCATCGACGGGCATGACACCAGTGCGCAGAACGAAGACACCGATACCTTTGCGGTGGGCGCGGTGTTGACCTGGTAAGCGGCCGGGGCCG‑CTAggcg                                                                                                                                                  >  1:1787670/1‑150 (MQ=255)              aTCAACGACAACCTGAAACTGGTGGCCGAGTACAACCAGTTTTCCATCGACGGGCATGACACCAGTGCGCAGAACGAAGACACCGATACCTTTGCGGTGGGCGCGGTGTTGACCTGGTAAGCGGCCGGGGCCG‑CTAGGCGGCCCCGAATg                                                                                                                                        >  2:947606/1‑150 (MQ=255)               tCAACGACAACCTGAAACTGGTGGCCGAGTACAACCAGTTTTCCATCGACGGGCATGACACCAGTGCGCAGAACGAAGACACCGATACCTTTGCGGTGGGCGCGGTGTTGACCTGGTAAGCGGCCGGGGCCG‑CTAGGCGGCCCCGAATGa                                                                                                                                       >  1:1332563/1‑150 (MQ=255)                cAACGACAACCTGAAACTGGTGGCCGAGTACAACCAGTTTTCCATCGACGGGCATGACAGCAGTGCGCAGAACGAAGACACCGAGGCCTTTGCGGTGGGCGCGGTGTTGACCTGGTAAGCGGCCGGGGCCG‑CTAGGCGGCCCCGAATGac                                                                                                                                      >  1:707002/1‑150 (MQ=255)                 aaCGACAACCTGAAACTGGTGGCCGAGTACAACCAGTTTTCCATCGACGGGCATGACACCAGTGCGCAGAACGAAGACACCGATACCTTTGCGGTGGGCGCGGTGTTGACCTGGTAAGCGGCCGGGGCCG‑CTAGGCGGCCCCGAATGaca                                                                                                                                     >  1:1178431/1‑150 (MQ=255)                 aaCGACAACCTGAAACTGGTGGCCGAGTACAACCAGTTTTCCATCGACGGGCATGACACCAGTGCGCAGAACGAAGACACCGATACCTTTGCGGTGGGCGCGGTGTTGACCTGGTAAGCGGCCGGGGCCG‑CTAGGCGGCCCCGAATGaca                                                                                                                                     >  2:1405647/1‑150 (MQ=255)                                  ggtggCCGAGTACAACCAGTTTTCCATCGACGGGCATGACACCAGTGCGCAGAACGAAGACACCGATACCTTTGCGGTGGGCGCGGTGTTGACCTGGTAAGCGGCCGGGGCCG‑CTAGGCGGCCCCGAATGACACGGTCTTTTACGCCCAt                                                                                                                    <  1:1386427/150‑1 (MQ=255)                                      gCCGAGTACAACCAGTTTTCCATCGACGGGCATGACACCAGTGCGCAGAACGAAGACACCGATACCTTTGCGGTGGGCGCGGTGTTGACCTGGTAAGCGGCCGGGGCCG‑CTAGGCGGCCCCGAATGACACGGTCTTTTACGCCCATGGaa                                                                                                                <  2:1178431/150‑1 (MQ=255)                                          aGTACAACCAGTTTTCCATCGACGGGCATGACACCAGTGCGCAGAACGAAGACACCGATACCTTTGCGGTGGGCGCGGTGTTGACCTGGTAAGCGGCCGGGGCCG‑CTAGGCGGCCCCGAATGACACGGTCTTTTACGCCCATGGAACCGc                                                                                                            <  1:1004113/150‑1 (MQ=255)                                               aaCCAGTTTTCCATCGACGGGCATGACACCAGTGCGCAGAACGAAGACACCGATACCTTTGCGGTGGGCGCGGTGTTGACCTGGTAAGCGGCCGGGGCCG‑CTAGGCGGCCCCGAAt                                                                                                                                         >  1:190998/1‑116 (MQ=255)                                               aaCCAGTTTTCCATCGACGGGCATGACACCAGTGCGCAGAACGAAGACACCGATACCTTTGCGGTGGGCGCGGTGTTGACCTGGTAAGCGGCCGGGGCCG‑CTAGGCGGCCCCGAAt                                                                                                                                         <  2:190998/116‑1 (MQ=255)                                               aaCCAGTTTTCCATCGACGGGCATGACACCAGTGCGCAGAACGAAGACACCGATACCTTTGCGGTGGGCGCGGTGTTGACCTGGTAAGCGGCCGGGGCCG‑CTAGGCGGCCCCGAATGACACGGTCTTTTACGCCCATGGaa                                                                                                                >  1:1644383/1‑141 (MQ=255)                                               aaCCAGTTTTCCATCGACGGGCATGACACCAGTGCGCAGAACGAAGACACCGATACCTTTGCGGTGGGCGCGGTGTTGACCTGGTAAGCGGCCGGGGCCG‑CTAGGCGGCCCCGAATGACACGGTCTTTTACGCCCATGGaa                                                                                                                <  2:1644383/141‑1 (MQ=255)                                               aaCCAGTTTTCCATCGACGGGCATGACACCAGTGCGCAGAACGAAGACACCGATACCTTTGCGGTGGGCGCGGTGTTGACCTGGTAAGCGGCCGGGGCCG‑CTAGGCGGCCCCGAATGACACGGTCTTTTACGCCCATGGAACCGCCCAcc                                                                                                       >  2:1627435/1‑150 (MQ=255)                                               aaCCAGTTTTCCATCGACGGGCATGACACCAGTGCGCAGAACGAAGACACCGATACCTTTGCGGTGGGCGCGGTGTTGACCTGGTAAGCGGCCGGGGCCG‑CTAGGCGGCCCCGAATGACACGGTCTTTTACGCCCATGGAACCGCCCAcc                                                                                                       >  2:753426/1‑150 (MQ=255)                                                aCCAGTTTTCCATCGACGGGCATGACACCAGTGCGCAGAACGAAGACACCGATACCTTTGCGGTGGGCGCGGTGTTGACCTGGTAAGCGGCCGGGGCCG‑CTAGGCGGCCCCGAATGACACGGTCTTTTACGCCCATGGAACCGCCCAccc                                                                                                      >  1:1918362/1‑150 (MQ=255)                                                 ccAGTTTTCCATCGACGGGCATGACACCAGTGCGCAGAACGAAGACACCGATACCTTTGCGGTGGGCGCGGTGTTGACCTGGTAAGCGGCCGGGGCCG‑CTAGGCGGCCCCGAATGACACGGTCTTTTACGCCCATGGAACCGCCCACCCg                                                                                                     >  1:772539/1‑150 (MQ=255)                                                    gTTTTCCATCGACGGGCATGACACCAGTGCGCAGAACGAAGACACCGATACCTTTGCGGTGGGCGCGGTGTTGACCTGGTAAGCGGCCGGGGCCG‑CTAGGCGGCCCCGAATGACACGGTCTTTTACGCCCATGGAACCGCCCACCCGCTa                                                                                                  <  2:984108/150‑1 (MQ=255)                                                     ttttCCATCGACGGGCATGACACCAGTGCGCAGAACGAAGACACCGATACCTTTGCGGTGGGCGCGGTGTTGACCTGGTAAGCGGCCGGGGCCG‑CTAGGCGGCCCCGAATGACACGGTCTTTTACGCCCATGGAACCGCCCACCCGCTAc                                                                                                 >  1:1073856/1‑150 (MQ=255)                                                          cATCGACGGGCATGACACCAGTGCGCAGAACGAAGACACCGATACCTTTGCGGTGGGCGCGGTGTTGACCTGGTAAGCGGCCGGGGCCGCCTA‑GCGGCCCCg                                                                                                                                            <  1:766535/102‑1 (MQ=255)                                                          cATCGACGGGCATGACACCAGTGCGCAGAACGAAGACACCGATACCTTTGCGGTGGGCGCGGTGTTGACCTGGTAAGCGGCCGGGGCCGCCTA‑GCGGCCCCg                                                                                                                                            >  2:766535/1‑102 (MQ=255)                                                              gACGGGCATGACACCAGTGCGCAGAACGAAGACACCGATACCTTTGCGGTGGGCGCGGTGTTGACCTGGTAAGCGGCCGGGGCCG‑CTAGGCGGCCCCGAATGACACGGTCTTTTACGCCCATGGAACCGCCCACCCGCTACCCAAGTAGc                                                                                        >  1:710078/1‑150 (MQ=255)                                                                  ggCATGACACCAGTGCGCAGAACGAAGACACCGATACCTTTGCGGTGGGCGCGGTGTTGACCTGGTAAGCGGCCGGGGCCG‑CTAGGCGGCCCCGAATGACACGGTCTTTTACGCCCATGGAACCGCCCACCCGCTACCCAAGTAGCATTc                                                                                    <  1:410241/150‑1 (MQ=255)                                                                     aTGACACCAGTGCGCAGAACGAAGACACCGATACCTTTGCGGTGGGCGCGGTGTTGACCTGGTAAGCGGCCGGGGCCG‑CTAGGCGGCCCCGAATGACACGGTCTTTTACGCCCATGGAACCGCCCACCCGCTACCCAAGTAGCATTCGTc                                                                                 >  2:771313/1‑150 (MQ=255)                                                                          acCAGTGCGCAGAACGAAGACACCGATACCTTTGCGGTGGGCGCGGTGTTGACCTGGTAAGCGGCCGGGGCCG‑CTAGGCGGCCCCGAATGACACGGTCTTTTACGCCCAt                                                                                                                    >  1:567269/1‑110 (MQ=255)                                                                          acCAGTGCGCAGAACGAAGACACCGATACCTTTGCGGTGGGCGCGGTGTTGACCTGGTAAGCGGCCGGGGCCG‑CTAGGCGGCCCCGAATGACACGGTCTTTTACGCCCAt                                                                                                                    <  2:567269/110‑1 (MQ=255)                                                                          acCAGTGCGCAGAACGAAGACACCGATACCTTTGCGGTGGGCGCGGTGTTGACCTGGTAAGCGGCCGGGGCCG‑CTAGGCGGCCCCGAATGACACGGTCTTTTACGCCCATGGAACCGCCCACCCGCTACCCAAGTAGCATTCGTcaccca                                                                            <  1:216716/150‑1 (MQ=255)                                                                             aGTGCGCAGAACGAAGACACCGATACCTTTGCGGTGGGCGCGGTGTTGACCTGGTAAGCGGCCGGGGCCG‑CTAGGCGGCCCCGAATGACACGGTCTTTTACGCCCATGGAACCGCCCACCCGCTACCCAAGTAGCATTCGTCACccacca                                                                         >  1:953098/1‑150 (MQ=255)                                                                                  gcAGAACGAAGACACCGATACCTTTGCGGTGGGCGCGGTGTTGACCTGGTAAGCGGCCGGGGCCG‑CTAGGCGGCCCCGAATGACACGGTCTTTTACGCCCATGGAACCGCCCACCCGCTACCCAAGTAGCATTCGTCACCCaccaccacc                                                                    >  1:1691860/1‑150 (MQ=255)                                                                                  gcAGAACGAAGACACCGATACCTTTGCGGTGGGCGCGGTGTTGACCTGGTAAGCGGCCGGGGCCG‑CTAGGCGGCCCCGAATGACACGGTCTTTTACGCCCATGGAACCGCCCACCCGCTACCCAAGTAGCATTCGTCACCCaccaccacc                                                                    >  2:1258591/1‑150 (MQ=255)                                                                                    aGAACGAAGACACCGATACCTTTGCGGTGGGCGCGGTGTTGACCTGGTAAGCGGCCGGGGCCG‑CTAGGCGGCCCCGAATGACACGGTCTTTTACGCCCATGGAACCGCCCACCCGCTAc                                                                                                 <  2:1695584/119‑1 (MQ=255)                                                                                    aGAACGAAGACACCGATACCTTTGCGGTGGGCGCGGTGTTGACCTGGTAAGCGGCCGGGGCCG‑CTAGGCGGCCCCGAATGACACGGTCTTTTACGCCCATGGAACCGCCCACCCGCTAc                                                                                                 >  1:1695584/1‑119 (MQ=255)                                                                                    aGAACGAAGACACCGATACCTTTGCGGTGGGCGCGGTGGTGACCTGGTAAACGGCAGGAAGAA‑AGCGGGGGCCCCGAATCATACGGTCTTTTACGCCCATGGAACCGCCCACCCGCTACCCAAGTAGCATTCGTCACCCACCACCAcctc                                                                  <  1:611089/150‑1 (MQ=255)                                                                                        cGAGGACACCGATACCTTTGCGGTGGGCGCGGTGTTGACCTGGTAAGCGGCCGGGGCCG‑CTAGGCGGCCCCGAATGACACGGTCTTTTACGCCCATGGAACCGCCCACCCGCTACCCAAGTAGCATTCGTCACCCACCACCACCTCCTTa                                                              <  1:1405647/150‑1 (MQ=255)                                                                                        agAAGACACCGATACCTTTGCGGTGGGCGCGGTGTTGACCTGGTAAGCGGCCGGGGCCG‑CTAGGCGGCCCCGAATGACACGGTCTTTTACGCCCATGGAACCGCCCACCCGCTACCCAAGTAGCATTCGTCACCCACCACCACCTCCTTa                                                              <  2:1735973/149‑1 (MQ=255)                                                                                           aGACACCGATACCTTTGCGGTGGGCGCGGTGTTGACCTGGTAAGCGGCCGGGGCCG‑CTAGGCGGCCCCGAATGACACGGTCTTTTACGCCCATGGAACCGCCCACCCGCTACCCAAGTAGCATTCGTCACCCACCACCACCTCCTTAAAc                                                           <  2:1332563/150‑1 (MQ=255)                                                                                               acCGATACCTTTGCGGTGGGCGCGGTGTTGACCTGGTAAGCGGCCGGGGCCG‑CTAGGCGGCCCCGAATGACACGGTCTTTTACGCCCATGGAACCGCCCACCCGCTACCCAAGTAGCATTCGTCACCCACCACCACCTCCTTAAACTGAt                                                       <  1:203582/150‑1 (MQ=255)                                                                                                  gATACCTTTGCGGTGGGCGCGGTGTTGACCTGGTAAGCGGCCGGGGCCG‑CTAGGCGGCCCCGAATGACACGGTCTTTTACGCCCATGGAACCGCCCACCCGCTACCCAAGTAGCATTCGTCACCCACCACCACCTCCTTAAACTGATCAc                                                    <  1:1695641/150‑1 (MQ=255)                                                                                                   aTACCTTTGCGGTGGGCGCGGTGTTGACCTGGTAAGCGGCCGGGGCCG‑CTAGGCGGCCCCGAATGACACGGTCTTTTACGCCCATGGAACCGCCCACCCGCTACCCAAGTAGCATTCGTCACCCACCACCACCTCCTTAAACTGATCACt                                                   >  2:1599021/1‑150 (MQ=255)                                                                                                   aTACCTTTGCGGTGGGCGCGGTGTTGACCTGGTAAGCGGCCGGGGCCG‑CTAGGCGGCCCCGAATGACACGGTCTTTTACGCCCATGGAACCGCCCACCCGCTACCCAAGTAGCATTCGTCACCCACCACCACCTCCTTAAACTGATCACt                                                   <  2:1787670/150‑1 (MQ=255)                                                                                                     aCCTTTTCGGTGGGCGCCGGGTTTACCTGGTCAGCGGCCGGGTCCG‑CTAGGCGGCCCCGAATGACACGGTCTTTTACGCCCATCGAACCGTCCACCTGCTAGCCAAGTAGCATTCATCACCCACCATCACCTCCTTAAACTGATCACTCa                                                 <  2:707002/150‑1 (MQ=255)                                                                                                     aCCTTTGCGGTGGGCGCGGTGTTGACCTGGTAAGCGGCCGGGGCCG‑CTAGGCGGCCCCGAATGACACGGTCTTTTACGCCCATGGAACCGCCCACCCGCTACCCAAGTAGCATTCGTCACCCACCACCACCTCCTTAAACTGATCACTCa                                                 <  1:1627435/150‑1 (MQ=255)                                                                                                     aCCTTTGCGGTGGGCGCGGTGTTGACCTGGTAAGCGGCCGGGGCCG‑CTAGGCGGCCCCGAATGACACGGTCTTTTACGCCCATGGAACCGCCCACCCGCTACCCAAGTAGCATTCGTCACCCACCACCACCTCCTTAAACTGATCACTCa                                                 <  2:772539/150‑1 (MQ=255)                                                                                                                       gtgtTGACCTGGTAAGCGGCCGGGGCCG‑CTAGGCGGCCCCGAATGACACGGTCTTTTACGCCCATGGAACCGCCCACCCGCTACCCACGTAGCATTCGTCACCCACCACCACCTCCTTAAACTGATCACTCAGCCACCCCGTGCCGgaga                               <  1:482965/150‑1 (MQ=255)                                                                                                                             aCCTGGTAAGCGGCCGGGGCCG‑CTAGGCGGCCCCGAATGACACGGTCTTTTCAGCCCAAGGAACCGCCCACCCGCTACCCAAGTAGCCTTCGTCACCCACCACCACCTCCTTAAACTGATCACTCAGCCACCCCGTGCCGGAGACCACGa                         >  2:464315/1‑150 (MQ=255)                                                                                                                             aCCTGGTAAGCGGCCGGGGCCG‑CTAGGCGGCCCCGAATGACACGGTCTTTTACGCCCATGGAACCGCCCACCCGCTACCCAAGTAGCATTCGTCACCCACCACCACCTCCTTAAACTGATCACTCAGCCACCCCGTGCCGGAGACCACGa                         <  1:771313/150‑1 (MQ=255)                                                                                                                             aCCTGGTAAGCGGACGGGGCCG‑CTAGGCGGCCCAGAATGACACGGTCTTTTACGCCCATGGAACCGCCCACCCGCTACCCAAGTAGCATTCGTCACCCACCACCACCTCCTTAAACTGATCACTCAGCCACCCCGTGCCGGAGACCACGa                         <  2:865220/150‑1 (MQ=255)                                                                                                                               cTGGTAAGCGGCCGGGGCCG‑CTAGGCGGCCCCGAATGACACGGTCTTTTACGCACATGGAACCGCCCACCCGCTACCCAAGTAGCATTCGTCACCCACCACCACCTCCTTAAACTGATCACTCAGCCACCCCGTGCCGGAGACCACGATg                       >  2:463915/1‑150 (MQ=255)                                                                                                                                tGGTAAGCGGCCGGGGCCG‑CTAGGCGGCCCCGAATGACACGGTCTTTTACGCCCATGGAACCGCCCACCCGCTACCCAAGTAGCATTCGTCACCCACCACCACCTCCTTAAACTGATCACTCAGCCACCCCGTGCCGGAGACCACGATgc                      >  1:173365/1‑150 (MQ=255)                                                                                                                                   tAAGCGGCCGGGGCCG‑CTAGGCGGCCCCGAATGACACGGTCTTTTACGCCCATGGAACCACCCACCCGCTACCCAAGTAGCATTCGTCACCCACCACCACCTCCTTAAACTGATCACTCAGCCACCCCGTGCCGga                                 <  1:1274505/136‑1 (MQ=255)                                                                                                                                   tAAGCGGCCGGGGCCG‑CTAGGCGGCCCCGAATGACACGGTCTTTTACGCCCATGGAACCACCCACCCGCTACCCAAGTAGCATTCGTCACCCACCACCACCTCCTTAAACTGATCACTCAGCCACCCCGTGCCGga                                 >  2:1274505/1‑136 (MQ=255)                                                                                                                                    taGCGGCCCCGGCCG‑CTAGGCGGCCCCGAATGACACGGTCTTTTACGCCCATGGAACCGCCCACCCGCTACCCAAGTAGCTTGCGTCACCCACCACCACCTCCTTAAACTGATCACTCAGCCACCCCGTGCCg                                   >  1:1554548/2‑133 (MQ=255)                                                                                                                                    taGCGGCCCCGGCCG‑CTAGGCGGCCCCGAATGACACGGTCTTTTACGCCCATGGAACCGCCCACCCGCTACCCAAGTAGCATTCGTCACCCACCACCACCTCCTTAAACTGATCACTCAGCCACCCCGTGCCg                                   <  2:1554548/132‑1 (MQ=255)                                                                                                                                     aGCGGCCGGGGCCG‑CTAGGCGGCCCCGAATGACACGGTCTTTTACGCCCATGGAACCGCCCACCCGCTACCCAAGTAGCATTCGTCACCCACCACCACCTCCTTAAACTGATCACTCAGCCACCCCGTGCCGGAGACCACGATGcagcag                 <  2:1691860/150‑1 (MQ=255)                                                                                                                                      gCGGCCGGGGCCG‑CTAGGCGGCCCCGAATGACAAGGTCTTTTCCGCCCATGGAAGCGCCCACCCGCTACCCCCCGACCATTCGTCACCCACCACCACCTCCTTAAACTGATCACTCAGCCACCCCGTGCCGGAGACCACGATGCAGCAgg                <  1:888761/150‑1 (MQ=255)                                                                                                                                        attcGGGGCCG‑CATAGCGGCACCAAATGACACGGTCTTTTACGCCCATGGAACCGCCCACCCGCTACCCAAGTAGCATTCGTCACCCACCACCACCTCCTTAAACTGATCACTCAGCCACCCCGTGCCGGAGACCACGATGCAGCAGGCt              >  2:168842/4‑150 (MQ=255)                                                                                                                                         ttcGGGGCCGCCTA‑GCGGCCCCGAATGACACGGTCTTTTACGCCCATGGAACCGCCCACCCGCTACCCAAGTAGCATTCGTCACCCACCACCACCTCCTTAAAc                                                           >  2:185924/3‑104 (MQ=255)                                                                                                                                         ttcGGGGCCGCCTA‑GCGGCCCCGAATGACACGGTCTTTTACGCCCATGGAACCGCCCACCCGCTACCCAAGTAGCATTCGTCACCCACCACCACCTCCTTAAAc                                                           <  1:185924/102‑1 (MQ=255)                                                                                                                                         ttcGGGGCCG‑CTAGGCGGCCCCGAATGACACGGTCTTTTACGCCCATGGAACCGCCCACCCGCTACCCAAGTAGCATTCGTCACCCACCACCACCTCCTTAAACTGATCACTCAGCCACCCCGTGCCGGAGACCACGAt                        <  2:1140040/137‑1 (MQ=255)                                                                                                                                         ttcGGGGCCG‑CTAGGCGGCCCCGAATGACACGGTCTTTTACGCCCATGGAACCGCCCACCCGCTACCCAAGTAGCATTCGTCACCCACCACCACCTCCTTAAACTGATCACTCAGCCACCCCGTGCCGGAGACCACGAt                        >  1:1140040/3‑139 (MQ=255)                                                                                                                                         gCCGGGGCCG‑CTAGGCGGCCCCGAATGACACGGTCTTTTACGCCCATGGAACCGCCCACCCGCTACCCAAGTAGCATTCGTCACCCACCACCACCTCCTTAAACTGATCACTCAGCCACCCCGTGCCGGAGACCACGATGCAGCAGGCTc             >  2:1775139/1‑150 (MQ=255)                                                                                                                                         gCCGGGGCCG‑CTAGGCGGCCCCGAATGACACGGTCTTTTACGCCCATGGAACCGCCCACCCGCTACCCAAGTAGCATTCGTCACCCACCACCACCTCCTTAAACTGATCACTCAGCCACCCCGTGCCGGAGACCACGATGCAGCAGGCTc             >  2:126188/1‑150 (MQ=255)                                                                                                                                          ccGGGGCCG‑CTAGGCGGCCCCGAATGACACGGTCTTTTACGCCCATGGAACCGCCCACCCGCTACCCAAGTAGCATTCGTca                                                                                >  2:1769217/1‑82 (MQ=255)                                                                                                                                          ccGGGGCCG‑CTAGGCGGCCCCGAATGACACGGTCTTTTACGCCCATGGAACCGCCCACCCGCTACCCAAGTAGCATTCGTca                                                                                <  1:1769217/82‑1 (MQ=255)                                                                                                                                            ggggCCGCCTA‑GCGGCCCCGAATGACACGGTCTTTTACGCCCATGGAACCGCCCACCCGCTACCCAAGTAGCATTCGTCACCCACCACCACCTCCTTAAACTGATCACTCAGCCAcc                                           >  1:1622906/1‑117 (MQ=255)                                                                                                                                            ggggCCGCCTA‑GCGGCCCCGAATGACACGGTCTTTTACGCCCATGGAACCGCCCACCCGCTACCCAAGTAGCATTCGTCACCCACCACCACCTCCTTAAACTGATCACTCAGCCAcc                                           <  2:1622906/117‑1 (MQ=255)                                                                                                                                             gggCCG‑CTAGGCGGCCCCGAATGACACGGTCTTTTACGCCCATGGAACCGCCCACCCGCTCCCCAAGTAGCATTCGTCACCCACCACCACCTCCTTAAACTGATCACTCAGCCACCCCGTGCCGGAGACCACGATGCAGCAGGCTCATaa         <  2:639357/150‑1 (MQ=255)                                                                                                                                                ccG‑CTAGGCGGCCCCGAATGACACGGTCTTTTACGCCCATGGAACCGCCCACCCGCTACCCAAGTAGCATTCGTCACCCACCACCACCTCCTTAAACTGATCACTCAGCCACCCCGTGCCGGAGACCACGATGCAGCAGGCTCATAacga      <  2:710078/150‑1 (MQ=255)                                                                                                                                                     tAGGCGGCCCCGAATGACACGGTCTTTTACGCCCATGGAACCGCCCACCCGCTACCCAAGTAGCATTCGTCACCCACCACCACCTCCTTAAACTGATCACTCAGCCACCCCGTGCCGGAGACCACGATGCAGCAGGCTCATAACGACGgt  <  2:173365/150‑1 (MQ=255)                                                                                                                                                     tAGGCGGCCCCGAATGACACGGTCTTTTACGCCCATGGAACCGCCCACCCGCTACCCAAGTAGCATTCGTCACCCACCACCACCTCCTTAAACTGATCACTCAGCCACCCCGTGCCGGAGACCACGATGCAGCAGGCTCATAACGACGgt  <  2:953098/150‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                    GCTGTTCCATGACATCAACGACAACCTGAAACTGGTGGCCGAGTACAACCAGTTTTCCATCGACGGGCATGACACCAGTGCGCAGAACGAAGACACCGATACCTTTGCGGTGGGCGCGGTGTTGACCTGGTAAGCGGCCGGGGCCG‑CTAGGCGGCCCCGAATGACACGGTCTTTTACGCCCATGGAACCGCCCACCCGCTACCCAAGTAGCATTCGTCACCCACCACCACCTCCTTAAACTGATCACTCAGCCACCCCGTGCCGGAGACCACGATGCAGCAGGCTCATAACGACGGT  >  NC_002947/3049939‑3050235

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