Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,040,433	C→A	A319D (GCC→GAC)	msrP →	protein‑L‑methionine sulfoxide reductase catalytic subunit MsrP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,040,433	0	C	A	100.0%	159.8 / NA	45	A319D (GCC→GAC)	msrP	protein‑L‑methionine sulfoxide reductase catalytic subunit MsrP
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (24/21);  total (24/21)

GTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGCCCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCGACCTTGTTAATCACCCCTCTGG  >  NC_000913/2040186‑2040675                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   gTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCAGCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCa                                                                                                                                                                                                                                                 <  2:200742/251‑1 (MQ=255)     gACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTg                                                                                                                                                                                                                                             <  1:331048/251‑1 (MQ=255)      aCTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTgg                                                                                                                                                                                                                                            >  1:462065/1‑251 (MQ=255)            tgtgCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCATCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGc                                                                                                                                                                                                                                      >  1:185631/1‑251 (MQ=255)              ggCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACAACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTg                                                                                                                                                                                                                                    <  1:164764/250‑1 (MQ=255)                        atatGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGcc                                                                                                                                                                                                                          >  2:270674/1‑251 (MQ=255)                                   aaGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCgg                                                                                                                                                                                                               >  1:33292/1‑251 (MQ=255)                                   aaGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCgg                                                                                                                                                                                                               >  2:402763/1‑251 (MQ=255)                                      ggATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCggg                                                                                                                                                                                                              >  2:450799/1‑249 (MQ=255)                                            aaTCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTc                                                                                                                                                                                                      >  1:21229/1‑251 (MQ=255)                                                     tccgatctAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATTCGt                                                                                                                                                                                             >  1:186748/8‑251 (MQ=255)                                                         tATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTg                                                                                                                                                                                          >  1:104682/1‑250 (MQ=255)                                                             aaGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAg                                                                                                                                                                                      <  1:270674/250‑1 (MQ=255)                                                               gCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCATCGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCaa                                                                                                                                                                                   <  2:389189/251‑1 (MQ=255)                                                                  gACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGACTACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGCATTTCTAAATGCGTCTGACAGCAAAAc                                                                                                                                                                                >  2:362357/1‑251 (MQ=255)                                                                    cccGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAg                                                                                                                                                                              <  2:280211/251‑1 (MQ=255)                                                                       gcgAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTg                                                                                                                                                                           <  1:181316/251‑1 (MQ=255)                                                                                                   gnCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCt                                                                                                                                                 <  2:468083/249‑1 (MQ=255)                                                                                                   ggCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCtt                                                                                                                                                <  1:402763/250‑1 (MQ=255)                                                                                                     cAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGc                                                                                                                                                                                                                 >  1:284394/1‑183 (MQ=255)                                                                                                     cAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGc                                                                                                                                                                                                                 <  2:284394/183‑1 (MQ=255)                                                                                                     cAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGcc                                                                                                                                             <  2:257398/251‑1 (MQ=255)                                                                                                          ccTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATc                                                                                                                                                                                                                                <  2:98001/163‑1 (MQ=255)                                                                                                          ccTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATc                                                                                                                                                                                                                                >  1:98001/1‑163 (MQ=255)                                                                                                               cGAATACGGTTTGTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGc                                                                                                                                   <  2:512345/251‑1 (MQ=255)                                                                                                                    aCGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGttt                                                                                                                              >  1:431040/1‑251 (MQ=255)                                                                                                                      ggtgtttACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGTCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTtct                                                                                                                            <  2:423827/247‑1 (MQ=255)                                                                                                                           ttACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTctgg                                                                                                                        <  1:355732/250‑1 (MQ=255)                                                                                                                                           ccGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGatca                                                                                                        >  1:434756/1‑250 (MQ=255)                                                                                                                                                    gATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTgg                                                                                              >  2:441011/1‑251 (MQ=255)                                                                                                                                                             cgcTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTCTCGTGGCCTGGATTCGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCACTCACGGTGGACTgggcgc                                                                                     >  1:353061/1‑248 (MQ=255)                                                                                                                                                              gcTGGTCACAGGCTACAGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGACTTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGcc                                                                                    >  2:394488/1‑251 (MQ=255)                                                                                                                                                               cTGGGCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGAGCCg                                                                                   >  2:488006/1‑251 (MQ=255)                                                                                                                                                                    cacaGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGAACAATCACGGTGGACTGGGTGCCGATCCg                                                                              >  2:195580/1‑251 (MQ=255)                                                                                                                                                                             ccGAACGATGTATTTGTGCAGGCGGCATCCTCGATGTACAGCGCCAGCCAATGTTACTGTTGAATGGTTACGCCGACCAGGTGGCGTCGCTGTAGCGTGTCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGTCAGCGAAACAGGTGCCATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCGGGTCTGGGCGATCAATCACGGTGGACTCGGTGCCGATCCGGTGAAAGat                                                                     <  1:103343/251‑1 (MQ=255)                                                                                                                                                                                 aCGATTTATTGGTTCAGGCGGCATCCTCGCTGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGGGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTc                                                                 >  1:355604/1‑251 (MQ=255)                                                                                                                                                                                    aTTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGc                                                              <  2:238464/251‑1 (MQ=255)                                                                                                                                                                                     tttATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGc                                                              >  1:378758/1‑250 (MQ=255)                                                                                                                                                                                        aTTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCAttt                                                          <  1:362357/251‑1 (MQ=255)                                                                                                                                                                                                 ggcggcATCATCGATGTACAGCGCCAGCCAACGCTACTGTTTACGGGTTACGCCGAACAGACGGCATCGCTGTATCGTGGTCTGGATTTGCGGGAGAATTTCTCAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAGCGTTTGCCTGCATCTTGCCGGGTTGATGCAGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGGTGTTCAGCATTTTACTGGTCg                                                 >  2:256010/1‑251 (MQ=255)                                                                                                                                                                                                      cATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCgcactg                                            <  2:449364/251‑1 (MQ=255)                                                                                                                                                                                                      cATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCgcactg                                            <  1:227365/251‑1 (MQ=255)                                                                                                                                                                                                                 tACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAAttt                                 >  2:267038/1‑251 (MQ=255)                                                                                                                                                                                                                                              ggTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCATTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCGACCTTGTTAATCACCCctct    <  2:260759/251‑1 (MQ=255)                                                                                                                                                                                                                                                ttACGCCGACCAGGTGGCATCGCTGTATCGTGGACTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCGACCTTGTTAATCACCCCTCTgg  >  1:546/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   GTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGCCCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCGACCTTGTTAATCACCCCTCTGG  >  NC_000913/2040186‑2040675

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