Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	447,554	G→A	A51V (GCG→GTG)	cyoE ←	protoheme IX farnesyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	447,554	0	G	A	100.0%	268.7 / NA	82	A51V (GCG→GTG)	cyoE	protoheme IX farnesyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (36/46);  total (36/46)

TACCCAGCAACGTGGCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACGCCACAACCAGTGACA‑CCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCA  >  NC_000913/447423‑447686                                                                                                                                    |                                                                                                                                      tACCCAGCAACGTGGCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACACCACAAc                                                                                                                                <  2:1574153/139‑1 (MQ=255) tACCCAGCAACGTGGCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACACCACAAc                                                                                                                                <  2:702887/139‑1 (MQ=255)  aCCCAGCAACGTGGCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACACCACAAcc                                                                                                                               <  1:1519472/139‑1 (MQ=255)      aGCAACGTGGCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACACCACAACCAGTg                                                                                                                           <  2:463128/139‑1 (MQ=255)        cAACGAGGCGTACAGCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACACCACAACCAGTGac                                                                                                                         <  1:1375780/139‑1 (MQ=255)             tGGCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACACCACAACCAGTGACA‑cccc                                                                                                                   >  2:519751/1‑139 (MQ=255)                  tACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACACCACAACCAGTGACA‑CCCCAACCa                                                                                                              <  2:1162800/139‑1 (MQ=255)                  tACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACACCACAACCAGTGACA‑CCCCAACCa                                                                                                              >  1:1162800/1‑139 (MQ=255)                     acCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACACCACAACCAGTGACA‑CCCCAACCAGCg                                                                                                           <  2:1464390/139‑1 (MQ=255)                       cAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACACCACAACCAGTGACA‑CCCCAACCAGCgtg                                                                                                         <  1:1202289/139‑1 (MQ=255)                       cAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACACCACAACCAGTGACA‑CCCCAACCAGCgtg                                                                                                         <  2:922801/139‑1 (MQ=255)                         gCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACACCACAACCAGTGACA‑CCCCAACCAGCGTGTa                                                                                                       >  1:1456398/1‑139 (MQ=255)                          cGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACACCACAACCAGTGACA‑CCCCAACCAGCGTGTAg                                                                                                      >  1:1529009/1‑139 (MQ=255)                                 gcagGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACACCACAACCAGTGACA‑CCCCAACCAGCGTGTAGATAAACa                                                                                               >  1:1002351/1‑139 (MQ=255)                                          aTCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACACCACAACCAGTGACA‑CCCCAACCAGCGTGTAGATAAACAGGGGATAAt                                                                                      >  2:685900/1‑139 (MQ=255)                                             aGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACACCACAACCAGTGACA‑CCCCAACCAGCGTGTAGATAAACAGGGGATAATCaa                                                                                   <  1:1231690/139‑1 (MQ=255)                                               gCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTACACACACAACCCGACACCACAACCAGTg                                                                                                                           >  1:806316/1‑98 (MQ=255)                                               gCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACACCACAACCAGTg                                                                                                                           <  2:806316/98‑1 (MQ=255)                                                    tCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACACCACAACCAGTGACA‑CCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGcc                                                                            >  1:972409/1‑139 (MQ=255)                                                     cACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACACCACAACCAGTGACA‑CCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGccc                                                                           >  2:936980/1‑139 (MQ=255)                                                           cACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACACCACAACCAGTGACA‑CCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAg                                                                     <  2:927343/139‑1 (MQ=255)                                                           cACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACACCACAACCAGTGACA‑CCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAg                                                                     <  1:519751/139‑1 (MQ=255)                                                           cACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATATAGTTGTTAAACACACAACCCGACACCACAACCAGTGACA‑CCCCAACCAGCg                                                                                                           <  2:1209127/101‑1 (MQ=255)                                                           cACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATATAGTTGTTAAACACACAACCCGACACCACAACCAGTGACA‑CCCCAACCAGCg                                                                                                           >  1:1209127/1‑101 (MQ=255)                                                             cccGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACACCACAAc                                                                                                                                <  2:1601166/79‑1 (MQ=255)                                                             cccGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACACCACAAc                                                                                                                                >  1:1601166/1‑79 (MQ=255)                                                             cccGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACACCACAACCAGTGACA‑CCCCAACCAGCGTGTAGATAAAc                                                                                                >  1:27851/1‑110 (MQ=255)                                                             cccGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACACCACAACCAGTGACA‑CCCCAACCAGCGTGTAGATAAAc                                                                                                <  2:27851/110‑1 (MQ=255)                                                             cccGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACACCACAACCAGTGACA‑CCCCAACCAGCGTGTAGATAAACAgggg                                                                                           <  2:600069/115‑1 (MQ=255)                                                             cccGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACACCACAACCAGTGACA‑CCCCAACCAGCGTGTAGATAAACAgggg                                                                                           >  1:600069/1‑115 (MQ=255)                                                             cccGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACACCACAACCAGTGACA‑CCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCttt                                                                        >  1:5439/1‑134 (MQ=255)                                                             cccGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACACCACAACCAGTGACA‑CCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCttt                                                                        <  2:5439/134‑1 (MQ=255)                                                              ccGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACACCACAACCAGTGACA‑CCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGcc                                                                  >  2:1297136/1‑139 (MQ=255)                                                                      tcGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACACCACAACCAGTGACA‑CCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGaa                                                          >  2:1024085/1‑139 (MQ=255)                                                                       cGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACACCACAACCAGTGACA‑CCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAAt                                                         <  1:1096944/139‑1 (MQ=255)                                                                         tCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACACCACAACCAGTGACA‑CCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATcc                                                       <  2:408763/139‑1 (MQ=255)                                                                              tCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACACCACAACCAGTGACA‑CCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCaa                                                  >  1:286553/1‑139 (MQ=255)                                                                               ccATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACACCACAACCAGTGACA‑CCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAAt                                                 >  1:320291/1‑139 (MQ=255)                                                                                 aTCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACACCACAACCAGTGACA‑CCCCAACCAGCGTGTAGATAAAc                                                                                                <  2:371944/90‑1 (MQ=255)                                                                                 aTCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACACCACAACCAGTGACA‑CCCCAACCAGCGTGTAGATAAAc                                                                                                >  1:371944/1‑90 (MQ=255)                                                                                 aTCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACACCACAACCAGTGACA‑CCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCa                                               <  2:293366/139‑1 (MQ=255)                                                                                 aTCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACACCACAACCAGTGACA‑CCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCa                                               <  1:685900/139‑1 (MQ=255)                                                                                   cTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACACCACAACCAGTGACA‑CCCCAACCAGCGTGTAGATAAACa                                                                                               <  2:179719/89‑1 (MQ=255)                                                                                   cTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACACCACAACCAGTGACA‑CCCCAACCAGCGTGTAGATAAACa                                                                                               >  1:179719/1‑89 (MQ=255)                                                                                       cTGTCGATATCCCTGTCGATGTAGTTGTTAAACACTCAACCCGACACCACAACCAGTGACA‑CCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCgaga                                         <  2:972409/139‑1 (MQ=255)                                                                                           cGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACACCACAACCAGTGACA‑CCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAg                                     <  2:1456398/139‑1 (MQ=255)                                                                                              tatCCCTGTCGATGTAGTTGTTAAACACACAACCCGACACCACAACCAGTGACA‑CCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGtt                                  >  2:1190111/1‑139 (MQ=255)                                                                                                  ccTGTCGATGTAGTTGTTAAACACACAACCCGACACCACAACCAGTGACA‑CCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCa                              <  1:1297136/139‑1 (MQ=255)                                                                                                    tGTCGATGTAGTTGTTAAACACACAACCCGACACCACAACCAGTGACA‑CCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCcag                                                                >  1:1024129/1‑103 (MQ=255)                                                                                                    tGTCGATGTAGTTGACAAACACACAACCCGACACCACAACCAGTGACA‑CCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCcag                                                                <  2:1024129/103‑1 (MQ=255)                                                                                                      tCGATGTAGTTGTTAAACACACAACCCGACACCACAACCAGTGACA‑CCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAga                          <  1:1555165/139‑1 (MQ=255)                                                                                                        gATGTAGTTGTTAAACACACAACCCGACACCACAACCAGTGACA‑CCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAgatg                        <  1:1024085/139‑1 (MQ=255)                                                                                                          tGTAGTTGTTAAACACACAACCCGACACCACAACCAGTGACA‑CCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCaa                                                  <  1:425009/111‑1 (MQ=255)                                                                                                          tGTAGTTGTTAAACACACAACCCGACACCACAACCAGTGACA‑CCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCaa                                                  >  2:425009/1‑111 (MQ=255)                                                                                                           gTAGTTGTTAAACACACAACCCGACACCACAACCAGTGACA‑CCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGatgatg                     >  1:586384/1‑139 (MQ=255)                                                                                                             aGTTGTTAAACACACAACCCGACACCACAACCAGTGACA‑CCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGcc                   <  2:320291/139‑1 (MQ=255)                                                                                                              gttgttAAACACACAACCCGACACCACAACCAGTGACA‑CCCCAACCAGCg                                                                                                           <  1:879967/50‑1 (MQ=255)                                                                                                              gttgttAAACACACAACCCGACACCACAACCAGTGACA‑CCCCAACCAGCg                                                                                                           >  2:879967/1‑50 (MQ=255)                                                                                                                tgttAAACACACAACCCGACACCACAACCAGTGACACCCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCa                                                   >  1:1412605/1‑105 (MQ=255)                                                                                                                tgttAAACACACAACCCGACACCACAACCAGTGACACCCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCa                                                   <  2:1412605/105‑1 (MQ=255)                                                                                                                tgttAAACACACAACCCGACACCACAACCAGTGACA‑CCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCa                              >  1:667020/1‑125 (MQ=255)                                                                                                                tgttAAACACACAACCCGACACCACAACCAGTGACA‑CCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCa                              <  2:667020/125‑1 (MQ=255)                                                                                                                tgttAAACACACAACCCGACACCACAACCAGTGACA‑CCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTgg                >  1:948610/1‑139 (MQ=255)                                                                                                                tgttAAACACACAACCCGACACCACAACCAGTGACA‑CCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTgg                >  2:1410638/1‑139 (MQ=255)                                                                                                                    aaaCACACAACCCGACACCACAACCAGTGACA‑CCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAg                                                                     >  2:252380/1‑82 (MQ=255)                                                                                                                    aaaCACACAACCCGACACCACAACCAGTGACA‑CCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAg                                                                     <  1:252380/82‑1 (MQ=255)                                                                                                                      acacacAACCCGACACCACAACCAGTGACA‑CCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGt          <  2:404396/139‑1 (MQ=255)                                                                                                                      acacacAACCCGACACCACAACCAGTGACA‑CCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGt          <  2:1529009/139‑1 (MQ=255)                                                                                                                       cacacaACCCGACACCACAACCAGTGACA‑CCCCAACCAGCGTGTAGATCAACAGGGGATAATCAATGCTGCCCTTTGAGGCCATCAGGAATCCCCCAACCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGtt         <  1:41898/139‑1 (MQ=255)                                                                                                                        acacaACCCGACACCACAACCAGTGACA‑CCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTa        <  2:286553/139‑1 (MQ=255)                                                                                                                        acacaACCCGACACCACAACCAGTGACA‑CCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTa        <  2:1002351/139‑1 (MQ=255)                                                                                                                         cacaACCCGACACCACAACCAGTGACA‑CCCCAACCAGCg                                                                                                           <  1:130476/39‑1 (MQ=255)                                                                                                                         cacaACCCGACACCACAACCAGTGACA‑CCCCAACCAGCg                                                                                                           >  2:130476/1‑39 (MQ=38)                                                                                                                         cacaACCCGACACCACAACCAGTGACA‑CCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGt          >  1:866072/1‑136 (MQ=255)                                                                                                                         cacaACCCGACACCACAACCAGTGACA‑CCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGt          <  2:866072/136‑1 (MQ=255)                                                                                                                            aaCCCGACACCACAACCAGTGACA‑CCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGt          >  2:583263/1‑133 (MQ=255)                                                                                                                            aaCCCGACACCACAACCAGTGACA‑CCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGt          <  1:583263/133‑1 (MQ=255)                                                                                                                            aaCCCGACACCACAACCAGTGACA‑CCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTg    >  1:1428244/1‑139 (MQ=255)                                                                                                                            aaCCCGACACCACAACCAGTGACA‑CCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTg    <  1:1190111/139‑1 (MQ=255)                                                                                                                             aCCCGACACCACAACCAGTGACA‑CCCCAACCAGCgtg                                                                                                         <  2:483413/37‑1 (MQ=38)                                                                                                                             aCCCGACACCACAACCAGTGACA‑CCCCAACCAGCgtg                                                                                                         >  1:483413/1‑37 (MQ=255)                                                                                                                              cccGACACCACAACCAGTGACA‑CCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCa  >  2:115870/1‑139 (MQ=255)                                                                                                                                    |                                                                                                                                      TACCCAGCAACGTGGCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACGCCACAACCAGTGACA‑CCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCA  >  NC_000913/447423‑447686

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