Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,722,554	A→T	5.3%	D352E (GAT→GAA)	htpG ←	chaperone protein HtpG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,722,554	0	A	T	5.3%	287.3 / 2.0	90	D352E (GAT→GAA)	htpG	chaperone protein HtpG
Reads supporting (aligned to +/- strand):  ref base A (42/43);  new base T (3/2);  total (45/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTGTTGGCGAAGTCTTCGGCCGGGCCCTCTTTCAGCACCTGACCGAAGTTTTTCCAGAAGCCCTTGTATTTCTCGGGTTCGTTCTTCGCCAGCTTCTCCAGCATGTCCAGCACACGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCAGCGACAGGTCATTGGAATCAACCACACCCTTGATGAAGCGCAGGTACAGCGGCAGGAACGACTCGGCCTGGTCCATGATGAACACACGCTGCACGTACAGCTTCA  >  NC_002947/4722413‑4722702                                                                                                                                              |                                                                                                                                                     cTTGTTGGCGAAGTCTTCGGCCGGGCCCTCTTTCAGCACCTGACCGAAGTTTTTCCAGAAGCCCTTGTATTTCTCGGGTTCGTTCTTCGCCAGCTTCTCCAGCATGTCCAGCACACGCTTGGTCAGCGCGGTCTTCATCGAATCgatgat                                                                                                                                              >  1:2102242/1‑150 (MQ=255)  ttgttgGCGAAGTCTTCGGCCGGGCCCTCTTTCAGCACCTGACCGAAGTTTTTCCAGAAGCCCTTGTATTTCTCGGGTTCGTTCTTCGCCAGCTTCTCCAGCATGTCCAGCACACGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATc                                                                                                                                             <  1:77270/150‑1 (MQ=255)  ttgttgGCGAAGTCTTCGGCCGGGCCCTCTTTCAGCACCTGACCGAAGTTTTTCCAGAAGCCCTTGTATTTCTCGGGTTCGTTCTTCGCCAGCTTCTCCAGCATGTCCAGCACACGCTTGGTCAGCGCGGTCTTCATCGAATCAATGATc                                                                                                                                             >  2:1942645/1‑150 (MQ=255)   tgttgGCGAAGTCTTCGGCCGGGCCCTCTTTCAGCACCTGACCGAAGTTTTTCCAGAAGCCCTTGTATTTCTCGGGTTCGTTCTTCGCCAGCTTCTCCAGCATGTCCAGCACACGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCg                                                                                                                                            <  1:1755854/150‑1 (MQ=255)      tgGCGAAGTCTTCGGCCGGGCCCTCTTTCAGCACCTGACCGAAGTTTTTCCAGAAGCCCTTGTATTTCTCGGGTTCGTTCTTCGCCAGCTTCTCCAGCATGTCCAGCACACGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGt                                                                                                                                         <  1:378188/150‑1 (MQ=255)        gCGAAGTCTTCGGCCGGGCCCTCTTTCAGCACCTGACCGAAGTTTTTCCAGAAGCCCTTGTATTTCTCGGGTTCGTTCTTCGCCAGCTTCTCCAGCATGTCCAGCACACGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCt                                                                                                                                       <  1:1695040/150‑1 (MQ=255)           aaGTCTTCGGCCGGGCCCTCTTTCAGCACCTGGCCGAAGTTTTTCCAGAAGCCCTTGTATTTCTCGGAGGCGTTCTTCGCCAGCTTCTCCAGCATGTCCAGCACACGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTc                                                                                                                                    <  1:1355521/150‑1 (MQ=255)                ttCGGCCGGGCCCTCTTTCAGCACCTGACCGAAGGTTTTCCAGAAGCCCTTGTATTTCTCGGGTTCGTTCTTCGCCAGCTTCTCCAGCATGTCCAGCACACGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAg                                                                                                                               >  1:921969/1‑150 (MQ=255)                 tCGGCCGGGCCCTCTTTCAGCACCTGACCGAAGTTTTTCCAGAAGCCCTTGTATTTCTCGGGTTCGTTCTTCGCCAGCTTCTCCAGCATGTCCAGCACACGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAgg                                                                                                                              <  2:2024560/150‑1 (MQ=255)                    gCCGGGCCCTCTTTCAGCACCTGACCGAAGTTTTTCCAGAAGCCCTTGTATTTCTCGGGTTCGTTCTTCGCCAGCTTCTCCAGCATGTCCAGCACACGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGAtt                                                                                                                           >  1:2033523/1‑150 (MQ=255)                    gCCGGACCCTCTTTCAGCACCTGACCGAGGTTTTTCCAGAAGCCCTTGTATTTCTCGGGTTCGTTCTTCGCCAGCTTCTCCAGCATGTCCAGCACACGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGAtt                                                                                                                           >  2:1843707/1‑150 (MQ=255)                       gggCCCTCTTTCAGCACCTGACCGAAGTTTTTCCAGAAGCCCTTGTATTTCTCGGGTTCGTTCTTCGCCAGCTTCTCCAGCATGTCCAGCACACGCTTGGTCAGCGCGGTCTTCATCGATTCgatgaaga                                                                                                                                            <  1:1796520/130‑4 (MQ=255)                       gggCCCTCTTTCAGCACCTGACCGAAGTTTTTCCAGAAGCCCTTGTATTTCTCGGGTTCGTTCTTCGCCAGCTTCTCCAGCATGTCCAGCACACGCTTGGTCAGCGCGGTCTTCATCGATTCgatgaaga                                                                                                                                            >  2:1796520/1‑127 (MQ=255)                       gggCCCTCGTGCAGCGCCTGACCGAAGTTTTTCCAGAAGCCCTTGTATTTCTCGGGTTCGTTCTTCGCCAGCTTCTCCAGCATGTCCAGCACACGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGATTTCa                                                                                                                        <  2:786031/150‑1 (MQ=255)                             tctTTCAGGACCTGACCGAGGGTTTGCCAGAAGCCCTTGTATTTCTCGGGTTCGTTCTTCGCCAGCTGCTCCAGCATGTCCAGCACACGTTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGATTTCACGCGac                                                                                                                  <  2:980833/150‑1 (MQ=255)                             tctTTCAGCACCTGACCGAAGTTTTTCCAGAAGCCCTTGTATTTCTCGGGTTCGTTCTTCGCCAGCTTCTCCAGCATGTCCAGCACACGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGa                                                                                                                             >  2:544866/1‑139 (MQ=255)                             tctTTCAGCACCTGACCGAAGTTTTTCCAGAAGCCCTTGTATTTCTCGGGTTCGTTCTTCGCCAGCTTCTCCAGCATGTCCAGCACACGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGa                                                                                                                             <  1:544866/139‑1 (MQ=255)                                ttCAGCACCTGACCGAAGTTTTTCCAGAAGCCCTTGTATTTCTCGGGTTCGTTCTTCGCCAGCTTCTCCAGCATGTCCAGCACACGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGATTTCACGCGACACg                                                                                                               >  2:752789/1‑150 (MQ=255)                                   aGCACCTGACCGAAGTTTTTCCAGAAGCCCTTGTATTTCTCGGGTTCGTTCTTCGCCAGCTTCTCCAGCATGTCCAGCACACGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTc                                                                                                            >  1:1973267/1‑150 (MQ=255)                                    gCACCTGACCGAAGTTTTTCCAGAAGCCCTTGTATTTCTCGGGTTCGTTCTTCGCCAGCTTCTCCAGCATGTCCAGCACACGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCa                                                                                                           >  1:1339060/1‑150 (MQ=255)                                       ccTGACCGAAGTTTTTCCAGAAGCCCTTGTATTTCTCGGGTTCGTTCTTCGCCAGCTTCTCCAGCATGTCCAGCACACGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCAGCg                                                                                                        >  1:2345629/1‑150 (MQ=255)                                        cTGACCGAAGTTTTTCCAGAAGCCCTTGTATTTCTCGGGTTCGTTCTTCGCCAGCTTCTCCAGCATGTCCAGCACACGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCgg                                                                                                                                           <  1:71838/114‑1 (MQ=255)                                        cTGACCGAAGTTTTTCCAGAAGCCCTTGTATTTCTCGGGTTCGTTCTTCGCCAGCTTCTCCAGCATGTCCAGCACACGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCgg                                                                                                                                           >  2:71838/1‑114 (MQ=255)                                          gACCGAAGTTTTTCCAGAAGCCCTTGTATTTCTCGGGTTCGTTCTTCGCCAGCTTCTCCAGCATGTCCAGCACACGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCAGCGACa                                                                                                     >  1:1134533/1‑150 (MQ=255)                                             cGAAGTTTTTCCAGAAGCCCTTGTATTTCTCGGGTTCGTTCTTCGCCAGCTTCTCCAGCATGTCCAGCACACGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGATTTCAc                                                                                                                       <  2:1903292/129‑1 (MQ=255)                                             cGAAGTTTTTCCAGAAGCCCTTGTATTTCTCGGGTTCGTTCTTCGCCAGCTTCTCCAGCATGTCCAGCACACGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGATTTCAc                                                                                                                       >  1:1903292/1‑129 (MQ=255)                                             cGAAGTTTTTCCAGAAGCCCTTGTATTTCTCGGGTTCGTTCTTCGCCAGCTTCTCCAGCATGTCCAGCACACGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCAGCGACAGGt                                                                                                  >  2:577700/1‑150 (MQ=255)                                             cGAAGTTTTTCCAGAAGCCCTTGTATTTCTCGGGTTCGTTCTTCGCCAGCTTCTCCAGCATGTCCAGCACACGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCAGCGACAGGt                                                                                                  >  2:1981631/1‑150 (MQ=255)                                              gAAGTTTTTCCAGAAGCCCTTGTATTTCTCGGGTTCGTTCTTCGCCAGCTTCTCCAGCATGTCCAGCACACGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCAGCGACAGGTc                                                                                                 >  1:904861/1‑150 (MQ=255)                                                 gTTTTTCCAGAAGCCCTTGTATTTCTCGGGTTCGTTCTTCGCCAGCTTCTCCAGCATGTCCAGCACACGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCAGCGACAGGTCAtt                                                                                              <  2:2345629/150‑1 (MQ=255)                                                  tttttCCAGAAGCCCTTGTATTTCTCGGGTTCGTTCTTCGCCAGCTTCTCCAGCATGTCCAGCACACGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCAGCGACAGGTCATTg                                                                                             <  1:1843707/150‑1 (MQ=255)                                                     ttCCAGAAGCCCTTGTATTTCTCGGGTTCGTTCTTCGCCAGCTTCTCCAGCATGTCCAGCACACGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCAGCGACAGGTCATTGGaa                                                                                          >  2:70514/1‑150 (MQ=255)                                                     ttCCAGAAGCCCTTGTATTTCTCGGGTTCGTTCTTCGCCAGCTTCTCCAGCATGTCCAGCACACGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCAGCGACAGGTCATTGGaa                                                                                          <  2:2102242/150‑1 (MQ=255)                                                         aGAAGCCCTTGTATTTCTCGGGTTCGTTCTTCGCCAGCTTCTCCAGCATGTCCAGCACACGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCAGCGACAGGTCATTGGAATCaa                                                                                      >  1:1409709/1‑150 (MQ=255)                                                           aaGCCCTTGTATTTCTCGGGTTCGTTCTTCGCCAGCTTCTCCAGCATGTCCAGCACACGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCAGCGACAGGTCATTGGAATCAAcc                                                                                    >  2:1987939/1‑150 (MQ=255)                                                                 ttGTATTTCTCGGGTTCGTTCTTCGCCAGCTTCTCCAGCATGTCCAGCACACGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCAGCGACAGGTCATTGGAATCAACCACAccc                                                                              >  2:2216052/1‑150 (MQ=255)                                                                 ttGTATTTCTCGGGTTCGTTCTTCGCCAGCTTCTCCAGCATGTCCAGCACACGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCAGCGACAGGTCATTGGAATCAACCACAccc                                                                              >  1:1567845/1‑150 (MQ=255)                                                                   gTATTTCTCGGGTTCGTTCTTCGCCAGCTTCTCCAGCATGTCCAGCACACGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCAGCGACAGGTCATTGGAATCAACCACACCCtt                                                                            <  1:1942645/150‑1 (MQ=255)                                                                   gTATTTCTCGGGTTCGTTCTTCGCCAGCTTCTCCAGCATGTCCAGCACACGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCAGCGACAGGTCATTGGAATCAACCACACCCtt                                                                            <  2:904861/150‑1 (MQ=255)                                                                      tttCTCGGGTTCGTTCTTCGCCAGCTTCTCCAGCATGTCCAGCACACGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCAGCGACAGGTCATTGGAATCAACCACACCCTtgat                                                                         <  2:637538/150‑1 (MQ=255)                                                                         ctcGGGTTCGTTCTTCGCCAGCTTCTCCAGCATGTCCAGCACACGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCAGCGACAGGTCATTGGaa                                                                                          >  1:1960386/1‑130 (MQ=255)                                                                         ctcGGGTTCGTTCTTCGCCAGCTTCTCCAGCATGTCCAGCACACGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCAGCGACAGGTCATTGGaa                                                                                          <  2:1960386/130‑1 (MQ=255)                                                                         ctcGGGTTCGTTCTTCGCCAGCTTCTCCAGCATGTCCAGCACACGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCAGCGACAGGTCATTGGAATCAACCACACCCTTGATGaa                                                                      <  2:1339060/150‑1 (MQ=255)                                                                           cGGGTTCGTTCTTCGCCAGCTTTTCCTGCATGTCCAGCACACGCTTTGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCAGCGACAGGTCATTGGAATCAACCACACCCTTGATGAAgc                                                                    <  2:921969/150‑1 (MQ=255)                                                                               ttcgttcTTCGCCAGCTTCTCCAGCATGTCCAGCACACGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCAGCGACAGGTCATTGGAATCAACCACACCCTTGATGAAGCGCAg                                                                <  2:184166/150‑1 (MQ=255)                                                                                   ttcttcGCCAGCTTCTCCAGCATGTCCAGCACACGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCAGCGACAGGTCATTGGAATCAACCACACCCTTGATGAAGCGCAGGTAc                                                            >  2:559905/1‑150 (MQ=255)                                                                                   ttcttcGCCAGCTTCTCCAGCATGTCCAGCACACGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCAGCGACAGGTCATTGGAATCAACCACACCCTTGATGAAGCGCAGGTAc                                                            >  2:1107013/1‑150 (MQ=255)                                                                                      ttcGCCAGCTTCTCCAGCATGTCCAGCACACGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCAGCGACAGGTCATTGGAATCAACCACACCCTTGATGAAGCGCAGGTACAGc                                                         >  1:1849106/1‑150 (MQ=255)                                                                                      ttcGCCAGCTTCTCCAGCATGTCCAGCACACGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCAGCGACAGGTCATTGGAATCAACCACACCCTTGATGAAGCGCAGGTACAGc                                                         >  1:1474581/1‑150 (MQ=255)                                                                                      ttcGCCAGCTTCTCCAGCATGTCCAGCACACGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTAAGCGACAGGTCATTGGAATCAACCACGCCCTTGATGATGCGCCGGTAGAGc                                                         <  2:1900376/150‑1 (MQ=255)                                                                                          ccAGCTTCTCCAGCATGTCCAGCACACGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCAGCGACAGGTCATTGGAATCAACCACACCCTTGATGAAGCGCAGGTACAGCGGCa                                                     >  2:992658/1‑150 (MQ=255)                                                                                          ccAGCTTCTCCAGCATGTCCAGCACACGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCAGCGACAGGTCATTGGAATCAACCACACCCTTGATGAAGCGCAGGTACAGCGGCa                                                     <  1:1107013/150‑1 (MQ=255)                                                                                           cAGCTTCTCCAGCATGTCCAGCACACGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCAGCGACAGGTCATTGGAATCAACCACACCCTTGATGAAGCGCAGGTACAGCGGCAg                                                    >  2:2270545/1‑150 (MQ=255)                                                                                            aGCTTCTCCAGCATGTCCAGCACCCGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGTGTCTTTCTGCAGGATTTCACGCGACACGTTCAGCGACAGGTCATTGGAATCAACCACACCCTTGATGAAGCGCAGGTACAGCGGCAgg                                                   <  1:10611/150‑1 (MQ=255)                                                                                            aGCTTCTCCAGCATGTCCAGCACACGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGTATTTCACGCGACACGTTCAGCGACAGGTCATTGGAATCAACCACACCCTTGATGAAGCGCAGGTACAGCGGc                                                      <  1:773218/147‑1 (MQ=255)                                                                                            aGCTTCTCCAGCATGTCCAGCACACGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCAGCGACAGGTCATTGGAATCAACCACACCCTTGATGAAGCGCAGGTACAGCGGc                                                      >  2:773218/1‑147 (MQ=255)                                                                                              cTTCTCCAGCATGTCCAGCACACGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCAGCGACAGGTCATTGGAATCAACCACACCCTTGATGAAGCGCAGGTACAGCGGCAGGaa                                                 <  2:1973267/150‑1 (MQ=255)                                                                                              cTTCTCCAGCATGTCCAGCACACGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCAGCGACAGGTCATTGGAATCAACCACACCCTTGATGAAGCGCAGGTACAGCGGCAGGaa                                                 <  1:1881046/150‑1 (MQ=255)                                                                                              cTTCTCCAGCATGTCCAGCACACGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCAGCGACAGGTCATTGGAATCAACCACACCCTTGATGAAGCGCAGGTACAGCGGCAGGaa                                                 <  1:2216052/150‑1 (MQ=255)                                                                                                  tcCAGCATGTCCAGCACACGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCAGc                                                                                                         >  2:1346998/1‑90 (MQ=255)                                                                                                  tcCAGCATGTCCAGCACACGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCAGc                                                                                                         <  1:1346998/90‑1 (MQ=255)                                                                                                    cAGCATGTCCAGCACACGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCAGCGACAGGTCATTGGAATCAACCACACCCTTGATGAAGCGCAGGTACAGCGGCAGGAACGACTc                                           <  1:1981631/150‑1 (MQ=255)                                                                                                      gCATGTCCAGCACACGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCAGCGACAGGTCATTGGAATCAACCACACCCTTGATGAAGCGCAGGTACAGCGGCAGGAACGACTCgg                                         >  2:154629/1‑150 (MQ=255)                                                                                                      gCATGTCCAGCACACGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCAGCGACAGGTCATTGGAATCAACCACACCCTTGATGAAGCGCAGGTACAGCGGCAGGAACGACTCgg                                         <  1:2270545/150‑1 (MQ=255)                                                                                                       cATGTCCAGCACACGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCAGCGACAGGTCATTGGAATCAACCACACCCTTGATGAAGCGCAGGTACAGCGGCAGGAACGACTCGGc                                        <  1:490413/150‑1 (MQ=255)                                                                                                           tCCAGCACACGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCAGCGACAGGTCATTGGAATCAACCACACCCTTGATGAAGCGCAGGTACAGCGGCAGGAACGACTCGGCCTgg                                    <  2:2033523/150‑1 (MQ=255)                                                                                                           tCCAGCACACGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCAGCGACAGGTCATTGGAATCAACCACACCCTTGATGAAGCGCAGGTACAGCGGCAGGAACGACTCGGCCTgg                                    >  1:1118679/1‑150 (MQ=255)                                                                                                           tCCAGCACACGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCAGCGACAGGTCATTGGAATCAACCACACCCTTGATGAAGCGCAGGTACAGCGGCAGGAACGACTCGGCCTgg                                    >  1:252481/1‑150 (MQ=255)                                                                                                            ccAGCACACGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCAGCGACAGGTCATTGGAATCAACCACACCCTTGATGAAGCGCAGGTACAg                                                          >  2:1014292/1‑127 (MQ=255)                                                                                                            ccAGCACACGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCAGCGACAGGTCATTGGAATCAACCACACCCTTGATGAAGCGCAGGTACAg                                                          <  1:1014292/127‑1 (MQ=255)                                                                                                             cAGCACACGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCAGCGACAGGTCATTGGAATCAACCACACCCTTGATGAAGCGCAGGTACAGCGGCAGGAACGACTCGGCCTGGTc                                  >  2:1820387/1‑150 (MQ=255)                                                                                                             cAGCACACGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCAGCGACAGGTCATTGGAATCAACCACACCCTTGATGAAGCGCAGGTACAGCGGCAGGAACGACTCGGCCTGGTc                                  <  2:1474581/150‑1 (MQ=255)                                                                                                               gCACACGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCAGCGACAGGTCATTGGAATCaa                                                                                      <  2:1355164/96‑1 (MQ=255)                                                                                                               gCACACGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCAGCGACAGGTCATTGGAATCaa                                                                                      >  1:1355164/1‑96 (MQ=255)                                                                                                                cacaCGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCAGCGACAGGTCATTGGAATCAACCACACCCTTGATGAAGCGCAGGTACAGCGGCAGGAACGACTCGGCCTGGTCCat                               <  1:577700/150‑1 (MQ=255)                                                                                                                cacaCGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCAGCGACAGGTCATTGGAATCAACCACACCCTTGATGAAGCGCAGGTACAGCGGCAGGAACGACTCGGCCTGGTCCat                               <  1:992658/150‑1 (MQ=255)                                                                                                                  cacGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCAGCGACAGGTCATTGGAATCAACCACACCCTTGATGAAGCGCAGGTACAGCGGCAGGAACGACTCGGCCTGGTCCatga                             <  2:1134533/150‑1 (MQ=255)                                                                                                                  cacGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCAGCGACAGGTCATTGGAATCAACCACACCCTTGATGAAGCGCAGGTACAGCGGCAGGAACGACTCGGCCTGGTCCatga                             <  1:1820387/150‑1 (MQ=255)                                                                                                                     gCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCAGCGACAGGTCATTGGAATCAACCACACCCTTGATGAAGCGCAGGTACAGCGGCAGGAACGACTCGGCCTGGTCCAtgatga                          >  2:2227539/1‑150 (MQ=255)                                                                                                                       ttGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCAGCGACAGGTCATTGGAATCAACCACACCCTTGATGAAGCGCAGGTACAGCGGCAGGAACGACTCGGCCTGGTCCATGATGaaa                        >  1:1843060/1‑149 (MQ=255)                                                                                                                          gTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCAGCGACAGGTCATTGGAATCAACCACACCCTTGATGAAGCGCAGGTACAGCGGCAGGAACGACTCGGCCTGGTCCATGATGAacaca                     >  1:50073/1‑150 (MQ=255)                                                                                                                          gTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCAGCGACAGGTCATTGGAATCAACCACACCCTTGATGAAGCGCAGGTACAGCGGCAGGAACGACTCGGCCTGGTCCATGATGAacaca                     <  1:2227539/150‑1 (MQ=255)                                                                                                                          gTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCAGCGACAGGTCATTGGAATCAACCACACCCTTGATGAAGCGCAGGTACAGCGGCAGGAACGACTCGGCCTGGTCCATGATGAacaca                     >  2:482287/1‑150 (MQ=255)                                                                                                                           tCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCAGCGACAGGTCATTGGAATCAACCACACCCTTGATGAAGCGCAGGTACAGCGGCAGGAACGACTCGGcc                                       >  2:179620/1‑131 (MQ=255)                                                                                                                           tCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCAGCGACAGGTCATTGGAATCAACCACACCCTTGATGAAGCGCAGGTACAGCGGCAGGAACGACTCGGcc                                       <  1:179620/131‑1 (MQ=255)                                                                                                                                 acgcggtCATCGATTCGAGGACCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCCGCGACAGGTCATTGGAATCAACCCCACCCTTGATCAAGCGCAGGTACAGCGGCAGGAACCACTCGGCCTGGGCCCTTATGAACACCCGCTGcc              >  2:1196226/7‑149 (MQ=255)                                                                                                                                   ccgatCATCGATTCGATGATCGGGTCTTTATGCAGGATTTCACGCGACACGTTCAGCGACAGGTCATTGGAATCAACCACACCCTTGATGAAGCGCAGGTACAGCGGCAGGAACGACTCGGCCTGGTCCATGATGAACACACGCTGCACg            >  1:317819/5‑150 (MQ=255)                                                                                                                                      atCATCGATTCGATGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCAGCGACAGGTCATTGGAATCAACCACACCCTTGATGAAGCGCAGGTACAGCGGCAGGAACGACTCGGCCTGGTCCATGATGAACACACGCTGCACGTAc         <  2:317819/149‑1 (MQ=255)                                                                                                                                        cATCGAATCGATGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCAGCGACAGGTCATTGGAATCAACCACACCCTTGATGAAGCGCAGGTACAGCGGCa                                                     >  1:409510/1‑104 (MQ=255)                                                                                                                                        cATCGAATCGATGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCAGCGACAGGTCATTGGAATCAACCACACCCTTGATGAAGCGCAGGTACAGCGGCa                                                     <  2:409510/104‑1 (MQ=255)                                                                                                                                          tcgaatcgaTGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCAGCGACAGGTCATTGGAATCAACCACACCCTTGATGAAGCGCAGGTACAGCGGCAGGAACGACTCGGCCTGGTCCATGATGAACACACGCTGCACGTACAGCt     >  2:1320448/1‑150 (MQ=255)                                                                                                                                          tcgaatcgaTGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCAGCGACAGGTCATTGGAATCAACCACACCCTTGATGAAGCGCAGGTACAGCGGCAGGAACGACTCGGCCTGGTCCATGATGAACACACGCTGCACGTACAGCt     >  1:763568/1‑150 (MQ=255)                                                                                                                                            gaatcgaTGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCAGCGACAGGTCATTGGAATCAACCACACCCTTGATGAAGCGCAGGTACAGCGGCAGGAACGACTCGGCCTGGTCCa                                <  1:1097556/121‑1 (MQ=255)                                                                                                                                            gaatcgaTGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCAGCGACAGGTCATTGGAATCAACCACACCCTTGATGAAGCGCAGGTACAGCGGCAGGAACGACTCGGCCTGGTCCa                                >  2:1097556/1‑121 (MQ=255)                                                                                                                                            gaatcgaTGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCAGCGACAGGTCATTGGAATCAACCACACCCTTGATGAAGCGCAGGTACAGCGGCAGGAACGACTCGGCCTGGTCCATGATGAACACACGCTGCACGTACAGCTTc   <  1:1320448/150‑1 (MQ=255)                                                                                                                                             aatcgaTGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCAGCGACAGGTCATTGGAATCAACCACACCCTTGATGAAGCGCAGGTACAGCGGCAGGAACGACTCGGCCTGGTCCATGATGAACACACGCTGCACGTACAGCTTCa  <  1:1987939/150‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                     CTTGTTGGCGAAGTCTTCGGCCGGGCCCTCTTTCAGCACCTGACCGAAGTTTTTCCAGAAGCCCTTGTATTTCTCGGGTTCGTTCTTCGCCAGCTTCTCCAGCATGTCCAGCACACGCTTGGTCAGCGCGGTCTTCATCGAATCGATGATCGGGTCTTTCTGCAGGATTTCACGCGACACGTTCAGCGACAGGTCATTGGAATCAACCACACCCTTGATGAAGCGCAGGTACAGCGGCAGGAACGACTCGGCCTGGTCCATGATGAACACACGCTGCACGTACAGCTTCA  >  NC_002947/4722413‑4722702

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