Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,926,900	T→G	6.1%	D647A (GAC→GCC)	ptrB ←	protease II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,926,900	0	T	G	6.1%	269.0 / 8.3	99	D647A (GAC→GCC)	ptrB	protease II
Reads supporting (aligned to +/- strand):  ref base T (41/52);  new base G (3/3);  total (44/55)
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.38e-01

TAACACTGAAACATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTA‑GCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCAT  >  NC_000913/1926751‑1927047                                                                                                                                                      |                                                                                                                                                     tAACACTGAAACATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGt                                                                                                                                                      >  1:550511/1‑150 (MQ=255)    cacTGAAACATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCgg                                                                                                                                                   <  1:769722/150‑1 (MQ=255)    cacTGAAACATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCgg                                                                                                                                                   >  1:1505154/1‑150 (MQ=255)      cTGAAACATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTa                                                                                                                                                 >  1:1022810/1‑150 (MQ=255)      cTGAAACATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTa                                                                                                                                                 >  1:1486721/1‑150 (MQ=255)      cTGAAACATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGACATGGACTCAGGCCATGGCGGTa                                                                                                                                                 <  1:727951/150‑1 (MQ=255)      cTGAAACATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGAGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTa                                                                                                                                                 >  1:1614561/1‑150 (MQ=255)       tGAAACATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTAc                                                                                                                                                <  2:1556256/150‑1 (MQ=255)        gAAACATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACa                                                                                                                                               >  1:338518/1‑150 (MQ=255)           aCATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGACATGGACTCAGGCCATGGCGGTaaatct                                                                                                                                            <  2:1583257/150‑6 (MQ=255)           aCATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGACATGGACTCAGGCCATGGCGGTaaatct                                                                                                                                            <  2:144015/150‑6 (MQ=255)            cATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGc                                                                                                                                           >  2:1101633/1‑150 (MQ=255)               tATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAAt                                                                                                                                        <  1:1376094/150‑1 (MQ=255)               tATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAACATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAAt                                                                                                                                        >  2:550615/1‑150 (MQ=255)                 tCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGCTTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATaa                                                                                                                                      >  1:172393/1‑150 (MQ=255)                   tGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATaaaa                                                                                                                                    <  2:1290769/150‑1 (MQ=255)                    ggAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAg                                                                                                                                   <  1:663643/150‑1 (MQ=255)                     gAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGa                                                                                                                                  >  2:443750/1‑150 (MQ=255)                       aaaTACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATg                                                                                                                                >  1:819101/1‑150 (MQ=255)                       aaaTACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATg                                                                                                                                <  1:1639604/150‑1 (MQ=255)                            cTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCAt                                                                                                                           >  2:912286/1‑150 (MQ=255)                               aGTCCGCAGGCGGAGCGGGTAATGTTCACTGCGCCCGCGCGCCCAGAAAAGCATATTCCATAGCAACCCCTGCGTACGATTTAAAGCGACCAGATTTACCGCCAAAGCCGGAGGCCATGTAGGTTCAGAGCAAGAAAACATGTTCATCgg                                                                                                                        >  2:1656725/1‑150 (MQ=255)                                    gCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTc                                                                                                                   <  2:1396287/150‑1 (MQ=255)                                     cAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATGGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCa                                                                                                                  <  1:239514/150‑1 (MQ=255)                                     cAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCa                                                                                                                  >  2:1677272/1‑150 (MQ=255)                                     cAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCATTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCa                                                                                                                  <  1:1440400/150‑1 (MQ=255)                                       ggCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGACATAGGTTTTCAGc                                                                                                                >  2:262355/1‑150 (MQ=255)                                            aGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTcgcg                                                                                                           <  2:771792/150‑1 (MQ=255)                                            aGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTcgcg                                                                                                           <  1:1333538/150‑1 (MQ=255)                                            aGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTcgcg                                                                                                           <  2:874156/150‑1 (MQ=255)                                               gggTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCaa                                                                                                        >  2:1344794/1‑150 (MQ=255)                                                 gTAATGCTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAAtt                                                                                                      <  2:209372/150‑1 (MQ=255)                                                  tAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAAttt                                                                                                     <  2:1234211/150‑1 (MQ=255)                                                    aTGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCCTGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTA‑g                                                                                                  <  2:301109/150‑1 (MQ=255)                                                     tGTTGCCTTCGCCACCGGGAACAGAAAAGCATATTCCCTCGCTACGCCGTCGTTCGATTTAAAGCGACCACATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTA‑Gc                                                                                                 <  1:472504/150‑1 (MQ=255)                                                      gTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTA‑GCg                                                                                                >  2:97532/1‑150 (MQ=255)                                                        tCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTA‑GCGAc                                                                                              >  2:1511688/1‑150 (MQ=255)                                                          ccTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATATAAGATGGTCTTCGGTTTTCAGCTCGCGCAATTTA‑GCGAccc                                                                                            >  2:265284/1‑150 (MQ=255)                                                          ccTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTA‑GCGAccc                                                                                            >  2:88580/1‑150 (MQ=255)                                                           cTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCTAAATTGC‑GCGagctg                                                                                           <  2:1100622/150‑5 (MQ=255)                                                               gcCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTA‑GCGACCCAtttt                                                                                       <  1:488136/150‑1 (MQ=255)                                                                ccAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTA‑GCGACCCATTTTg                                                                                      >  1:616360/1‑150 (MQ=255)                                                                  aGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTA‑GCGACCCATTTTGcc                                                                                    <  2:1229903/150‑1 (MQ=255)                                                                     gcgACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTA‑GCGACCCATTTTGCCGGt                                                                                 <  1:1388533/150‑1 (MQ=255)                                                                     gcgACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTA‑GCGACCCATTTTGCCGGt                                                                                 <  2:910807/150‑1 (MQ=255)                                                                         ccAGAAAAGCATATTCCATCGCTACGCCTTAGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTA‑GCGACCCATTTTGCCGGTTccc                                                                             >  1:292732/1‑150 (MQ=255)                                                                              aaaGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTA‑GCGACCCATTTTGCCGGTTCCCAatat                                                                        <  1:1197166/150‑1 (MQ=255)                                                                               aatcATATTCCATCGCTACGCCTTCGTAAGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTA‑GCGACCCATTTTGCCGGTTCCCAATAtt                                                                       <  1:1373880/147‑1 (MQ=255)                                                                                aacATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTA‑GCGACCCATTTTGCCGGTTCCCAATATTg                                                                      <  1:550615/148‑1 (MQ=255)                                                                                   atatTGAATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTA‑GCGACCCATTTTGCCGGTTCCCAATATTGCAc                                                                   <  1:1109091/150‑1 (MQ=255)                                                                                     atTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTA‑GCGACCCATTTTGCCGGTTCCCAATATTGCACCt                                                                 >  1:648746/1‑150 (MQ=255)                                                                                     atTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTA‑GCGACCCATTTTGCCGGTTCCCAATATTGCACCt                                                                 <  1:912286/150‑1 (MQ=255)                                                                                      ttCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTA‑GCGACCCATTTTGCCGGTTCCCAATATTGCACCTg                                                                <  1:1029304/150‑1 (MQ=255)                                                                                         cATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTA‑GCGACCCATTTTGCCGGTTCCCAATATTGCACCTgaga                                                             >  2:632735/1‑150 (MQ=255)                                                                                             gCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTA‑GCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCg                                                         >  2:769028/1‑150 (MQ=255)                                                                                              cTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTA‑GCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGt                                                        >  2:53890/1‑150 (MQ=255)                                                                                                 cGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTA‑GCGACCCAGTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCa                                                     >  2:1536372/1‑150 (MQ=255)                                                                                                    cTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTA‑GCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAAc                                                  <  1:672318/150‑1 (MQ=255)                                                                                                     ttCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTA‑GCGACCCATTTTGCCGGTTCACAATATTGCACCTGAGAATCGTGCAAAcc                                                 <  2:304908/150‑1 (MQ=255)                                                                                                       cGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTA‑GCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCgg                                               >  1:1257203/1‑150 (MQ=255)                                                                                                         tACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTA‑GCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTc                                             <  1:1101633/150‑1 (MQ=255)                                                                                                          aCGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTA‑GCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCg                                            >  2:725763/1‑150 (MQ=255)                                                                                                            gATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTA‑GCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGtt                                          >  1:928792/1‑150 (MQ=255)                                                                                                            gATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTA‑GCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGtt                                          >  1:455180/1‑150 (MQ=255)                                                                                                               ttAAAGCGACCAGATTTACCGCCTTGGCCTGAGTCCATGTCGGTAGAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTA‑GCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTAcc                                       <  2:1395017/150‑1 (MQ=255)                                                                                                               ttAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTA‑GCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTAcc                                       <  1:787292/150‑1 (MQ=255)                                                                                                               ttAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTA‑GCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTAcc                                       <  2:1150089/150‑1 (MQ=255)                                                                                                                tAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGATCGCGCAATTTC‑GCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGGGCAACCCGGGCGTTAccc                                      >  1:340713/1‑149 (MQ=255)                                                                                                                 aaaGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTA‑GCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAg                                     <  2:118463/150‑1 (MQ=255)                                                                                                                 aaaGAGACCTGATTTACCGCCATGGCCGGAGTCCATGTCGGTACAGAGCAATACAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTA‑GCTACCCATTTTGCCGTTTCCCAATATTTCACATGAGAATCGTGCAAACCGGTCGTTAACAg                                     <  2:292615/150‑1 (MQ=255)                                                                                                                  aaGCGACCAGAATTACCGCCATGGCCTGAGTCCATGTCTGTACAGAGCAATAAAAGATGGGCATCGGTTTTCAGCTCTCGCAATTTA‑GCGACCCATTTTGCCTGATCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGt                                    <  1:559336/150‑1 (MQ=255)                                                                                                                   aGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTA‑GCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTa                                   <  1:1087498/150‑1 (MQ=255)                                                                                                                   aGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTA‑GCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTa                                   <  1:1048427/150‑1 (MQ=255)                                                                                                                       aCCAGTTTTACCGCCATTGCCTGAGTCCATGTCGGTACAGAGAAATAAAAGATGGTCTTCGGTTTTCAGCTCGCGCAATTTA‑GCGACCCATTTTGCCGGTTCCCAATATTGCACCTTAGAATCGTGCAAACAGGTCGTTACAAGTAAATg                               <  2:72634/150‑1 (MQ=255)                                                                                                                       aCCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTA‑GCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAATCCTGTCGTTACCAGTAAAtt                               >  1:395938/1‑149 (MQ=255)                                                                                                                       aCCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTA‑GCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATg                               >  2:1679909/1‑150 (MQ=255)                                                                                                                        ccAGATTTACCGCCATGGCCTGCGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTA‑GCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACGGGTCGTTACCAGTAAATGc                              >  1:863969/1‑150 (MQ=255)                                                                                                                          aGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTA‑GCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTTAATGCgg                            >  2:1517086/1‑150 (MQ=255)                                                                                                                           gATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTA‑GCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGa                           >  2:1387726/1‑150 (MQ=255)                                                                                                                             tttACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTA‑GCGACCCATTTTGCCGGTTCCCAATATTGCACCTAAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATa                         <  2:338518/150‑1 (MQ=255)                                                                                                                                aCCGACATGGACTCAGGCCATGGCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTA‑GCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGc                      >  2:675376/1‑150 (MQ=255)                                                                                                                                aCCGACATGGACTCAGGCCATGGCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTA‑GCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGc                      >  1:824566/1‑150 (MQ=255)                                                                                                                                 ccGCCATGGCCTGGGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTA‑GCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGcc                     <  1:434918/150‑1 (MQ=255)                                                                                                                                 ccGACATGGACTCAGGCCATGGCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCATCTCGCGCAATTTA‑GCGACCCATTTTGCCGGTTCCCAATATTGAACATGAGACTCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGcc                     >  2:174743/1‑150 (MQ=255)                                                                                                                                   gCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTA‑GCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCtg                   >  1:614077/1‑150 (MQ=255)                                                                                                                                      ctGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTGTTCAGCTCGCGCAATTTA‑GCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGGCGTTACCAGTAAATGCGGATAAGCCTGGGc                <  2:1411804/149‑1 (MQ=255)                                                                                                                                        ggCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTA‑GCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCgg              >  1:730011/1‑150 (MQ=255)                                                                                                                                         gCCGGAGTCCATGTCGGTTCAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTA‑GCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGt             <  2:1077634/150‑1 (MQ=255)                                                                                                                                          ccTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTA‑GCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTg            <  1:705273/150‑1 (MQ=255)                                                                                                                                          ccTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTA‑GCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTg            >  2:840995/1‑150 (MQ=255)                                                                                                                                          ccTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTATCAGCTCGCGCAATTTA‑GCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTg            >  1:980600/1‑150 (MQ=255)                                                                                                                                             gAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTA‑GCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACg         <  2:1614561/150‑1 (MQ=255)                                                                                                                                                tCCTTGTCGGTACAGAGCAAAAAAAGATGGGCATCGGTTTTCAGCTCGCGCAATTTA‑GCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGGTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTg      <  1:1497865/150‑1 (MQ=255)                                                                                                                                                tCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTA‑GCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTg      <  2:111943/150‑1 (MQ=255)                                                                                                                                                 ccATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTA‑GCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGt     <  1:1472206/150‑1 (MQ=255)                                                                                                                                                  cATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTAAGCTCGCGCAATTTA‑GCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGACGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTc    >  1:626233/1‑150 (MQ=255)                                                                                                                                                  cATGTCGGTAAATTTCAATAAAACGAGGGCCTCCGTTTTCCCCCCCGGGAATTTACGCTACCAATTTT‑CCGGGTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTc    <  2:611798/150‑1 (MQ=255)                                                                                                                                                   aTGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTA‑GCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCa   <  1:1656725/150‑1 (MQ=255)                                                                                                                                                    tGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTA‑GCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCAt  <  2:1486721/150‑1 (MQ=255)                                                                                                                                                    tGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTA‑GCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCAt  <  2:1022810/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                     TAACACTGAAACATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTA‑GCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCAT  >  NC_000913/1926751‑1927047

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