Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,577,194	Δ1 bp	7.4%	coding (1244/1776 nt)	PP_2259 ←	sigma‑54 dependent transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,577,194	0	A	.	7.4%	216.0 / 10.3	68	coding (1244/1776 nt)	PP_2259	sigma‑54 dependent transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (27/36);  new base . (3/2);  total (30/38)
Fisher's exact test for biased strand distribution p-value = 6.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

GCCGGCGATGCGGTAGTACAGGTCCTGGCGAAAGCGCCCGGCGCGCACTTCGTCGGCCAGGTTGCGATTGCTCGCCGACACCACGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCG‑ACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGCCGTGGTCGGCCTTGCGGATCAGTCCCGGGTTGCCCTTGTGGTGGGCGCCGGTGAACGCGCCTTTGTCATAGCCGAACAGTTCCGA  >  NC_002947/2577049‑2577342                                                                                                                                                   |                                                                                                                                                     gCCGGCGATGCGGTAGTACAGGTCCTGGCGAAAGCGCCCGGCGCGCACTTCGTCGGCCAGGTTGCGATTGCTCGCCGACACCACGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCG‑ACTGCAGTAcgcg                                                                                                                                                  >  2:753940/1‑150 (MQ=255)   cGGCGATGCGGTAGTACAGGTCCTGGCGAAAGCGCCCGGCGCGCACTTCGTCGGCCAGGTTGCGATTGCTCGCCGACACCACGCGGATATCCACCGGCACCGGTTCGCCACTGGCAAGCGGCTGGATGCTGCGCG‑ACTGCAGTACGCgca                                                                                                                                                >  1:728413/1‑150 (MQ=255)   cGGCGATGCGGTAGTACAGGTCCTGGCGAAAGCGCCCGGCGCGCACTTCGTCGGCCAGGTTGCGATTGCTCGCCGACACCACGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCG‑ACTGCAGTACGCgca                                                                                                                                                >  1:1017507/1‑150 (MQ=255)            ggTAGTACAGGTCCTGGCGAAAGCGCCCGGCGCGCACTTCGTCGGCCAGGTTGCGATTGCTCGCCGACACCACGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCGTACTGCAGT‑CGCgcagcatccagc                                                                                                                                       <  1:664126/150‑7 (MQ=255)                 taCAGGTCCTGGCGAAAGCGCCCGGCGCGCACTTCGTCGGCCAGGTTGCGATTGCTCGCCGACACCACGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCG‑ACTGCAGTACGCGCAGCAGGCGGGCCtgg                                                                                                                                  <  2:697600/150‑1 (MQ=255)                          tGGCGAAAGCGCCCGGCGCGCACTTCGTCGGCCAGGTTGCGATTGCTCGCCGACACCACGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCG‑ACTGCAGTACGCGCAGCAGGc                                                                                                                                          <  2:2619/133‑1 (MQ=255)                          tGGCGAAAGCGCCCGGCGCGCACTTCGTCGGCCAGGTTGCGATTGCTCGCCGACACCACGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCG‑ACTGCAGTACGCGCAGCAGGc                                                                                                                                          >  1:2619/1‑133 (MQ=255)                          gggCGAAAGCGCCCGGCGCGCACTTCGTCGGCCAGGTTGCGATTGCTCGCCGACACCACGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCG‑ACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGt                                                                                                                         <  2:426467/149‑1 (MQ=255)                           ggCGAAAGCGACCGGCGCGCACTTCGTCGGCCAGGTTGCGATTGCTCGCCGACACCACGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCG‑ACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTg                                                                                                                        >  2:951883/1‑150 (MQ=255)                            gCGAAAGCGCCCGGCGCGCACTTCGTCGGCCAGGTTGCGATTGCTCGCCGACACCACGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCG‑ACTGCAGTACGCGCAGCAGGCgg                                                                                                                                        >  1:1171005/1‑133 (MQ=255)                            gCGAAAGCGCCCGGCGCGCACTTCGTCGGCCAGGTTGCGATTGCTCGCCGACACCACGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCG‑ACTGCAGTACGCGCAGCAGGCgg                                                                                                                                        <  2:1171005/133‑1 (MQ=255)                               aaaGCGCCCGGCGCGCACTTCGTCGGCCAGGTTGCGATTGCTCGCCGACACCACGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCG‑ACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCAt                                                                                                                    <  1:753940/150‑1 (MQ=255)                                  gcgcCCGGCGCGCACTTCGTCGGCCAGGTTGCGATTGCTCGCCGACACCACGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCG‑ACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTc                                                                                                                 <  2:1026242/150‑1 (MQ=255)                                    gcCCGGCGCGCACTTCGTCGGCCAGGTTGCGATTGCTCGCCGACACCACGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCG‑ACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGc                                                                                                               >  1:1066050/1‑150 (MQ=255)                                         gcgcgcACTTCGTCGGCCAGGTTGCGATTGCTCGCCGACACCACGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCG‑ACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAAtc                                                                                                          <  2:359797/150‑1 (MQ=255)                                              cACTTCGTCGGCCAGGTTGCGATTGCTCGCCGACACCACGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCG‑ACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATc                                                                                                     <  2:890287/150‑1 (MQ=255)                                                 ttCGTCGGCCAGGTTGCGATTGCTCGCCGACACCACGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCG‑ACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAg                                                                                                  >  1:558793/1‑150 (MQ=255)                                                      cgGCCAGGTTGCGATTGCTGGCCGACACTGCGCGGATATCCACCGGCACCGGTTCGCCACTGCACAGAAACATGATGCTGCGCG‑ACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACa                                                                                             <  1:1169110/150‑1 (MQ=255)                                                      cgGCCAGGTTGCGATTGCTCGCCGACACCACGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCG‑ACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACa                                                                                             <  2:1122788/150‑1 (MQ=255)                                                       ggCCAGGTTGCGATTGCTCGCCGAAACCACGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCG‑ACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAg                                                                                            <  2:761139/150‑1 (MQ=255)                                                             gTTGCGATTGCTCGCCGACACCACGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCG‑ACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTg                                                                                                                        <  1:1147008/116‑1 (MQ=255)                                                             gTTGCGATTGCTCGCCGACACCACGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCG‑ACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTg                                                                                                                        >  2:1147008/1‑116 (MQ=255)                                                             gTTGCGATTGCTCGCCGACACCACGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCG‑ACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGcc                                                                                      >  1:187581/1‑150 (MQ=255)                                                              gtGCGATGGCTCGCCGACACCACGCGGATATCCCCCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCG‑ACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTAATCCAGAAACAGGATGCCg                                                                                     <  2:108160/149‑1 (MQ=255)                                                                    ttGCTCGCCGACACCACGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCG‑ACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGc                                                                                       >  2:783274/1‑142 (MQ=255)                                                                    ttGCTCGCCGACACCACGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCG‑ACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGc                                                                                       <  1:783274/142‑1 (MQ=255)                                                                        tCGCCGACACCACGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCG‑ACTGCAGTACGCGcagcag                                                                                                                                            <  1:101440/85‑1 (MQ=255)                                                                        tCGCCGACACCACGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCG‑ACTGCAGTACGCGcagcag                                                                                                                                            >  2:101440/1‑85 (MQ=255)                                                                        tCGCCGACACCACGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCG‑ACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATc                                                                                                     >  2:68920/1‑124 (MQ=255)                                                                        tCGCCGACACCACGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCG‑ACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATc                                                                                                     <  1:68920/124‑1 (MQ=255)                                                                          gCCGACACCACGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCG‑ACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGCCGTGGTCGGCCTTg                                                                         <  2:1017507/150‑1 (MQ=255)                                                                           ccGACACCACGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCG‑ACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGCCGTGGTCGGCCTTGc                                                                        <  2:1066050/150‑1 (MQ=255)                                                                             gACACCACGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCG‑ACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGCCGTGGTCGGCCTTGCgg                                                                      >  1:842955/1‑150 (MQ=255)                                                                              acacCACGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCG‑ACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGCCGTGGTCGGCCTTGCGGa                                                                     >  2:851850/1‑150 (MQ=255)                                                                                accacGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCG‑ACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGCCGTGGTCGGCCTTGCGGATc                                                                   >  1:711492/1‑150 (MQ=255)                                                                                accacGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCG‑ACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGCCGTGGTCGGCCTTGCGGATc                                                                   >  2:107636/1‑150 (MQ=255)                                                                                 ccacGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCG‑ACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGCCGTGGTCGGCCTTGCGGATCa                                                                  <  2:711492/150‑1 (MQ=255)                                                                                   acGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCG‑ACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCGTGTCGCCAATCTCATCCAGAAACAGGATGCCGTGGTCGGCCTTGCGGATCAGt                                                                >  1:1150965/1‑150 (MQ=255)                                                                                     gcgGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCG‑ACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGCCGTGGTCGGCCTTGCGGATCAGTcc                                                              <  2:187581/150‑1 (MQ=255)                                                                                        gATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCG‑ACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGCCGTGGTCGGCCTCGCGGATCAGTCCCgg                                                           <  1:804337/150‑1 (MQ=255)                                                                                            tCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCG‑ACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGCCGTGGTCGGCCTTGCGGATCAGTCCCGGGTTg                                                       <  1:107636/150‑1 (MQ=255)                                                                                             ccACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCG‑ACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATc                                                                                                     <  2:1097748/103‑1 (MQ=255)                                                                                             ccACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCG‑ACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATc                                                                                                     >  1:1097748/1‑103 (MQ=255)                                                                                              gaccgtgaccggATCCCCACTGCCCAGCGGCTGGATGCTGCGCG‑ACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGCCGTGGTCGGCCTTGCGGATCAGTCCCGGGTTGcc                                                     <  2:558793/143‑1 (MQ=255)                                                                                              caccggcaccggTTCGCCACTGCCCAGCGGCTGGATGCTGCGCG‑ACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGCCGTGGTCGGCCTTg                                                                         <  1:52320/130‑1 (MQ=255)                                                                                              caccggcaccggTTCGCCACTGCCCAGCGGCTGGATGCTGCGCG‑ACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGCCGTGGTCGGCCTTg                                                                         >  2:52320/1‑130 (MQ=255)                                                                                               accggcaccggTTCGCCACTGCCCAGCGGCTGGATGCTGCGCG‑ACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGCCGTGGTCGGCCTTGCGGATCAGTCCCGGGTTGccc                                                    <  1:851850/150‑1 (MQ=255)                                                                                               accggcaccggTTCGCCACTGCCCAGCGGCTGGATGCTGCGCG‑ACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGCCGTGGTCGGCCTTGCGGATCAGTCCCGGGTTGccc                                                    >  1:1198846/1‑150 (MQ=255)                                                                                                   gcaccggTTCGCCACTGCCCAGCGGCTGGATGCTGCGCG‑ACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGCCGTGGTCGGCCTTGCGGATCAGTCCCGGGTTGCCCTtgt                                                <  1:1172651/150‑1 (MQ=255)                                                                                                    caccggTTCGCCACTGCCCAGCGGCTGGATGCTGCGCG‑ACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGCCGTGGTCGGCCTTGCGGATCAGTCCCGGGTTGCCCTtgtg                                               <  1:951883/150‑1 (MQ=255)                                                                                                       cggTTCGCCACTGCCCAGCGGCTGGATGCTGCGCG‑ACTGCAGTACGCGCAGCAGGCGGGCCTTGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGCCGTGGTCGGCCTTGCGGATCAGTCCCGGGTTGCCCTTgtggtg                                            <  2:461693/150‑1 (MQ=255)                                                                                                              ccACTGCCCAGCGGCTGGATGCTGCGCG‑ACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCaa                                                                                                            >  2:84759/1‑79 (MQ=255)                                                                                                              ccACTGCCCAGCGGCTGGATGCTGCGCG‑ACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCaa                                                                                                            <  1:84759/79‑1 (MQ=255)                                                                                                                    cccAGCGGCTGGATGCTGCGCG‑ACTGCAGTACGCGCAGCAGGCGGCCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGCCGTGGTCGGCCTTGCGGATCAGTCCCGGGTTGCCCTTGTGGTGGGCGCCGGTGAAc                               <  1:131069/150‑1 (MQ=255)                                                                                                                      cAGCGGCTGGATGCTGCGCG‑ACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGCCGTGGTCGGCCTTGCGGATCAGTCCCGGGTTGCCCTTGTGGTGGGCGCCGGTGAAcgc                             >  1:322794/1‑150 (MQ=255)                                                                                                                        gCGGCTGGATGCTGCGCG‑ACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGCCGTGGTCGGCCTTGCGGATCAGTCCCGGGTTGCCCTTGTGGTGGGCGCCGGTGAACgcgc                           >  1:1076468/1‑150 (MQ=255)                                                                                                                              ggATGCTGCGCG‑ACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGCCGTGGTCGGCCTTGCGGATCAGTCCCGGGTTGCCCTTGTGGTGGGCGCCGGTGAACGCGCCTTTGt                     <  2:752364/150‑1 (MQ=255)                                                                                                                                  gctgcGCG‑AATGAAGTAAGTCCAGAAAGCGGACATGGGTGGGAAGTGGCATGTCGCAAATCTCATCCAGAAATAGTATGCTGTGGTCGGCCTGTCGGATCAGTCCCGGTTTGCCCTTGTGGTGGGCGCCGGTGAACGCGCCTTTGTcctt                 <  1:777078/150‑4 (MQ=255)                                                                                                                                    tgcGCGTACTGCAGT‑CGCGCAGCATGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGCCGTGGTCGGCCTTGCGGATCAGTCCCGGGTTGCCCTTGTGGTGGGCGCCGGTGAACgcgc                           >  2:597108/1‑138 (MQ=255)                                                                                                                                    tgcGCGTACTGCAGT‑CGCGCAGCATGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGCCGTGGTCGGCCTTGCGGATCAGTCCCGGGTTGCCCTTGTGGTGGGCGCCGGTGAACgcgc                           <  1:597108/138‑1 (MQ=255)                                                                                                                                    tgcGCGTACTGCAGT‑CGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGCCGTGGTCGGCCTTGCGGATCAGTCCCGGGTTGCCCTTGTGGTGGGCGCCGGTGAACGCGCCTTTGTCATAGc               >  1:1143903/1‑150 (MQ=255)                                                                                                                                     gcgcgtaCTGCAGT‑CGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGCCGTGGTCGGCCTTGCGGATCAGTCCCGGGTTGCCCTTGTGGTGGGCGCCGGTGAACGCGCCTTTGTCATAGcc              >  1:817147/7‑150 (MQ=255)                                                                                                                                           aCTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGCCGTGGTCGGCCTTGCGGATCAGTCCCGGGTTGCCCTTGTGGTGGGCGCCGGTGAACGCGCCTTTGTCATAGCCGAACa         <  2:1076468/150‑1 (MQ=255)                                                                                                                                             tGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGCCGTGGTCGGCCTTGCGGATCAGTCCCGGGTTGCCCTTgtggtg                                            >  2:130295/1‑113 (MQ=255)                                                                                                                                             tGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGCCGTGGTCGGCCTTGCGGATCAGTCCCGGGTTGCCCTTgtggtg                                            <  1:130295/113‑1 (MQ=255)                                                                                                                                                 gTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGCCGTGGTCGGCCTTGCGGATCAGTCCCGGGTTGCCCTtg                                                 >  2:482069/1‑104 (MQ=255)                                                                                                                                                 gTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGCCGTGGTCGGCCTTGCGGATCAGTCCCGGGTTGCCCTtg                                                 <  1:482069/104‑1 (MQ=255)                                                                                                                                                  tACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGCCGTGGTCGGCCTTGCGGATCAGTCCCGGGTTGCCCTTGTGGTGGGCGCCGGTGAACGCGCCTTTGTCATAGCCGAACAGTTCCGa  >  1:1072564/1‑150 (MQ=255)                                                                                                                                                   |                                                                                                                                                     GCCGGCGATGCGGTAGTACAGGTCCTGGCGAAAGCGCCCGGCGCGCACTTCGTCGGCCAGGTTGCGATTGCTCGCCGACACCACGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCG‑ACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGCCGTGGTCGGCCTTGCGGATCAGTCCCGGGTTGCCCTTGTGGTGGGCGCCGGTGAACGCGCCTTTGTCATAGCCGAACAGTTCCGA  >  NC_002947/2577049‑2577342

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