Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,859,421	T→C	23.8%	A777A (GCT→GCC)	mutS →	methyl‑directed mismatch repair protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,859,421	0	T	C	23.8%	66.4 / 3.0	44	A777A (GCT→GCC)	mutS	methyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base T (4/25);  new base C (8/3);  total (15/29)
Fisher's exact test for biased strand distribution p-value = 7.38e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGACCCAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGG‑TTGCAGCTC‑TGGCAGGCGTGCCAAAA‑GAGGTTATTAAGCG‑‑CGCACGGCAAAAG‑CTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACG‑CAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGCCTGG  >  NC_000913/2859203‑2859647                                                                                                                                                                                                                            |                                                                                                                                                                                                                                         cTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGACCCAGTTACCGGAGAAAATGGAAAGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCCAGCAAAAGCTACGGCCTGGCGG‑TTGGAGCCC‑CGGAAgga                                                                                                                                                                                                                                >  1:224364/1‑227 (MQ=255)   gTCGCTGGCGTGGGCGTGCGCGGAAAAGCTGGCGAATAAGATTAAGGAATTGACGTTATTTGCTACCCACTATTTCGAGCTGACCCAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGG‑TTGCAGCTC‑TGGCAGGCGt                                                                                                                                                                                                                              <  2:649536/228‑1 (MQ=255)    tCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGACCCAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGCGCAAAAGCTACGGCCTGGCGG‑GTGCACCCC‑TTGCCGGCggg                                                                                                                                                                                                                             >  2:74219/1‑226 (MQ=255)    tCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGACCCAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCCCGGCGAGCAAAAGCTACGGCCTGGCGG‑CTGTCGACC‑TGCCAGGCggg                                                                                                                                                                                                                             >  2:641842/1‑226 (MQ=255)         ggcgtgggcgtgCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGACCCAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGG‑TTGCAGCTC‑CTGCAGGCGTGCCaaa                                                                                                                                                                                                                        >  2:32365/1‑228 (MQ=255)            gtggggtgggccgGAAAATTGGGGCAAAAAGGTTAAGGGATTGAAGTTTTTTGCTGCCCATTATTGCGAGCTGACCCAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGG‑TTGCAGCTC‑TGGCAGGCGTGCCAAAa‑g                                                                                                                                                                                                                     <  1:186266/216‑1 (MQ=255)            gtgggcgtgCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGACCCAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGG‑TTGCAGCTC‑TGGCAGGCGTGCCAAAa‑g                                                                                                                                                                                                                     <  2:197878/227‑1 (MQ=255)               ggcgtgCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTCTTTGCTACCCACTATTTCGAGCTGACCCAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGAGCATCTCGATGCACTGGAGCACGGCGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGG‑TTTTAGCTC‑TGGCAGGCGGGCCcaaaagggg                                                                                                                                                                                                                  >  2:431778/1‑223 (MQ=255)                                               aGGCTTTGACGTTATTTGCTACCCCCTATTTCGAGCTGACCCATGTACCGGAGAAAATGAAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGG‑TTGCAGCTC‑TGGCAGGCGTGCCAAAA‑GAGGTTATTAAGCG‑‑CGCACGGCAAAAG‑CTGCGTGAg                                                                                                                                                                               <  1:412/227‑1 (MQ=255)                                                 gCATTGAAGTTATTTGCTACCCACTATTTCGAGCTGACCCAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGG‑TTGCAGCTC‑TGGCAGGCGTGCCAAAA‑GAGGTTATTAAGCG‑‑CGCACGGCAAAAG‑CTGCGTGAGCTg                                                                                                                                                                            <  2:483818/228‑1 (MQ=255)                                                    ttGACGTTATTTGCTACCCACTATTTCGAGCTGACCCAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGG‑GTGCAGCCC‑TGGCAGGCGTGGCCAAA‑GAGGGTATTAAGCG‑‑CGCCCGGCCAAAG‑CTGCCGGAGCCGGaa                                                                                                                                                                         >  2:142241/1‑228 (MQ=255)                                                     tGACGTTATTTGCTACCCACTATTTCGAGCTGACCCAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTTGCGGC‑C‑TGGGCGGCGTGCCCAAA‑GAGGTTCTTAAACG‑‑GGCACGGCACAAG‑CTTCGTGAGGTGGaa                                                                                                                                                                         >  2:408478/1‑227 (MQ=255)                                                     tGACGTTATTTGCTACCCACTATTTCGAGCTGACCCAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGG‑TTGCAGCCC‑TGGCCGGGGTGCCAAAA‑GAGGGTATTAACCG‑‑GGCACGGCCAAAG‑AGGCGGGGGCTGGgaa                                                                                                                                                                        >  2:138380/1‑228 (MQ=255)                                                     tGACGTTATTTGCTACCCACTATTTCGAGCTGACCCAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGG‑ATGCCGGTCTTGCCGGGGGGGCCAAAA‑G‑GGGTATTAAGAG‑‑GGCCCGGCAAAAA‑Ccgcggggggtgggaa                                                                                                                                                                        >  1:614693/1‑217 (MQ=255)                                                     tGACGTTATTTGCTACCCACTATTTCGAGCTGACCCAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGACACCATTGCC‑TTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGG‑TTGCAGCCC‑TGGCAGGCGTGCCAAAA‑ACGGGTATTAAAAG‑‑GGCCCGGCAAAAG‑ATGCGGGGGGCGGAAAg                                                                                                                                                                       >  2:7257/1‑228 (MQ=255)                                                                        cTATGTCGAGCTGACCCAGTTTCCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGG‑TTGCAGCTC‑TGGCAGGCGTGCCAAAA‑GAGGTTATTAAGCG‑‑CGCACGGCAAAAG‑CTGc                                                                                                                                                                                    <  2:292518/197‑1 (MQ=255)                                                                             tCGAGCTGACCCAGGTACCGGAGAAAATGGAAGGCGTCGCTAACGTTCATCTCGATGCACTGGAGCACGGCGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGG‑TTGCAGCTC‑TGGCAGGCGTGCCAAAA‑GAGGTTATTAAGCG‑‑CGCACGGCAAAAG‑CTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTa                                                                                                                                                <  2:130502/228‑1 (MQ=255)                                                                                aGCTGACCCAGTTACCGGAGAAAATGAAAGGCGTCGCTAACGCGCATCTCGATGAACTGGAGCACGGCGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGGAAAAGCTACGGCCTAGCGG‑TTTCCGCCC‑TCGCAGGCGTGGCAAAA‑AAGGTGATTAAACGGCCGCCTGGCAAGA‑‑TTGC‑AGAGCTGGAAAACAATTCGCCGagaccaccccgaccag                                                                                                                                             >  1:217824/1‑213 (MQ=255)                                                                                     aCCCAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGCGCAAAAGCTACGGCCTGGCGG‑TTTCAGCTC‑TGGCCGGCGGGCCCAAA‑GGGGGTATTAAGCG‑‑CGCACGGCCAAAGACTGTGGGCGCT‑GAAAGCACTTCGCCGAACGCGGCCGCTACG‑Cagggg                                                                                                                                       >  2:3790/1‑224 (MQ=255)                                                                                     aCCCAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGG‑GTGCAGCGC‑TGGCGAGCGCGGCCAAA‑AAGGGTATTAAGCG‑‑CGCCCGGCGAAAA‑CGTCGTCGGCTGGAAAAGACTTTGCCGAGcaccgccccacaccaagag                                                                                                                                        >  1:413736/1‑216 (MQ=255)                                                                                       ccAGTGACCGGAGAAAATGGAAGGCTCGGCTTACGGGCAGCGCGATGCACGGGGTCAGGGCGCAACCGTTGCCTTGATGGCCAGCGGGCAGGATGGGGCTGCGAGCAAAAGCTTCGGCCTGGCGG‑TTGCAGCTC‑TGGCAGGCGTGCCAAAA‑GAGGTTATTAAGCG‑‑CGCACGGCAAAAG‑CTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACG‑CAAGTGGa                                                                                                                                     <  1:590296/228‑1 (MQ=255)                                                                                       ccAGTGACCGGAGAAAATGGAAGGCGTAGCTAACGTGCATCCAGATGCACTGGAGCCCGGCGACCCCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGG‑TTGCAGCTC‑TGGCAGGCGTGCCAAAA‑GAGGTTATTAAGCG‑‑CGCACGGCAAAAG‑CTGCGTGAGATGGAAAGCATTTCGCCGAACGCCGCCGCTACG‑CAAGTGGa                                                                                                                                     <  2:130340/228‑1 (MQ=255)                                                                                         aGTGACCGGAAAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGACACCATTGCCTTTATGCCCAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGG‑TTGCAGCTC‑TGGCAGGCGTGCCAAAA‑GAGGTTATTAAGCG‑‑CGCACGGCAAAAG‑CTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACG‑CAAGTGGATg                                                                                                                                   <  2:429132/228‑1 (MQ=255)                                                                                           ttACCGGAAAAAATGGAGGGCGTCGCTAAAGTGCCTCCCGATGCCCTGGGGCAGGGCGACAACATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGG‑TTGCAGCTC‑TGGCAGGCGTGCCAAAA‑GAGTTTATTAAGCG‑‑CGCACGGCAAAAG‑CTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACG‑CAAGTGGATGGt                                                                                                                                 <  1:525350/228‑1 (MQ=255)                                                                                              cgggAGAAAATGGAGGGCGTGCCTAACGGTCATCTCGAGTCACTGGGCCACGGCGACCCCATTGCCTTTATGCACAGCGTGCAGGATGGGGGGGCGAGCAAAAGCTACGGCCTGGCGG‑TTGCAGCTC‑TGGCAGGCGTGCCAAAA‑GAGGTTATTAAGCG‑‑CGCACGGCAAAAG‑CTGGGTGAGCTGGAAAGCATTTTGCCGAACGCCGCCGCTACG‑CAAGTGGATGGTACg                                                                                                                              <  1:87347/226‑1 (MQ=255)                                                                                                  agaAAATGGAAGGGGGCTCTAACGTGCATCTCGATGCACTGGAGCACGGGGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGG‑TTGCAGCTC‑TGGCAGGCGTGCCAAAA‑GAGGTTATTAAGCG‑‑CGCACGGCAAAAG‑CTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACG‑CAAGTGGATGGTACGCaaa                                                                                                                          <  2:158986/228‑1 (MQ=255)                                                                                                  aaaaaaTGGAAGGGGTCGCTAAGGGGCATCGGGATGCACTGGAGCCCGGCGACACCATTGCCTTTATGCACGGCGTGCAGGATGGCGCGGGGAGCAAAAGCTACGGCCTGGCGG‑TTGCAGCTC‑TGGCAGGCGTGCCAAAA‑GAGGTTATTAAGCG‑‑CGCACGGCAAAAG‑CTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACG‑CAAGTGGATGGTACGCaaa                                                                                                                          <  1:142241/226‑1 (MQ=255)                                                                                                                gTCGCTAATGGGCCTCGCGGGGCCCTGGGGCAGGGGCACCCCATTGCCTTTAGGCACAGCGTGCAGGATGGGGGGGCGAGCAAAAGCTCCGGCCTGGCGG‑TTGCAGCTC‑GGGCAGGCGTGCCAAAA‑GAGGTTATTAAGCG‑‑CGCACGGCAAAAG‑CTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACG‑CAAGTGGATGGTACGCAAATGTCTTTGCTGTc                                                                                                             <  1:641842/227‑1 (MQ=255)                                                                                                                          tGCATCTCGATGCACTGGAGCACGGCGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGG‑TTGCAGCTC‑TGGCAGGCGTGCCAAAA‑GAGGTTATTAAGCG‑‑CGCACGGCAAAAG‑CTGCGTGAGCTGGAAAGCATTTCGCCGAAcgccgc                                                                                                                                                     <  2:453792/178‑1 (MQ=255)                                                                                                                            cATCTCGATGCACTGGAGCACGGCGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGG‑TTGCAGCGC‑TTGCAGGCGGGGCCAAA‑GAGGTTATTAAGAG‑‑CGCCCGGGCAAAG‑CGGCGGGAGCTGGAAAACAATTGGCCCAACGCCGCCCCGACC‑ACAGTGGATGGGAAGCGAAAGTCTGTGCTGTGAGTACCAagagaa                                                                                                >  1:415261/1‑228 (MQ=255)                                                                                                                            cATCTCGATGCACTGGAGCACGGCGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGG‑TTGCAGCCC‑TGGCAGGCGTGCCAAAA‑GAGGGGATTAAGCG‑‑CGCCCGGCCAAAA‑GTGCGTGAGCTGGAAAGGATTTTTCCCAAAGGCGCCCGCACGCCAGGGGGGTGGAAGGCAAAGGTCTTGGTTGCCAGTACCaaaaga                                                                                                 >  1:128543/1‑228 (MQ=255)                                                                                                                                      gatgggggCACGGCGACACCATTGCCTTTATGCAAAGTGTGCAGGATGGCGGGGCGAGCAAAAGCTACGGCCTGGCGG‑TTGCAGCTC‑TGGCAGGCGTGCCAAAA‑GAGGTTATTAAGCG‑‑CGCACGGCAAAAG‑CTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACG‑CAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCt                                                                                       <  2:95010/220‑1 (MQ=255)                                                                                                                                          ggggtagggCGACCCCATTGCCTTGTGGCAAGGGGGGCAGGAGGGGGGGGCGAGCAAAAGCTGCGGCCTGGCGG‑TGGCAGCCC‑TGGAGGGCGTGCCAAAA‑GAGGTTTTTAAGCG‑‑CGAACGGCAAAAG‑ATGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACG‑CAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGt                                                                                  <  1:623027/221‑1 (MQ=255)                                                                                                                                                  gCGACACCTTTGCCTTTATGGACGGTGTTCAGGGTGGCGCGGCGAGAAAAAGCTACGGCCGGGCGG‑TTGCAGCCC‑CGGCAGGGGTGCAAAAA‑GAGGTTATTAAGCG‑‑CGCTCGGCAAAAG‑CTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACG‑CAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAGCTGTGCCTGCGGTCGAAGct                                                                           <  1:65802/227‑1 (MQ=255)                                                                                                                                                  gCGACACCATTGCCTTTAGGCACAGCGTGCAGGATGGCCCGGCGAGCAAAAGCTACGGCCTGGCGG‑TTGCAGCCC‑TGGCAGGCGTGCCAAAA‑GAGGTTATTAAGCG‑‑CGCACGGCAAAAG‑CTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACG‑CAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGct                                                                           <  2:508564/227‑1 (MQ=255)                                                                                                                                                                              gCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGG‑TTGCAGCTC‑TGGCAGGCTTGCCAAAA‑GAGGTTATTAAGGG‑‑CGCACGGCAAAAG‑CTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACG‑CAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATtcactcac                                              <  2:636269/228‑1 (MQ=255)                                                                                                                                                                                            gatgCAAAAGATACGGCTGGGCGG‑TTGCGGCTT‑TGGCAGG‑GGGCCAAAACGAGGTTATTAGGGC‑‑GGCACGGCAAAGG‑CTGGGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACG‑CAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGgcg                                 <  1:32365/224‑1 (MQ=255)                                                                                                                                                                                             taGCAAAACCTACTGCCTGGCGT‑TTGCAGCGC‑TGGCAGGCGTGCCAAAA‑GAGGTTGTTAAGCC‑‑CGCACGGCAAAAG‑CTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACG‑CAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGgcgc                                <  2:81282/226‑1 (MQ=255)                                                                                                                                                                                                 aaaaGCTTCGGCCGGGCGG‑TTGCAGCTC‑TGGCGGGCGGGCAAAAA‑GAGGTTATTAAGCG‑‑CCCACGGCAAAAG‑CTGCGTGAGCTGGAAAGCTTTTCGCCGAACGCCGCCGCTACG‑CAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGa                            <  1:408478/227‑1 (MQ=255)                                                                                                                                                                                                         cGGCTTGGCGG‑TTGCAGCTC‑TGGCGGGCGTTCCAAAG‑GAGGTTAATAAGCG‑‑CGCACGGCAAAAG‑CTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACG‑CAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGAttt                    <  2:523813/227‑1 (MQ=255)                                                                                                                                                                                                                  gg‑TGGCAGCTC‑GGGCAGGCGTGCCAAAA‑GAGGTTATTAAGCG‑‑CGCACGGCAAAAG‑CTGCGTGAGCTGAAAAGCATTTCGCCGAACGCAGCCGCTACG‑CAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGaa          <  1:619422/228‑1 (MQ=255)                                                                                                                                                                                                                   g‑TTGCAGCTC‑TGGCAGGCGTGCCAAAA‑GAGGTTATTAAGCG‑‑CGCACGGCAAAAG‑CTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACG‑CAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGaa          <  2:436153/227‑1 (MQ=255)                                                                                                                                                                                                                        agcgtc‑tTGCAGGGGTTCCAAAA‑CGGGTTATGAACCG‑‑GCCACGGAAAAAG‑CGGCGTGAGGGGTAAAGCATTTCGCCGAACGCCGCCGCTTCG‑CAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGcc     <  1:431778/224‑1 (MQ=255)                                                                                                                                                                                                                         agttc‑tGGCAGGCGTGCAAAAA‑GAGGTTATTAAGCG‑‑CGCACGGCAAAAG‑CTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACG‑CAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGcc     <  2:314402/224‑1 (MQ=255)                                                                                                                                                                                                                         aGCTC‑TGGCAGGCGTGCCAAAA‑GAGGTTATTAAGCT‑‑GGCACGGCAAAAG‑GTGCGTGAGCTGGAAAGCATTGCGCCGAACGCCGCCGCTACG‑CAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGCCt    <  2:527785/228‑1 (MQ=255)                                                                                                                                                                                                                          gCTC‑TGGCAGGCGTGCCAAAA‑GAGGTTATTAAGCG‑‑CGCACGGCAAAAG‑CTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACG‑CAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGCCt    <  2:415261/227‑1 (MQ=255)                                                                                                                                                                                                                           ctc‑tGGCAGGCGTGCCAAAA‑GAGGTTATTAAGCG‑‑CGCACGGCAAAAG‑CTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACG‑CAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGAATTCACTCCCCCCGCGCCAGGCGCTGGGGAGGATTTATCGCTTTaaagacctgg  >  1:314402/1‑220 (MQ=255)                                                                                                                                                                                                                            |                                                                                                                                                                                                                                         CTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGACCCAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGG‑TTGCAGCTC‑TGGCAGGCGTGCCAAAA‑GAGGTTATTAAGCG‑‑CGCACGGCAAAAG‑CTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACG‑CAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGCCTGG  >  NC_000913/2859203‑2859647

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