Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,229,201:1	+C	9.8%	intergenic (‑143/+33)	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,201	1	.	C	9.8%	226.4 / 18.9	61	intergenic (‑143/+33)	PP_4606/PP_4607	ferric siderophore receptor/transmembrane sensor
Reads supporting (aligned to +/- strand):  ref base . (27/28);  new base C (4/2);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.63e-01

TCATGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAAT‑CCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGCTGCCGCTGACC  >  NC_002947/5229055‑5229345                                                                                                                                                    |                                                                                                                                                  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)         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                                                                                                             <  3:204605/137‑1 (MQ=255)                                                 tttCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCg                                                                                                             >  4:204605/1‑137 (MQ=255)                                                 tttCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCa                                                                                                  >  4:203282/1‑148 (MQ=255)                                                     aCTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCaaa                                                                                                <  8:144613/146‑1 (MQ=255)                                                     aCTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCaaa                                                                                                >  7:144613/1‑146 (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                                                                                             <  8:112526/142‑1 (MQ=255)                                                            acacGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGc                                                                                             >  7:112526/1‑142 (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                                                    >  2:40208/1‑149 (MQ=255)                                                                                              cTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAg                                                    >  6:17749/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                   >  5:116458/1‑146 (MQ=255)                                                                                                                                 gTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAAT‑CCAGCGCCAGCAGCGGCTCGGCCACCg                   <  6:116458/146‑1 (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)                                                                                                                                           cccTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAAT‑CCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGCTGCCGCt      >  2:280139/1‑149 (MQ=255)                                                                                                                                               aCGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCAt                                                                <  4:143458/88‑1 (MQ=255)                                                                                                                                               aCGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCAt                                                                >  3:143458/1‑88 (MQ=255)                                                                                                                                                cGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAAT‑CCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGCTGCCGCTGAcc  <  5:224545/148‑1 (MQ=255)                                                                                                                                                cGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAAT‑CCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGCTGCCGCTGAc   <  3:222240/147‑1 (MQ=255)                                                                                                                                                cGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAAT‑CCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGCTGCCGCTGAc   <  2:123741/147‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                  TCATGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAAT‑CCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGCTGCCGCTGACC  >  NC_002947/5229055‑5229345

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