Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,229,170	G→A	8.7%	intergenic (‑112/+64)	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,170	0	G	A	8.7%	157.8 / 4.6	58	intergenic (‑112/+64)	PP_4606/PP_4607	ferric siderophore receptor/transmembrane sensor
Reads supporting (aligned to +/- strand):  ref base G (29/23);  new base A (3/2);  total (32/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.68e-01

AACCAGAGAGGGACGCAACGCGGGGTGAACTTCATGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGT‑GCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCG  >  NC_002947/5229024‑5229316                                                                                                                                                    |                                                                                                                                                    aaCCAGAGAGGGACGCAACGCGGGGTGAACTTCATGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGT‑GCCTGGAAATATTTTTGCTGTggcggc                                                                                                                                                   <  8:337970/149‑1 (MQ=255)  aCCAGAGAGGGACGCAACGCGGGGTGAACTTCATGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTCGCCTGGAAATATTTTTGCTGTggcggc                                                                                                                                                   <  1:16953/149‑1 (MQ=255)  aCCAGAGAGGGACGCAACGCGGGGTGAACTTCATGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGT‑GCCTGGAAATATTTTTGCTGTGCCAGcc                                                                                                                                                  >  2:186837/1‑149 (MQ=255)    cAGAGAGGGACGCAACGCGGGGTGAACTTCATGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGT‑GCCTGGAAATATTTTTGCTGTGGCGGCctc                                                                                                                                                >  2:333660/1‑149 (MQ=255)     agagagGGACGCAACGCGGGGTGAACTTCATGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGT‑GCCTGGAAATATTTTTGCTGAGGCGGCctct                                                                                                                                               >  6:329059/1‑149 (MQ=255)           ggACGCAACGCGTGTTGAACTTCATGTACGAAGACCTGTGAGGGTGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGT‑GCCTTGAAATATTTTTGCTGTGGAGGCCTCTTCGCgg                                                                                                                                         <  8:353955/149‑1 (MQ=255)           ggACGCAACGCGGGGTGAACTTCATGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGT‑GCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCgg                                                                                                                                         >  5:348071/1‑149 (MQ=255)                 aaCGCGGGGTGAACTTCATGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGT‑GCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAAc                                                                                                                                   >  4:296101/1‑149 (MQ=255)                 aaCGCGGGGTGAACTTCATGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGT‑GCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAAc                                                                                                                                   >  3:211486/1‑149 (MQ=255)                 aaCGCGGGGTGAACTTCATGGACGCAGACCTGTCAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGT‑GCCTGGAAATATTTTTGATGTGGCGGCCTCTTCGCGGGTAAAc                                                                                                                                   >  6:167208/1‑149 (MQ=255)                  aCGCGGGGTGAACTTCATGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGT‑GCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAAcc                                                                                                                                  >  1:229228/1‑149 (MQ=255)                    gcgGGGTGAACTTCATGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGT‑GCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCg                                                                                                                                >  5:276285/1‑149 (MQ=255)                      ggggTGAACTTCATGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGT‑GCCTGGAAATATTTTTGCTGTGGCGGCCTCTTcgcg                                                                                                                                          >  7:140392/1‑137 (MQ=255)                      ggggTGAACTTCATGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGT‑GCCTGGAAATATTTTTGCTGTGGCGGCCTCTTcgcg                                                                                                                                          <  8:140392/137‑1 (MQ=255)                       gggTGAACTTCATGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGT‑GCCTGGAAATATTTTTGCTGTGGAGGCCTCTTCGCGGGTAACCCCGccc                                                                                                                             >  8:217975/1‑149 (MQ=255)                         gTGAACTTCATGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGT‑GCCTGGAAATATTTTTGCTGTGCCAGCCTCTTCACGGGTAAAccccccgta                                                                                                                           <  1:186837/149‑4 (MQ=255)                                tCATGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGT‑GCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACggg‑gg                                                                                                                   >  3:300295/1‑149 (MQ=255)                                 cATGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGT‑GCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACggg‑gg                                                                                                                   <  6:276285/148‑1 (MQ=255)                                     gACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGT‑GCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGttt                                                                                                               <  4:276451/148‑1 (MQ=255)                                     gACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGT‑GCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTa                                                                                                              >  3:290307/1‑149 (MQ=255)                                     gACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGT‑GCCTGGAAATATTTTTGCTGTGGCAGCCTCTTCACGGGTAAACCC‑CCCGTAGGGGCGGGTTTa                                                                                                              >  8:171838/1‑149 (MQ=38)                                        gCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTTCACGGATGCCACGTGCTGTTACCTCACCTGT‑GCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGccc                                                                                                                             >  6:235185/1‑132 (MQ=255)                                        gCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTTCACGGATGCCACGTGCTGTTACCTCACCTGT‑GCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGccc                                                                                                                             <  5:235185/132‑1 (MQ=255)                                        gCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGT‑GCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCC‑CCCGTAGGGGCGGGTTTAccc                                                                                                           <  8:31724/149‑1 (MQ=39)                                         cAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGT‑GCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCg                                                                                                          <  8:123565/149‑1 (MQ=255)                                            aCCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGT‑GCCTGGAAATATTTTTGCTGTGGCAGCCTCTTCACGGGTAAACCC‑CCCGTAGGGGCGGGTTTACCCGCGa                                                                                                       <  7:171838/149‑1 (MQ=37)                                                  gAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGT‑GCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAgg                                                                                                  <  4:211486/148‑1 (MQ=255)                                                    gggggTGATTATTGCGAATCAATCGCACTTGCACTCATTACACGGATGCCACGTGCTGTTACCTCACGTGT‑GCCTTGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTg                                                                                               <  5:167208/149‑1 (MQ=255)                                                     ggggTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGT‑GCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTCCCCGTGAAGAGGCTgg                                                                                              >  4:222240/1‑149 (MQ=255)                                                      gggTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGT‑GCCTGGAAATATTTTTGCTGTggcggc                                                                                                                                                   >  7:207428/1‑96 (MQ=255)                                                      gggTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGT‑GCCTGGAAATATTTTTGCTGTggcggc                                                                                                                                                   <  8:207428/96‑1 (MQ=255)                                                      gattGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGT‑GCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGc                                                                                             <  6:175426/146‑1 (MQ=255)                                                        gTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGT‑GCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAAccc                                                                                                                                 <  2:130464/112‑1 (MQ=255)                                                           attattGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGT‑GCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAgg                                                                                        >  5:203484/1‑149 (MQ=255)                                                               ttGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGT‑GCCTGGAAATATTTTTGATGTGGCAGCCTCTTCACGGGTAAACCC‑CCCCTAGGGGCGGGTTTACACTCGAAGAGGCCGCCACAgcaaa                                                                                     >  1:267498/1‑144 (MQ=21)                                                                  cGAATCAATCGCACTTTCACTCATTCCCCGGATGCCACGTGCTGTTACCTCACCTGT‑CCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTa                                                                                  <  7:217975/148‑1 (MQ=255)                                                                             cACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGT‑GCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGt                                                                      <  5:329059/149‑1 (MQ=255)                                                                                tttCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGT‑GCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCg                                                                               >  4:204605/1‑137 (MQ=255)                                                                                tttCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGT‑GCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCg                                                                               <  3:204605/137‑1 (MQ=255)                                                                                tttCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGT‑GCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCa                                                                    >  4:203282/1‑148 (MQ=255)                                                                                    aCTCATTACACGGATGCCACGTGCTGTTACCTCACCTGT‑GCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCaaa                                                                  >  7:144613/1‑146 (MQ=255)                                                                                    aCTCATTACACGGATGCCACGTGCTGTTACCTCACCTGT‑GCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCaaa                                                                  <  8:144613/146‑1 (MQ=255)                                                                                    aCTCATTACACGGATGCCACGTGCTGTTACCTCACCTGT‑GCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGc                                                               <  6:203484/149‑1 (MQ=255)                                                                                      tCATTACACGGATGCCACGTGCTGTTACCTCACCTGT‑GCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGt                                                                                                                 >  1:149052/1‑97 (MQ=255)                                                                                      tCATTACACGGATGCCACGTGCTGTTACCTCACCTGT‑GCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGt                                                                                                                 <  2:149052/97‑1 (MQ=255)                                                                                        aTTACACGGATGCCACGTGCTGTTACCTCACCTGT‑GCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCa                                                           <  5:35593/149‑1 (MQ=255)                                                                                           acacGGATGCCACGTGCTGTTACCTCACCTGT‑GCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGc                                                               >  7:112526/1‑142 (MQ=255)                                                                                           acacGGATGCCACGTGCTGTTACCTCACCTGT‑GCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGc                                                               <  8:112526/142‑1 (MQ=255)                                                                                             acGGATGCCACGTGCTGTTACCTCACCTGT‑GCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACg                                                      >  8:215265/1‑149 (MQ=255)                                                                                                     cACGTGCTGTTACCTCACCTGT‑GCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCg                                              >  1:123741/1‑149 (MQ=255)                                                                                                       cGTGCTGTTACCTCACCTGT‑GCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTg                                            <  3:296101/149‑1 (MQ=255)                                                                                                        gTGCTGTTACCTCACCTGT‑GCCTGGAAATATTTTTGCTGTGCCTGCCTCTTCGCGGGTAAACCCGCCCCTACGGT‑GGGTTTACCCGTGAAGAGGCTGGCACATGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGa                                           <  2:229228/149‑1 (MQ=255)                                                                                                           cTGTTACCTCACCTGT‑GCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGaaa                                         <  7:215265/148‑1 (MQ=255)                                                                                                            tGTTACCTCACCTGT‑GCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGt                                             >  1:7616/1‑143 (MQ=255)                                                                                                            tGTTACCTCACCTGT‑GCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGt                                             <  2:7616/143‑1 (MQ=255)                                                                                                                    cACCTGT‑GCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAgg                                                                                        >  3:17037/1‑92 (MQ=255)                                                                                                                    cACCTGT‑GCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAgg                                                                                        <  4:17037/92‑1 (MQ=255)                                                                                                                       cTGT‑GCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCtt                            >  6:310036/1‑149 (MQ=255)                                                                                                                              cTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAg                      >  2:40208/1‑149 (MQ=255)                                                                                                                              cTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAg                      >  6:17749/1‑149 (MQ=255)                                                                                                                               tGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGc                     >  8:66680/1‑149 (MQ=255)                                                                                                                                    atatTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATc                >  6:224545/1‑149 (MQ=255)                                                                                                                                                 ggcggcCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCagcagc   >  1:7574/1‑149 (MQ=255)                                                                                                                                                  gcggcCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCg  >  4:55521/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                    AACCAGAGAGGGACGCAACGCGGGGTGAACTTCATGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGT‑GCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCG  >  NC_002947/5229024‑5229316

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