Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,336,951	A→C	22.1%	G27G (GGT→GGG)	adiC ←	arginine:agmatine antiporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,336,951	0	A	C	22.1%	161.2 / 13.0	76	G27G (GGT→GGG)	adiC	arginine:agmatine antiporter
Reads supporting (aligned to +/- strand):  ref base A (40/16);  new base C (2/18);  total (42/34)
Fisher's exact test for biased strand distribution p-value = 2.29e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATGGCGATATTGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTG‑CACCGATAATCGTCACCAACCATCCATAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAG‑AAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGA‑CGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTTATTAATATAGCTCTTCGTTAAGGTTAAAAATCAGGCGACAGTGCGTTCT  >  NC_000913/4336720‑4337186                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                             caTGGCGATATTGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGGG‑CACCGATAATCGGCACCAACCATCCATAAATGG‑CAATCCCGCCAGGAGAGGCCAGGGTTGCAGGTAACAG‑AAAAACACCTGAc                                                                                                                                                                                                                                        >  1:348237/1‑238 (MQ=255)    ggCGATATTGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTTGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTG‑CACCGATAATCGTCACCAAACATCCATAAATGG‑CAATCCCGCCAGTAAAGGCCAGGCTTTCAGGGAAAAG‑AAAAACACCTGAccccc                                                                                                                                                                                                                                    >  1:329967/1‑238 (MQ=255)           ttGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGGAAGAAACAACAGGACTTGGTTCGCGGAACGACCTTTTGGCGGAAACCAACTAAGGCCCCCGTG‑CAGCGATAATCAGCAACAACCATCCATAAAAGA‑CAAACCCCCCCACAGAGGCCCAGCTGGCCAGGAAAAA‑AAAAACCCCAGCCCCCataaaaat                                                                                                                                                                                                                             >  2:654909/1‑235 (MQ=255)                 aTCCAGCAGGCCAGCCAGTAGGGGACGTTGGTTTGATAACCGAGAACCGGGCCAAAGCAGCGGCGGGCGTAAGGGTAAGAACCCCCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTG‑CACCGATAATCGTCACCAACCATCCATAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAG‑AAAAACACCTGACCCCATAATATTCCCCGa                                                                                                                                                                                                                       <  1:713290/239‑1 (MQ=255)                       cagGCCAGCCAGTAGAGGACGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTG‑CACCGATAATCGTCACCAACCATCCATAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAG‑AAAAACACCTGACCCCCTAATATTCCCCGccacca                                                                                                                                                                                                                  >  1:204820/1‑238 (MQ=255)                        agGCCAGCCAGTAGAGGACGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTG‑CACCGATAATCGTCACCAACCATCCATAAATGG‑CAATCCCGCCCGTAGAGGCCAGGATTTCAGGTAACAG‑AAAAACACCTGCCCCCATAAAATTCCCCGACACCATc                                                                                                                                                                                                                >  1:404907/1‑239 (MQ=255)                            cagccagTAGAGGACGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTG‑CACCGATAATCGTCACCAACCATCCATAAATGG‑CAATCCCGCCCGTAGAGGCCAGATTTGCCGGGAACAG‑AAAAACACCTGACCCCCTAATATTCCCCGACACCCTCAggg                                                                                                                                                                                                            >  1:382324/1‑239 (MQ=255)                                    aggggACGTTGGTTTGAAAACCGAGAAGCGGGCAAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACTTTTTGGCGTATCCCATCGAGAGCCCCAGTG‑CACCGATAATCGTCCCCAACCATCCATAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAG‑AAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGA‑Cg                                                                                                                                                                                                      <  1:263308/233‑1 (MQ=255)                                    agagGACGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTG‑CACCGATAATCGTCACCAACCATCCATAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAG‑AAAAACACCTGACCCCATAAAATTCCCCGACACCATCAGGGTGA‑Cgggg                                                                                                                                                                                                   >  1:47268/1‑239 (MQ=255)                                           gTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTG‑CACCGGTAATCGTCCCCAACCATCCATAAAAGG‑CCATCCCCCCAGTAGAGGCCAGGTTTGCAGGGAACAGAAAAAACCCCTAACCCCA‑AATATTTCCCGACACCACCAGGGGGACCGGGGATTaa                                                                                                                                                                                              >  2:400227/1‑238 (MQ=255)                                             tGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTG‑CACCGATAATCGTCACCAACCATCCATAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGAAACAG‑AAAAAAACCTGACCCCATAATATTCCCCGAACCCATCAGGGGGA‑CGGGGATTAAGccc                                                                                                                                                                                          >  1:232845/1‑239 (MQ=255)                                             tGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACGTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTG‑CACCGATAATCGTCACCAACCATCCATAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTTCAGGGAACAG‑AAAAACACCTGCCCCCCTAAAATTCCCCGAAACAATAAGGGGGA‑AGGGGATAAAAccc                                                                                                                                                                                          >  2:432373/1‑239 (MQ=255)                                                  tGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAAAGCCCCAGTG‑CACCGGTAATCGTCACCAACCATCCATAAATGG‑CAATCCCCCCACTAGAGGCCAGGTTTGCAGGGAACAA‑AAAAACACCTCACCCCAAAATATTCCCCGCCACCATCAGGGGAA‑CGGGGAGTAAACCCACttt                                                                                                                                                                                     >  2:883174/1‑239 (MQ=255)                                                         cGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTG‑CACCGATAATCGTCACCAACCATCCATAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAG‑AAAAACACCTGACCCCCTAATATTCCCCGACCCCATCAGGGGGA‑CGGGGATTAAACCCACttttgtggca                                                                                                                                                                              >  2:100297/1‑232 (MQ=255)                                                         cGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTG‑CACCGATAATCGTCACCAACCATCCATAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAG‑AAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGA‑CGGGGATTAAGCCCACTTTGTGagca                                                                                                                                                                              >  1:626282/1‑239 (MQ=255)                                                          gagaAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTG‑CACCGATAATCGTCACCAACCATCCATAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAG‑AAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGGGA‑CGGGGATTAAACCCCCTTTGTGagcct                                                                                                                                                                             >  1:231020/1‑237 (MQ=255)                                                                       aaGCAGGGGCGGGCGAAAGCGTAAGACCCCCCGGGATTTGGGTCGAGGAACGACATTTTGGCTTATACCATCGAGAGCCCCAGTC‑CACCGATAACCGTCACCACACATCAAAAAAGGC‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAG‑AAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGA‑CGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGa                                                                                                                                                                <  1:608389/239‑1 (MQ=255)                                                                           aggggggggcgTAACCGTAAGACCCACAGGGTCGTGGGTGGAGGACCGACATTTTGGCGTATCCCATCGAGGCCCCCAGGG‑CACCGATAATCGTCACCAACCACCAATAAATGC‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAG‑AAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGA‑CGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACa                                                                                                                                                              <  1:344500/231‑1 (MQ=255)                                                                                  ggcgTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTG‑CACCGATAATCGTCACCAACCATCCATAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAG‑AAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGA‑CGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTaaa                                                                                                                                                      >  1:779965/1‑238 (MQ=255)                                                                                             aagAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTG‑CACCGATAATCGTCACCAACCATCCATAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAG‑AAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGA‑CGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCGCCGAAga                                                                                                                                          >  1:249850/1‑238 (MQ=255)                                                                                             aaccaccacgAGGACTGGGTCCGAGGAACGATTTTTGGGGTTATTCCATCGAGAGCCCCATTA‑CACCGAGAATCGACCCCAAACATCAATAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAG‑AAAAACACCTGACCCCATAATATTCCCCGACCCCATCAGGGTGA‑CGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAgg                                                                                                                                          <  1:654909/235‑1 (MQ=255)                                                                                                    accaGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTG‑CACCGATAATCGTCACAAACCATCCATAAATGG‑CAATCCCGCCAGTAGAGCCCAGGTATGCAGTTATGAG‑AAAAACAACTGACCCCCTAATATTTCCCGACACCATCAGGGTGA‑CGGGGATCAAGCCCACTTTCTGCGCATCAGCGTCCGAAGACCTAATAAAACTCCTCCAAAgacaaa                                                                                                                                      >  2:62170/1‑231 (MQ=255)                                                                                                     ccaGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTG‑CACCGATAATCGTCACCAACCATCCATAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAG‑AAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGGGA‑CGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTc                                                                                                                                  >  1:404362/1‑239 (MQ=255)                                                                                                     ccaGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTG‑CACCGATAATCGTCACCAACCATCCATAAATGG‑CAATCCCGCCAGTAGAAGCCAGGTTTGCAGGTAACAG‑AAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGA‑CGGGGCTTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTc                                                                                                                                  >  1:311016/1‑239 (MQ=255)                                                                                                                   gAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTG‑CACCGATAATCGTCACCAACCATCCATAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAG‑AAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGA‑CGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAAtt                                                                                                                     >  1:789492/1‑238 (MQ=255)                                                                                                                    aGGCACGACTTTTTGGCGTATCCAATCGAGAGCCCC‑GTGCCCCCGATAATCGTCACCAACCAGCCATAAAGGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAG‑AAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGA‑CGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTa                                                                                                                    <  1:84868/238‑1 (MQ=255)                                                                                                                    aGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTG‑CACCGATAATCGTCACCAACCATCCATAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAG‑AAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGA‑CGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTAc                                                                                                                   >  1:640546/1‑239 (MQ=255)                                                                                                                    aGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTG‑CACCGATAATCGTCACCAACCATCCATAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAG‑AAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGA‑CGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAAACATAATTAAACCCCTGCGAAGGGGAGCTTCGTGACAAAAAATAAc                                                                                                                   >  2:519276/1‑239 (MQ=255)                                                                                                                      gaccGACTTTTGGGCATATACCATCGAGAGCCCGAGCG‑CACGGATAATCGTCCCCAACCACCCAAAAATGG‑CAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAG‑AAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGA‑CGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACg                                                                                                                  <  1:519276/235‑1 (MQ=255)                                                                                                                                      tataCCATCGAGAGCCCCAGTG‑CACCGATAATCGTCACCAACCATCCATAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAG‑AAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGA‑CGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCt                                                                                                   >  1:804804/1‑237 (MQ=255)                                                                                                                                       atcccATCGGGAGCCCGAGTG‑CACCGATAATCGTCCCCAACCATCCATAAATGG‑CAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAG‑AAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGA‑CGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCa                                                                                                <  1:883174/236‑1 (MQ=255)                                                                                                                                            tgggaagaccccATGG‑CACGAATAAACGTCAACAACCATCGATAAATGG‑CAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAG‑AAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGA‑CGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGa                                                                                           <  1:365163/231‑1 (MQ=255)                                                                                                                                             tCGAGAGCCCCAGTG‑CACCGCTAATCGTCACCAACCATCCATAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAG‑AAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGA‑CGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTCCGAAGGCGGACTTCGGGACAATAAATTACGGAATAACAAAAGCCTGCATTCgct                                                                                          >  2:165566/1‑237 (MQ=255)                                                                                                                                             tCGAGAGCCCCAGTG‑CACCGATAATCGTCACCAACCATCCATAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAG‑AAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGA‑CGGGGATTAAGCCCACTTTGTGAGCATCAGCACACGAAGACATATTTAAACTCCTGCGAAGCCGAGCTCCGCGACAAAAAATTACGTCATATCAAACCCCTTCATTCGac                                                                                          >  2:737706/1‑238 (MQ=255)                                                                                                                                               gaacacccaAGTG‑CCCCGATAACCGACCCCAACCATCCA‑AAATGGCCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAG‑AAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGA‑CGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGa                                                                                           <  1:100297/231‑1 (MQ=255)                                                                                                                                                agagCCCCAGTG‑CACCGATAATCGTCACCAACCATCCATAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAG‑AAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGAGA‑CGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGGGACAATAAATTACGTCATATAAAACGCCTGCATTCGAAAGt                                                                                       >  2:825083/1‑239 (MQ=255)                                                                                                                                                  agCCCCAGTG‑CACCGATAATCGTCACCAACCATCCATAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAG‑AAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGA‑CGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGGGACAATAAATTACGTCATATCATACGCCTGCATTCAATAGTTg                                                                                     >  1:614420/1‑239 (MQ=255)                                                                                                                                                  agCCCCAGTG‑CACCGATAATCGTCACCAACCATCCATAAATGG‑CAATCCCGCCAGTAGAGGCCAGGCTTGCAGGTAACAG‑AAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGA‑CGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCAAGCTTCGTGACAATAAATTACGTCAAATCATAACCCTGCATTCGAGAGTTg                                                                                     >  1:248091/1‑239 (MQ=255)                                                                                                                                                    ccccAGTG‑AACCGATAATCGTCACCAACCATCCATAAATGG‑CAATCCGGCCAGTAGAGACCAGTGTTGCAGGTAACCG‑AAAAACACCTGACCCCATAATATTCCCCGACACCATCAGAGCGA‑CGGGGATTAAGCCCCCTTGGTGAGCATCAGCATACGAAGACATAATTAAACTCCTGCGAAGGCGAGCTACGAGACAATAAATTACGTCATATCATACCCCCGCATTCGAGAGGTGaa                                                                                   >  2:707050/1‑239 (MQ=255)                                                                                                                                                         gTG‑CACCGATAATCGTCACCAACCATCCATAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAG‑AAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGA‑CGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTGGAattatt                                                                              >  1:189721/1‑239 (MQ=255)                                                                                                                                                             cACCGATAATCGTCACCAACCATCCATAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAG‑AAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGA‑CGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGCCAAATCATACGCCTGCATTCGATAGTTGAATTATtctc                                                                           >  1:285995/1‑239 (MQ=255)                                                                                                                                                               acaaTACTCCTCACCAA‑CAGCCATAAAAGGCCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAG‑AAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGA‑CGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGa                                                                         <  1:436378/236‑1 (MQ=255)                                                                                                                                                                ggatacgcGTCCCCAACCATCAAAAAAGGA‑CAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAG‑AAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGA‑CGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGAc                                                                        <  1:781683/232‑1 (MQ=255)                                                                                                                                                                  aTAATCGTCACCAACCATCCATAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAG‑AAAACCACCTGACCCCATAATATTCCCCGACACCATCAGGGTGA‑CGGGGATTAAGCCCACTTTTTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGCGCTTCGTGACAATAAATTACGTCATATCAAACGCCTGCATTTGATAGTTGAATTATTATCAACCg                                                                      >  2:320419/1‑239 (MQ=255)                                                                                                                                                                  aTAATCGTCACCAACCATCCATAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAG‑AAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGA‑CGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCCTACGCCTGCATTCGTTAGTTGAATTATTTTCGACCg                                                                      >  1:569294/1‑239 (MQ=255)                                                                                                                                                                  aTAATCGTCACCAACCATCCATAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAG‑AAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGA‑CGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCg                                                                      >  1:652141/1‑239 (MQ=255)                                                                                                                                                                     cactgcatcaaccaacAAAAAAGGG‑CAAGCCCTCCAGTAGAAGCCAAGGTTGCAGGTAACAG‑AAAAACCCCTGACCCCATAATATTCCCCGACACCATCAGGGTGA‑CGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCggtt                                                                   <  2:874935/231‑3 (MQ=255)                                                                                                                                                                       cGTCACCAACCACCAAAAAATGG‑CAACCCCGACAGTAGAGGCCAGGTTTGCAGGTAACAG‑AAAAACCCCTGACCCCATAATTATCCCCGGAACCATCAGGGTGA‑CGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTa                                                                 <  1:737706/239‑1 (MQ=255)                                                                                                                                                                          cACCAACCATCCATAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAG‑AAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGA‑CGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATc                                                               >  1:835739/1‑238 (MQ=255)                                                                                                                                                                          cACCAACCATCCATAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAG‑AAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGA‑CGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATc                                                               >  1:874202/1‑238 (MQ=255)                                                                                                                                                                                ccacccatAAAGGG‑CACTCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAG‑AAAAACCCCTGACCCCATAATATTCCCCGACACCATCAGGGTGA‑CGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATg                                                        <  1:880271/239‑1 (MQ=255)                                                                                                                                                                                  actcatAAAGGG‑GAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAG‑AAAAACCCCTGACCCCATAATATTCCCCGACACCATGAGGGTGA‑CGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGcc                                                      <  1:62170/236‑1 (MQ=255)                                                                                                                                                                                   tccatAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAG‑AAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGA‑CGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCAAACGCCTGCATTCGATAGTTGAATTATTATAAACCGGGTTAACGCGGATGcgg                                                     >  2:778663/1‑237 (MQ=255)                                                                                                                                                                                    ccatAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAG‑AAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGA‑CGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACACTAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGt                                                    >  1:467017/1‑239 (MQ=255)                                                                                                                                                                                     aaaaaaGGG‑CAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAG‑AAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGA‑CGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGtt                                                   <  1:129380/236‑1 (MQ=255)                                                                                                                                                                                       tAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAG‑AAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGA‑CGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGc                                                                                                     <  1:656882/187‑1 (MQ=255)                                                                                                                                                                                       tAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAG‑AAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGA‑CGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGc                                                                                                     >  2:656882/1‑187 (MQ=255)                                                                                                                                                                                       aaagtGG‑CAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAG‑AAAAACCCCTGACCCCATAATATTCCCCGACACCATCAGGGTGA‑CGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGttat                                                 <  1:432373/235‑1 (MQ=255)                                                                                                                                                                                        aaagtgcaaaTCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAG‑AAAAACCCCTGGCCCCATAATATTCCCCGACACCATCAGGGTGA‑CGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGttat                                                 <  2:789492/231‑1 (MQ=255)                                                                                                                                                                                        aaaacgcccacccccccAGTTGAGGCCAGGTTTGCAGGTAACAG‑AAAAACCCCTGACCCCATAATATTCCCCGACCCCATCAGGGTGA‑CGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGttat                                                 <  2:626282/228‑1 (MQ=255)                                                                                                                                                                                                  ccccGCCAGTAGAGGCCAGGTTTGCAGGTAACAG‑AAAAACCCCTGACCCCATAATATTCCCCGACACCAGCAGGGTGA‑CGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTTATTAATATAGct                                       <  1:400227/238‑1 (MQ=255)                                                                                                                                                                                                  ccccGCCAGTAGAGGCCAGGTTTGCAGGTAACAG‑AAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGA‑CGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTTATTAATATAGct                                       <  2:47268/238‑1 (MQ=255)                                                                                                                                                                                                            tagggCCAGGTTTGGAGGGAACAA‑AAAAAACCCTGACCCCATAATATTCCCCGACACCATCAGGGTGA‑GGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACCTAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTTATTAATATAGCTCTTCGtt                                <  2:231020/233‑1 (MQ=255)                                                                                                                                                                                                            agagGCCAGGTTTGCAGGTAACAG‑AAAAACCCCTGACCCCATAATATTCCCCGACACCATCAGGGTGA‑CGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTTATTAATATAGCTCTTCGTTaa                              <  1:825083/238‑1 (MQ=255)                                                                                                                                                                                                            aaagGCCAGGTTTGAAGGGAACAG‑AAAAACCCCTTACCCCATAATATTCCCCGACACCATCAGGGTGA‑AGGGGATTAAGCCCACTTTGTGAGCATCAGCAACCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTTATTAATATAGCTCTTCGtt                                <  2:640546/234‑1 (MQ=255)                                                                                                                                                                                                                gCCAGGTTTGCAGGTAACAG‑AAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGA‑CGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTTATTAATATAGCTCTTCGTTAAGGTTa                         >  1:783170/1‑239 (MQ=255)                                                                                                                                                                                                                gCCAGGTTTGCAGGTAACAG‑AAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGA‑CGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTTATTAATATAGCTCTTCGTTAAGGTTa                         >  1:536211/1‑239 (MQ=255)                                                                                                                                                                                                                     gTTTGCAGGTAACAG‑AAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGA‑CGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTAATTAATATAGTTCTTCGTTAAGGTTAAAAAt                    >  2:297319/1‑239 (MQ=255)                                                                                                                                                                                                                             gTAACAG‑AAAAACCCCTGCCCCCATAATTTTCCCCGACACCATCAGGGTGA‑CGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTTATTAATATAGCTCTTCGTTAAGGTTAAAAATCAGGCGAc            <  2:372009/239‑1 (MQ=255)                                                                                                                                                                                                                              tAACAG‑AAAAAACCCTGTCCCCATAATATTCCCCGACACCATCAGGGTGA‑CGGGGATTAAGCCCACTTTGTCAGCATCAGCATCCGAAGACATAATTAAACTCCTTCGAAGGTGAGATTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTTATTAATATAGCTCTTCGTTAAGGTTAAAAATCAGGCGAc            <  1:803486/238‑1 (MQ=255)                                                                                                                                                                                                                              tAACAA‑AAAACCCCCTGCCCCCATAATATTCCCCGACCCCATCAGGGTGA‑CGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTTATTAATATAGCTCTTCGTTAAGGTTAAAAATCAGGCg              <  2:404362/236‑1 (MQ=255)                                                                                                                                                                                                                                  cg‑AAAAACACCTGACCCCATAATAGTCCCCGACACCATCGGGGGTA‑CGGGGATTAAGCCCCCTTTGGGAGCATCATAATCATAAGACATAATTAAACACCTGGAAAGGCGAGCTTCTTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTTATTAATATAGCTCTTCGTTAAGGTTAAAAATCAGGCGACAGTGc       <  2:534440/238‑1 (MQ=255)                                                                                                                                                                                                                                  aG‑AAAACCCCCGGGCCCCAAAATATTCCCCCCCACCAACCGGGTGA‑AGGGGATTAAGCCCACTTTGTGGGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGGCAGTTTCGTGACAATTAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTTATTAATATAGCTCTTCGTTAAGGTTAAAAATCAGGCGACAGTGc       <  2:547296/239‑1 (MQ=255)                                                                                                                                                                                                                                  aG‑AAAAACCCCTGACCCCATAATATTCCCCGACACCATCAGGGTGA‑CGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTTATTAATATAGCTCTTCGTTAAGGTTAAAAATCAGGCGACAGTGc       <  1:153396/239‑1 (MQ=255)                                                                                                                                                                                                                                   aaaaaaaCCCCTGACCCCGTAAGATGTCCCGACGCCATCAGGGGGA‑AGGGGAGTAAGCCCACTTTGTGAGCCACAGCAGCCGAAAACATAATGAAACTCCTGCGAAGGCGAGCTTCGTGGCAATAAATTACGTCATATCATACGCCTGCATTCGATAGGTGAATTATTATCGACCGGGTTATCGCTGATGCCGTTATTAATATAGCTCTTCGTTAAGGTTAAAAATCAGGCGACAGTGc       <  2:329967/237‑1 (MQ=255)                                                                                                                                                                                                                                        aaCCCCTGACCCCAGAATTTTCCCCGACACCATCAGGGTGA‑CGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTTATTAATATAGCTCTTCGTTAAGGTTAAAAATCAGGCGACAGTGCGTtct  <  1:274779/239‑1 (MQ=255)                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                             CATGGCGATATTGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTG‑CACCGATAATCGTCACCAACCATCCATAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAG‑AAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGA‑CGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTTATTAATATAGCTCTTCGTTAAGGTTAAAAATCAGGCGACAGTGCGTTCT  >  NC_000913/4336720‑4337186

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