Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,229,191	Δ1 bp	10.0%	intergenic (‑133/+43)	PP_4606 ← / ← PP_4607	ferric siderophore receptor/transmembrane sensor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,229,191	0	G	.	10.0%	183.3 / 14.8	60	intergenic (‑133/+43)	PP_4606/PP_4607	ferric siderophore receptor/transmembrane sensor
Reads supporting (aligned to +/- strand):  ref base G (28/26);  new base . (4/2);  total (32/28)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.35e-01

GCGGGGTGAACTTCATGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAAT‑CCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGCTGC  >  NC_002947/5229043‑5229337                                                                                                                                                     |                                                                                                                                                     gcgGGGTGAACTTCATGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCg                                                                                                                                                      >  5:276285/1‑149 (MQ=255)    gggTGAACTTCATGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGAGGCCTCTTCGCGGGTAACCCCGccc                                                                                                                                                   >  8:217975/1‑149 (MQ=255)      gTGAACTTCATGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGCCAGCCTCTTCACGGGTAAAccccccgta                                                                                                                                                 <  1:186837/149‑4 (MQ=255)             tCATGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACggg‑gg                                                                                                                                         >  3:300295/1‑149 (MQ=255)              cATGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACggg‑gg                                                                                                                                         <  6:276285/148‑1 (MQ=255)                  gACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGttt                                                                                                                                     <  4:276451/148‑1 (MQ=255)                  gACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTa                                                                                                                                    >  3:290307/1‑149 (MQ=255)                  gACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCAGCCTCTTCACGGGTAAACCC‑CCCGTAGGGGCGGGTTTa                                                                                                                                    >  8:171838/1‑149 (MQ=38)                     gCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTTCACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGccc                                                                                                                                                   <  5:235185/132‑1 (MQ=255)                     gCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTTCACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGccc                                                                                                                                                   >  6:235185/1‑132 (MQ=255)                     gCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCC‑CCCGTAGGGGCGGGTTTAccc                                                                                                                                 <  8:31724/149‑1 (MQ=39)                      cAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCg                                                                                                                                <  8:123565/149‑1 (MQ=255)                         aCCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCAGCCTCTTCACGGGTAAACCC‑CCCGTAGGGGCGGGTTTACCCGCGa                                                                                                                             <  7:171838/149‑1 (MQ=37)                               gAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAgg                                                                                                                        <  4:211486/148‑1 (MQ=255)                                 gggggTGATTATTGCGAATCAATCGCACTTGCACTCATTACACGGATGCCACGTGCTGTTACCTCACGTGTGCCTTGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTg                                                                                                                     <  5:167208/149‑1 (MQ=255)                                  ggggTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTCCCCGTGAAGAGGCTgg                                                                                                                    >  4:222240/1‑149 (MQ=255)                                   gattGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGc                                                                                                                   <  6:175426/146‑1 (MQ=255)                                        attattGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAgg                                                                                                              >  5:203484/1‑149 (MQ=255)                                            ttGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGATGTGGCAGCCTCTTCACGGGTAAACCC‑CCCCTAGGGGCGGGTTTACACTCGAAGAGGCCGCCACAgcaaa                                                                                                           >  1:267498/1‑144 (MQ=21)                                               cGAATCAATCGCACTTTCACTCATTCCCCGGATGCCACGTGCTGTTACCTCACCTGTCCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTa                                                                                                        <  7:217975/148‑1 (MQ=255)                                                          cACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGt                                                                                            <  5:329059/149‑1 (MQ=255)                                                             tttCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCg                                                                                                     >  4:204605/1‑137 (MQ=255)                                                             tttCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCg                                                                                                     <  3:204605/137‑1 (MQ=255)                                                             tttCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCa                                                                                          >  4:203282/1‑148 (MQ=255)                                                                 aCTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCaaa                                                                                        >  7:144613/1‑146 (MQ=255)                                                                 aCTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCaaa                                                                                        <  8:144613/146‑1 (MQ=255)                                                                 aCTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGc                                                                                     <  6:203484/149‑1 (MQ=255)                                                                   tCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGt                                                                                                                                       >  1:149052/1‑97 (MQ=255)                                                                   tCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGt                                                                                                                                       <  2:149052/97‑1 (MQ=255)                                                                     aTTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCa                                                                                 <  5:35593/149‑1 (MQ=255)                                                                        acacGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGc                                                                                     >  7:112526/1‑142 (MQ=255)                                                                        acacGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGc                                                                                     <  8:112526/142‑1 (MQ=255)                                                                          acGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACg                                                                            >  8:215265/1‑149 (MQ=255)                                                                                  cACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCg                                                                    >  1:123741/1‑149 (MQ=255)                                                                                    cGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTg                                                                  <  3:296101/149‑1 (MQ=255)                                                                                     gTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGCCTGCCTCTTCGCGGGTAAACCCGCCCCTACGGT‑GGGTTTACCCGTGAAGAGGCTGGCACATGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGa                                                                 <  2:229228/149‑1 (MQ=255)                                                                                        cTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGaaa                                                               <  7:215265/148‑1 (MQ=255)                                                                                         tGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGt                                                                   >  1:7616/1‑143 (MQ=255)                                                                                         tGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGt                                                                   <  2:7616/143‑1 (MQ=255)                                                                                                 cACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAgg                                                                                                              <  4:17037/92‑1 (MQ=255)                                                                                                 cACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAgg                                                                                                              >  3:17037/1‑92 (MQ=255)                                                                                                    cTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCtt                                                  >  6:310036/1‑149 (MQ=255)                                                                                                          cTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAg                                            >  6:17749/1‑149 (MQ=255)                                                                                                          cTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAg                                            >  2:40208/1‑149 (MQ=255)                                                                                                           tGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGc                                           >  8:66680/1‑149 (MQ=255)                                                                                                                atatTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAAT‑c                                     >  6:224545/1‑149 (MQ=255)                                                                                                                             ggcggcCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAAT‑CCAGCGCCagcagc                        >  1:7574/1‑149 (MQ=255)                                                                                                                              gcggcCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAAT‑CCAGCGCCAGCAGCg                       >  4:55521/1‑149 (MQ=255)                                                                                                                              ccagcCTCTTCACGGGTAAACCC‑CCCGTAGGGGCGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAAT‑CCAGCGCCAGCAGCg                       >  1:87120/4‑149 (MQ=37)                                                                                                                                   ctctTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAAT‑CCAGCGCCAGCAGCGGCTCg                  <  3:55521/149‑1 (MQ=255)                                                                                                                                   ctctTCACGGGTAAACCC‑CCCGTAGGGGCGGGTTTACCCGCGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAAT‑CCAGCGCCAGCAGCGGCTCg                  >  5:352410/1‑149 (MQ=37)                                                                                                                                      ttCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAAT‑CCAGCGCCAGCAGCGGCTCGGc                >  3:298147/1‑148 (MQ=255)                                                                                                                                         gcgGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAAT‑CCAGCGCCAGCAGCGGCTCGGCCAcc            >  3:57289/1‑149 (MQ=255)                                                                                                                                           gggTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGcc                                                                     <  6:100243/91‑1 (MQ=255)                                                                                                                                           gggTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGcc                                                                     >  5:100243/1‑91 (MQ=255)                                                                                                                                           gggTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAAT‑CCAGCGCCAGCAGCGGCTCGGCCACCgg          >  8:92676/1‑149 (MQ=255)                                                                                                                                            ggTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAAT‑CCAGCGCCAGCAGCGGCTCGGCCACCgg          >  6:370751/1‑148 (MQ=255)                                                                                                                                             gTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAAT‑CCAGCGCCAGCAGCGGCTCGGCCACCg           <  6:116458/146‑1 (MQ=255)                                                                                                                                             gTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAAT‑CCAGCGCCAGCAGCGGCTCGGCCACCg           >  5:116458/1‑146 (MQ=255)                                                                                                                                                aaCCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATACCAGCGCCAGCAGCGGCTCGGCCACCGGGaa       <  1:40208/148‑1 (MQ=255)                                                                                                                                                 aCCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAAT‑CCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGCt    <  5:17749/149‑1 (MQ=255)                                                                                                                                                   ccGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAAT‑CCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGCTGc  <  4:290307/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     GCGGGGTGAACTTCATGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAAT‑CCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGCTGC  >  NC_002947/5229043‑5229337

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