Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,612,525	C→T	5.4%	Q221* (CAG→TAG)	yhhT →	putative transporter YhhT

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,612,525	0	C	T	5.4%	232.1 / 2.1	74	Q221* (CAG→TAG)	yhhT	putative transporter YhhT
Reads supporting (aligned to +/- strand):  ref base C (37/33);  new base T (2/2);  total (39/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.73e-01

TTTTGCGCTGAATAATCCACAGATTCACATCGCGGGATTACACCGCGCACTTAAAGGCGTTTCGCACTATCTTGCATTGAAGACGCTACTCAGTTTATGGACAGGTGTCATCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGTGGGCAGTACTGGCGTTTTTGCTCAACTACGTGCCCAATATCGGCGCGGTAATTTCCGCCGTACCGCCAATGATTCAGGTGCTGCTGTTTAATGGTGTTTACGAATGTATTCTGGTCGGCGCATTGTTTTTAGTG  >  NC_000913/3612383‑3612674                                                                                                                                               |                                                                                                                                                      ttttGCGCTGAATAATCCACAGATTCACATCGCGGGATTACACCGCGCACTTAAAGGCGTTTCGCACTATCTTGCATTGAAGACGCTACTCAGTTTATGGACAGGTGTCATCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTgc                                                                                                                                                >  2:541177/1‑150 (MQ=255)   ttGCGCTGAATAATCCACAGATTCACATCGCGGGATTACACCGCGCACTTAAAGGCGTTTCGCACTATCTTGCATTGAAGACGCTACTCAGTTTATGGACAGGTGTCATCGTCTGGCTGGGGCTGGAGCTGATGCGGGTGCAGTTTgcgc                                                                                                                                              >  2:1169010/1‑150 (MQ=255)     gcgcTGAATAATCCACAGATTCACATCGCGGGATTACACCGCGCACTTAAAGGCGTTTCGCACTATCTTGCATTGAAGACGCTACTCAGTTTATTGACAGGTGTCATCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTg                                                                                                                                            >  1:26908/1‑150 (MQ=255)      cgcTGAATAATCCACAGATTCACATCGCGGGATTACACCGCGCACTTAAAGGCGTTTCGCACTATCTTGCATTGAAGACGCTACTCAGTTTATGGACAGGTGTCATCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGa                                                                                                                                           >  1:1299918/1‑150 (MQ=255)        cTGAATAATCCACAGATTCACATCGCGGGATTACACCGCGCACTTAAAGGCGTTTCGCACTATCTTGCATTGAAGACGCTACTCAGTTTATGGACAGGTGTCATCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGAtg                                                                                                                                         <  1:335559/150‑1 (MQ=255)             taatCCACAGATTCACATCGCGGGATTACACCGCGCACTTAAAGGCGTTTCGCACTATCTTGCATTGAAGACGCTACTCAGTTTATGGACAGGTGTCATCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGTGGGc                                                                                                                                    <  2:840198/150‑1 (MQ=255)              aatCGACAGATTCACATCGCGGGATTACACCGCGCACTTAAAGGCGTTTCGCACTATCTTGCATTGAAGACGCTACTCAGTTTATGGACAGGTGTCATCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGTGGGCa                                                                                                                                   >  1:418136/1‑150 (MQ=255)                tCCACAGATTCACATCGCGGGATTACACCGCGCACTTAAAGGCGTTTCGCACTATCTTGCATTGAAGACGCTACTCAGTTTATGGACAGGTGTCATCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGTGGGCAGt                                                                                                                                 <  2:198315/150‑1 (MQ=255)                  cacaGATTCACATCGCGGGATTACACCGCGCACTTAAAGGCGTTTCGCACTATCTTGCATTGAAGACGCTACTCAGTTTATGGACAGGTGTCATCGTGTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGTGGGCAGTAc                                                                                                                               <  2:4331/150‑1 (MQ=255)                  cacaGATTCACATCGCGGGATTACACCGCGCACTTAAAGGCGTTTCGCACTATCTTGCATTGAAGACGCTACTCAGTTTATGGACAGGTGTCATAGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGTGGGCAGTAc                                                                                                                               >  2:440141/1‑150 (MQ=255)                  cacaGATTCACATCGCGGGATTACACCGCGCACTTAAAGGCGTTTCGCACTATCTTGCATTGAAGACGCTACTCAGTTTATCGACAGGTGTCATCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGTGGGCAGTAc                                                                                                                               >  1:515661/1‑150 (MQ=255)                    caGATTCACATCGCGGGATTACACCGCGCACTTAAAGGCGTTTCGCACTATCTTGCATTGAAGACGCTACTCAGTTTATGGACAGGTGTCATCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGTGGACAGTACTg                                                                                                                             >  1:523707/1‑150 (MQ=255)                      gATTCACATCGCGGGATTACACCGCGCACTTAAAGGCGTTTCGCACTATCTTGCATTGAAGACGCTACTCAGTTTATGGACAGGTGTCATCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGTGGGCAGTACTGGc                                                                                                                           >  1:219384/1‑150 (MQ=255)                       aTTCACATCGCGGGATTACACCGCGCACTTAAAGGCGTTTCGCACTATCTTGCATTGAAGACGCTACTCAGTTTATGGACAGGTGTCATCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGTGGGCAGTACTGGCg                                                                                                                          <  2:700528/150‑1 (MQ=255)                             aTCGCGGGATTACACCGCGCACTTAAAGGCGTTTCGCACTATCTTGCATTGAAGACGCTACTCAGTTTATGGACAGGTGTCATCGTTTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGTGGGCAGTACTGGCGTTTTTg                                                                                                                    >  1:1141145/1‑150 (MQ=255)                             aTCGCGGGATTACACCGCGCACTTAAAGGCGTTTCGCACTATCTTGCATTGAAGACGCTACTCAGTTTATGGACAGGTGTCATCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGTGGGCAGTACTGGCGTTTTTg                                                                                                                    <  1:800842/150‑1 (MQ=255)                              tCGCGGGATTACACCGCGCACTTAAAGGCGTTTCGCACTATCTTGCATTGAAGACGCTACTCAGTTTATGGACAGGTGTCATCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGTGGGCAGTACTGGCGTTTTTGc                                                                                                                   <  1:1171002/150‑1 (MQ=255)                                gcgGGATTACACCGCGCACTTAAAGGCGTTTCGCACTATCTTGCATTGAAGACGCTACTCAGTTTATGGACAGGTGTCATCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGTGGGCAGTACTGGCGTTTTTGCTc                                                                                                                 >  1:105334/1‑150 (MQ=255)                                gcgGGATTACACCGCGCACTTAAAGGCGTTTCGCACTATCTTGCATTGAAGACGCTACTCAGTTTATGGACAGGTGTCATCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGTGGGCAGTACTGGCGTTTTTGCTc                                                                                                                 >  1:1132420/1‑150 (MQ=255)                                 cgGGATTACACCGCGCACTTAAAGGCGTTTCGCACTATCTTGCATTGAAGACGCTACTCAGTTTATGGACAGGTGTCATCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGTGGGCAGTACTGGCGTTTTTGCTCa                                                                                                                <  1:1336112/150‑1 (MQ=255)                                   ggATTACACCGCGCACTTAAAGGCGTTTCGCACTATCTTGCATTGAAGACGCTACTCAGTTTATGGACAGGTGTCATCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGTGGGCAGTACTgg                                                                                                                            <  1:818541/136‑1 (MQ=255)                                   ggATTACACCGCGCACTTAAAGGCGTTTCGCACTATCTTGCATTGAAGACGCTACTCAGTTTATGGACAGGTGTCATCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGTGGGCAGTACTgg                                                                                                                            >  2:818541/1‑136 (MQ=255)                                      ttACACCGCGCACTTAAAGGCGTTTCGCACTATCTTGCATTGAAGACGCTACTCAGTTTATGGACAGGTGTCATCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGTGGGCAGTACTGGCGTTTTTGCTCAACTAc                                                                                                           >  1:1279775/1‑150 (MQ=255)                                                cACTTAAAGGCGTTTCGCACTATCTTGCATTGAAGACGCTACTCAGTTTATGGACAGGTGTCATCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGTGGGCAGTACTGGCGTTTTTGCTCAACTACGTGCCCAata                                                                                                 >  1:165636/1‑150 (MQ=255)                                                    tAAAGGCGTTTCGCACTATCTTGCATTGAAGACGCTACTCAGTTTATGGACAGGTGTCATCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGTGGGCAGTACTGGCGTTTTTGCTCAACTACGTGCCCAATATCgg                                                                                             >  2:396900/1‑150 (MQ=255)                                                     aaaGGCGTTTCGCACTATCTTGCATTGAAGACGCTACTCAGTTTATGGACAGGTGTCATCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGTGGGCAGTACTGGCGTTTTTGCTCAACTACGTGCCCAATATCGgc                                                                                            >  2:468126/1‑150 (MQ=255)                                                     aaaGGCGTTTCGCACTATCTTGCATTGAAGACGCTACTCAGTTTATGGACAGGTGTCATCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGTGGGCAATACTGGCGTTTTTGCTCAACTACGTGCCCAATATCGgc                                                                                            >  2:1296743/1‑150 (MQ=255)                                                      aaGGCGTTTCGCACTATCTTGCATTGAAGACGCTACTCAGTTTATGGACAGGTGTCATCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGTGGGCAGTACTGGCGTTTTTGCTCAACTACGTGCCCAATATCGgcg                                                                                           <  1:893409/150‑1 (MQ=255)                                                         gCGTTTCGCACTATCTTGCATTGAAGACGCTACTCAGTTTATGGACAGGTGGCATCGTCTGTCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGTGGGCAGTACTGGCGTTTTTGCTCAACTACGTGCCCAATATCGGCGCgg                                                                                        <  2:300048/150‑1 (MQ=255)                                                           gTTTCGCACTATCTTGCATTGAAGACGCTACTCAGTTTATGGACAGGTGTCATCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGTGGGCAGTACTGGCGTTTTTGCTCAACTACGTGCCCAATATCGGCGCGGTa                                                                                      >  2:567184/1‑150 (MQ=255)                                                            tttCGCACTATCTTGCATTGAAGACGCTACTCAGTTTATGGACAGGTGTCATCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGTAGTTTGCGCTGATGTGGGCAGTACTGGCGTTTTTGCTCAACTACGTGCCCAATATCGGCGCGGTaa                                                                                     >  2:1388253/1‑150 (MQ=255)                                                            tttCGCACTATCTTGCATTGAAGACGCTACTCAGTTTATGGACAGGTGTCATCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGTAGTTTGCGCTGATGTGGGCAGTACTGGCGTTTTTGCTAAACTACGTGCCCAATATCGGCGCGGTaa                                                                                     >  2:412908/1‑150 (MQ=255)                                                             ttCGCACTATCTTGCATTGAAGACGCTACTCAGTTTATGGACAGGTGTCATCGTCTGGCTGGGGCGGGAGCTGATGGGGGTGCAGTTTGCGCTGATGTGGGCAGTACTGGCGTTTTTGCTCAACTACGTGCCCAATATCGGCGCGGTAAt                                                                                    <  1:691594/150‑1 (MQ=255)                                                                  aCTATCTTGCATTGAAGACGCTACTCAGTTTATGGACAGGTGTCATCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGTGGGCAGTACTGGCGTTTTTGCTCAACTACGTGCCCAATATCGGCGCGGTAATTTccg                                                                               <  2:219384/150‑1 (MQ=255)                                                                  aCTATCTTGCATTGAAGACGCTACTCAGTTTATGGACAGGTGTCATCGGCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGTGGGCAGTACTGGCGTTTTTGCTCAACTACGTGCCCAATATCGGCGCGGTAATTTccg                                                                               <  1:284652/150‑1 (MQ=255)                                                                      tCTTGCATTGAAGACGCTACTCAGTTTATGGACAGGTGTCATCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGTGGGCAGTACTGGCGTTTTTGCTCAACTACGTGCCCAATATCGGCGCGGTAATTTCCGCCGt                                                                           <  1:7552/150‑1 (MQ=255)                                                                          gaaTTGAAGACGCTACTCAGTTTATGGACAGGTGTCATCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTCCGCAGATGTGGGCAGTACTGGCGTTTTTGCTCAACTAAGTGCCCAAAATCGTCGCGGTAATTTCCGCCGTAACg                                                                       >  2:41812/3‑150 (MQ=255)                                                                                aaGACGCTACTCAGTTTATGGACAGGTGTCATCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGTGGGCAGTACTGGCGTTTTTGCTCAACTACGTGCCCAATATCGGCGCGGTAATTTCCGCCGTACCGCCAATg                                                                 >  1:1056956/1‑150 (MQ=255)                                                                                      cTACTCAGTTTATGGACAGGTGTCATCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGTGGGCAGTACTGGCGTTTTTGCTCAACTACGTGCCCAATATCGGCGCGGTAATTTCCGCCGTACCGCCAATGATTCAg                                                           <  1:718555/150‑1 (MQ=255)                                                                                       tACTCAGTTTATGGACAGGTGTCATCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGTGGGCAGTACTGGCGTTTTTGCTCAACTACGTGCCCAATATCGGCGCGGTAATTTCCGCCGTACCGCCAATGATTCAgg                                                          >  2:1100629/1‑150 (MQ=255)                                                                                       tACTCAGTTTATGGACAGGTGTCATCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGTGGGCAGTACTGGCGTTTTTGCTCAACTACGTGCCCAATATCGGCGCGGTAATTTCCGCCGTACCGCCAATGATTCAgg                                                          <  2:1132420/150‑1 (MQ=255)                                                                                            aGTTTATGGACAGGTGTCATCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGTGGGCAGTACTGGCGTTTTTGCTCAACTACGTGCCCAATATCGGCGCGGTAATTTCCGCCGTACCGCCAATGATTCAGGtgctg                                                     >  2:468120/1‑150 (MQ=255)                                                                                             gTTTATGGACAGGTGTCATCGTCTGGCTGGGGCTGGAGCTGATGGGGTTGCAGTTTGCGCTGATGTGGGCAGTACTGGCGTTTTTGCTCAACTACGTGCCCAATATCGGCGCGGTAATTTCCGCCGTACCGCCAATGATTCAGGtgctgc                                                    <  1:687848/150‑1 (MQ=255)                                                                                                tATGGACAGGTGTCATCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGTGGGCAGTACTGGCGTTTTTGCTCAACTACGTGc                                                                                                       >  1:1363081/1‑96 (MQ=255)                                                                                                tATGGACAGGTGTCATCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGTGGGCAGTACTGGCGTTTTTGCTCAACTACGTGc                                                                                                       <  2:1363081/96‑1 (MQ=255)                                                                                                 aTGGACAGGTGTCATCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGTGGGCAGTACTGGCGTTTTTGCTCAACTACGTGCCCAATATCGGCGCGGTAATTTCCGCCGTACCGCCAATGATTCAGGTGCTGCTGtt                                                <  1:440141/150‑1 (MQ=255)                                                                                                   ggACAGGTGTCATCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGTGGGCAGTACTGGCGTTTTTGCTCAACTACGTGCCCAATATCGGCGCGGTAATTTCCGCCGTACCGCCAATGATTCAGGTGCTGCTGTTTa                                              <  2:418136/150‑1 (MQ=255)                                                                                                   ggACAGGTGTCATCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGTGGGCAGTACTGGCGTTTTTGCTCAACTACGTGCCCAATATCGGCGCGGTAATTTCCGCCGTACCGCCAATGATTCAGGTGCTGCTGTTTa                                              >  2:237058/1‑150 (MQ=255)                                                                                                    gACAGGTGTCATCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGTGGGCAGTACTGGCGTTTTTGCTCAACTACGTGCCCAATATCGGCGCGGTAATTTCCGCCGTACCGCCAATGATTCAGGTGCTGCTGTTTaa                                             >  1:1076705/1‑150 (MQ=255)                                                                                                         gtgtCATCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGTGGGCAATACTGGCGTTTTTGCTCAACTACGTGCCCAATATCGGCGCGGTAATTTCCGCCGTACCGCCAATGATTCAGGTGCTGCTGTTTAATGgtg                                        <  1:1296743/150‑1 (MQ=255)                                                                                                          tgtCATCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGTGGGCAGTACTGGCGTTTTTGCTCAACTACGTGCCCAATATCGGCGCGGTAATTTCCGCCGTACCGCCAATGATTCAGGTGCTGCTGTTTAATGgtgt                                       <  2:515661/150‑1 (MQ=255)                                                                                                              aTCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGTGGGCAGTACTGGCGTTTTTGCTCAACTACGTGCCCAATATCGGCGCGGTAATTTCCGCCGTACCGCCAATGATTCAGGTGCTGCTGTTTAATGGTGTTTAc                                   >  2:209462/1‑150 (MQ=255)                                                                                                              aTCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGTGGGCAGTACTGGCGTTTTTGCTCAACTACGTGCCCAATATCGGCGCGGTAATTTCCGCCGTACCGCCAATGATTCAGGTGCTGCTGTTTAATGGTGTTTAc                                   >  1:755710/1‑150 (MQ=255)                                                                                                               tCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGTGGGCAGTACTGGCGTTTTTGCTCAACTACGTGCCCAATATCgg                                                                                             >  2:573257/1‑91 (MQ=255)                                                                                                               tCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGTGGGCAGTACTGGCGTTTTTGCTCAACTACGTGCCCAATATCgg                                                                                             <  1:573257/91‑1 (MQ=255)                                                                                                               tCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGTGGACAGTACTGGCGTTTTTGCTCAACTACGTGCCCAATATCGGCGCGGTAATTTCCGCCGTACCGCCAATGATTCAGGTGCTGCTGTTTAATGGTGTTTACg                                  <  2:523707/150‑1 (MQ=255)                                                                                                                cGTCTGTCTGGGGCTGGAGCTTATGGGTGTGCAGTTTGCGCTGATGTGGGCAGTACTGGCGTTTTTGCTCAACTACGTGCCCAATATCTGCGCTGTAATTTCCGCCGTACCGCCAATGATTCAGGTGCTGCTGTTTAATGGTGTTTACGa                                 <  1:41812/150‑1 (MQ=255)                                                                                                                 gTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGTGGGCAGTACTGGCGTTTTTGCTCAACTACGTGCCCAATATCGGCGCGGTAATTTCCGCCGTACCGCCAATGATTCAGGTGCTGCTGTTTAATGGTGTTTACGaa                                >  1:754336/1‑150 (MQ=255)                                                                                                                    tggctggGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGTGGGCAGTACTGGCGTTTTTGCTCAACTACGTGCCCAATATCGGCGCGGTAATTTCCGCCGTACCGCCAATGATTCAGGTGCTGCTGTTTAATGGTGTTTACGAATGt                             <  1:889222/150‑1 (MQ=255)                                                                                                                     ggctggggcgggAGCTGATGGGGGTGCAGTTTGCGCGGATGTGGGCAGTACTGGCGTTTTTGCTCAACTACGTGCCCAATATCGGCGCGGTAATTTCCGCCGTACCGCCAATGATTCAGGTGCTGCTGTTTAATGGTGTTTACGAATGTa                            <  2:921641/150‑1 (MQ=255)                                                                                                                      gctggggctggAGCTGATGGGGGTGTAGTTTGCGCTGATGTGGGCAGTACTGGCGTTTTTGCTCAACTACGTGCCCAATATCGGCGCGGTAATTTCCGCCGTACCGCCAATGATTCAGGTGCTGCTGTTTAATGGTGTTTACGAATGTAt                           <  1:412908/150‑1 (MQ=255)                                                                                                                      gctggggctggAGCTGATGGGGGTGTAGTTTGCGCTGATGTGGGCAGTACTGGCGTTTTTGCTCAACTACGTGCCCAATATCGGCGCGGTAATTTCCGCCGTACCGCCAATGATTCAGGTGCTGCTGTTTAATGGTGTTTACGAATGTAt                           <  1:1388253/150‑1 (MQ=255)                                                                                                                        tggggctggAGCTGATGGGGGTGCAGTTTGCGCTGATGTGGGCAGTACTGGCGTTTTTGCTCAACTACGTGCCCAATATCGGCGCGGTAATTTCCGCCGTACCGCCAATGATTCAGGTGCTGCTGTTTAATGGTGTTTACGAATGTATTc                         >  1:779543/1‑150 (MQ=255)                                                                                                                        tggggctggAGCTGATGGGGGTGCAGTTTGCGCTGATGTGGGCAGTACTGGCGTTTTTGCTCAACTACGTGCCCAATATCGGCGCGGTAATTTCCGCCGTACCGCCAATGATTCAGGTGCTGCTGTTTAATGGTGTTTACGAATGTATTc                         >  1:17396/1‑150 (MQ=255)                                                                                                                           ggctggAGCTGATGGGGGTGCAGTTTGCGCTGATGTGGGCAGTACTGGCGTTTTTGCTCAACTACGTGCCCAATATCGGCGCGGTAATTTCCGCCGTACCGCCAATGATTCAGGTGCTGCTGTTTAATGGTGTTTACGAATGTATTCTgg                      >  1:638359/1‑150 (MQ=255)                                                                                                                           ggctggAGCTGATGGGGGTGCAGTTTGCGCTGATGTGGGCAGTACTGGCGTTTTTGCTCAACTACGTGCCCAATATCGGCGCGGTAATTTCCGCCGTACCGCCAATGATTCAGGTGCTGCTGTTTAATGGTGTTTACGAATGTATTCTgg                      <  2:1279775/150‑1 (MQ=255)                                                                                                                               ggAGCTTATGGGGGTGCAGTTTGCGCTGATGTGGGCAGTACTGGCGTTTTTGCTCAACTACGTGCCCAATATCGGCGCGGTAATTTCCGCCGTACCGCCAATGATTCAGGTGCTGCTGTTTAATGGTGTTTACGAATGTATTCTGGTCgg                  <  1:216949/150‑1 (MQ=255)                                                                                                                                   cTGATGGGGGTGCAGTTTGCGCTGATGTGGGCAGTACTGGCGTTTTTGCTCAACTACGTGCCCAATATCGGCGCGGTAATTTCCGCCGTACCGCCAATGATTCAGGTGCTGCTGTTTAATGGTGTTTACGAATGTATTCTGGTCGGCGCa              <  1:1100629/150‑1 (MQ=255)                                                                                                                                     gATGGGGGTGCAGTTTGCGCTGATGTGGGCAGTACTGGCGTTTTTGCTCAACTACGTGCCCAATATCGGCGCGGTAATTTCCGCCGTACCGCCAATGATTCAGGTGCTGCTGTTTAATGGTGTTTACGAATGTATTCTGGTCGGCGCAtt            >  1:382629/1‑150 (MQ=255)                                                                                                                                     gATGGGGGTGCAGTTTGCGCTGATGTGGGCAGTACTGGCGTTTTTGCTCAACTACGTGCCCAATATCGGCGCGGTAATTTCCGCCGTACCGCCAATGATTCAGGTGCTGCTGTTTAATGGTGTTTACGAATGTATTCTGGTCGGCGCAtt            >  2:102516/1‑150 (MQ=255)                                                                                                                                          gggTGCAGTTTGCGCTGATGTGGGCAGTACTGGCGTTTTTGCTCAACTACGTGCCCAATATCGGCGCGGTAATTTCCGCCGTACCGCCAATGATTCAGGTGCTGCTGTTTAATGGTGTTTACGAATGTATTCTGGTCGGCGCATTGtttt       >  2:1099170/1‑150 (MQ=255)                                                                                                                                           ggTGCAGTTTGCGCTGATGTGGGCAGTACTGGCGTTTTTGCTCATCTACGTGCCCAATATCGGCGCGGTAATTTCCGCCGTACCGCCAATGATTCAGGTGCTGCTGTTTAATGGTGTTTACGAATGTATTCTGGTCGGCGCATTGttttt      <  1:142616/150‑1 (MQ=255)                                                                                                                                             tGCAGTTTGCGCTGATGTGGGCAGTACTGGCGTTTTTGCTCAACTACGTGCCCAATATCGGCGCGGTAATTTCCGCCGTACCGCCAATGATTCAGGTGCTGCTGTTTAATGGTGTTTACGAATGTATTCTGGTCGGCGCATTGTTTTTAg    >  1:595179/1‑150 (MQ=255)                                                                                                                                              gCAGTTTGCGCTGATGTGGGCAGTACTGGCGTTTTTGCTCAACTACGTGCCCAATATCGGCGCGGTAATTTCCGCCGTACCGCCAATGATTCAGGTGCTGCTGTTTAATGGTGTTTACGAATGTATTCTGGTCGGCGCATTGTTTTTAGt   <  1:396900/150‑1 (MQ=255)                                                                                                                                               cAGTTTGCGCTGATGTGGGCAGTACTGGCGTTTTTGCTCAACTACGTGCCCAATATCGGCGCGGTAATTTCCGCCGTACCGCCAATGATTCAGGTGCTGCTGTTTAATGGTGTTTACGAATGTATTCTGGTCGGCGCATTGTTTTTAGTg  >  2:152665/1‑150 (MQ=255)                                                                                                                                               |                                                                                                                                                      TTTTGCGCTGAATAATCCACAGATTCACATCGCGGGATTACACCGCGCACTTAAAGGCGTTTCGCACTATCTTGCATTGAAGACGCTACTCAGTTTATGGACAGGTGTCATCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGTGGGCAGTACTGGCGTTTTTGCTCAACTACGTGCCCAATATCGGCGCGGTAATTTCCGCCGTACCGCCAATGATTCAGGTGCTGCTGTTTAATGGTGTTTACGAATGTATTCTGGTCGGCGCATTGTTTTTAGTG  >  NC_000913/3612383‑3612674

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