Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	228,066	G→A	6.0%	G298S (GGC→AGC)	PP_0176 →	putative efflux transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	228,066	0	G	A	6.0%	242.1 / 3.2	82	G298S (GGC→AGC)	PP_0176	putative efflux transporter
Reads supporting (aligned to +/- strand):  ref base G (45/32);  new base A (2/3);  total (47/35)
Fisher's exact test for biased strand distribution p-value = 6.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

CCTCGACGACGCACTGGTAGCCCTGGAAGGCCGCGCCCCGGCCCAGTTCGCCCAGCTTCAACCCGTGCTGGAAGCCGCTCGTGAATGGCTCAAAGGCACCGCCGACAGCGCTTCGGTAGCGCACTGGGCCA‑CGTTGCACGAGCAGATCGGCCGCCTGCAGCCCGCTACAGCCGCCCTCGACCAGCGCGCCGAGCTGCTGTTGTCCAACGCCCTGTACCGCCTGACCGAATGGGCCGACCTGTGGCAGGACTGCTGCACCCTGCAGCACGCCCTGCGTACCGATGATGCCAAGCCTTGG  >  NC_002947/227918‑228215                                                                                                                                                      |                                                                                                                                                      ccTCGACGACGCACTGGTAGCCCTGGAAGGCCGCGCCCCGGCCCAGTTCGCCCAGCTTCAACCCGTGCTGGAAGCCGCTCGTGAATGGCTCAAAGGCACCGCCGACAGCGCTTCGGTAGCGCACTGGGCCA‑CGTTGCACGAGCAGATCgg                                                                                                                                                      >  1:1536949/1‑150 (MQ=255)  cgcgacgaCGCACTGGTCGAGCTGGAAGGCCGCGCCCCGTCCGAGTTCGCCCATCTGCAACCCGTGCTGGAAGCCGCTCGTGAATGGCTCCACGGCACCGCCGACAGCGCTTCGGTAGCGCACGGGGACC‑CGTTGCACGAGCAGCTCGgc                                                                                                                                                     <  2:1647860/148‑1 (MQ=255)  cTCGACGACGCACTGGTAGCCCTGGAAGGCCGCGCCCCGGCCCAGTTCGCCCAGCTTCAACCCGTGCTGGAAGCCGCTCGTGAATGGCTCAAAGGCACCGCCGACAGCGCTTCGGTAGCGCACTGGGCCA‑CGTTGCACGAGCAGATCGgc                                                                                                                                                     <  2:1813846/150‑1 (MQ=255)   tCGACGACGCACTGGTAGCCCTGGAAGGCCGCGCCCCGGCCCAGTTCGCCCAGCTTCAACCCGTGCTGGAAGCCGCTCGTGAATGGCTCAAAGGCACCGCCGACAGCGCTTCGGTAGCGCACTGGGCCA‑CGTTGCACGAGCAGATCGgcc                                                                                                                                                    >  2:876717/1‑150 (MQ=255)           gtaCTGGTAGCCCTGGAAGGCCGCGCCCCGGCCCAGTTCGCCCAGCTTCAACCCGTGCTGGAAGCCGCTCGTGAATGGCTCAAAGGCACCGCCGACAGCGCTTCGGTAGCGCACTGGGCCA‑CGTTGCACGAGCAGATCGGCCGCCTGCAg                                                                                                                                            <  2:922498/148‑1 (MQ=255)           gCACTGGTAGCCCTGGAAGGCCGCGCCCCGGCCCAGTTCGCCCAGCTTCAACCCGTGCTGGAAGCCGCTCGTGAATGGCTCAAAGGCACCGCCGACAGCGCTTCGGTAGCGCACTGGGCCA‑CGTTGCACGCGCAGATCGGCCGCCTGCAg                                                                                                                                            >  1:1193404/1‑150 (MQ=255)           gCACTGGTAGCCCTGGAAGGCCGCGCCCCGGCCCAGTTCGCCCAGCTTCAACCCGTGCTGGAAGCCGCTCGTGAATGGCTCAAAGGCACCGCCGACAGCGCTTCGGTAGCGCACTGGGCCA‑CGTTGCACGAGCAGATCGGCCGCCTGCAg                                                                                                                                            >  2:1243724/1‑150 (MQ=255)           gCACTGGTAGCCCTGGAAGGCCGCGCCCCGGCCCAGTTCGCCCAGCTTCAACCCGTGCTGGAAGCCGCTCGTGAATGGCTCAAAGGCACCGCCGACAGCGCTTCGGTAGCGCACTGGGCCA‑CGTTGCACGAGCAGATCGGCCGCCTGCAg                                                                                                                                            >  2:1615298/1‑150 (MQ=255)                  tAGCCCTGGAAGGCCGCGCCCCGGCCCAGTTCGCCCAGCTTCAACCCGTGCTGGAAGCCGCTCGTGAATGGCTCAAAGGCACCGCCGACAGCGCTTCGGTAGCGCACTGGGCCA‑CGTTGCACGAGCAGATCGGCCGCCTGCAGCCCGCTa                                                                                                                                     >  1:713846/1‑150 (MQ=255)                   aGCCCTGGAAGGCCGCGCCCCGGCCCAGTTCGCCCAGCTTCAACCCGTGCTGGAAGCCGCTCGTGAATGGCTCAAAGGCACCGCCGACAGCGCTTCGGTAGCGCACTGGGCCA‑CGTTGCACGAGCAGATCGGCCGCCTGCAGCCCGCTAc                                                                                                                                    >  1:1401676/1‑150 (MQ=255)                     cccTGGAAGGCCGCGCCCCGGCCCAGTTCGCCCAGCTTCAACCCGTGCTGGAAGCCGCTCGTGAATGGCTCAAAGGCACCGCCGACAGCGCTTCAGTAGCGCACTGGGCCA‑CGCCTCACGAGCAGATCGGCCGCCTGCAGCCCGCTACAg                                                                                                                                  >  2:1400745/1‑150 (MQ=255)                       cTGGAAGGCCGCGCCCCGGCCCAGTTCGCCCAGCTTCAACCCGTGCTGGAAGCCGCTCGTGAATGGCTCAAATGCACCGCCGACAGCGCTTCGGTAGCGCACTGGGCCA‑CGTTGCACGAGCAGATCGGCCGCCTGCAGCCCGc                                                                                                                                       >  1:1076268/1‑143 (MQ=255)                       cTGGAAGGCCGCGCCCCGGCCCAGTTCGCCCAGCTTCAACCCGTGCTGGAAGCCGCTCGTGAATGGCTCAAAGGCACCGCCGACAGCGCTTCGGTAGCGCACTGGGCCA‑CGTTGCACGAGCAGATCGGCCGCCTGCAGCCCGc                                                                                                                                       <  2:1076268/143‑1 (MQ=255)                       cTGGAAGGCCGCGCCCCGGCCCAGTTCGCCCAGCTTCAACCCGTGCTGGAAGCCGCTCGTGAATGGCTCAAAGGCACCGCCGACAGCGCTTCGGTAGCGCACTGGGCCA‑CGTTGCACGAGCAGATCAGCCGCCTGCAGCCCGCTACAgcc                                                                                                                                >  1:911135/1‑150 (MQ=255)                       cTGGAAGGCCGCGCCCCGGCCCAGTTCGCCCAGCTTCAACCCGTGCTGGAAGCCGCTCGTGAATGGCTCAAAGGCACCGCCGACAGCGCTTCGGTAGCGCACTGGGCCA‑CGTTGCACGAGCAGATCAGCCGCCTGCAGCCCGCTACAgcc                                                                                                                                >  1:1897859/1‑150 (MQ=255)                        tGGAAGGCCGCGCCCCGGCCCAGTTCGCCCAGCTTCAACCCGTGCTGGAAGCCGCTCGTGAATGGCTCAAAGGCACCGCCGACAGCGCTTCGGTAGCGCACTGGGCCA‑CGTTGCACGAGCAGATCGGCCGCCTGCAGCCCGCTACAgccg                                                                                                                               >  1:30233/1‑150 (MQ=255)                           aaGGCCGCGCCCCGGCCCAGTTCGCCCAGCTTCAACCCGTGCTGGAAGCCGCTCGTGAATGGCTCAAAGGCACCGCCGACAGCGCTTCGGTAGCGCACTGGGCCA‑CGTTGCACGAGCAGATCGGCCGCCTGCAGCCCGCTACAGCCGccc                                                                                                                            >  2:426601/1‑150 (MQ=255)                             tgCCGCGCCCCGGCCCAGTTCGCCCAGCTTCAACCCGTGCTGGAAGCCGCTCGTGAATGGCTCAAAGGCACCGCCGACAGCGCTTCGGTAGCGCACTGGGCCA‑CGTTGCACGAGCAGATCGGCTGCCTGCAGCCCGCTACAGCCGCCCTc                                                                                                                          <  1:1790521/149‑1 (MQ=255)                               ccGCGCCCCGGCCCAGTTCGCCCAGCTTCAACCCGTGCTGGAAGCCGCTCGTGAATGGCTCAAAGGCACCGCCGACAGCGCTTCGGTAGCGCACTGGGGCA‑CGTGGCACGAGCAGATCGGCCGCCTGCAGCCCGCTACAGCCGCCCTCGa                                                                                                                        >  2:1790521/1‑150 (MQ=255)                               ccGCGCCCCGGCCCAGTTCGCCCAGCTTCAACCCGTGCTGGAAGCCGCTCGTGAATGGCTCAAAGGCACCGCCGACAGCGCTTCGGTAGCGCACTGGGCCA‑CGTTGCACGAGCAGATCAGCCGCCTGCAGCCCGCTACAGCCGCCCTCGa                                                                                                                        <  2:1897859/150‑1 (MQ=255)                               ccGCGCCCCGGCCCAGTTCGCCCAGCTTCAACCCGTGCTGGAAGCCGCTCGTGAATGGCTCAAAGGCACCGCCGACAGCGCTTCGGTAGCGCACTGGGCCA‑CGTTGCACGAGCAGATCAGCCGCCTGCAGCCCGCTACAGCCGCCCTCGa                                                                                                                        <  2:911135/150‑1 (MQ=255)                                 gcgcCCCGGCCCAGTTCGCCCAGCTTCAACCCGTGCTGGAAGCCGCTCGTGAATGGCTCAAAGGCACCGCCGACAGCGCTTCGGTAGCGCACTGGGCCA‑CGTTGCACGAGCAGATCGGCCGCCTGCAGCCCGCTACAGCCGCCCTCGAcc                                                                                                                      >  1:930817/1‑150 (MQ=255)                                   gcCCCGGCCCAGTTCGCCCAGCTTCAACCCGTGCTGGAAGCCGCTCGTGAATGGCTCAAAGGCACCGCCGACAGCGCTTCGGTAGCGCACTGGGCCA‑CGTTGCACGAGCAGATCGGCCGCCTGCAGCCCGCTACAGCCGCCCTCGACCAg                                                                                                                    <  1:524795/150‑1 (MQ=255)                                    ccccGGCCCAGTTCGCCCAGCTTCAACCCGTGCTGGAAGCCGCTCGTGAATGGCTCAAAGGCACCGCCGACAGCGCTTCGGTAGCGCACTGGGCCA‑CGTTGCACGAGCAGATCGGCCGCCTGCAGCCCGCTACAGCCGCCCTCGACCAgc                                                                                                                   >  1:1073751/1‑150 (MQ=255)                                     cccGGCCCAGTTCGCCCAGCGTCAACCCGTGCTCGAAGCCGCTCGTGAATGGCTCAAAGGCACCGCCGACAGCGCTTCGGTAGCGCACTGGGGCA‑CGTTGCACGAGCAGATCGGCCGCCTGCAGCCCGCTACAGCCGCCCTCGACCAgcg                                                                                                                  >  1:961798/1‑150 (MQ=255)                                          cccAGTTCGCCCAGCTTCAACCCGTGCTGGAAGCCGCTCGTGAATGGCTCAAAGGCACCGCCGACAGCGCTTCGGTAGCGCACTGGGCCA‑CGTTGCACGAGCAGATCGGCCGCCTGCAGCCCGCTACAGCCGCCCTCGACCAGCGCGCCg                                                                                                             >  1:1304839/1‑150 (MQ=255)                                            cAGTTCGCCCAGCTTCAACCCGTGCTGGAAGCCGCTCGTGAATGGCTAAAAGGCACCGCCGACAGCGCTTCGGTAGCGCACTGGGCCA‑CGTTGCACGAGCAGATCGGCCGCCTGCAGCCCGCTACAGCCGCCCTCGACCAGCGCGCCGAg                                                                                                           >  2:1633006/1‑150 (MQ=255)                                             aGTTCGCCCAGCTTCAACCCGTGCTGGAAGCCGCTCGTGAATGGCTCAAAGGCACCGCCGACAGCGCTTCGGTAGCGCACTGGGCCA‑CGTTGCACGAGCAGATCGGCCGCCTGCAGCCCGCTACAGCCGCCCTCGACCAGCGCGCCGAgc                                                                                                          <  2:1304839/150‑1 (MQ=255)                                             aGTTCGCCCAGCTTCAACCCGTGCTGGAAGCCGCTCGTGAATGGCTCAAAGGCACCGCCGACAGCGCTTCGGTAGCGCACTGGGCCA‑CGTTGCACGAGCAGATCGGCCGCCTGCAGCCCGCTACAGCCGCCCTCGACCAGCGCGCCGAgc                                                                                                          <  2:713846/150‑1 (MQ=255)                                               ttCGCCCAGCTTCAACCCGTGCTGGAAGCCGCTCGTGAATGGCTCAAAGGCACCGCCGACAGCGCTTCGGTAGCGCACTGGGCCA‑CGTTGCACGAGCAGATCGGCCGCCTGCAGCCCGCTACAGCCGCCCTCGAcc                                                                                                                      <  1:103933/136‑1 (MQ=255)                                               ttCGCCCAGCTTCAACCCGTGCTGGAAGCCGCTCGTGAATGGCTCAAAGGCACCGCCGACAGCGCTTCGGTAGCGCACTGGGCCA‑CGTTGCACGAGCAGATCGGCCGCCTGCAGCCCGCTACAGCCGCCCTCGAcc                                                                                                                      >  2:103933/1‑136 (MQ=255)                                                 cGCCCAGCTTCAACCCGTGCTGGAAGCCGCTCGTGAATGGCTCAAAGGCACCGCCGACAGCGCTTCGGTAGCGCACTGGGCCA‑CGTTGCACGAGCAGATCGGCCGCCTGCAGCCCGCTACAGCCGCCCTCGACCAGCGCGCCGAgctgct                                                                                                      <  1:1725926/150‑1 (MQ=255)                                                    ccAGCTTCAACCCGTGCTGGAAGCCGCTCGTGAATGGCTCAAAGGCACCGCCGACAGCGCTTCGGTAGCGCACTGGGCCA‑CGTTGCACGAGCAGATCGGCCGCCTGCAGCCCGCTACAGCCGCCCTCGACCAgc                                                                                                                   >  2:1763863/1‑134 (MQ=255)                                                    ccAGCTTCAACCCGTGCTGGAAGCCGCTCGTGAATGGCTCAAAGGCACCGCCGACAGCGCTTCGGTAGCGCACTGGGCCA‑CGTTGCACGAGCAGATCGGCCGCCTGCAGCCCGCTACAGCCGCCCTCGACCAgc                                                                                                                   <  1:1763863/134‑1 (MQ=255)                                                    ccAGCTTCAACCCGTGCTGGAAGCCGCTCGTGAATGGCTCAAAGGCACCGCCGACAGCGCTTCGGTAGCGCACTGGGCCA‑CGTTGCACGAGCAGATCGGCCGCCTGCAGCCCGCTACAGCCGCCCTCGACCAGCGCGCCGAGCTGCtgtt                                                                                                   >  1:1004440/1‑150 (MQ=255)                                                    ccAGCTTCAACCCGTGCTGGAAGCCGCTCGTGAATGGCTCAAAGGCACCGCCGACAGCGCTTCGGTAGCGCACTGGAAAATTGATTAATG‑GCAGATCGTCCGCCTGGAGCACCCTACACATGCCCTCGTTCATCGCGCCGAGCTGCtgtt                                                                                                   >  1:1468715/1‑150 (MQ=255)                                                    ccAGCTTCAACCCGTGCTGGAAGCCGCTCGTGAATGGCTCAAAGGCACCGCCGACAGCGCTTCGGAAGCGCACTGGGCCA‑CGTTGCACGAGCAGATCGGCCGCCTGCAGCCCGCTACAGCCGCCCTCGACCAGCGCGCCGAGCTGCtgtt                                                                                                   <  2:1907778/150‑1 (MQ=255)                                                      aGCTTCAACCCGTGCTGGAAGCCGCTCGTGAATGGCTAAAAGGCACCGCCGACAGCGCTTCGGTAGCGCACTGGGCCA‑CGTTGCACGAGCAGATCGGCCGCCTGCAGCCCGCTACAGCCGCCCTCGACCAGCGCGCCGAGCTGCtgttgt                                                                                                 <  1:1633006/150‑1 (MQ=255)                                                             aCCCGTGCTGGAAGCCGCTCGTGAATGGCTCAAAGGCACCGCCGACAGCGCTTCGGTAGCGCACTGGGCCA‑CGTTGCACGAGCAGATCGGCCGCCTGCAGCCCGCTACAGCCGCCCTCGACCAGCGCGCCGAGCTGCTGTTGTCCAACGc                                                                                          >  2:1104584/1‑150 (MQ=255)                                                             aCCCGTGCTGGAAGCCGCTCGTGAATGGCTCAAAGGCACCGCCGACAGCGCTTCGGTAGCGCACTGGGCCA‑CGTTGCACGAGCAGATCGGCCGCCTGCAGCCCGCTACAGCCGCCCTCGACCAGCGCGCCGAGCTGCTGTTGTCCAACGc                                                                                          >  2:1104716/1‑150 (MQ=255)                                                              cccGTGCTGGAAGCCGCTCGTGAATGGCTCAAAGGCACCGCCGAAAGAGCTTCGGGAGCGCACTGGGCCA‑CGTTGCACGAGCAGATCGGCCGCCTGCAGCCCGCTACAGCCGCCCTCGACCAGCGCGCCGAGCTGCTGTTGTCCAACGcc                                                                                         <  2:1297960/150‑1 (MQ=255)                                                                cGTGCTGGAAGCCGCTCGTGAATGGCTCAAAGGCACCGCCGACAGCGCTTCGGTAGCGCACTGGGCCA‑CGTTGCACGAGCAGATCGGCCGCCTGCAGCCCGCTACAGCCGCCCTCGACCAGCGCGCCGAGCTGCTGTTGTCCAACGCCCt                                                                                       >  2:727348/1‑150 (MQ=255)                                                                  tGCTGGAAGCCGCTCGTGAATGGCTCAAAGGAACCGCCTACATCGCTTCGGTAGCGCACTGGGCCA‑CGTTGCACGAGCAGATCGGCCGCCTGCAGCCCGCTACAGCCGCCCTCGACCAGCGCGCCGAGCTGCTGTTGTCCAACGCCCTGt                                                                                     <  2:1589563/150‑1 (MQ=255)                                                                     tGGAAGCCGCTCGTGAATGGCTCAAAGGCACCGCCGACAGCGCTTCGGTAGCGCACTGGGCCA‑CGTTGCACGAGCAGATCGGCCGCCTGCAGCCCGCTACAGCCGCCCTCGACCAGCGCGCCGAGCTGCTGTTGTCCAACGCCCTGTAcc                                                                                  >  2:1592674/1‑150 (MQ=255)                                                                     tGGAAGCCGCTCGTGAATGGCTCAAAGGCACCGCCGACAGCGCTTCGGTAGCGCACTGGGCCA‑CGTTGCACGAGCAGATCGGCCGCCTGCAGCCCGCTACAGCCGCCATCGACAAGCGCGCCGAGCTGCTGTTGTCCAACGCCCTGTAcc                                                                                  >  2:1590109/1‑150 (MQ=255)                                                                       gAAGCCGCTCGTGAATGGCTCAAAGGCACCGCCGACAGCGCTTCGGTAGCGCACTGGGCCA‑CGTTGCACGAGCAGATCGGCCGCCTGCAGCCCGCTACAGCCGCCCTCGACCAGCGCGCCGAGCTGCTGTTGTCCAACGCCCTGTACCGc                                                                                <  1:1321931/150‑1 (MQ=255)                                                                          gCCGCTCGTGAATGGCTCAAAGGCACCGCCGACAGCGCTTCGGTAGCGCACTGGGCCA‑CGTTGCACGAGCAGATCGGCCGCCTGCAGCCCGCTACAGCCGCCCTCGACCAGCGCGCCGAGCTGCTGTTGTCCAACGCCCTGTACCGCCTg                                                                             >  2:1872640/1‑150 (MQ=255)                                                                            cGCTCGTGAATGGCTCAAAGGCACCGCCGACAGCGCTTCGGTAGCGCACTGGGCCA‑CGTTGCACGAGCAGATCGGCCGCCTGCAGCCCGCTACAGCCGCCCTCGACCAGCGCGCCGAGCTGCTGTTGTCCAACGCCCTGTACCGCCTGAc                                                                           >  2:745481/1‑150 (MQ=255)                                                                             gCTCGTGAATGGCTCAAAGGCACCGCCGACAGCGCTTCGGTAGCGCACTGGGCCA‑CGTTGCACGAGCAGATCGGCCGCCTGCAGCCCGCTACAGCCGCCCTCGACCAGCGCGCCGAGCTGCTGTTGTCCAACGCCCTGTACCGCCTGAcc                                                                          >  1:653755/1‑150 (MQ=255)                                                                             gCTCGTGAATGGCTCAAAGGCACCGCCGACAGCGCTTCGGTAGCGCACTGGGCCA‑CGTTGCACGAGCAGATCGGCCGCCTGCAGCCCGCTACAGCCGCCCTCGACCAGCGCGCCGAGCTGCTGTTGTCCAACGCCCTGTACCGCCTGAcc                                                                          >  1:170153/1‑150 (MQ=255)                                                                             gCTCGTGAATGGCTCAAAGGCACCGCCGACAGCGCTTCGGTAGCGCACTGGGCCA‑CGTTGCACGAGCAGATCAGCCGCCTGCAGCCCGCTACAGCCGCCCTCGACCAGCGCGCCGAGCTGCTGTTGTCCAACGCCCTGTACCGCCTGAcc                                                                          <  2:1401676/150‑1 (MQ=255)                                                                                cGTGAATGGCTCAAAGGCACCGCCGACAGCGCTTCGGTAGCGCACTGGGCCA‑CGTTGCACGAGCAGATCGGCCGCCTGCAGCCCGCTACAGCCGCCCTCGACCAGCGCGCCGAGCTGCTGTTGTCCAACGCCCTGTACCGCCTGACCGaa                                                                       <  1:540107/150‑1 (MQ=255)                                                                                 gTGAATGGCTCAAAGGCACCGCCGACAGCGCTTCGGTAGCGCACTGGGCCA‑CGTTGCACGAGCAGATCGGCCGCCTGCAGCCCGCTACAGCCGCCCTCGACCAGCGCGCCGAGCTGCTGTTGTCCAACGCCCTGTACCGCCTGACCGAAt                                                                      >  2:169737/1‑150 (MQ=255)                                                                                  tGAATGGCTCAAAGGCACCGCCGACAGCGCTTCGGTAGCGCACTGGGCCA‑CGTTGCACGAGCAGATCGGCCGCCTGCAGCCCGCTACAGCCGCCCTCGACCAGCGCGCCGAGCTGCTGTTGTCCAACGCCCTGTACCGCCTGACCGAATg                                                                     >  1:1560035/1‑150 (MQ=255)                                                                                    aaTGGCTCAAAGGCACCGCCGACAGCGCTTCGGTAGCGCACTGGGCCA‑CGTTGCACGAGCAGATCGGCCGCCTGCAGCCCGCTACAGCCGCCCTCGACCAGCGCGCCGAGCTGCTGTTGTCCAACGCCCTGTACCGCCTGACCGAATggg                                                                   <  1:455211/150‑1 (MQ=255)                                                                                    aaTGGCTCAAAGGCACCGCCGACAGCGCTTCGGTAGCGCACTGGGCCA‑CGTTGCACGAGCAGATCGGCCGCCTGCAGCCCGCTACAGCCGCCCTCGACCAGCGCGCCGAGCTGCTGTTGTCCAACGCCCTGTACCGCCTGACCGAATggg                                                                   <  1:727348/150‑1 (MQ=255)                                                                                     aTGGCTCAAAGGCACCGCCGACAGCGCTTCGGTAGCGCACTGGGCCA‑CGTTGCACGAGCAGATCGGCCGCCTGCAGCCCGCTACAGCCGCCCTCGACCAGCGCGCCGAGCTGCTGTTGTCCAACGCCCCGTACCGCCTGACCGAATGGGc                                                                  <  1:876717/150‑1 (MQ=255)                                                                                      tGGCTCAAAGGCACCGCCGACAGCGCTTCGGTAGCGCACTGGGCCA‑CGTTGCACGAGCAGATCGGCCGCCTGCAGCCCGCTACAGCCGCCCTCGACCAGCGCGCCGAGCTGCTGTTGTCCAACGCCCTGTACCGCCTGACCGAATGGGcc                                                                 <  2:127982/150‑1 (MQ=255)                                                                                          tCAAAGGCACCGCCGACAGCGCTTCGGTAGCGCACTGGGCCA‑CGTTGCACGAGCAGATCGGCCGCCTGCAGCCCGCTACAGCCGCCCTCGACCAGCGCGCCGAGCTGCTGTTGTCCAACGCCCTGTACCGCCTGACCGAATGGGCCGAcc                                                             >  1:1229226/1‑150 (MQ=255)                                                                                          tCAAAGGCACCGCCGACAGCGCTTCGGTAGCGCACTGGGCCA‑CGTTGCACGAGCAGATCGGCCGCCTGCAGCCCGCTACAGCCGCCCTCGACCAGCGCGCCGAGCTGCTGTTGTCCAACGCCCTGTACCGCCTGACCGAATGGGCCGAcc                                                             >  2:1352035/1‑150 (MQ=255)                                                                                          tCAAAGGCACCGCCGACAGCGCTTCGGTAGCGCACTGGGCCA‑CGTTGCACGAGCAGATCGGCCGCCTGCAGCCCGCTACAGCCGCCCTCGACCAGCGCGCCGAGCTGCTGTTGTCCAACGCCCTGTACCGCCTGACCGAATGGGCCGAcc                                                             >  2:1909105/1‑150 (MQ=255)                                                                                          tCAAAGGCACCGCCGACAGCGCTTCGGTAGCGCACTGGGCCA‑CGTTGCACGAGCAGATCGGCCGCCTGCAGCCCGCTACAGCCGCCCTCGACCAGCGCGCCGAGCTGCTGTTGTCCAACGCCCTGTACCGCCTGACCGAAGGGGCCGAcc                                                             >  1:1734743/1‑150 (MQ=255)                                                                                             aaGGCACCGCCGACAGCGCTTCGGTAGCGCACTGGGCCA‑CGTTGCACGAGCAGATCGGCCGCCTGCAGCCCGCTACAGCCGCCCTCGACCAGCGCGCCGAGCTGCTGTTGTCCAACGCCCTGTACCGCCTGACCGAATGGGCCGACCtgt                                                          <  1:1851666/150‑1 (MQ=255)                                                                                                  aCCGCCGACAGCGCTTCGGTAGCGCACTGGGCCA‑CGTTGCACGAGCAGATCGGCCGCCTGCAGCCCGCTACAGCCGCCCTCGACCAGCGCGCCGAGCTGCTGTTGTCCAACGCCCTGTACCGCCTGACCGAATGGGCCGACCTGTGGCAg                                                     >  2:443389/1‑150 (MQ=255)                                                                                                         aCAGCGCTTCGGTAGCGCACTGGGCCA‑CGTTGCACGAGCAGATCGGCCGCCTGCAGCCCGCTACAGCCGCCCTCGACCAGCGCGCCGAGCTGCTGTTGTCCAACGCCCTGTACCGCCTGACCGAATGGGCCGACCTGTGGCAGGActgct                                              <  2:1536949/150‑1 (MQ=255)                                                                                                             cgcTTCGGTAGCGCACTGGGCCA‑CGTTGCACGAGCAGATCGGCCGCCTGCAGCCCGCTACAGCCGCCCTCGACCAGCGCGCCGAGCTGCTGTTGTCCAACGCCCTGTACCGCCTGACCGAATGGGCCGACCTGTGGCAGGACTGCTGCAc                                          <  2:1430140/150‑1 (MQ=255)                                                                                                             cgcTTCGGTAGCGCACTGGGCCA‑CGTTGCACGAGCAGATCGGCCGCCTGCAGCCCGCTACAGCCGCCCTCGACCAGCGCGCCGAGCTGCTGTTGTCCAACGCCCTGTACCGCCTGACCGAATGGGCCGACCTGTGGCAGGACTGCTGCAc                                          >  2:637054/1‑150 (MQ=255)                                                                                                              gcTTCGGTAGCGCACTGGGCCA‑CGTTGCACGAGCAGATCGGCAGCCTGCAGCCCGCTACAGCCGCCCTCGACCAGCGCGCCGAGCTGCTGTTGTCCAACGCCCTGTACCGCCTGACCGTATGGGCCGACCTGTGCCAGGACTGCCGCAcc                                         >  2:436498/1‑150 (MQ=255)                                                                                                                 tCGGTAGCGCACTGGGCCA‑CGTTGCACGAGCAGATCGGCCGCCTGCAGCCCGCTACAGCCGCCCTCGACCAGCGCGCCGAGCTGCTGTTGTCCAACGCCCTGTACCGCCTGACCGAATGGGCCGACCTGTGGCAGGACTGCTGCACCCTg                                      <  1:169737/150‑1 (MQ=255)                                                                                                                   ggTAGCGCACTGGGCCA‑CGTTGCACGAGCAGATCGGCCGCCTGCAGCCCGCTACAGCCGCCCTCGACCAGCGCGCCGAGCTGCTGTTGTCCAACGCCCTGTACCGCCTGACCGAATGGGCCGACCTGTGGCAGGACTGCTGCACCCTgca                                    <  1:1352035/150‑1 (MQ=255)                                                                                                                   ggTAGCGCACTGGGCCA‑CGTTGCACGAGCAGATCGGCCGCCTGCAGCCCGCTACAGCCGCCCTCGACCAGCGCGCCGAGCTGCTGTTGTCCAACGCCCTGTACCGCCTGACCGAATGGGCCGACCTGTGGCAGGACTGCTGCACCCTgca                                    <  2:1004440/150‑1 (MQ=255)                                                                                                                    gTAGCGCACTGGGCCA‑CGTTGCACGAGCAGATCGGCCNCCTGCAGCCNGCTANAGCCGCCCTCGACCAGCGCGCCGAGCCGCTGTTGTCCAACGCCCTGTACCGCCTGACCGAATGGGCCGACCTGTGGCAGGACTGCTGCACCCTgnag                                   >  1:46936/1‑150 (MQ=255)                                                                                                                          cACTGGGCCA‑CGTTGCACGAGCAGATCGGCCGCCTGCAGCCCGCTACAGCCGCCCTCGACCAGCGCGCCGAGCTGCTGTTGTCCAACGCCCTGTACCGCCTGACCGAATGGGCCGACCTGTGGCAGGACTGCTGCACCCTGCAGCANGcc                             <  2:1193404/150‑1 (MQ=255)                                                                                                                                 ccA‑CGTTGCACGAGCAGATCGGCCGCCTGCAGCCCGCTACAGCCGCCCTCGACCAGCGCGCCGAGCTGCTGTTGTCCAACGCCCTGTACCGCCTGACCGAATGGGCCGACCTGTGGCAGGACTGCTGCACCCTGCAGCACGCCCTGCGTa                      <  2:30233/150‑1 (MQ=255)                                                                                                                                  cA‑CGTTGCACGAGCAGATCGGCCGCCTGCAGCCCGCTACACACGAACTCGACCAGCGCGCCGAGCTGCTGTTGTCCAACGCCCTGTACCGCCTGACCGAATGGGCCGACCTGTGGCAGGACTGCTGCACCCTGCAGCACGCCCTGCGTAc                     >  1:976850/1‑150 (MQ=255)                                                                                                                                      gTTGCACGAGCAGATCGGCCGCCTGCAGCCCGCTACAGCCGCCCTCGACCAGCGCGCCGAGCTGCTGTTGTCCAACGCCCTGTACCGCCTGACCGAATGGGCCGACCTGTgg                                                        <  1:550141/112‑1 (MQ=255)                                                                                                                                      gTTGCACGAGCAGATCGGCCGCCTGCAGCCCGCTACAGCCGCCCTCGACCAGCGCGCCGAGCTGCTGTTGTCCAACGCCCTGTACCGCCTGACCGAATGGGCCGACCTGTgg                                                        >  2:550141/1‑112 (MQ=255)                                                                                                                                          cACGAGCAGATCGGCCGCCTGCAGCCCGCTACAGCCGCCCTCGACCAGCGCGCCGAGCTGCTGTTGTCCAACGCCCTGTACCGCCTGACCGAATGGGCCGACCTGTGGCAGGACTGCTGCACCCTGCAGCACGCCCTGCGTACCgatgat              <  2:81820/150‑1 (MQ=255)                                                                                                                                                  gATCGGCCGCCTGCAGCCCGCTACAGCCGCCCTCGACCAGCGCGCCGAGCTGCTGTTGTCCAACGCCCTGTACCGCCTGACCGAATGGGCCGACCTGTGGCAGGACTGCTGCACCCTGCAGCACGCCCTGCGTACCgatgat              <  1:502364/142‑1 (MQ=255)                                                                                                                                                  gATCGGCCGCCTGCAGCCCGCTACAGCCGCCCTCGACCAGCGCGCCGAGCTGCTGTTGTCCAACGCCCTGTACCGCCTGACCGAATGGGCCGACCTGTGGCAGGACTGCTGCACCCTGCAGCACGCCCTGCGTACCgatgat              >  2:502364/1‑142 (MQ=255)                                                                                                                                                  gATCGGCCGCCTGCAGCCCGCTACAGCCGCCCTCGACCAGCGCGCCGAGCTGCTGTTGTCCAACGCCCTGTACCGCCTGACCGAATGGGCCGACCTGTGGCAGGACTGCTGCACCCTGCAGCACGCCCTGCGTACCGATGATGCCAAGcc      >  1:209785/1‑150 (MQ=255)                                                                                                                                                   aTCGGCCGCCTGCAGCCCGCTACAGCCGCCCTCGACCAGCGCGCCGAGCTGCTGTTGTCCAACGCCCTGTACCGCCTGACCGAATGGGCCGACCTGTGGCAGGACTGCTGCACCCTGCAGCACGCCCTGCGTACCGATGATGCCAAGCCt     >  1:667554/1‑150 (MQ=255)                                                                                                                                                   aTCGGCCGCCTGCAGCCCGCTACAGCCGCCCTCGACCAGCGCGCCGAGCTGCTGTTGTCCAACGCCCTGTACCGCCTGACCGAATGGGCCGACCTGTGGCAGGACTGCTGCACCCTGCAGCACGCCCTGCGTACCGATGATGCCAAGCCt     >  2:1437220/1‑150 (MQ=255)                                                                                                                                                      ggCCGCCTGCAGCCCGCTACAGCCGCCCTCGACCAGCGCGCCGAGCCGCTGTTGTCCAACGCCCTGTACCGCCTGACCGAATGGGCCGACCTGTGGCAGGACTGCTGCACCCTGCAGCACGCCCTGCGTACCGATGATGCCAAGCCTTgg  <  2:46936/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                      CCTCGACGACGCACTGGTAGCCCTGGAAGGCCGCGCCCCGGCCCAGTTCGCCCAGCTTCAACCCGTGCTGGAAGCCGCTCGTGAATGGCTCAAAGGCACCGCCGACAGCGCTTCGGTAGCGCACTGGGCCA‑CGTTGCACGAGCAGATCGGCCGCCTGCAGCCCGCTACAGCCGCCCTCGACCAGCGCGCCGAGCTGCTGTTGTCCAACGCCCTGTACCGCCTGACCGAATGGGCCGACCTGTGGCAGGACTGCTGCACCCTGCAGCACGCCCTGCGTACCGATGATGCCAAGCCTTGG  >  NC_002947/227918‑228215

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