Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,878,707	A→C	10.1%	T331T (ACA→ACC)	yeaP →	diguanylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,878,707	0	A	C	10.1%	205.0 / 2.0	77	T331T (ACA→ACC)	yeaP	diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base A (36/29);  new base C (2/9);  total (38/39)
Fisher's exact test for biased strand distribution p-value = 4.67e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTTGTTTCACTAAACAACGAGAATGCGGATATTTCGTCGCTGCGAGAACGCATTCAGCAGCAAATACGTGGAGAATATCACTTAGGTGATGTTGATTTGTATTATCCCGGTGCCAGTCTTGGCATAGTAGAAGTCGATCCTGAAA‑CAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAA‑AACACAAACA‑GAAAACACCTTTTGTC‑GCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGCTGCCCCGGGAGCATAAAGATGAAAAAAACAACGATTATTATGATGGGTGTGGCGATTATTGTCGTACTCGGCACTGAGCTGGGA  >  NC_000913/1878488‑1878939                                                                                                                                                                                                                               |                                                                                                                                                                                                                                          gatgttttaaagggacagggagaAGGGGGATATTTCGTCGATGGGAGAACGCATGCAGCAGCAAATACGTGGAGAATTTCACTTAGGTTATGTTGATTTGTATTATCCCGGTGCCAGTCTTGGCATAGTAGAAGTCGATCCTGAAA‑CAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAA‑AACACAAACA‑GAAAACACCTTTTGTC‑GCGCATCCAg                                                                                                                                                                                                                       <  1:175889/220‑1 (MQ=255)      agtcccttaacaacGAGAATGCGGAATTTTCGGCGCTGCGAGAACGCATTCAGCAGCAAATACGTGGAGAATATCACTTAGGTGATGTTGATTTGTATTATCCCGGTGCCAGTCTTGGCATAGTAGAAGTCGATCCTGAAA‑CAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAA‑AACACAAACA‑GAAAACACCTTTTGTC‑GCGCATCCAGCGCTa                                                                                                                                                                                                                  <  2:831260/231‑1 (MQ=255)                   cGAGAATGCGGATATTTCGTCGCTGCGAGAACGCATTCAGCAGCAAATACGTGGAGAATATCACTTAGGTGATGTTGATTTGTATTATCCCGGTGCCAGTCTTGGCATAGTAGAAGTCGATCCTGAAA‑CAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAA‑AACACAAACA‑GAAAAAACCTTTTGTC‑GCGCATCCAGCGCAACATACCTGAgggg                                                                                                                                                                                                     >  1:878092/1‑237 (MQ=255)                   cGAGAATGCGGATATTTCGTCGCTGCGAGAACGCATTCAGCAGCAAATACGCGGAGAATATCACTTAGGTGATGTTGATTTGTATTATCCCGGTGCCAGTCTTGCCATCGGAGAAGTCGCTTCTGAAC‑CAACCGATGCTGACAGTCCCCTGCATCTTGCCGATATTGCGAT‑TCTCCGGGGAAAACAACCACAACA‑GAAAACCCTTTTTGCC‑GCGATTCCCACCCTCAATTCCTAAGGCg                                                                                                                                                                                                     >  1:831260/1‑239 (MQ=255)                        aTGGGGATATTTTGTCGGTGCGAGAACGCCTTCAGCAGGAAATACGTGGAGAATATCACTTAGGTGATGTTGATTTGTATTATCCCGGTGCCAGTCTTGGCATAGTAGAAGTCGATCCTGAAA‑CAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAA‑AACACAAACA‑GAAAACACCTTTTGTC‑GCGCATCCAGCGCTACATTCCTGAGGCGTATTc                                                                                                                                                                                                <  1:602924/239‑1 (MQ=255)                           cGGATATTTCGTCGGTGCGAGAACGCATTCAGCAGCAAATACGTGGAGAATATCACTTAGGTGATGTTGATTTGTATTATCCCGGTGCCAGTCTTGGCATAGTAGAAGTCGATCCTGAAA‑CAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAA‑AACACAAACA‑GAAAACACCTTTTGTC‑GCGCATCCAGCGCTACATTCCTGAGGCGTATTcaca                                                                                                                                                                                             <  2:507623/239‑1 (MQ=255)                              atatTTCGTCGCTGCGAGAACGCATTCAGCAGCAAATACGTGGAGAATATCACTTAGGTGATGTTGATTTGTATTATCCCGGTGCCAGTCTTGGCATAGTAGAAGTCGATCCTGAAA‑CAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAA‑AACACAAACA‑GAAAACACCTTTTGTC‑GCGCATCCAGCGCTACATTc                                                                                                                                                                                                             >  2:175889/1‑220 (MQ=255)                                    cgtcgCTGCGAGAACGCATTCAGCAGCAAATACGTGGAGAATATCACTTAGGTGATGTTGATTTGTATTATCCCGGTGCCAGTCTTGGCATAGTAGAAGTCGATCCTGAAA‑CAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAA‑AACACAAACA‑GAAAACACCTTTTGTC‑GCGCACCCAGCGCTACATTCCTGAGGCGTATTAACATCCTTTTGa                                                                                                                                                                                    >  2:342348/1‑239 (MQ=255)                                     gtcgCTTCGAGAACCGTTTCAGCATCAAATTCGTGGAGAATATCCATTAGGTGATGTTGATTTGGATTATCCCGGTTCCAGTCTTGGCATAGTAGAAGTCGATCCTGAAA‑CAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAA‑AACCCAAACA‑GAAAACACCTTTTGTC‑GCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGAt                                                                                                                                                                                   <  2:996033/239‑1 (MQ=255)                                     gtcgCTGCGAGAACGCATTCAGCAGCAAATACGTGGAGAATATCACTTAGGTGATGTTGATTTGTATTATCCCGGTGCCAGTCTTGGCATAGTAGAAGTCGATCCTGAAA‑CAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAA‑AACACAAACA‑GAAAACACCTTTTGTC‑GCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGAt                                                                                                                                                                                   <  2:365806/239‑1 (MQ=255)                                              ataaCGCATTCAGCAGCAAATCCGTGGAGAATATCACTTAGGTGATTTTGATTTGTATTATCCCGGTGCCAGTCTTGGCATAGTAGAAGTCGATCCTGAAA‑CAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAA‑AACACAAACA‑GAAAACACCTTTTGTC‑GCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATa                                                                                                                                                                           <  2:52674/236‑1 (MQ=255)                                                           cagcaAATACGTGGAGAATATCACTTAGGTGATGTTGATTTGTATTATCCCGGTGCCAGTCTTGGCATAGTAGAAGTCGATCCTGAAA‑CAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAA‑AACACAAACA‑GAAAACACCTTTTGTC‑GCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGt                                                                                                                                                             >  2:544826/1‑239 (MQ=255)                                                              caAATACGTGGAGAATATCACTTAGGTGATGTTGATTTGTATTATCCCGGTGCCAGTCTTGGCATAGTAGAAGTCGATCCTGAAA‑CAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAA‑AACACAAACA‑GAAAACACCTTTTGTC‑GCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGtatt                                                                                                                                                          >  2:607726/1‑239 (MQ=255)                                                                     ggggAGAATATCAATTAGGTGATGTTGATTTGTATTATCCCGGTGCCAGTCTTGGCATAGTAGAAGTCGATCCTGAAA‑CAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAA‑AACACAAACA‑GAAAACACCTTTTGTC‑GCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTc                                                                                                                                                   <  2:748394/237‑1 (MQ=255)                                                                          gtatatCACTTAGGTGATGTTGATTTGTATTATCCCGGTACCAGTCTTGGCATAGTAGAAGTCGATCCTGAAA‑CAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAA‑AACACAAACA‑GAAACAACCTTTTGTC‑GCGCATCCAGCGCTACATTCCTGAGGCGAATTTACATCCTTTTGATTGGGGAAAAAATGGGAATCGGTTTTATTTTTcgtgtg                                                                                                                                              >  1:988144/3‑234 (MQ=255)                                                                           aaTATCGCTTAGGTGATGTTGATTTGTATTATCCCGGTGCCAGTCTTGGCATAGTAGAAGTCGATCCTGAAA‑CAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAA‑AACACAAACA‑GAAAACACCTTTTGTC‑GCGCATc                                                                                                                                                                                                                          <  2:829404/162‑1 (MQ=255)                                                                             tatCCCTTAGGTGATGTTGATTTGTATTATCCCGGTGCCAGTCTTGGCATAGTAGAAGTCGATCCTGAAA‑CAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAA‑AACACAAACA‑GAAAACACCTTTTGTC‑GCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTa                                                                                                                                           <  2:680927/239‑1 (MQ=255)                                                                                     aGGTGATGTTGATTTGTATTATCCCGGTGCCAGTCTTGGCATAGTAGAAGTCTGTCCTGTAA‑CAACAGATGCAGACAGTGCCCTGCATGATTCCGGTATTGCGATGTATCAAGAGAAAA‑AACACAAACA‑GAGAACACCTTTTGTC‑GCGCATCCAGCGCTACATTCCTGACGCGTAGTCACATCATTTTGATTGTTGATAACATCCGCATCGGTATTATTTTTCCGGTTATAATCTTcct                                                                                                                                   >  2:418786/1‑237 (MQ=255)                                                                                     aGGTGATGTTGATTTGTATTATCCCGGTGCCAGTCTTGGCATAGTAGAAGTCGATCCTGAAA‑CAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAA‑AACACAAACA‑GAAAACACCTTTTGTC‑GCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTcct                                                                                                                                   >  2:317044/1‑237 (MQ=255)                                                                                     aGGTGATGTTGATTTGTATTATCCCGGTGCCAGTCTTGGCATAGTAGAAGTCGATCCTGAAA‑CAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAA‑AACACAAACA‑GAAAAACCCTTTTGTC‑GCGCAGCCAGCGCTCCa                                                                                                                                                                                                                >  1:829404/1‑162 (MQ=255)                                                                                     aGGTGATGTTGATTTGTATTATCCCGGTGCCAGTCTTGGCATAGTAGAAGTCGATCCTGAAA‑CAACCGATACAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAA‑AACACAAACA‑GAAAACACCTTTTGTC‑GCGCATCCAGCGCTACATTCCTGAGGCGTATTCCCATCCTTTTGTTTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCttatt                                                                                                                                   >  1:390856/1‑236 (MQ=255)                                                                                       gTGATGTTGATTTGTATTATCCCGGTGCCAGTCTTGGCATAGTAGAAGTCGATCCTGAAA‑CAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAA‑AACACAAACA‑GAAAACACCTTTTGTC‑GCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAAt                                                                                                                                         >  2:201322/1‑231 (MQ=255)                                                                                       gTGATGTTGATTTGTATTATCCCGGTGCCAGTCTTGGCATAGTAGAAGTCGATCCTGAAA‑CAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAA‑AACACAAACA‑GAAAACACCTTTTGTC‑GCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAAt                                                                                                                                         <  1:201322/231‑1 (MQ=255)                                                                                       gTGATGTTGATTTGTATTATCCCGGTGCCAGTCTTGGCATAGTAGAAGTCGATCCTGAAA‑CAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAA‑AACACAAACA‑GAAAACACCTTTTGTC‑GCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTa                                                                                                                                 >  2:574161/1‑239 (MQ=255)                                                                                        ttgttttttttttTATTTTCCCGGGTCCGGTCGTGGCCTAGGAGAAGTGGATCCTGAAACCACCCGAGCCAG‑CAGGGCCCTGCATGCTGCCTATATTGCGATGTATCAGGAGAAAA‑AACACAAACA‑GAAAACACCTTTTGTC‑GCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTa                                                                                                                                 <  1:544826/229‑1 (MQ=255)                                                                                          aTGTTGATTTGTATTATCCCGGTGCCAGTCTTGGCATAGTAGAAGTCGATCCTGAAA‑CAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAA‑AACACAAACA‑GAAAACACCTTTTGTC‑GCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAg                                                                                                                              >  1:224969/1‑239 (MQ=255)                                                                                               gATTTGTATTATCCCGGTGCCAGTCTTGGCATAGTAGAAGTCGATCCTGAAA‑CAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAA‑AACACAAACA‑GAAAACACCTTTTGTC‑GCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACACCGGTc                                                                                                                         >  2:399529/1‑239 (MQ=255)                                                                                                           cccGGTGCCAGTCTTGGCATAGTAGAAGTCGATCCTGAAA‑CAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAA‑AACACAAACA‑GAAAACACCTTTTGTC‑GCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCa                                                                                                             <  2:825546/239‑1 (MQ=255)                                                                                                               gTGCCAGTCTTGGCATAGTAGAAGTCGCTCCTGAAA‑CAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGCAAA‑CACA‑AAACACGAACACACCATTTTTC‑GCGCATCAAGCGCCACATACCTGAGGAGTATTCACATCCTTTTGCTTGGTGATAACCTGCGAATCGGTATTATTTTTCCGGTTATAATCTTCATTCCAGCGGTGATATTTTTAACATGGt                                                                                                         >  1:784737/1‑239 (MQ=255)                                                                                                               gTGCCAGTCTTGGCATAGTAGAAGTCGATCCTGAAA‑CAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGCTGTATCAGGAGAAAA‑AACACAAACA‑GAAAACACCTTTTGTC‑GCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCGTTCAACCGGTCGTTTTTTTAGCATGGt                                                                                                         >  1:679133/1‑239 (MQ=255)                                                                                                               gTGCCAGTCTTGGCATAGTAGAAGTCGATCCTGAAA‑CAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAA‑AACACAAACA‑GAAAACACCTTTTGTC‑GCGCATCCAGCGCTACATTCCTGAGGCGTATTcac                                                                                                                                                                                              <  1:998811/154‑1 (MQ=255)                                                                                                               gTGCCAGTCTTGGCATAGTAGAAGTCGATCCTGAAA‑CAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAA‑AACACAAACA‑GAAAACACCTTTTGTC‑GCGCATCCAGCGCTACATTCCTGAGGCGTATTcac                                                                                                                                                                                              >  2:998811/1‑154 (MQ=255)                                                                                                               gTGCCAGTCTTGGCATAGTAGAAGTCGATCCTGAAA‑CAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAA‑AACACAAACA‑GAAAACACCTTTTGTC‑GCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGTTAACATGCGAATCGGTATTTTTTTTCCGGTTGTAATCTTATTTTCAGCGGTGGAATTTTTAgatgggt                                                                                                         >  1:263222/1‑233 (MQ=255)                                                                                                               gTGCCAGTCTTGGCATAGTAGAAGTCGATCCTGAAA‑CAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAA‑AACACAAACA‑GAAAACACCTTTTGTC‑GCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAAt                                                                                                                                                                 <  2:704660/183‑1 (MQ=255)                                                                                                               gTGCCAGTCTTGGCATAGTAGAAGTCGATCCTGAAA‑CAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAA‑AACACAAACA‑GAAAACACCTTTTGTC‑GCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGGATTATTTTTCCGGTTGAAATCTTCGTTACAGCGGTCGGATTTTTAACATGGt                                                                                                         >  1:901813/1‑239 (MQ=255)                                                                                                               gTGCCAGTCTTGGCATAGTAGAAGTCGATCCTGAAA‑CAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAA‑AACACAAACA‑GAAAACACCTTTTGTC‑GCGCATCCAGCGCTACATTCCTGAGGCGGATTCCAACCATTTTGATTGGTGAAAACATGCGAATCGGTATTATTTTTCCGGTTTTAATCTTCTTTAAAGCGGTGGGATTTTAAGCATgtt                                                                                                         >  1:580973/1‑237 (MQ=255)                                                                                                               gTGCCAGTCTTGGCATAGTAGAAGTCGATCCTGAAA‑CAACCGATGCAGACACTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAA‑AACACAAACA‑GAAAACACCTTTTGTC‑TCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAAt                                                                                                                                                                 >  1:704660/1‑183 (MQ=255)                                                                                                               gTGCCAGTCTTGGCATAGTAGAAGTCGATCCTGAAA‑CAACCCATGCAGACAGTGCACTGCATGCTGTCGATATTGCGATGTATCAGGAGAAAA‑AACACAAACA‑GAAAACACCTTTTGGC‑GCGCATCCAGAGCTACATTCCTGAGGCGCATTTACATCCTTTTGCTTGGTGACAACATGCGAATCGGGATTATTTTTCCGGTTGTCATCTTCATTATAGCGGTCGAATTTTTAGCATGGt                                                                                                         >  2:738702/1‑239 (MQ=255)                                                                                                                tGCCAGTCTTGGCATAGTAGAAGTCGATCCTGAAA‑CAACCGATGCAGACAGTGACCTGCATGCGGCCGATATTGCGATGTATCAGGAGAAAA‑AACACAACCA‑GAAAACACCTTTTGTC‑GCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAAt                                                                                                                                                                 <  2:885558/182‑1 (MQ=255)                                                                                                                             atatTAGAATGCGGTCCTGAAA‑CAACCAGTGCAAACAGGGCCCGGCATGGTGCCGTTATTGCGATGTATCAGGGAAAAA‑AACACAAACA‑GAAAACACCTTTTGTC‑GCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTggcgg                                                                                            <  1:317044/238‑1 (MQ=255)                                                                                                                                gtagAAGTCGATCCTGAAA‑CAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAA‑AACACAAACA‑GAAAACACCTTTTGTC‑GCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCg                                                                                                                            >  1:834864/1‑203 (MQ=255)                                                                                                                                gtagAAGTCGATCCTGAAA‑CAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAA‑AACACAAACA‑GAAAACACCTTTTGTC‑GCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCg                                                                                                                            <  2:834864/203‑1 (MQ=255)                                                                                                                                  agAAGTCGATCCTGAAA‑CAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAA‑AACACAAACA‑GAAAACACCTTTTGTC‑GCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTAtgc                                                                                      <  2:224969/239‑1 (MQ=255)                                                                                                                                         gATCCTGAAA‑CAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAA‑AACACAAACA‑GAAAACACCTTTTGTC‑GCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTCACGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCttattgtccccg                                                                               >  1:571296/1‑228 (MQ=255)                                                                                                                                           tCCTGAAA‑CAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAA‑AACACAAACA‑GAAAACACCTTTTGTC‑GCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGACTCGGTATTATTTTTCCGGTTGTAATCTTCATTTCAGCGGTCTTATTTTTAGCATGGTTTTTTTTTGGCGGGTTAGCCGCCCCggg                                                                             >  1:14494/1‑239 (MQ=255)                                                                                                                                           tCCTGAAA‑CAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAA‑AACACAAACA‑GAAAACACCTTTTGTC‑GCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCAt                                                                                                                                   >  2:186982/1‑185 (MQ=255)                                                                                                                                           tCCTGAAA‑CAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAA‑AACACAAACA‑GAAAACACCTTTTGTC‑GCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTAATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGTCGGCTTTTCTTCCCCtgg                                                                             >  2:998848/1‑239 (MQ=255)                                                                                                                                           tCCTGAAA‑CAACCGATGAAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAA‑AACACAAACA‑GAAAACACCTTTTGTC‑GCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCAt                                                                                                                                   <  1:186982/185‑1 (MQ=255)                                                                                                                                                 aa‑caacCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAA‑AACACAAACA‑GAAAACACCTTTTGTC‑GCGCATCCAGCGCTACATTCCTGAGGCGGATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTggcggccttattctcccgggggcgca                                                                       >  1:937511/1‑219 (MQ=255)                                                                                                                                                      acCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAA‑AACACAAACA‑GAAACCACCTTTTGTC‑GCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTc                                                                                                                         <  2:561338/185‑1 (MQ=255)                                                                                                                                                                aaagggCCCTGCATGGTGCCGATTTTGCGATGTATCAGGAGAAAA‑AACACAAACA‑GAAAACACCTTTTGTC‑GCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTTATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGCTGCCCCGGGAGCATAAAGATGAAAAAaac                                                         <  2:263222/234‑1 (MQ=255)                                                                                                                                                                aCAGTTCCCTTGATGCTGCCGATATTGCGATTTATCCGGAGAAAA‑AACACAAACA‑GAAAACACCTTTTGTC‑GCGCATCCAGCGCTTCATTCCTGAGGCCTATTCACATCCTTTTGATTTGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGCTGCCCCGGGAGCATAAAGATGAAAAAaac                                                         <  2:878092/239‑1 (MQ=255)                                                                                                                                                                  aGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAA‑AACACAAACA‑GAAAACACCTTTTGTC‑GCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGCTGCCCCGGGAGCATAAAGATGAAAAAaacaa                                                       <  2:1000688/239‑1 (MQ=255)                                                                                                                                                                   tggCCCTGACTTCTGACGATATTGCGATGTATCAGGAGAAAA‑AACACAAACA‑GAAAACCCCTTTTGTC‑GCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGCTGCCCCGGGAGCATAAAGATGAAAAAaacaac                                                      <  1:607726/237‑1 (MQ=255)                                                                                                                                                                     gCCCTGCATGCTGCCGATTTTGCGATGTATCAGGAGAAAA‑AACACAAACA‑GAAAACCCCTTTTGTC‑GCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCAGTGTTTTTTATTGGCGGCTATGCTGCCCCGGGAGCATAAAGATGAAAAAAACAACGa                                                    <  2:988144/239‑1 (MQ=255)                                                                                                                                                                      cccTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAA‑AACACAAACA‑GAAAACCCCTTTTGTC‑GCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGCTGCCCCGGGAGCATAAAGATGAAAAAAACAACGa                                                    <  2:937511/238‑1 (MQ=255)                                                                                                                                                                       ccTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAA‑AACACAAACA‑GAAAACACCTTTTGTC‑GCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGCTGCCCCGGGAGCATAAAGATGAAAAAAACAACGatt                                                  <  1:574161/239‑1 (MQ=255)                                                                                                                                                                          gCATGCTGCCGATATTGCGATGTATCAGGAGAAAA‑AACACAAACA‑GAAAACACCTTTTGTC‑GCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCCTTACAGCGGGCGCATTTTTAGCATGGTTTTTTATTTGCGGCTATTCGTCCCCCGGGGCCTAAAGAAGaaaaaaaaaaaaaaaatt                                               >  1:102949/1‑231 (MQ=255)                                                                                                                                                                               ctgcCGATATTGCGATGTATCAGGAGAAAA‑AACACAAACA‑GAAAACCCCTTTTGTC‑GCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGCTGCCCCGGGAGCATAAAGATGAAAAAAACAACGATTATTatgat                                          <  2:464980/239‑1 (MQ=255)                                                                                                                                                                                tgcCGATATTGCGATGTATCAGGAGAAAA‑AACACAAACA‑GAAAACACCTTTTGTC‑GCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAAt                                                                                                                                                                 >  1:116970/1‑119 (MQ=255)                                                                                                                                                                                tgcCGATATTGCGATGTATCAGGAGAAAA‑AACACAAACA‑GAAAACACCTTTTGTC‑GCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAAt                                                                                                                                                                 <  2:116970/119‑1 (MQ=255)                                                                                                                                                                                   cctttttGCGATGTAACAGGAAAAAA‑AACACAAACA‑GAAAACCCCTTTTGTC‑GCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGCTGCCCCGGGAGCATAAAGATGAAAAAAACAACGATTATTATGATGGgt                                      <  1:858223/234‑1 (MQ=255)                                                                                                                                                                                     atatTGCGATGTATCAGGAGAAAA‑AACACAAACA‑GAAAACACCTTTTGTC‑GCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGa                                                                                                                                                                                    <  2:392043/95‑1 (MQ=255)                                                                                                                                                                                     atatTGCGATGTATCAGGAGAAAA‑AACACAAACA‑GAAAACACCTTTTGTC‑GCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGa                                                                                                                                                                                    >  1:392043/1‑95 (MQ=255)                                                                                                                                                                                     atatTGCGATGTATCAGGAGAAAA‑AACACAAACA‑GAAAACACCTTTTGTC‑GCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGTTGCCCCCGGGGCTTAAAAATGTAAAAAAAAAAAATATTTATAATggggtt                                    >  2:858223/1‑236 (MQ=255)                                                                                                                                                                                           cGATGTATCAGGAGAAAA‑AACACAAACA‑GAAAACACCTTTTGTC‑GCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCCTTCATTCCCCCGGGGCATAAAAAAAAAAAAAACAAAGaatatttataaggggtggggggg                              >  1:1018295/1‑222 (MQ=255)                                                                                                                                                                                            gATGTATCAGGAGAAAA‑AACACAAACA‑GAAAACACCTTTTGTC‑GCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTA‑TATTTTTCCGGTTGTAAt                                                                                                                                         >  1:531009/1‑130 (MQ=255)                                                                                                                                                                                            gATGTATCAGGAGAAAA‑AACACAAACA‑GAAAACACCTTTTGTC‑GCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTA‑TATTTTTCCGGTTGTAAt                                                                                                                                         <  2:531009/130‑1 (MQ=255)                                                                                                                                                                                              ttttttagggaaaaa‑aaCCCAAACC‑AAAAACCCCTTTTGTCGGCGTATCAAG‑GCTTCCTTCCTGGGGCTTATTCCCATCCTTTTGATTGGTTATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGCTGCCCCGGGAGCATAAAGATGAAAAAAACAACGATTATTATGATGGGTGTGGCGatta                            <  1:547404/227‑1 (MQ=255)                                                                                                                                                                                                       ggaAAA‑AACACAAACA‑GAAAACTCCTTTTTTC‑GCGCATCCCGGGCTACATTCCTGAGGCGTATTCACATCCTTTTGATGGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGCTGCCCCGGGAGCATAAAGATGAAAAAAACAACGATTATTATGATGGGTGTGGCGATTATTGTCGTAc                   <  1:342348/237‑1 (MQ=255)                                                                                                                                                                                                        gaAAA‑AACACAAACA‑GAAAACACCTTTTGTC‑GCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCCTTGCTGCCCCGGGGGCATAAAAATGAAAAAAAAAACGATTATTTTGATGGgtggtggggtttttgttggact                  >  2:547404/1‑219 (MQ=255)                                                                                                                                                                                                         gaca‑aaCCCAAACA‑GAAAACCCCTTTTGTC‑GGCCATCCAGCGCTACATTCCTGGGGCTTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGCTGCCCCGGGAGCATAAAGATGAAAAAAACAACGATTATTATGATGGGTGTGGCGATTATTGTCGTACt                  <  1:544649/234‑1 (MQ=255)                                                                                                                                                                                                           aa‑aaCACAAACA‑GAAAACACCTTTTGTC‑GCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGCTGCCCCGGGAGGATAAAGATTGAAAAAACAACCATTATTTTGATGGGTGTTGCGGTTATTGTCGCACTccg               >  2:363128/1‑236 (MQ=255)                                                                                                                                                                                                            a‑aaCACAAACA‑GAAAACCCCTTTTGTC‑GCGCATCCAGCGCTTCATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGCTGCCCCGGGAGCATAAAGATGAAAAAAACAACGATTATTATGATGGGTGTGGCGATTATTGTCGTACTCGGCa             <  2:780474/239‑1 (MQ=255)                                                                                                                                                                                                                cacaAACA‑GAAAACACCTTTTGTC‑GCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGCTGCCCCGGGAGCATAAAGATGAAAAAAACAACGATTATTATGATGGGTGTGGCGATTATTGTCGTACTCGGCACTg          <  2:679133/239‑1 (MQ=255)                                                                                                                                                                                                                           aaaCCCCTTTTGTC‑GCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGCTGCCCCGGGAGCATAAAGATGAAAAAAACAACGATTATTATGATGGGTGTGGCGATTATTGTCGTACTCGGCACTGAGCTGGGa  <  2:571296/237‑1 (MQ=255)                                                                                                                                                                                                                           aaaCACCTTTTGTC‑GCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGCTGCCCCGGGGGCATAAAAATGAAAAAAACAACGATTATTATTATGGGGGTGGCGGTTATTGTCGTACTCGGCACTGAGcccgtc  >  2:544649/1‑232 (MQ=255)                                                                                                                                                                                                                               |                                                                                                                                                                                                                                          AGTTGTTTCACTAAACAACGAGAATGCGGATATTTCGTCGCTGCGAGAACGCATTCAGCAGCAAATACGTGGAGAATATCACTTAGGTGATGTTGATTTGTATTATCCCGGTGCCAGTCTTGGCATAGTAGAAGTCGATCCTGAAA‑CAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAA‑AACACAAACA‑GAAAACACCTTTTGTC‑GCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGCTGCCCCGGGAGCATAAAGATGAAAAAAACAACGATTATTATGATGGGTGTGGCGATTATTGTCGTACTCGGCACTGAGCTGGGA  >  NC_000913/1878488‑1878939

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