Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,104,042	T→G	L344V (TTA→GTA)	mutY →	adenine DNA glycosylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,104,042	0	T	G	100.0%	326.3 / NA	85	L344V (TTA→GTA)	mutY	adenine DNA glycosylase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (36/49);  total (36/49)

CCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGATTTTCAGCTGTACCCCGGCGAGCTGGGAAAACGCATCTATAAC  >  NC_000913/3103903‑3104188                                                                                                                                            |                                                                                                                                                    ccACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGGTACGCACTGgc                                                                                                                                          <  1:1065214/151‑1 (MQ=255)   aCTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGGTACGCACTGgcgc                                                                                                                                        <  1:428502/151‑1 (MQ=255)       aGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGGTACGCACTGGCGCGCCg                                                                                                                                    <  1:675193/151‑1 (MQ=255)          tatTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGGTACGCACTGGCGCGCCGGtt                                                                                                                                 >  1:35327/1‑151 (MQ=255)          tatTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGGTACGCACTGGCGCGCCGGtt                                                                                                                                 >  1:112180/1‑151 (MQ=255)              gtgCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGGTACGCACTGGCGCGCCGGTTTAg                                                                                                                              <  1:1110948/150‑1 (MQ=255)                 ccTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGGTACGCACTGGCGCGCCGGTTTAGcgcg                                                                                                                          >  1:460467/1‑151 (MQ=255)                  cTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGGTACGCACTGGCGCGCCGGTTTAGcgcg                                                                                                                          >  2:1004083/1‑150 (MQ=255)                     tgtgGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGGTACGCACTGGCGCGCCGGTTTAGCGCGTg                                                                                                                        >  2:187054/1‑149 (MQ=255)                      gtgGCTTCCCGTGTCGTCATTCCCCGGCTGCTTGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGGTACGCACTGGCGCGCCGGTTTAgcg                                                                                                                            <  1:541089/144‑1 (MQ=255)                      gtgGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGGTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑g                                                                                                                     <  1:760182/150‑1 (MQ=255)                       tgGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGGTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTc                                                                                                                   >  1:967971/1‑151 (MQ=255)                        ggCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGGTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTc                                                                                                                   <  1:543394/150‑1 (MQ=255)                        ggCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGGGACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTCg                                                                                                                  >  1:536811/1‑151 (MQ=255)                         gCTTCCCGTGTCGTCATTCACCGGCTGCAGGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGGTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTc                                                                                                                   <  1:798976/149‑1 (MQ=255)                             cccGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGGTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTCGATa                                                                                                               >  2:104556/1‑149 (MQ=255)                             cccGTGTCGTCATTCACCGGCTGCAGGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGGTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTCGAt                                                                                                                <  2:845614/148‑1 (MQ=255)                                gtgtCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGGTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTCGATAAaga                                                                                                           <  2:366566/150‑1 (MQ=255)                                gtgtCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGGTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTCGATAAag                                                                                                            >  2:869293/1‑149 (MQ=255)                                   tcgtcATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGGTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTCGATAAAGAGgat                                                                                                       >  2:708582/1‑151 (MQ=255)                                    cgtcATTCACCGGCTGCATGGGTGAAGGCAAGGCGTTGTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGGTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTCGATAAAGAGgat                                                                                                       <  1:652626/150‑1 (MQ=255)                                    cgtcATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGGTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTCGATAAAGAgg                                                                                                         <  2:253538/148‑1 (MQ=255)                                       cATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGGTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTCGATAAAGAGGATGAtt                                                                                                   >  1:20464/1‑151 (MQ=255)                                          tCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGGTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTCGATAAAGAGGATGATTTAt                                                                                                >  1:363128/1‑151 (MQ=255)                                           cACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGGTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTCGATAAAGAGGATGATTTATg                                                                                               >  1:242701/1‑151 (MQ=255)                                             ccGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGGTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTCGATAAAGAGGATGATTTAt                                                                                                <  2:1148171/148‑1 (MQ=255)                                             ccGCCTGCATGAATGAAGGCATTGCGCTCTGTTAAAACTTAGCGCATCCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGGTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTCGATAAAGAGGATGATTTATGAg                                                                                             <  2:977186/151‑1 (MQ=255)                                                gCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGGTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTCGATAAAGAGGATGATTTATGAGCAg                                                                                          <  1:537305/151‑1 (MQ=255)                                                     aTGGATGAAGGCAATGCGCTCTGGTATACCTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGGTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTCGATAAAGAGGATGATTTATGAGCAGAACg                                                                                      <  2:853507/150‑1 (MQ=255)                                                          tgcaGGCAATGAGCTCTGGTATAACATAGCGTAACCGCCGTCAGTTGGCCTAGCTGCTCCCGTGGAGCGTTTGTTACAGCAGGTACGCACTGGCGCGCCGGTATAGCGCGTGA‑GTCGATAAAGAGGCTGATTAATGAGCAGAACGAttttt                                                                                >  2:780957/4‑151 (MQ=255)                                                          tGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGGTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTCGATAAAGAGGATGATTTATGAGCAGAACGAtttt                                                                                 <  2:573214/150‑1 (MQ=255)                                                           gAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGGTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTCGATAAAGAGGATGATTTATGAGCAGAACGAttttt                                                                                <  2:86918/150‑1 (MQ=255)                                                            aaGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGGTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTCGATAAAGAGGATGATTTATGAGCAGAACGAtttttt                                                                               <  1:508809/150‑1 (MQ=255)                                                             aGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGGTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGt                                                                             >  1:289953/1‑151 (MQ=255)                                                               gCAATGCCCTCTGGTATAACTTAGCGCAGCCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGGTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGt                                                                             <  1:1155887/149‑1 (MQ=255)                                                                cAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGGTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACt                                                                          >  2:449994/1‑151 (MQ=255)                                                                 aaTGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGGTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACt                                                                          <  1:837217/150‑1 (MQ=255)                                                                  aTGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGGTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACttt                                                                        <  1:585678/151‑1 (MQ=255)                                                                  aTGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGGTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACttt                                                                        >  1:963417/1‑151 (MQ=255)                                                                  aTGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCAAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGGTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACttt                                                                        >  1:885744/1‑151 (MQ=255)                                                                    gcgcTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGGTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACttt                                                                        <  2:314141/149‑1 (MQ=255)                                                                     cgcTCTGGTCTAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGGTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTcc                                                                      >  2:567382/1‑150 (MQ=255)                                                                     cgcTCGGTTATACCTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGGTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCt                                                                     <  1:457446/151‑1 (MQ=255)                                                                       ctctGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGGTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGc                                                                   >  2:404120/1‑151 (MQ=255)                                                                        tctGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGGTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTg                                                                    <  1:301485/149‑1 (MQ=255)                                                                           ggTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGGTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAAc                                                                <  1:1090787/150‑1 (MQ=255)                                                                              ataACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGGTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTg                                                             <  1:239728/150‑1 (MQ=255)                                                                              ataACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGGTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTg                                                             >  2:771003/1‑150 (MQ=255)                                                                                  cTTAGCGCATCCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGGTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGc                                                         <  1:527984/150‑1 (MQ=255)                                                                                   ttAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGGTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGc                                                         >  2:875149/1‑149 (MQ=255)                                                                                    tAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGGTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAg                                                          >  1:701931/1‑147 (MQ=255)                                                                                    tAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGGTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGc                                                         >  1:124784/1‑148 (MQ=255)                                                                                    tAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGGTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCa                                                        <  1:1116750/149‑1 (MQ=255)                                                                                    tAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGGTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAg                                                       <  2:408457/150‑1 (MQ=255)                                                                                      gcgcAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGGTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAg                                                    <  1:1210921/151‑1 (MQ=255)                                                                                       cgcAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGGTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGa                                                      <  1:1046839/148‑1 (MQ=255)                                                                                         cAACCGCCGTCAGTTGGCCTAGCGGCTGCCGTGGAGCGTTTGTTACAGCAGGTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAg                                                    <  2:457044/148‑1 (MQ=255)                                                                                           aCCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGGTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTc                                                 >  2:1161187/1‑149 (MQ=255)                                                                                            ccgccgTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGGTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAg                                               <  2:289415/150‑1 (MQ=255)                                                                                            ccgccgTCAGTTGGCCGAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGGTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAgg                                              <  2:624057/151‑1 (MQ=255)                                                                                             cgccgTCAGTTGGCCTAGCTGCTCCCGTGGAGCGTTTGTTACAGCAGGTACGCTCTGGCGCGCCGGTTTAGCGCGTGC‑GTCGATAGAGAGGATGAGTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAg                                               <  2:62229/149‑1 (MQ=255)                                                                                             cgccgTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGGTACGCTCTGGCGCGCCGGTTTAGCGCGTGA‑GTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTGCTGCAACGTGAAGCAGAAGGTCa                                                <  2:860641/148‑1 (MQ=255)                                                                                                cgTCAGTTGGCCTAGCGGCTCCCGTGGGGCGTTTGTTACAGCAGGTACGGACTGGCGCGCCGGTTTTGCGCGTTAGGTCG‑TAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGAtt                                           <  2:385809/150‑1 (MQ=255)                                                                                                cgTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGGTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGAttt                                          >  1:299407/1‑151 (MQ=255)                                                                                                cgTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGGTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGAttt                                          >  1:1084418/1‑151 (MQ=255)                                                                                                cgTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGGTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGAttt                                          <  1:1196746/151‑1 (MQ=255)                                                                                                cgTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGGTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGAttt                                          <  1:1175628/151‑1 (MQ=255)                                                                                                   cAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGGTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGATTTTc                                        >  2:502197/1‑150 (MQ=255)                                                                                                    aGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGGTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGATTTTCAg                                      <  2:878969/151‑1 (MQ=255)                                                                                                          ccTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGGAACGCACTCGCGCGCCGGTTTAGAGCGTGA‑GTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGATTTTCAGCTGTAc                                >  2:880969/1‑151 (MQ=255)                                                                                                             aGCGGCTCCCGTGGAGCGTTTGTTACAGCAGGTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGATTTTCAGCTGTAcc                               <  2:115985/149‑1 (MQ=255)                                                                                                                 gCTCCCGTGGAGCGTTTGTTACAGCAGGTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGATTTTCAGCTGTACCCCGGCg                         >  1:1079834/1‑151 (MQ=255)                                                                                                                  cTCCCGTGGAGCGTTTGTTACAGCAGGTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGATTTTCAGCTGTACCCCGGCg                         >  2:1073122/1‑150 (MQ=255)                                                                                                                        tGGAGCGTTTGTTACAGCAGGTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGATTTTCAGCTGTACCCCGGCGAGCTggg                  >  1:1157281/1‑151 (MQ=255)                                                                                                                        tGGAGCGTTTGTTACAGCAGGTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGATTTTCAGCTGTACCCCGGCGAGCTggg                  >  1:240753/1‑151 (MQ=255)                                                                                                                         tgAGCGTTTGTTACAGCAGGTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGATTTTCAGCTGTACCCCGGCGAGCTggg                  >  2:131790/2‑150 (MQ=255)                                                                                                                            gCGTTTGTTACAGCAGGTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGATTTTCAGCTGTACCCCGGCGAGCTGGGaaaa              <  2:1142749/151‑1 (MQ=255)                                                                                                                             cGTTTGTTACAGCAGGTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGATTTTCAGCTGTACCCCGGCGAGCTGGGAAAAc             <  1:444541/151‑1 (MQ=255)                                                                                                                               tttGTTACAGCAGGTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGATTTTCAGCTGTACCCCGGCGAGCTGGGAAAACGc           >  1:715061/1‑151 (MQ=255)                                                                                                                                 tGTTACAGCAGGTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGATTTTCAGCTGTACCCCGGCGAGCTGGGAAAACGCa          <  2:280592/150‑1 (MQ=255)                                                                                                                                  gTTACAGCAGGTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGATTTTCAGCTGTACCCCGGCGAGCTGGGAAAACGCATc        >  2:655751/1‑151 (MQ=255)                                                                                                                                     aCAGCAGGTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTCGATAAAGAGGATGATTTAGGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGATTTTCAGCTGTACCCCGGCGAGCTGGGAAAACGCATCta      <  2:1018465/150‑1 (MQ=255)                                                                                                                                      cagcagGTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGATTTTCAGCTGTACCCCGGCGAGCTGGGAAAACGCATCtata    <  2:509108/151‑1 (MQ=255)                                                                                                                                       agcagGTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGATTTTCAGCTGTACCCCGGCGAGCTGGGAAAACGCATCTATaa   <  2:460445/151‑1 (MQ=255)                                                                                                                                        gcagGTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGATTTTCAGCTGTACCCCGGCGAGCTGGGAAAACGCATCTATAAc  <  2:337018/151‑1 (MQ=255)                                                                                                                                            |                                                                                                                                                    CCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAGCGCGTGA‑GTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGATTTTCAGCTGTACCCCGGCGAGCTGGGAAAACGCATCTATAAC  >  NC_000913/3103903‑3104188

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