Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,077,754	G→A	A62A (GCG→GCA)	cbtA →	CP4‑44 prophage; toxin of the YeeV‑YeeU toxin‑antitoxin system

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,077,754	0	G	A	100.0%	235.1 / NA	62	A62A (GCG→GCA)	cbtA	CP4‑44 prophage; toxin of the YeeV‑YeeU toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (31/31);  total (31/31)

CGCCCGTCTCCTGTTGAAATCTGGCAGATACTGCTGTCCCGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGA  >  NC_000913/2077611‑2077898                                                                                                                                                |                                                                                                                                                 cGCCCGTCTCCTGTTGAAATCTGGCAGATACTGCTGTCCCGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGa                                                                                                                                              <  2:536449/148‑1 (MQ=255)   cccGTCTCCTGTTGAAATCTGGCAGATACTGCTGTCCCGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACttt                                                                                                                                         <  2:497270/151‑1 (MQ=255)       tctcCTGTTGAAATCTGGCAGATACTGCTGTCCCGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCg                                                                                                                                     >  1:564729/1‑151 (MQ=255)          ccTGTTGAAATCTGGCAGATACTGCTGTCCCGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTgg                                                                                                                                  <  2:229798/151‑1 (MQ=255)              ttGAAATCTGGCAGATACTGCTGTCCCGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGaaa                                                                                                                               >  2:313698/1‑150 (MQ=255)               tGAAATCTGGCAGATACTGCTGTCCCGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGaaaa                                                                                                                              <  2:86030/150‑1 (MQ=255)                gAAATCTGGCAGATACTGCTGTCCCGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGaaaaa                                                                                                                             <  2:669216/150‑1 (MQ=255)                  aaTCTGGCAGATACTGCTGTCCCGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATAc                                                                                                                          >  1:317774/1‑151 (MQ=255)                      tGGCAGATACTGCTGTCCCGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACgcgc                                                                                                                      >  1:627064/1‑151 (MQ=255)                       ggCAGATACTGCTGTCCCGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATAcg                                                                                                                         <  2:296613/147‑1 (MQ=255)                          aGATACTGCTGTCCCGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGt                                                                                                                  <  1:653160/151‑1 (MQ=255)                          aGATACTGCTGTCCCGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGt                                                                                                                  >  2:306064/1‑151 (MQ=255)                              aCTGCTGTCCCGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCg                                                                                                               <  2:441547/150‑1 (MQ=255)                                 gctgTCCCGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTa                                                                                                             <  2:702733/149‑1 (MQ=255)                                 gctgTCCCGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTAcc                                                                                                           >  2:463391/1‑151 (MQ=255)                                          aCTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCg                                                                                                  >  2:357856/1‑151 (MQ=255)                                             gctgGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCgg                                                                                                 <  2:50933/149‑1 (MQ=255)                                                ggACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGAt                                                                                              <  1:306039/149‑1 (MQ=255)                                                ggACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTc                                                                                            >  1:821428/1‑151 (MQ=255)                                                 gACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCa                                                                                           <  2:566592/151‑1 (MQ=255)                                                 gACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCa                                                                                           <  1:236208/151‑1 (MQ=255)                                                        aCTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTg                                                                                    >  1:255525/1‑151 (MQ=255)                                                         cTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGt                                                                                   >  2:403900/1‑151 (MQ=255)                                                         cTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGt                                                                                   >  1:555033/1‑151 (MQ=255)                                                          tATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGt                                                                                   >  2:675265/1‑150 (MQ=255)                                                                cTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCa                                                                                           <  1:806070/136‑1 (MQ=255)                                                                cTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGAGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTc                                                                                            >  2:806126/1‑135 (MQ=255)                                                                      cTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGt                                                                      >  2:108753/1‑151 (MQ=255)                                                                                cacCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACa                                                             <  2:670786/150‑1 (MQ=255)                                                                                  ccTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGc                                                           >  2:227435/1‑150 (MQ=255)                                                                                     tttGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCAt                                                         <  1:694990/149‑1 (MQ=255)                                                                                         ccGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGatat                                                   >  2:123770/1‑151 (MQ=255)                                                                                         ccGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGatat                                                   <  2:227657/151‑1 (MQ=255)                                                                                             tGAACGTGTGATTGAGCAGCATATTGGGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATAtcctc                                               >  1:623389/1‑151 (MQ=255)                                                                                              gAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATAtcctcc                                              >  2:617775/1‑151 (MQ=255)                                                                                                   tgtgATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCt                                         >  1:17016/1‑151 (MQ=255)                                                                                                    gtgATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTGCCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTc                                        <  1:815498/151‑1 (MQ=255)                                                                                                    gtgATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCt                                         <  1:357833/150‑1 (MQ=255)                                                                                                    gtgATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTc                                        <  1:636955/151‑1 (MQ=255)                                                                                                    gtgATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTc                                        <  1:721585/151‑1 (MQ=255)                                                                                                       aTTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGc                                      >  1:688689/1‑150 (MQ=255)                                                                                                        ttGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGc                                      >  2:458639/1‑149 (MQ=255)                                                                                                         tGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAgg                                   >  2:398770/1‑151 (MQ=255)                                                                                                               gcaTATTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGAcc                             >  1:187919/1‑151 (MQ=255)                                                                                                                 atatTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCg                            >  1:421635/1‑150 (MQ=255)                                                                                                                   atTGAGGCAGGCATTTCTCTGTGTGATGCAGTGACCTTTCTCGTGGAAAAATACGCGCTGGTGCGTAGCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCATTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGc                          <  1:769504/150‑1 (MQ=255)                                                                                                                   atTGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGcc                         <  1:460217/151‑1 (MQ=255)                                                                                                                     tGAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCtg                       >  2:4848/1‑151 (MQ=255)                                                                                                                      gAGGCAGGCATTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCtg                       <  2:7167/150‑1 (MQ=255)                                                                                                                           aggcaTTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGAcc                 >  2:112384/1‑151 (MQ=255)                                                                                                                            ggcaTTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGAccc                <  2:679540/151‑1 (MQ=255)                                                                                                                            ggcaTTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGAccc                >  2:147301/1‑151 (MQ=255)                                                                                                                               aTTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCcgcg             >  2:471987/1‑151 (MQ=255)                                                                                                                               aTTTCACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCcgcg             <  1:523102/151‑1 (MQ=255)                                                                                                                                   cACTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACa          >  1:704002/1‑150 (MQ=255)                                                                                                                                     cTGTGTGATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAAt        <  1:182507/150‑1 (MQ=255)                                                                                                                                       gtgtgATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAAtt       <  1:670325/149‑1 (MQ=255)                                                                                                                                       gtgtgATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTAc     <  1:333651/151‑1 (MQ=255)                                                                                                                                        tgtgATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTAGCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGACTGATGACCTGCGACAATTAc     >  2:527132/1‑150 (MQ=255)                                                                                                                                         gtgATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAg   <  1:403705/151‑1 (MQ=255)                                                                                                                                          tgAGGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAg   <  2:619070/150‑1 (MQ=255)                                                                                                                                           gATGCAGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGa  >  2:20896/1‑150 (MQ=255)                                                                                                                                                |                                                                                                                                                 CGCCCGTCTCCTGTTGAAATCTGGCAGATACTGCTGTCCCGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGA  >  NC_000913/2077611‑2077898

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