Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,303,241	A→G	16.7%	S85S (TCT→TCC)	napF ←	ferredoxin‑type protein, role in electron transfer to periplasmic nitrate reductase NapA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,303,241	0	A	G	16.7%	175.0 / 5.3	78	S85S (TCT→TCC)	napF	ferredoxin‑type protein, role in electron transfer to periplasmic nitrate reductase NapA
Reads supporting (aligned to +/- strand):  ref base A (21/41);  new base G (12/3);  total (34/44)
Fisher's exact test for biased strand distribution p-value = 2.64e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGAGATACTCCGCAGTTATGGCTGATACCGGACAGCTGGCTGCACATGCTCCACATCCGTTGCAAAGTTGGCTATTAAGTTGCGGCTGGTAGATCCCGGAGAGCGTCGGGCGAAAGATAATTGCCATTGGCTCACAGCTATCCTGACAGCGGCGGCATTCAACTGACTGATACGCCAGGCACGCGTCCCCGATGGTGAACTGC‑AAATCCCAGG‑CCCTGGT‑GTGGCGCGGAG‑AAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCA‑GCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGG  >  NC_000913/2303011‑2303478                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                  tGGAGATCCTCCGGAGTTATGGCTGATACCGGACAGCTGGCTGAACATGCTCCACATCCGTTGCAAAGTTGGCTATTAAGTTGCGGCTGGTAGATCCCGGAGAGCGTCGGGCGAAAGATAATTGCCATTGGCTCACAGCTATCCTGACAGCGGCGGCATTCAACTGACTGATACGCCAGGCACGCGTCCCCGATGGTGAACTGC‑AAATCCCAGG‑CCCTGGT‑GTGGCGCGGAG‑AAAATAg                                                                                                                                                                                                                                          <  1:110317/239‑1 (MQ=255) tGGAGATACTCCGCAGTTATGGCTGATACCGGACAGCTGGCTGCACATGCTCCACATCCGTTGCAAAGTTGGCTATTAAGTTGCGGCTGGTAGATCCCGGAGAGCGTCGGGCGAAAGATAATTGCCATTGGCTCACAGCTATCCTGACAGCGGCGGCATTCAACTGACTGATACCCCAGGCCCCCGTCCCCCATGGCGAACTGC‑CAATCCCAGG‑CCCTGGG‑GTGGCGCGGGG‑gaaaat                                                                                                                                                                                                                                           >  1:195014/1‑236 (MQ=255)     gaTACTCCGCAGTTATGGCTGATACCGGACAGCTGGCTGCACATGCTCCACATCCGTTGCAAAGTTGGCTATTAAGTTGCGGCTGGTAGATCCCGGAGAGCGTCGGGCGAAAGATAATTGCCATTGGCTCACAGCTATCCTGACAGCGGCGGCATTCAACCGACTGATACCCCAGGGACCGGTCCCCGATGGTGAACTGC‑CAACCCCAGG‑CCCCTGG‑GTGGCGCGGCG‑CAAATAACGAt                                                                                                                                                                                                                                      >  1:166785/1‑239 (MQ=255)       tACTCCCCAGTTATGGCTGATACCGGACAGCTGGCTGCACATGCTCCACATCCGTTGCAGAGTTGGCTATTAAGTTGCGGCTGGTAGATCCCGGAGAGCGTCGGGCGAAAGATAATTGCCATTGGCTCACAGCTATCCTGACAGCGGCGGCATTCAACTGACTGATACGCCAGGCACGCGTCCCCGATGGTGAACTGC‑AAATCCCAGG‑CCCTGGT‑GTGGCGCGGAG‑AAAATAGCGAtt                                                                                                                                                                                                                                     <  1:143460/238‑1 (MQ=255)       tACCCCGCAGTTAGGGTGGTTCCCGGACGGCTGGCGGCCTTAGCTCCGACGCCGGTGCAAAGTTGGCTATTAATGTGCGGCTGGTAGATCCTGGAGAGCGTCGGGCGAAAGATAATTGCCATTGGCTCACAGCTATCCTGACAGCGGCGGCATTCAACTGACTGATACGCCAGGCACGCGTCCCCGATGGTGAACTGC‑AAATCCCGGG‑CACTGGT‑GTGGCGCGGAG‑AAAATAGCGAtt                                                                                                                                                                                                                                     <  2:517099/238‑1 (MQ=255)         cTCCGCAGTTATGGCTGATACCGGACAGCTGGCTTCACATGCTCCACATCCGTTGCAAAGTTGGCTATTAAGTTGCGGCTGGTAGATACCTGAGAGCGTCCGGCGCAAAATAAATGCCATTTGCTCCCAGCTATCCTGGCAGCGGCGGGCATCCACCGGCTTGTAACCCCGGCACGCGTTCCCCATGGTGAACTGC‑ACATCCCCGGCCCCGGGT‑GTGGCGGGGGC‑aaaaaacggatgcg                                                                                                                                                                                                                                   >  2:508492/1‑229 (MQ=255)          tCCGCAGTTATGGCTGATACCGGACAGCTGGCTGCACATGCTCCACATCCGTTGCAAAGTTGGCTATTAAGTTGCGGCTGGTAGATCCCGGAGAGCGTCGGGCGAAAGATAATTGCCATTGGCTCACAGCTATCCTGACAGCGGCGGCATTCAACTGACTGATACGCCAGGCACGCGTCCCCGATGGTGAACTGC‑AAATCCCAGG‑CCCTCGG‑GTGGGCCGGGG‑AAAAAACGGCTTCggg                                                                                                                                                                                                                                 >  1:378055/1‑239 (MQ=255)            cgcgtttttggcttaacccgggaaGCTGGCTGCATAGGCCCCACATCCGTTGTAAAGTTGGATATTAAGTTGCGGCTGGGAGATCCCGGAGAGCGTCGGGCGAAAGATAATTGCCATTGGCTCACAGCTATCCTGACAGCGGCGGCATTCAACTGACTGATACGCCAGGCACGCGTCCCCGATGGTGAACTGC‑AAATCCCAGG‑CCCTGGT‑GTGGCGCGGAG‑AAAATAGCGACTCGGGGc                                                                                                                                                                                                                               <  2:758203/216‑1 (MQ=255)            cGCAGTTATGGCTGATACCGGACAGCTGGCTGCACATGCTCCACATCCGTTGCAAAGTTGGCTATTAAGTTGCGGCTGGTAGATCCCGGAGAGCGTCGGGCGAAAGATAATTGCCATTGGCTCACAGCTATCCTGCCAGCGGCGGCATTCAACTGACTGATACGCCAGGCACACCTCCCCGATGGTGAACTGC‑AAATCCCCGA‑CCCTGGT‑TTGGGGCGGAG‑AAAATAAGGGTTCggggt                                                                                                                                                                                                                               >  1:117598/1‑238 (MQ=255)              cAGTTATGGCTGATACCGGACCGCTGGCCGCACATGCTCCACCTCCGTTGCAAAGTTGGCTATTAAGTTGCGGCTGGTAGATCCCGAAGCGCGTCGGGCGAAAGATAATTGCCATTGGCCCCCAGCCATCCTCACAGCGGCGGCACTTCACTGACCTATACGCCCGGCACCCCTTCCCGACTGCGAGATGT‑CAATCCCCTG‑CCCTGCG‑GTGCGGCGGGG‑aaaacaacagtctcgggca                                                                                                                                                                                                                              >  2:145569/1‑223 (MQ=255)              cAGTTATGGCTGATACCGGACAGCTGGCTGCACATGCTCCACATCCGTTGCAAAGTTGGCTATTAAGTTGCGGCTGGTAGATTCCGGAGAGCCGCGGGCGGAAGATACTTGCCCTTGGCTCACAGCCATCCTGACAGCGGCGCGATTTAAATGACGCATACGCCAGGCCACCGTCCCCCCCGGGTAACGCC‑AAAACCCCGG‑ACCTGGTGGGGGCGGGGGG‑CAAA‑AGCGAGTCgggggca                                                                                                                                                                                                                             >  2:687935/1‑236 (MQ=255)                   aTGGCTGATACCGGACAGCTGGCTGCACATGCTCCACCTCCGTTGCAAAGTTGGCTATTAAGTTGCGGCTGGTAGATCCCGGAGAGCGTCGGGCGCAAGATAATTTCCATTGGCTCCCAGCTATCCTGACAGCGGCGGCATTCAACTGCCTGGTACCCCCGGCCCCCGTCCCCGCTGGGGAACTGC‑AAATCCCCGC‑CCCTGTG‑GTGGCGCGGGG‑caaaaaccggtctcgggcacaccc                                                                                                                                                                                                                         >  2:506929/1‑218 (MQ=255)                   aTGGCTGATACCGGACAGCTGGCTGCACATGCTCCACATCCGTTGCAAAGTTGGCTATTAAGTTGCGGCTGGTAGATCCCGGAGAGCGTCGGGCGAAAGATAATTGCCATTGGCTCACAGCCATCCTGACAGCGGGGGGCTTCAAATGACAGGTACGCCCGGGACCCCTCCCCCAGGGGGAAACGC‑CAATCCACGG‑CCCCGGT‑TGGGGGCGGGG‑GAAAAAGGGAATCGGGGGaaac                                                                                                                                                                                                                           >  2:519615/1‑234 (MQ=255)                    tGGCTGATACCGGACAGCTGGCTGCACATGCTCCACATCCGTTGCAAAGTTGGCTATTAAGTTGCGGCTGGTAGATCCCGGAGAGCGTCGGGCGAAAGATAATTGCCATTGGCGCACAGCTATCCTGACCACGGCGGCATTCAACTGCACGATACGGCCGGCAAGCGTCCCCCCTGGGGGACTGA‑AAAATCCAAG‑CCCCGGT‑GCGGCGCGGGG‑gaaaaacagattaggaggacacccgc                                                                                                                                                                                                                       >  2:73738/1‑217 (MQ=255)                         gATTCCGGACAGCTGGCTGCACATGCTCCACATCCGTTGCAAAGTTGGCTATTAAGTTGCGGCTGGTAGATCCCGGAGAGCGTCGGGCGAAAGATAATTGCCATTGGCTCACAGCTATCCTGACAGCGGCGGCATTCAACTGACTGATACGCCAGGCACGCGTCCCCGATGGTGAACTGC‑AAATCCCAGG‑CCCTGGT‑GTGGCGCGGAG‑AAAATAGCGATTCGGGGCAAGCCTgcgcgc                                                                                                                                                                                                                   <  1:832265/238‑1 (MQ=255)                         gATACCGGACAGCTGGCTGCACATGCTCCACATCCGTTGCAAAGTTGGCTATTAAGTTGCGGCTGGTAGATCCCGGAGAGCGTCGGGCGAAAGATAATTGCCATTGGCTCACAGCTATCCTGACAGCGGCGGCATTCAACTGACTGATACGCCAGGCACGCGTCCCCGATGGTGAACTGC‑AAATCCCAGG‑CCCTGGT‑GTGGCGCGGAG‑AAAATAGCGATTCGGGGCAAGCCTgcgcgc                                                                                                                                                                                                                   <  1:835611/238‑1 (MQ=255)                           tACCGGACAGCTGGCTGCACATGCTCCACATCCGTTGCAAAGTTGGCTATTAAGTTGCGGCTGGTAGATCCCGGGGAGCGCCGGGCGAAAAATAAATGCCATTGGCTCCCAGCTATCCTGCCAGCGGCGGCAGTCAACGGACTGATACGCCCCGGCCGCGGCCCCGAGGGGGGACTGC‑AAAACCCAGC‑CCCCGGG‑GGGGGGCGGAG‑CAAATAACGAATCGGGGGAAgacggccgggcgg                                                                                                                                                                                                                >  2:846262/1‑227 (MQ=255)                                gACAGCTGGCTGCACATGCTCCACATCCGTTGCAAAGTTGGCTATTAAGTTGCGGCTGGTAGATCCCGGAGAGCGTCGGGCGAAAGATAATTGCCATTGGCTCACAGCTATCCTGACAGCGGCGGCATTCAACTGACTGATACGCCAGGCACGCGTCCCCGATGGTGAACTGC‑AAATCCCAGG‑CCCTGGG‑GTGGCGCGGAG‑CAAATAGCGGTTCGGGGCAAGCCTGCGCCCAGGCGGaa                                                                                                                                                                                                           >  1:903733/1‑239 (MQ=255)                                     ctggctgCACATGCTCCACATCCGTTGCAAAGTTGGCTATTAAGTTGCGGCTGGTAGATCCCGGAGAGCGTCGGGCGAAAGATAATTGCCATTGGCTCACAGCTATCCTGACAGCGGCGGCATTCAACTGACTGATACGCCAGGCACGCGTCCCCGATGGTGAACTGC‑AAATCCCAGG‑CCCTGGT‑GTGGCGCGGGG‑AAAATAGCGATTCGGGGCAAGCCTGCCCGCAGACGTAACAGaa                                                                                                                                                                                                      >  1:853008/1‑239 (MQ=255)                                     ctggctgCACATGCTCCACATCCGCTGCAAAGTTGGCTATTAAGTTGCGGCTGGTAGATCCCGGAGAGCGTCGGGCGAAAGATAATTGCCATTGGGTAACAGCTATCCTGACAGCGGCGGCATTCAACTGACTGATACGCCAGGCCCGCGTCCCCGATGGGGAACTGC‑AAAACCCAAG‑CCCCGGG‑GTGGCGCGGGG‑AAAAAAACGGATCTGGGGACGCCCGCGCGCggggggaaacaaa                                                                                                                                                                                                      >  2:886012/1‑229 (MQ=255)                                       ggttggccaTGCGCCCCATTCGTTGGAAGGTTGGCTTATAAGTGGCGGTGGGTAGATCCCGGAGAGCGTCGGGCGAAAGATAATTGCCAGTGGCTCACAGCTAGCCTGACAGCGGCGGCATTCAACTGACTGATACGCCAGGCACGCGTCCCCGATGGTGAACTGC‑AAATCCCAGG‑CCCTGGT‑GTGGCGCGGAG‑AAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGc                                                                                                                                                                                                    <  2:68497/232‑1 (MQ=255)                                       cggatatgtgtataagggacagTTGCAAAGTTGGCTATTAAGTTGCGGCTGGTAGATCCCGGAGAGCGTCGGGCGAAAGATAATTGCCATTGGCTCACAGCTATCCTGACAGCGGCGGCATTCAACTGACTGATACGCCAGGCACGCGTCCCCGATGGTGAACTGC‑AAATCCCAGG‑CCCTGGT‑GTGGCGCGGAG‑AAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAg                                                                                                                                                                                                     <  1:585377/217‑1 (MQ=255)                                               aTGCTCCACATCCGTTGCAAAGTTGGCTATTAAGTTGCGGCTGGTAGATCCCGGAGAGCGTCGGGCGAAAGATAATTGCCATTGGCTCACACCTATCCTGACAGCGGGGGCCTTCAACTGACTTATAACCCATGCCACCGTCCCCGCTGGTGAACTGCAAAACCCCAGC‑CCCGGCT‑GTGCCCGGGGC‑AAAATCGCG‑TTTGGGGCAAACCCGCGCGCAGGGGTAACACAAGccggacacg                                                                                                                                                                                            >  2:669711/1‑231 (MQ=255)                                                      acaTCCGTTGCAAAGTTGGCTATTAAGTTGCGGCTGGTAGATCCCGGAGAGCGTCGGGCGAAAGATAATTGCCATTGGCTCACAGCTATCCTGACAGCGGCGGCATTCAACTGACTGATACGCCAGGCACGCGTCCCCGATGGTGAACTGC‑AAATCCCAGG‑CCCTGGT‑GTGGCGCGGAG‑AAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACCCGTACtttt                                                                                                                                                                                     >  1:32815/1‑239 (MQ=255)                                                       caTCCGTTGCAAAGTTGGCTATTAAGTTGCGGCTGGTAGATCCCGGAGAGCGTCGGGCGAAAGATAATTGCCATTGGCTCACAGCTATCCTGACAGCGGCGGCATTCAACTGACTGATACGCCAGGCACGCGTCCCCGATGGTGAACTGC‑AAATCCCAGG‑CCCTGGT‑GTGGCGCGGGG‑AACATAGCGCTTCGGGGCAAGCCTGCGCGCAGGCGTAACCGAAACTGCAATCGTTAttttt                                                                                                                                                                                    >  1:434838/1‑239 (MQ=255)                                                             ttGCAAAGTGTGCTATTGAGTTGCGGCTGGTAAATCCCGGAGAGCTGCGGGCGAAAGATAATTGCCATTGGCTCACAGCCATCCTGACAGCGGCGGCATTCAACTTACTGATACGCCAGGCACGCGTCCCCGATGGTGAACTGC‑AAATCCCAGG‑CCCTGGT‑GTGGCGCGGAG‑AAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAAtt                                                                                                                                                                              <  2:34418/239‑1 (MQ=255)                                                                  agtgtGGCTATTAAGTGGCGGGTGGAACATCCCGGGGAGCGTCGGGCGAAAGAAAATTGCAGTTGGCTACCAGGTACCCTGACAGCGGCGGCATTCAACTGACTGATACGCCAGGCCCGCGTCCCCGATGGTGAACTGC‑AAATCCCAGG‑CCCTGGT‑GTGGCGCGGAG‑AAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGc                                                                                                                                                                       <  2:32815/235‑1 (MQ=255)                                                                    tttGGGTTTTTAGTGGCGGGTGGTAGATCCCGGGGAGGGTGGGGCGAAAAAGAATTGCCATTGGCTCACAGCTATCCTGACAGCGGCGGCATTCAACCGACTGAAACGCCAGGCACGCGTCCCCGATGGTGAACTGC‑AAATCCCGGG‑CCCTGGT‑GTGGGGCGGAG‑AAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTc                                                                                                                                                                     <  2:378055/238‑1 (MQ=255)                                                                        gCTATTAAGTTGCGGCTGGTAGATCCCGGAGAGCGTCGGGCGAAAGATAATTGCCATTGGCTCACAGCTATCCTGACAGCGGCGGCATTCAACTGACTGATACGCCAGGCACGCGTCCCCGATGGTGAACTGC‑AAATCCCAGG‑CCCTGGT‑GTGGCGCGGAG‑AAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCCCTCGTTATTTTTGAAATTAA‑‑CGCTCGGaa                                                                                                                                                                 >  1:767991/1‑238 (MQ=255)                                                                              aaGTTGCGGCTGTTAGATCCCGGAGGGCGTGGGGCGAAAGATAATTGCCATGGGCCCACAGCTATCCGGACAGCGGCGGCCTTCAACTGACTGATACGCCAGGCACGCGTCCCCGATGGTGAACTGC‑AAATCCCAGG‑CCCTGGT‑GTGGCGCGGAG‑AAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAgccgc                                                                                                                                                           <  2:903733/239‑1 (MQ=255)                                                                                      gCTGGTAGATCCCGGAGAGCGTCGGGCGAAAGATAATTGCCATTGGCTCACAGCTATCCTGACAGCGGCGGCATTCAACTGACTGATACGCCAGGCACGCGTCCCCGATGGGGCACTTC‑AAATCCCAGG‑CCCTGGT‑GTGGCGCGGAG‑AAAATAGCGCTTCGGGGCAAGCCTGCGCGCCGGCGTGACAGAAACTGCACTCGTTATTTTTTAAATTAA‑‑CGCCCGGATAACCTCccgcgccc                                                                                                                                                   >  2:98723/1‑238 (MQ=255)                                                                                             agtCCCGGAGACCGTCGGGCGAAAGTTAATTGCTATTGGCTCACAGCTCTCCTGACAGCGGCGGCATTCAACTGACTGATACGGCAGGCACGGGTCCCCGATGGTGAACTGC‑AAATCCCAGG‑CCCTGGT‑GTGGCGCGGAG‑AAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCa                                                                                                                                            <  2:166785/237‑1 (MQ=255)                                                                                                 ccGGCGAGGGTCGGGCGAAAAATAAGTGCCATTGGCTCAAAGCTTCCCTGACAGGGGCGGCATTCAACTGACTGATCCGCCAGGCACGCGTCCCCGAGTGTGAACTGC‑AAATCCCAGG‑CCCTGGT‑GTGGCGCGGAG‑AAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCccgcgccacaccagggcct                                                                                                                                        <  2:767991/239‑13 (MQ=255)                                                                                                         cGTCGGGCGAAAGATAATTGCCATTGGCTCACAGCTATCCTGACAGCGGCGGCATTCAACTGACTGATACGCCAGGCACGCGTCCCCGATGGTGAACTGC‑AAATCCCAGG‑CCCTGGT‑GTGGCGCGGAG‑AAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTCCAGAATGTTGttt                                                                                                                                 >  1:288956/1‑238 (MQ=255)                                                                                                         cGTCGGGCGAAAGATAATTGCCATTGGCTCACAGCTATCCTGACAGCGGCGGCATTCAACTGACTGATACGCCAGGCACGCGTCCCCGATGGTGAACTGC‑AAATCCCAGG‑CCCTGGG‑GTGGCGCGGAG‑AAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGttt                                                                                                                                 >  1:308166/1‑238 (MQ=255)                                                                                                            cGGGCGAAAGATAATTGCCTTTGGCTCACAGCTATCCGGACAGCGGCGGCATTCAACTGACTGATACGCCAGGCACGGGTCCCCGATGGTGAACTGC‑AAATCCCAGG‑CCCTGGT‑GTGGCGCGGAG‑AAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTca                                                                                                                              <  2:279962/238‑1 (MQ=255)                                                                                                            cGGGCGAAAGATAATTGCCATTGGCTCACAGCTATCCTGACAGCGGCGGCATTCAACTGACTGATACGCCAGGCACGCGTCCCCGATGGTGAACTGC‑AAATCCCAGG‑CCCTGGT‑GTGGCGCGGAG‑AAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTca                                                                                                                              <  1:519615/238‑1 (MQ=255)                                                                                                             ggggGAAAAATAATTGCCATGGGCCAACAGGTATCCGGCAAGCGGGGGCATTAAACTGACGGAGCCGCCGGGCACGCGTCCCCGATGGTGAACTGC‑AAATCCCAGG‑CCCTGGT‑GGGGCGCGGAG‑AAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTcaca                                                                                                                            <  2:6872/239‑1 (MQ=255)                                                                                                             gggCGAAAAATTATTGCCTTTGGCTCCCAGGTTCCCTGGCAGGGGCGGCATTCAACTGACTGATCCTCCCGGCAGGGGGCCCCGATGGTGAACTGC‑AAATCCCAGG‑CCCTGGT‑GGGGGGAGGAG‑AAAATAGCGATTCGGGGCAAGACTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTcaca                                                                                                                            <  2:541953/239‑1 (MQ=255)                                                                                                              ggCGAAAGATAATTGCCATTGGCTCACAGCTATCCTGACAGCGGCGGCATTCAACTGACTGATACGCCAGGCACGCGTCCCCGATGGTGAACTGC‑AAATCCCAGG‑CCCTGGT‑GTGGCGCGGAG‑AAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTcaca                                                                                                                            <  1:111253/238‑1 (MQ=255)                                                                                                              ggCGAAAGATAATTGCCATTGGCTCACAGCTATCCTGACAGCGGCGGCATTCAACTGACTGATACGCCAGGCACGCGTCCCCGATGGTGAACTGC‑AAATCCCAGG‑CCCTGGT‑GTGGCGCGGAG‑AAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTcaca                                                                                                                            <  1:846262/238‑1 (MQ=255)                                                                                                              ggCGAAAGATAAGTGCCATTGGCTCACAGCTATCCTGACAGCGGCGGCATTCAACTGACTGATACGCCAGGCACGCGTCCCCGATGGTGAACTGC‑AAATCCCAGG‑CCCTGGT‑GTGGCGCGGAG‑AAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTcaca                                                                                                                            <  1:277876/238‑1 (MQ=255)                                                                                                                  aaaGATAATTTCCATTGGCCCCCCGCTATCCTGACGGGGGCGGCATTTAACTGACTTATACGCCAGGCACGCGTCCCCGATGGTGAACTGC‑AAATCCCAGG‑CCCTGGT‑GTGGCGCGGAG‑AAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCAt                                                                                                                       <  2:770288/239‑1 (MQ=255)                                                                                                                   aaGATAAGTGCCATTGGCTCAAAGCTATCCTGACAGGGGCGGCATTCAACTGACTGATACGCCAGGCACGCGTCCCCGATGGTGAACTGC‑AAATCCCAGG‑CCCTGGT‑GTGGCGCGGAG‑AAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCAt                                                                                                                       <  1:99701/238‑1 (MQ=255)                                                                                                                             ccATTGGCTCACAGCTATCCTGACAGCGGCGGCATTCAACTGACTGATACGCCAGGCACTCGTCCCCGATGGTGAACTGC‑AAATCCCAGG‑CCCTGGT‑GTGGCGCGGAG‑AAAATAGCGATTCGGGGCAAGCCTGCCCGCAGGCGTAACAGAAGCTGCACTCTGTATTTTTGGAATTAA‑‑CGCTCGGATAGCCCCCCCAGCCACGATCCCGAATGGTGTTTTTACAACCAtacaaaaagcgc                                                                                                            >  2:587453/1‑228 (MQ=255)                                                                                                                                                 tgcaggcggcgGCTTTAAACGTACTGAGCCGCCAGGCACGCGTCCCCGATGGTGAACTGC‑AAAACCCAGG‑CCCGGGT‑GTGGCGCGGAG‑AAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCA‑GCGGGTGCAATGGGt                                                                                       <  2:434838/234‑1 (MQ=255)                                                                                                                                                 cgACAGCGGCGGCAATCAACTGACTGATACGCAAGGCACGCGTCCCCGATGGTGAACTGC‑AAATCCCAGG‑CCCTGGT‑GTGGCGCGGAG‑AAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCA‑GCGGGTGCAATGGGt                                                                                       <  1:669711/238‑1 (MQ=255)                                                                                                                                                   aCAGCGGCGGCATTCAACTGACTGATACGCCAGGCACGCGTCCCCGATGGTGAACTGC‑AAATCCCAGG‑CCCTGGT‑GTGGCGCGga                                                                                                                                                                                                                                                   >  2:873745/1‑85 (MQ=255)                                                                                                                                                   aCAGCGGCGGCATTCAACTGACTGATACGCCAGGCACGCGTCCCCGATGGTGAACTGC‑AAATCCCAGG‑CCCTGGT‑GTGGCGCGga                                                                                                                                                                                                                                                   <  1:873745/85‑1 (MQ=255)                                                                                                                                                   aCAGCGGCGGCATTCAACTGACTGATACGCCAGGCACGCGTCCCCGATGGTGAACTGC‑AAATCCCAGG‑CCCTGGT‑GTGGCGCGGAG‑AAAATACCGATTCGGGGCAAGCCTGCGCGCCGGCGCAACAGAAGTTGACTTCGTTATTTTTGCACTTAA‑‑CGCCCGGATAGCCGCCCCCGCCGCATTGCAGAATGGTGGTTTCACAAGCATAAATAAAGCCCTCCCA‑CCGGGGGAAAgggggaa                                                                                     >  2:692273/1‑236 (MQ=255)                                                                                                                                                       cggcggCATTCAACTGACTGATACGCCAGGCACGCGTCCCCGATGGTGAACTGC‑AAATCCCAGG‑CCCTGGT‑GTGGCGCGGAG‑AAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCA‑GCGGGTGCAATGGGTCAGaaa                                                                                 >  1:756085/1‑239 (MQ=255)                                                                                                                                                       cggcggCATTCAACTGACTGATACGCCAGGCACGCGTCCCCGATGGTGAACTGC‑AAATCCCAGG‑CCCTGGT‑GTGGCGCGGAG‑AAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGCCGCA‑GCGGGTGCAATTGGTCAGaaa                                                                                 >  1:541504/1‑239 (MQ=255)                                                                                                                                                               ttCAACTGACTGATACGCCAGGCACGCGTCCCCGATGGTGAACTGC‑AAATCCCAGG‑CCCTGGT‑GTGGCGCGGAG‑AAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCA‑GCGGGTGCAATGGGTCAGAAAATGAGAtt                                                                         <  1:3023/239‑1 (MQ=255)                                                                                                                                                                tCAACTGACTGATGCGCAAGGCACGCGTCCCCGATGGTGAACTGC‑AAATCCCAGG‑CCCTGGT‑GTGGCGCGGAG‑AAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCA‑GCGGGTGCAATGGGTCAGAAAATGAGATtc                                                                        <  1:98723/239‑1 (MQ=255)                                                                                                                                                                     tggctgaTCCGCCAGGCCCGCTTCCCCAAGGGGGAACGGC‑AAATCCCAGG‑CCCGGGT‑GTGGCGCGGAG‑AAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCA‑GCGGGTGCAATGGGTCAGAAAATGAGATTCATCAc                                                                   <  2:853008/236‑1 (MQ=255)                                                                                                                                                                     tgacggaTACGCCAGGCACGCGTCCCCGATGGTGAACTGC‑AAATCCCAGG‑CCCTGGT‑GTGGCGCGGAG‑AAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCA‑GCGGGTGCAATGGGTCAGAAAATGAGATTCATCAc                                                                   <  1:22775/239‑1 (MQ=255)                                                                                                                                                                           gtccGCCAGGCCCGCGTCCCCGAGGGTGAACTGC‑AAATCCCAGG‑CCCGGGG‑GTGGCGCGGGG‑AAAATAGCGGTTCGGGGCAAGCCTGCGCGCAGGGGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCA‑GCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGAcc                                                             <  2:288956/236‑1 (MQ=255)                                                                                                                                                                            tACGCCAGGCACGGGTCCCCGATGGTGAACTGC‑AAATCCCAGG‑CCCTGGT‑GGGGCGCGGAG‑AAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCA‑GCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGAcc                                                             <  1:178636/238‑1 (MQ=255)                                                                                                                                                                            tACCCCGGGCACGGCGCCCCGATGGTGAACTGG‑AAATCCCAGG‑CCCTGGT‑GTGGCGCGGAG‑AAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCA‑GCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGAcc                                                             <  2:308166/238‑1 (MQ=255)                                                                                                                                                                               gCCAGGCACGGGTCCCCGGTGGTGAACGGC‑AAATCCCAGG‑CCCTGGT‑GTGGCGCGGAG‑AAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCA‑GCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAggg                                                         <  1:145569/239‑1 (MQ=255)                                                                                                                                                                               ccctggccaggctCCCCGAGTGTGACCCGG‑AAACCCCAGG‑CCCGGGG‑GGGGCGCGGAG‑AAAAGAGCGTTTTGGGGCAAGCCGGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCA‑GCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAggg                                                         <  2:195014/227‑1 (MQ=255)                                                                                                                                                                                 cAGGCACGCGTCCCCGATGGTGAACTGC‑AAATCCCAGG‑CCCTGGT‑GTGGCGCGGAG‑AAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAACGCGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCA‑GCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAggg                                                         >  1:24863/1‑239 (MQ=255)                                                                                                                                                                                   ggCCCGGGTCCACGGTGGTGAACTGA‑AAATCCCGGG‑CCCTGGG‑TTGGCGCGGAG‑AAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCA‑GCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGcgga                                                     <  2:756085/239‑1 (MQ=255)                                                                                                                                                                                   ggCACGCGTCCCCGATGGTGAACTGC‑AAATCCCAGG‑CCCTGGT‑GTGGCGCGGAG‑AAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCA‑GCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGcgga                                                     <  1:886012/239‑1 (MQ=255)                                                                                                                                                                                   ggCACGCGTCCCCGATGGTGAACTGC‑AAATCCCAGG‑CCCTGGT‑GTGGCGCGGAG‑AAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCCCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACACGCATTAATACCGGCCTCGCACGCGGG‑GCAAAGGGGGAGAAAATGAGATTAAACAACGGACCAggggggg                                                     >  2:847033/1‑238 (MQ=255)                                                                                                                                                                                      ccgggcccccGATGGGGAACTTC‑AAATCCCAGG‑CCCTGGG‑GTGGCCGGGAG‑AAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCA‑GCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGcggacg                                                   <  2:916735/232‑1 (MQ=255)                                                                                                                                                                                      cacgttgcctgATGGTTAACCGA‑AACGCCCAGG‑CCCTGGT‑GTGGGGGGGGG‑AAAATGGCGATTCGTGTCAAGCCGGCGCGCAGGCGTAACGGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCA‑GCGGGTGCAATGGGTCAGAAAATGAGAATCATCGCCGGACCAGGGcggacg                                                   <  2:117598/228‑1 (MQ=255)                                                                                                                                                                                         cgTCCCCGATGGTGAACTGC‑AAATCCCAGG‑CCCTGGT‑GTGGCGCGGAG‑AAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGAATTAATACAGGCGTCGCA‑GCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGAt                                                >  1:916735/1‑238 (MQ=255)                                                                                                                                                                                               cGATGGTGAACTGC‑AAATCCCAGG‑CCCTGGT‑GTGGCGCGGAG‑AAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCCCA‑GCGGGGGCAATGGGTCAGAAAAATAGATTCATCACCGGCACAGGGCGGGCGGAAACCGtac                                         >  2:8446/1‑237 (MQ=255)                                                                                                                                                                                                                  cccAGG‑CCCTGGG‑GTGGCGCGGAG‑AAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCA‑GCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGa                       <  1:506929/239‑1 (MQ=255)                                                                                                                                                                                                                       g‑CCCTGGT‑GTGGCGCGGAG‑AAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGCTAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCA‑GCGGGCGCAATGGGTCAGAAAACGAGATTCATCACCGGACAAGGGGGGCAGGATACCGGTACGGGCTTTTGCCAAGCGACCAgg                  >  2:321830/1‑238 (MQ=255)                                                                                                                                                                                                                        gacagggT‑GTGG‑GCGGGGAAAAATTGCGATGCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCCTCACTTGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCA‑GCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCa                    <  2:424958/232‑1 (MQ=255)                                                                                                                                                                                                                            tGGT‑GTGGCGCGGAG‑AAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCA‑GCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGGGAGt              >  1:672530/1‑239 (MQ=255)                                                                                                                                                                                                                                 gtgGCGCGGAG‑AAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCA‑GCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATg           <  1:692273/238‑1 (MQ=255)                                                                                                                                                                                                                                   ggCGGGGAG‑AAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCA‑GCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGcc         <  1:226961/238‑1 (MQ=255)                                                                                                                                                                                                                                    gcgcGGAG‑AAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCA‑GCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGcccc       >  1:44138/1‑239 (MQ=255)                                                                                                                                                                                                                                    gcgcGGAG‑AAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCA‑GCGGGGGCAATGGGTCAGAAAATGAGATTCATAACCGGACCAGGGCGGACGAATACCGTTACTGGCTTTGCGCCACCGACCAGGGAGGAGGcccc       >  2:226961/1‑239 (MQ=255)                                                                                                                                                                                                                                      gcgGAG‑AAAATTGCTGTTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCATTCGTTATTTTTGAAATTAC‑‑ACCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCA‑GCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGa     <  2:617152/239‑1 (MQ=255)                                                                                                                                                                                                                                       cgGAG‑AAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCA‑GCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGAc    >  1:8231/1‑239 (MQ=255)                                                                                                                                                                                                                                       aggAG‑AAAATAGCCGTTGGGGGCAAGGCTGGGCGCAGGCGTAAGAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGGTTTCACAAGCATTAATACAGGCGTCGCA‑GCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGa     <  1:172287/237‑1 (MQ=255)                                                                                                                                                                                                                                         gag‑aAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGc                                                                                                                         <  2:416264/121‑1 (MQ=255)                                                                                                                                                                                                                                         gag‑aAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGc                                                                                                                         >  1:416264/1‑121 (MQ=255)                                                                                                                                                                                                                                         gag‑aAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACCGGCGTCGCA‑GCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGc          >  1:424958/1‑231 (MQ=255)                                                                                                                                                                                                                                         gag‑aAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCA‑GCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGCCGGATACCGTTACTGGCTTTGCGCCAGCGCCCAGTGAGTATGCCCCGACgg  >  1:84616/1‑239 (MQ=255)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                  TGGAGATACTCCGCAGTTATGGCTGATACCGGACAGCTGGCTGCACATGCTCCACATCCGTTGCAAAGTTGGCTATTAAGTTGCGGCTGGTAGATCCCGGAGAGCGTCGGGCGAAAGATAATTGCCATTGGCTCACAGCTATCCTGACAGCGGCGGCATTCAACTGACTGATACGCCAGGCACGCGTCCCCGATGGTGAACTGC‑AAATCCCAGG‑CCCTGGT‑GTGGCGCGGAG‑AAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCA‑GCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGG  >  NC_000913/2303011‑2303478

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