Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,491,905	A→T	5.9%	intergenic (+55/+17)	ydcA → / ← hokB	protein YdcA/toxin HokB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,491,905	0	A	T	5.9%	262.4 / 3.7	83	intergenic (+55/+17)	ydcA/hokB	protein YdcA/toxin HokB
Reads supporting (aligned to +/- strand):  ref base A (34/44);  new base T (3/2);  total (37/46)
Fisher's exact test for biased strand distribution p-value = 6.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

GTTCGAAAGGGGGGATCTCACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATGCAGATAATGAGCAGACACAC  >  NC_000913/1491756‑1492050                                                                                                                                                      |                                                                                                                                                  gTTCGAAAGGGGGGATCTCACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGaaa                                                                                                                                                   >  1:964293/1‑150 (MQ=255)  ttCGAAAGGGGGGATCTCACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGaaaa                                                                                                                                                  >  2:1070122/1‑150 (MQ=255)           ggggATCTCACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGTTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGATTTTCCCCgcgc                                                                                                                                         <  1:519112/150‑1 (MQ=255)            gggATCTCACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCt                                                                                                                                        <  1:218326/150‑1 (MQ=255)              gATCTCACACTGTACGGCAGGCGGCAAATTTGTCTGTAATTATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTg                                                                                                                                      >  2:994205/1‑150 (MQ=255)              gATCTCACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTg                                                                                                                                      <  1:778664/150‑1 (MQ=255)                tctcACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAGGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGATTTTCCCCgcgcatttg                                                                                                                                    <  1:199889/150‑6 (MQ=255)                  tcACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCt                                                                                                                                  >  1:1273807/1‑150 (MQ=255)                      acTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGGGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTAcc                                                                                                                              <  1:60711/150‑1 (MQ=255)                       cTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCt                                                                                                                             >  2:1269120/1‑150 (MQ=255)                         gtacttcAGGCTGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCTTGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTgg                                                                                                                           <  1:610013/144‑1 (MQ=255)                             ggCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGt                                                                                                                       >  2:403833/1‑150 (MQ=255)                                 ggcggcAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAg                                                                                                                   >  1:665022/1‑150 (MQ=255)                                   cggcAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGc                                                                                                                 >  1:73837/1‑150 (MQ=255)                                   cggcAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGc                                                                                                                 <  1:211777/150‑1 (MQ=255)                                    ggcAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGcc                                                                                                                >  1:807129/1‑150 (MQ=255)                                       aaaTTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATg                                                                                                             <  2:964293/150‑1 (MQ=255)                                          tttGTCTGTAATGATGGTTCTATTAGTGAATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATATGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAgc                                                                                                          >  2:213172/1‑150 (MQ=255)                                           ttGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAgcg                                                                                                         <  1:177367/150‑1 (MQ=255)                                             gTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCgca                                                                                                       <  2:1049363/150‑1 (MQ=255)                                              tCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCgcag                                                                                                      >  2:282901/1‑150 (MQ=255)                                               cTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCgcagc                                                                                                     <  2:199771/150‑1 (MQ=255)                                                 gTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGc                                                                                                                 <  1:279669/136‑1 (MQ=255)                                                 gTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGc                                                                                                                 >  2:279669/1‑136 (MQ=255)                                                  tAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAAc                                                                                                  >  1:903029/1‑150 (MQ=255)                                                     tgatgGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAAcctc                                                                                               >  2:573173/1‑150 (MQ=255)                                                     tgatgGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAAcctc                                                                                               >  1:1245193/1‑150 (MQ=255)                                                     tgatgGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAAcctc                                                                                               >  1:804540/1‑150 (MQ=255)                                                      gatgGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTAcc                                                                                                                              >  1:1363936/1‑118 (MQ=255)                                                      gatgGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTAcc                                                                                                                              <  2:1363936/118‑1 (MQ=255)                                                       ttgGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAAcctcct                                                                                             <  1:713940/149‑1 (MQ=255)                                                       atgGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGGGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAAcctcct                                                                                             <  1:369134/150‑1 (MQ=255)                                                         ggTTCTATTTGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTa                                                                                           <  1:1394196/150‑1 (MQ=255)                                                         ggTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTa                                                                                           >  2:232189/1‑150 (MQ=255)                                                               aTTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCg                                                                                     >  1:268321/1‑150 (MQ=255)                                                                    tGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCg                                                                                <  1:1070122/150‑1 (MQ=255)                                                                     gCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCgg                                                                               <  2:665022/150‑1 (MQ=255)                                                                     gCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCgg                                                                               >  2:684270/1‑150 (MQ=255)                                                                        tCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAAc                                                                            >  2:92293/1‑150 (MQ=255)                                                                         cGAAAAAAACATGCACTAACTGAAGGGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAAcc                                                                           <  1:1282404/150‑1 (MQ=255)                                                                          gAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCg                                                                          >  2:58398/1‑150 (MQ=255)                                                                              aaaaCATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGt                                                                      <  1:1205234/150‑1 (MQ=255)                                                                               aaaCATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGtt                                                                     <  1:213172/150‑1 (MQ=255)                                                                               aaaCATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGtt                                                                     >  2:668991/1‑150 (MQ=255)                                                                                 aCATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAAGCTGCTACTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCg                                                                   <  2:256658/150‑1 (MQ=255)                                                                                  cATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGtt                                                                     <  2:624647/147‑1 (MQ=255)                                                                                  cATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGtt                                                                     >  1:624647/1‑147 (MQ=255)                                                                                    tGCACTAACTGAAGTGTAAAAGGGGGGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTag                                                                <  1:931969/150‑1 (MQ=255)                                                                                       aCTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGc                                                             >  2:621477/1‑150 (MQ=255)                                                                                            cTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAg                                                                       <  1:1177786/135‑1 (MQ=255)                                                                                            cTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAg                                                                       >  2:1177786/1‑135 (MQ=255)                                                                                            cTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTACTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTg                                                        >  1:1292444/1‑150 (MQ=255)                                                                                             tGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGt                                                       <  2:848817/150‑1 (MQ=255)                                                                                             tGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGt                                                       <  2:843487/150‑1 (MQ=255)                                                                                              gAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGTGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGATCGCAGCAACCTCCTTTTCACCGTCCCGTAACCGCAGTTCGTAGAGCGTTTGTc                                                      <  2:198824/150‑1 (MQ=255)                                                                                               aaGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCg                                                     <  1:1002784/150‑1 (MQ=255)                                                                                                   gtAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTc                                                 >  2:574650/1‑150 (MQ=255)                                                                                                   gtAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTc                                                 <  1:574650/150‑1 (MQ=255)                                                                                                   gtAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTc                                                 >  1:396635/1‑150 (MQ=255)                                                                                                    tAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCa                                                <  1:573173/150‑1 (MQ=255)                                                                                                     aaaaGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAg                                               <  1:1269120/150‑1 (MQ=255)                                                                                                            gTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGtgtg                                        <  2:1152762/150‑1 (MQ=255)                                                                                                                cATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAAtt                                    <  2:1188514/150‑2 (MQ=255)                                                                                                                cATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATg                                    <  1:58398/150‑1 (MQ=255)                                                                                                                cATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATg                                    >  2:1221179/1‑150 (MQ=255)                                                                                                                        aaTGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATc                            <  2:251594/150‑1 (MQ=255)                                                                                                                            aaTCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTaa                        <  1:1166746/150‑1 (MQ=255)                                                                                                                             aTCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAAt                       >  1:979019/1‑150 (MQ=255)                                                                                                                                 gCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATGCAg                   >  2:1102873/1‑150 (MQ=255)                                                                                                                                   tAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATGCAGAt                 <  1:305540/150‑1 (MQ=255)                                                                                                                                   tAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATGCAGAt                 <  1:1221179/150‑1 (MQ=255)                                                                                                                                     agggcaagcgcGGGGATTTTCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATGCAGATaa               >  2:283956/8‑150 (MQ=255)                                                                                                                                       cAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATGCAGATAATg             <  1:621477/150‑1 (MQ=255)                                                                                                                                       cAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATGCAGATAATg             >  1:33068/1‑150 (MQ=255)                                                                                                                                       cAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATGCAGATAATg             >  1:322104/1‑150 (MQ=255)                                                                                                                                        aaaTGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGTGTCAGGAGTGTGAATGTCAGAATCGTAATGCAGATAATGa            <  1:45536/150‑1 (MQ=255)                                                                                                                                        aaaTGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATGCAGATAATGa            <  2:60480/150‑1 (MQ=255)                                                                                                                                        aaaTGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATGCAGATAATGa            >  2:736391/1‑150 (MQ=255)                                                                                                                                            gcgcGGGGATTTTCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATGCAGATAATGAGCAg        >  1:434755/1‑150 (MQ=255)                                                                                                                                             cgcgGGGATTTTCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATGCAGATAATGAGCAGa       >  2:1083856/1‑150 (MQ=255)                                                                                                                                             cgcgGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATGCAGATAATGAGCAGa       <  2:73837/150‑1 (MQ=255)                                                                                                                                             cgcgGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATGCAGATAATGAGCAGa       <  1:684270/150‑1 (MQ=255)                                                                                                                                              gcgGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATGCAGATAATGAGCAGac      <  1:994205/150‑1 (MQ=255)                                                                                                                                                ggggAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATGCAGATAATGAGCAGacac    <  1:282901/150‑1 (MQ=255)                                                                                                                                                  ggAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATGCAGATAATGAGCAGacacac  >  1:894064/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                  GTTCGAAAGGGGGGATCTCACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATGCAGATAATGAGCAGACACAC  >  NC_000913/1491756‑1492050

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