Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,980,111	C→G	9.4%	S120C (TCC→TGC) ‡	uspC →	universal stress protein C

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,980,111	0	C	G	9.4%	262.1 / 15.8	85	S120C (TCC→TGC) ‡	uspC	universal stress protein C
Reads supporting (aligned to +/- strand):  ref base C (39/37);  new base G (3/5);  total (43/42)
Fisher's exact test for biased strand distribution p-value = 7.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.12e-01

GTTAATTCAGGATGCGGGGTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTT  >  NC_000913/1979965‑1980260                                                                                                                                                    |                                                                                                                                                      gTTAATTCAGGATGCGGGGTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTg                                                                                                                                                    >  1:580162/1‑150 (MQ=255)  ttAATTCAGGATGCGGGGTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGc                                                                                                                                                   >  2:450673/1‑150 (MQ=255)      ttCAGGATGCGGGGTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCgg                                                                                                                                               >  2:85045/1‑150 (MQ=255)       tCAGGATGCGGGGTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCCGAGCAGGATGCTCgc                                                                                                                                               >  2:1593335/1‑149 (MQ=255)       tCAGGATGCGGGGTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGc                                                                                                                                              >  1:1114010/1‑150 (MQ=255)       tCAGGATGCGGGGTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGc                                                                                                                                              >  2:1503071/1‑150 (MQ=255)         aGGATGCGGGGTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCCGAGCAGGATGCTCgcga                                                                                                                                             <  2:297650/150‑4 (MQ=255)          ggATGCGGGGTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCCGAGCAGGATGCTCgcgaa                                                                                                                                            <  2:379822/150‑5 (MQ=255)          ggATGCGGGGTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGaa                                                                                                                                           <  2:141400/150‑1 (MQ=255)            aTGCGGGGTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCTGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGaaaa                                                                                                                                         <  2:720818/150‑1 (MQ=255)            aTGCGGGGTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCCGAGCAGGATGCTC‑GCGaaaa                                                                                                                                         <  2:1020791/150‑1 (MQ=255)            aTGCGGGGTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGaaaa                                                                                                                                         >  1:1392907/1‑150 (MQ=255)             tGCGGGGTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAag                                                                                                                                        >  2:1229464/1‑150 (MQ=255)               cGGGGTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCCGAGCAGGATGCTC‑GCGAAAagaa                                                                                                                                      <  1:274348/150‑2 (MQ=255)               cGGGGTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCCGAGCAGGATGCTC‑GCGAAAagaa                                                                                                                                      <  1:1593335/150‑2 (MQ=255)                ggggTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGt                                                                                                                                     <  2:1214926/150‑1 (MQ=255)                 gtgTATCCCGTTGACAAGACGTTTATTGACTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGtt                                                                                                                                    >  2:102388/3‑150 (MQ=255)                     aTCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTg                                                                                                                                <  2:807015/150‑1 (MQ=255)                     aTCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTg                                                                                                                                >  1:79025/1‑150 (MQ=255)                       cccGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGcc                                                                                                                              >  1:1203836/1‑150 (MQ=255)                               cAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGa                                                                                                                      <  1:1504741/150‑1 (MQ=255)                                 aGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAgg                                                                                                                    <  1:771950/150‑1 (MQ=255)                                 aGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAgg                                                                                                                    <  1:766209/150‑1 (MQ=255)                                 aGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAgg                                                                                                                    <  2:271239/150‑1 (MQ=255)                                      tttATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTACGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGAt                                                                                                               <  1:395711/150‑1 (MQ=255)                                       ttATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGAtg                                                                                                              <  1:714919/150‑1 (MQ=255)                                         aTTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATTATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGAtgtg                                                                                                            <  2:1531476/150‑1 (MQ=255)                                         aTTGCCTACGGTGAGTTAAGCGAACATATTCTAGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGAtgtg                                                                                                            <  1:942483/150‑1 (MQ=255)                                          ttGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCAACCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGc                                                                                                           <  1:521923/150‑1 (MQ=255)                                            gCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTg                                                                                                         >  1:1439689/1‑150 (MQ=255)                                             ccTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGt                                                                                                        <  2:248226/150‑1 (MQ=255)                                              cTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGtt                                                                                                       <  2:1392907/150‑1 (MQ=255)                                               tACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCTGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTa                                                                                                      <  2:1483360/150‑1 (MQ=255)                                                  ggTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGtt                                                                                                   <  1:691587/150‑1 (MQ=255)                                                      aGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCAc                                                                                               <  1:85045/150‑1 (MQ=255)                                                        ttAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACtt                                                                                             >  1:270055/1‑150 (MQ=255)                                                          aaGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTAc                                                                                           >  1:1511854/1‑150 (MQ=255)                                                             cGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACggg                                                                                        >  2:665713/1‑150 (MQ=255)                                                             cGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACggg                                                                                        >  1:173245/1‑150 (MQ=255)                                                              gAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGga                                                                                       <  2:1439689/150‑1 (MQ=255)                                                              gAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGga                                                                                       >  2:1201231/1‑150 (MQ=255)                                                                  atatTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGAtt                                                                                   >  2:308818/1‑150 (MQ=255)                                                                   tatTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTa                                                                                  >  2:1629110/1‑150 (MQ=255)                                                                    atTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTaa                                                                                 >  2:295038/1‑150 (MQ=255)                                                                      tCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAATAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAAcc                                                                               <  2:580162/150‑1 (MQ=255)                                                                      tCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAAcc                                                                               >  2:989608/1‑150 (MQ=255)                                                                      tCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAAcc                                                                               <  1:615927/150‑1 (MQ=255)                                                                       cTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCg                                                                              <  1:230065/150‑1 (MQ=255)                                                                        tGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGc                                                                             >  2:95670/1‑150 (MQ=255)                                                                          gAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTTCTGTTAGTTCCACTTACGGGAGATTAACCGCTc                                                                           >  2:294506/1‑150 (MQ=255)                                                                            aGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCt                                                                         <  2:1064946/150‑1 (MQ=255)                                                                                 gTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCa                                                                    >  1:1306439/1‑150 (MQ=255)                                                                                 gTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCa                                                                    >  1:1312177/1‑150 (MQ=255)                                                                                  tCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCaa                                                                   <  1:1503071/150‑1 (MQ=255)                                                                                    acaAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGc                                                                 >  1:805580/1‑150 (MQ=255)                                                                                      aaGCATCATTTCGATTTGGTTATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGATCATCATTCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCtt                                                               >  2:665377/1‑150 (MQ=255)                                                                                        acatcatTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTg                                                             >  2:1533027/2‑150 (MQ=255)                                                                                         catcatTTCGATTTGGTGATTTGCGGTAATCACAATAACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTgg                                                            >  1:10323/1‑150 (MQ=255)                                                                                           tcatTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGaa                                                          <  1:1189655/150‑1 (MQ=255)                                                                                            catTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGaaa                                                         >  1:151863/1‑150 (MQ=255)                                                                                             atTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAg                                                        >  1:1050042/1‑150 (MQ=255)                                                                                              tttCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTTCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTCCGCAAGCTTTGGAAAgg                                                       <  1:294506/150‑1 (MQ=255)                                                                                              tttCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAgg                                                       <  2:1314233/150‑1 (MQ=255)                                                                                               ttCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGt                                                      >  2:1272995/1‑150 (MQ=255)                                                                                                    tttGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCaa                                                 <  1:1028558/150‑1 (MQ=255)                                                                                                       ggTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACtt                                              >  2:543143/1‑150 (MQ=255)                                                                                                          gATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTAt                                           >  1:257545/1‑150 (MQ=255)                                                                                                            tttGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATcg                                         <  1:1229464/150‑1 (MQ=255)                                                                                                               gCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCgcgc                                      >  2:656499/1‑150 (MQ=255)                                                                                                                  gTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCgctgc                                   <  2:1114010/150‑1 (MQ=255)                                                                                                                   tAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCgctgct                                  >  2:1262170/1‑150 (MQ=255)                                                                                                                       cacaatcacaGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCt                              <  2:532453/150‑1 (MQ=255)                                                                                                                       cacaatcacaGTTTCTTTTCG‑CGAGCATACTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCAAATTACTGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCt                              >  2:1187761/1‑150 (MQ=255)                                                                                                                            tcacaGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCcg                         <  1:546311/150‑1 (MQ=255)                                                                                                                            tcacaGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTAACGCGCTGCTGGCTTTCcg                         >  1:1550537/1‑150 (MQ=255)                                                                                                                             cagaGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAATAGTTATTTCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCcgc                        <  2:493476/147‑1 (MQ=255)                                                                                                                              acaGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCcgcg                       <  1:857619/150‑1 (MQ=255)                                                                                                                                 gTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCtt                    >  1:1498376/1‑150 (MQ=255)                                                                                                                                   ctCTTTTCGCCGAGCAGGATGCTC‑GCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTcg                  >  1:622745/2‑150 (MQ=255)                                                                                                                                     cTTTTCGCCGAGCAGGATGCTC‑GCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTcgcg                >  1:1162768/1‑150 (MQ=255)                                                                                                                                     cTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTcgcg                <  1:989608/150‑1 (MQ=255)                                                                                                                                          cg‑cgAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGAAGCGCGGaaag           <  1:1262170/150‑3 (MQ=255)                                                                                                                                           g‑cgAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTa          >  1:91597/1‑150 (MQ=255)                                                                                                                                                 cATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGc     <  2:977968/150‑1 (MQ=255)                                                                                                                                                   tCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCtt   >  1:410239/1‑150 (MQ=255)                                                                                                                                                    ccTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCttt  <  2:1203836/150‑1 (MQ=255)                                                                                                                                                    ccTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCttt  <  1:450673/150‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                      GTTAATTCAGGATGCGGGGTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTT  >  NC_000913/1979965‑1980260

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