Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,397,557	G→T	G49V (GGT→GTT)	mutL →	methyl‑directed mismatch repair protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,397,557	0	G	T	100.0%	235.4 / NA	64	G49V (GGT→GTT)	mutL	methyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (37/27);  total (37/27)

ATGCCAATTCAGGTCTTACCGCCACAACTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGC  >  NC_000913/4397412‑4397696                                                                                                                                                  |                                                                                                                                            aTGCCAATTCAGGTCTTACCGCCACAACTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTgggg                                                                                                                                        <  2:814196/151‑1 (MQ=255)      aaTTCAGGTCTTACCGCCACAACTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGaaa                                                                                                                                   >  1:1231623/1‑151 (MQ=255)       aTTCAGGTCTTACCGCCACAACTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGaga                                                                                                                                   >  1:974606/1‑148 (MQ=255)       aTTCAGGTCTTACCGCCACAACTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGaaa                                                                                                                                   >  2:28914/1‑150 (MQ=255)       aTTCAGGTCTTACCGCCACAACTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGaaa                                                                                                                                   >  2:764557/1‑150 (MQ=255)            ggTCTTACCGCCACAACTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTAt                                                                                                                              >  2:909880/1‑150 (MQ=255)             gttttcccGCCACAACTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATc                                                                                                                             <  2:173343/144‑1 (MQ=255)                      ccACAACTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTg                                                                                                                   >  1:493322/1‑151 (MQ=255)                           aCTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAAc                                                                                                              <  2:86032/151‑1 (MQ=255)                           aCTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAAc                                                                                                              <  2:757172/151‑1 (MQ=255)                                      aGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTAt                                                                                                   >  1:580518/1‑151 (MQ=255)                                      aGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTAt                                                                                                   >  2:792549/1‑151 (MQ=255)                                           gCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCaaaa                                                                                              >  1:928468/1‑151 (MQ=255)                                            ccGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCaaaaa                                                                                             <  1:1142812/151‑1 (MQ=255)                                            ccGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCaaaaa                                                                                             >  1:22086/1‑151 (MQ=255)                                               cAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGa                                                                                          >  1:868710/1‑151 (MQ=255)                                                  gtgaggtgGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGa                                                                                          <  1:1205620/148‑1 (MQ=255)                                                     aggtgGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGCCGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAg                                                                                      >  2:730609/1‑149 (MQ=255)                                                       gtggtCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctgg                                                                                  >  1:1036497/1‑151 (MQ=255)                                                         ggtCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcg                                                                                >  1:1029413/1‑151 (MQ=255)                                                          gtCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcg                                                                                <  2:812660/150‑1 (MQ=255)                                                          gtCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcg                                                                                <  1:1080946/150‑1 (MQ=255)                                                          gtCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGGGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGAATCAAAAAAGATGAgctggcg                                                                                >  1:960837/1‑150 (MQ=255)                                                           tCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgct                                                                              <  1:415263/151‑1 (MQ=255)                                                                  aCCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcgct                                                                        <  2:132754/150‑1 (MQ=255)                                                                  aCCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcg                                                                          >  1:950761/1‑148 (MQ=255)                                                                          cGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATTGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCa                                                               <  2:1167478/151‑1 (MQ=255)                                                                          cGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCa                                                               <  1:104471/151‑1 (MQ=255)                                                                                tCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGcc                                                           <  1:722031/149‑1 (MQ=255)                                                                                tCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGcc                                                           <  2:816254/149‑1 (MQ=255)                                                                                tCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGAGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGcca                                                          >  2:892117/1‑150 (MQ=255)                                                                                 cAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGccacc                                                        <  1:415403/151‑1 (MQ=255)                                                                                 cAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGccac                                                         >  2:1225568/1‑150 (MQ=255)                                                                                   aaGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCCCCAg                                                      >  2:1244911/1‑151 (MQ=255)                                                                                     gAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTa                                                    >  1:930249/1‑151 (MQ=255)                                                                                      aaCTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGt                                                     <  2:1264797/149‑1 (MQ=255)                                                                                        cTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTaaaa                                                 >  2:389396/1‑151 (MQ=255)                                                                                         tAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAAt                                                <  1:1019959/151‑1 (MQ=255)                                                                                           gTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCCCGTCATGCCACCCGTAAAATc                                               >  2:566243/1‑150 (MQ=255)                                                                                           gTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGCTATCAAAAAAGATGACCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCg                                              >  2:1241785/1‑151 (MQ=255)                                                                                             ggAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGc                                             >  1:1120534/1‑150 (MQ=255)                                                                                             ggAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGc                                             >  1:381150/1‑150 (MQ=255)                                                                                               aaagCAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTGTCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGGGCTGGCGCTGGCGCTGGCTCGTCATGCCGCCAGTAAAATCGCctc                                          <  1:222183/151‑1 (MQ=255)                                                                                                         cGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGAt                                <  1:273922/151‑1 (MQ=255)                                                                                                           aTGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGAtct                              <  2:1222060/151‑1 (MQ=255)                                                                                                           aTGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGAt                                >  1:1181986/1‑149 (MQ=255)                                                                                                                ggTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTTTCTGGACGATCTCGAAg                         >  2:886394/1‑151 (MQ=255)                                                                                                                 gTGCGACGCGTATCGATATTGATATCGAACGGGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAg                         <  1:1211730/150‑1 (MQ=255)                                                                                                                  tGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGc                        <  2:1231103/150‑1 (MQ=255)                                                                                                                         cgTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATc                 <  1:28919/150‑1 (MQ=255)                                                                                                                         cgTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCa                >  1:960004/1‑151 (MQ=255)                                                                                                                           tATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAg               >  1:1173523/1‑150 (MQ=255)                                                                                                                           tATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAg               >  1:937299/1‑150 (MQ=255)                                                                                                                           tATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAg               >  1:780232/1‑150 (MQ=255)                                                                                                                           tATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGc              <  1:428627/151‑1 (MQ=255)                                                                                                                            aTCGATATTGATATCGACCGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGcc             <  2:867577/151‑1 (MQ=255)                                                                                                                                atatTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTggg         >  1:259952/1‑151 (MQ=255)                                                                                                                                atatTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTggg         >  2:452688/1‑151 (MQ=255)                                                                                                                                atatTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGCCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGACTggg         >  2:1075803/1‑151 (MQ=255)                                                                                                                                 tatTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGc        <  1:228693/151‑1 (MQ=255)                                                                                                                                   ttGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCtt      >  2:1230792/1‑151 (MQ=255)                                                                                                                                     gATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCttt     >  2:841677/1‑150 (MQ=255)                                                                                                                                      atatCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTcg   >  1:268442/1‑151 (MQ=255)                                                                                                                                       tatCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATACCGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCATTcgc  <  2:650046/151‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                            ATGCCAATTCAGGTCTTACCGCCACAACTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGC  >  NC_000913/4397412‑4397696

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