Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,649,376	C→A	7.4%	E576* (GAG→TAG) ‡	PP_1446 ←	TonB‑dependent receptor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,649,376	0	C	A	7.4%	286.1 / 4.5	91	E576* (GAG→TAG) ‡	PP_1446	TonB‑dependent receptor
Reads supporting (aligned to +/- strand):  ref base C (34/50);  new base A (4/3);  total (38/53)
Fisher's exact test for biased strand distribution p-value = 4.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

GTAGCCATCGGTGCTGGTTTCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGG  >  NC_002947/1649236‑1649521                                                                                                                                              |                                                                                                                                                  gtAGCCATCGGTGCTGGTTTCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTGTc                                                                                                                                          >  1:1975617/1‑150 (MQ=255)  tAGCCATCGGTGCTGGTTTCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTGTCg                                                                                                                                         >  2:482752/1‑150 (MQ=255)  tAGCCATCGGTGCTGGTTTCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGCGAGCCACTgccc                                                                                                                                         <  1:764065/150‑4 (MQ=255)   aGCCATCGGTGCTGGTTTCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTgcccc                                                                                                                                        <  1:390777/150‑5 (MQ=255)   aGCCATCGGTGCTGGTTTCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTGTCGa                                                                                                                                        <  1:482752/150‑1 (MQ=255)   aGCCATCGGTGCTGGTTTCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTGTCGa                                                                                                                                        <  2:1157745/150‑1 (MQ=255)   aGCCATCGGTGATGGTTTCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTGTCGa                                                                                                                                        <  1:1516256/150‑1 (MQ=255)    gCCATCGGTGCTGGTTTCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGCTCg                                                                                                                                                  >  2:1180829/1‑139 (MQ=255)    gCCATCGGTGCTGGTTTCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGCTCg                                                                                                                                                  <  1:1180829/139‑1 (MQ=255)    gCCATCGGTGCTGGTTTCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTGTCGAg                                                                                                                                       <  2:632670/150‑1 (MQ=255)    gCCATCGGTGCTGGTTTCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTGTCGAg                                                                                                                                       <  1:2399426/150‑1 (MQ=255)    gCCATCGGTGCTGGTTTCGTTGGCGGGTTTGCGGTGCTACGACGCTTCGTGCGGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGACAGTGGCGAGCCACTgcccc                                                                                                                                        <  1:1141563/150‑5 (MQ=255)       aTCGGTGCTGGTTTCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTGAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTGTCGAGGtt                                                                                                                                    >  1:305252/1‑150 (MQ=255)       aTCGGTGCTGGTTTCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTGTCGAGGtt                                                                                                                                    >  1:1507986/1‑150 (MQ=255)        tCGGTGCTGGTTTCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTGTCGAGGTTc                                                                                                                                   <  1:980810/150‑1 (MQ=255)           gTGCTGGTTTCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTGTCGAGGTTCTTg                                                                                                                                <  1:2104380/150‑1 (MQ=255)               tGGTTTCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTGTCGAGGTTCTTGGCTc                                                                                                                            >  1:155201/1‑150 (MQ=255)                 gTTTCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCgg                                                                                                                          <  2:2251794/150‑1 (MQ=255)                 gTTTCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCgg                                                                                                                          <  2:856167/150‑1 (MQ=255)                 gTTTCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCgg                                                                                                                          <  1:359916/150‑1 (MQ=255)                 gTTTCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCgg                                                                                                                          <  1:1472931/150‑1 (MQ=255)                 gTTTCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTCCAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCgg                                                                                                                          <  2:1349416/150‑1 (MQ=255)                  tttCGTTGGCGGGTCTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCggg                                                                                                                         >  2:1612802/1‑150 (MQ=255)                    tCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCACGGACCAGGCCGCTGTTCGTGCGCAGCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGta                                                                                                                       <  1:2099642/150‑1 (MQ=255)                      gTTAGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTCCACGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGtata                                                                                                                     <  1:1468731/150‑1 (MQ=255)                          gCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCa                                                                                                                 >  2:1149700/1‑150 (MQ=255)                          gCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCa                                                                                                                 >  2:1452688/1‑150 (MQ=255)                           cGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCAc                                                                                                                <  2:1314737/150‑1 (MQ=255)                                  gCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACaa                                                                                                         <  1:1149700/150‑1 (MQ=255)                                  gCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACaa                                                                                                         <  2:155201/150‑1 (MQ=255)                                    ggTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACt                                                                                                       >  2:1981308/1‑150 (MQ=255)                                     gTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTGTCGAGGTTCTTg                                                                                                                                >  2:723384/1‑124 (MQ=255)                                     gTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTGTCGAGGTTCTTg                                                                                                                                <  1:723384/124‑1 (MQ=255)                                         tgcGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTGTCGAg                                                                                                                                       <  1:1658765/113‑1 (MQ=255)                                         tgcGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTGTCGAg                                                                                                                                       >  2:1658765/1‑113 (MQ=255)                                          gcGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTGTCGAg                                                                                                                                       >  1:1340359/1‑112 (MQ=255)                                          gcGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTGTCGAg                                                                                                                                       <  2:1340359/112‑1 (MQ=255)                                              cGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGaa                                                                                             <  2:1975617/150‑1 (MQ=255)                                                 aGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCt                                                                                          >  1:1157245/1‑150 (MQ=255)                                                    gTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGcc                                                                                       <  1:2080078/150‑1 (MQ=255)                                                         gcACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTGTCGAGGTTCtt                                                                                                                                 <  2:493176/103‑1 (MQ=255)                                                         gcACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTGTCGAGGTTCtt                                                                                                                                 >  1:493176/1‑103 (MQ=255)                                                          cACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGa                                                                                                            <  1:1775067/123‑1 (MQ=255)                                                          cACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGa                                                                                                            >  2:1775067/1‑123 (MQ=255)                                                          cACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACg                                                                                 <  2:1157245/150‑1 (MQ=255)                                                             gTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGtt                                                                              >  2:1308794/1‑150 (MQ=255)                                                                 aCCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCg                                                                          <  1:1612802/150‑1 (MQ=255)                                                                  cccGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGtata                                                                                                                     <  2:2066032/106‑1 (MQ=255)                                                                  cccGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGtata                                                                                                                     >  1:2066032/1‑106 (MQ=255)                                                                   ccGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGc                                                                        <  1:1523660/150‑1 (MQ=255)                                                                        cttGCCAGCGGTCCAGCACCCAGAGCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGG‑CGGTGGCTCGCCACTGTCGAGGGTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACt                                                                   <  2:128722/148‑1 (MQ=255)                                                                          tGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGc                                                                 >  2:449408/1‑150 (MQ=255)                                                                           gCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGaa                                                                                             <  1:645750/121‑1 (MQ=255)                                                                           gCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGaa                                                                                             >  2:645750/1‑121 (MQ=255)                                                                              aGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCg                                                                                                              <  2:239800/101‑1 (MQ=255)                                                                              aGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCg                                                                                                              >  1:239800/1‑101 (MQ=255)                                                                                   tCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTcc                                                        >  2:1749474/1‑150 (MQ=255)                                                                                      aGCTCCCAGACCCGGGCGCTGTTCAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGATCGGGCCTGCCTGGGCAGATAGCGCGGGTATAGTCACCGGCCAAATCCAACGCGAATCTCCGGCGACGGTGCTCGGCCAACTGCCAGCGGTCCTgg                                                     <  1:1749474/150‑1 (MQ=255)                                                                                      aGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTgg                                                     >  2:525033/1‑150 (MQ=255)                                                                                           ccAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAg                                                              <  1:1592211/136‑1 (MQ=255)                                                                                           ccAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAg                                                              >  2:1592211/1‑136 (MQ=255)                                                                                           ccAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTc                                                >  2:511270/1‑150 (MQ=255)                                                                                              gACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGAt                                             <  1:449408/150‑1 (MQ=255)                                                                                                ccAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGc                                           <  1:525033/150‑1 (MQ=255)                                                                                                 cAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAgc                                                                                                  <  2:2224486/94‑1 (MQ=255)                                                                                                 cAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGCTAGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAgc                                                                                                  >  1:2224486/1‑94 (MQ=255)                                                                                                  aGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCg                                         <  2:1673927/150‑1 (MQ=255)                                                                                                    gCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTa                                       <  1:1981308/150‑1 (MQ=255)                                                                                                       gCTGTTCAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCg                                         >  2:1215441/1‑145 (MQ=255)                                                                                                       gCTGTTCAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGaa                                    >  1:1357420/1‑150 (MQ=255)                                                                                                       gCTGTTCAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCCCGATGCCg                                         <  1:1215441/145‑1 (MQ=255)                                                                                                         tGTTCAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAgc                                  <  1:1658816/150‑1 (MQ=255)                                                                                                                gcgcAGCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGTGCCTCGATGCCGTAGAAGCGCGCCCg                           >  1:2106354/1‑150 (MQ=255)                                                                                                                   cggCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCAc                        <  1:955674/148‑1 (MQ=255)                                                                                                                     gCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACgg                                                                                <  2:2140712/92‑1 (MQ=255)                                                                                                                     gCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACgg                                                                                >  1:2140712/1‑92 (MQ=255)                                                                                                                         ggCGATGCGGGG‑CAGAGGCTCGCCACTGTCGAGATTCTTGGCTCGGGTATAGTCGCCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTg                  >  2:637671/1‑150 (MQ=255)                                                                                                                             aTGCGGGG‑CAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCg                                                                                                              >  2:345077/1‑54 (MQ=255)                                                                                                                             aTGCGGGG‑CAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCg                                                                                                              <  1:345077/54‑1 (MQ=255)                                                                                                                               gCGGGG‑CAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTg            <  2:1507986/150‑1 (MQ=255)                                                                                                                               gCGGGG‑CAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTg            <  1:511270/150‑1 (MQ=255)                                                                                                                               gCGGGG‑CAGTCGCTCGCCAATGAAGAGGTTCTTGGCTTGGGTATAGTCACCGGAAAACTCCAGCGCGAAGCTGCGGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGTGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTg            <  2:305252/150‑1 (MQ=255)                                                                                                                                cctcgacAGTGGCTAGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTgg                 >  1:2384318/7‑145 (MQ=255)                                                                                                                                cctcgacAGTGGCGAGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTgg                 <  2:2384318/139‑1 (MQ=255)                                                                                                                                 gggg‑CAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGt          <  1:1452688/150‑1 (MQ=255)                                                                                                                                 ctcgacAGTGGCGAGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTCACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGGACTgg           >  2:363526/6‑150 (MQ=255)                                                                                                                                 ctcgacAGTGGCGAGCCACTGCCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTgg           >  1:1433365/6‑150 (MQ=255)                                                                                                                                    g‑CAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTTCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTAtt       >  1:2077123/1‑150 (MQ=255)                                                                                                                                          ggCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCggg  >  1:793003/1‑150 (MQ=255)                                                                                                                                             tCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTg            >  2:1095157/1‑137 (MQ=255)                                                                                                                                             tCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTg            <  1:1095157/137‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                  GTAGCCATCGGTGCTGGTTTCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGG‑CAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGG  >  NC_002947/1649236‑1649521

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