Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,165,497	G→A	V24M (GTG→ATG)	murI →	glutamate racemase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,165,497	0	G	A	100.0%	227.1 / NA	66	V24M (GTG→ATG)	murI	glutamate racemase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (33/33);  total (33/33)

CACCTCTCACCTGACAGTTCGTGGTAAAATAGCCAACCTGTTCGACAAAGATTATGAGACAGTCTATGGCTACCAAACTGCAGGACGGGAATACACCTTGTCTGGCAGCTACACCTTCTGAACCACGTCCCACCGTGCTGGTGTTTGACTCCGGCGTCGGTGGGTTGTCGGTCTATGACGAGATCCGGCATCTCTTACCGGATCTCCATTACATTTATGCTTTCGATAACGTCGCTTTCCCGTATGGCGAAAAAAGCGAAGCGTTTATTGT  >  NC_000913/4165363‑4165633                                                                                                                                       |                                                                                                                                         cACCTCTCACCTGACAGTTCGTGGTAAAATAGCCAACCTGTTCGACAAAGATTATGAGACAGTCTATGGCTACCAAACTGCAGGACGGGAATACACCTTGTCTGGCAGCTACACCTTCTGAACCACGTCCCACCATGCTg                                                                                                                                     >  1:21157/1‑140 (MQ=255)            tGACAGTTCCTGGTAAAATAGCCAACCTGTTCGACAAAGATTATGAGACAGTCTATGGCTACCAAACTGCAGGACGGGAATACACCTTGTCTGGCAGCTACACCTTCTGAACCACGTCCCACCATGCTGGTGTTTGACTc                                                                                                                          <  2:455243/140‑1 (MQ=255)              aCAGTTCGTGGTAAAATAGCCAACCTGTTCGACAAAGATTATGAGACAGTCTATGGCTACCAAACTGCAGGACGGGAATACGCCTTGTCTGGCAGCTACACCTTCTGAACCACGTCCCACCATGCTGGTGTTTGACTCCg                                                                                                                        >  2:698031/1‑140 (MQ=255)              aCAGTTCGTGGTAAAATAGCCAACCTGTTCGACAAAGATTATGAGACAGTCTATGGCTACCAAACTGCAGGACGGGAATACACCTTGTCTGGCAGCTACACCTTCTGAACCACGTCCCACCATGCTGGTGTTTGACTCCg                                                                                                                        >  1:464834/1‑140 (MQ=255)                  ttCGTGGTAAAATAGCCAACCTGTTCGACAAAGATTATGAGACAGTCTATGGCTACCAAACTGCAGGACGGGAATACACCTTGTCTGGCAGCTACACCTTCTGAACCACGTCCCACCATGCTGGTGTTTGACTCCGGCGt                                                                                                                    >  1:554797/1‑140 (MQ=255)                  ttCGTGGTAAAATAGCCAACCTGTTCGACAAAGATTATGAGACAGTCTATGGCTACCAAACTGCAGGACGGGAA‑ACACCTTGTCTGGCAGCTACACCTTCTGAACCACGTCCCACCATGCTGGTGTTTGACTCCGGCGTc                                                                                                                   >  2:183435/1‑140 (MQ=255)                      tGGTAAAATAGCCAACTTGTTCGACAAAGATTATGGGACCGTCTATGGTTACCAAACTGCAGGACGGGAATACACCTTGTCTGGCAGCTACCCCTTCTGAACCACGTCCCACCATGCTGGTGTTTGACTCCGGCGTCGGt                                                                                                                <  2:357716/140‑1 (MQ=255)                           aaaTAGCCAACCTGTTCGACAAAGATTATGAGACAGTCTATGGCTACCAAACTGCAGGACGGGAATACACCTTGTCTGGCAGCTACACCTTCTGAACCACGTCCCACCATGCTGGTGTTTGACTCCGGCGTCGGTGGGtt                                                                                                           >  1:434934/1‑140 (MQ=255)                             aTAGCCAACCTGTTCGACAAAGATTATGAGACAGTCTATGGCTACCAAACTGCAGGACGGGAATACACCTTGTCTGGCAGCTACACCTTCTGAACCACGTCCCACCATGCTGGTGTTTGACTCCGGCGTCGGTGGGTTGt                                                                                                         <  2:554797/140‑1 (MQ=255)                                 ccAACCTGTTCGACAAAGATTATGAGACAGTCTATGGCTACCAAACTGCAGGACGGGAATACACCTTGTCTGGCAGCTACACCTTCTGAACCACGTCCCACCATGCTGGTGTTTGACTCCGGCGTCGGTGGGTTGTCGGt                                                                                                     <  1:124040/140‑1 (MQ=255)                                 ccAACCTGTTCGACAAAGATTATGAGACAGTCTATGGCTACCAAACTGCAGGACGGGAATACACCTTGTCTGGCAGCTACACCTTCTGAACCACGTCCCACCATGCTGGTGTTTGACTCCGGCGTCGGTGGGTTGTCGGt                                                                                                     <  2:672637/140‑1 (MQ=255)                                      cTGTTCGACAAAGATTATGAGACAGTCTATGGCTACCAAACTGCAGGACGGGAATACACCTTGTCTGGCAGCTACACCTTCTGAACCACGTCCCACCATGCTGGTGTTTGACTCAGGCGTCGGTGGGTTGTCGGTCTATg                                                                                                <  2:496056/140‑1 (MQ=255)                                         ttCGACAAAGATTATGAGACAGTCTATGGCTACCAAACTGCAGGACGGGAATACACCTTGTCTGGCAGCTACACCTTCTGAACCACGTCCCACCATGCTGGTGTTTGACTCCGGCGTCGGTGGGTTGTCTGTCTATGACg                                                                                             >  1:64571/1‑140 (MQ=255)                                                 aGATTATGAGACAGTCTATGGCTACCAAACTGCAGGACGGGAATACACCTTGTCTGGCAGCTACACCTTCTGAACCACGTCCCACCATGCTGGTGTTTGAATCCGGCGTCGGTGGGTTGTCGCTCTATGACGAGATCCgg                                                                                     >  2:281807/1‑140 (MQ=255)                                                  gATTATGAGACAGTCTATGGCTACCAAACTGCAGGACGGGAATACACCTTGTCTGGCAGCTACACCTTCTGAACCACGTCCCACCATGCTGGTGTTTGACTCCGGCGTCGGTGGGTTGTCGGTCTATGACGAGATCCGGc                                                                                    <  1:131406/140‑1 (MQ=255)                                                    ttATGAGACAGTCTATGGCTACCAAACTGTAGGACGGGAATACACCTTGTCTGGCAGCTACACCTTCTGAACCACGTCCCACCATGCTGGTGTTTGACTCCGGCGTCGGTGGGTTGTCGGTCTATGACGAGATCCGGCAt                                                                                  <  1:323390/140‑1 (MQ=255)                                                    ttATGAGACAGTCTATGGCTACCAAACTGCAGGACGGGAATACACCTTGTCTGGCAGCTACACCTTCTGAACCACGTCCCACCATGCTGGTGTTTGACTCCGGCGTCGGTGGGTTGTCGGTCTATGACGAGATCCGGCAt                                                                                  <  1:671354/140‑1 (MQ=255)                                                     tATGAGACAGTCTATGGCTACCAAACTGCAGGACGGGAATACACCTTGTCTGGCAGCTACACCTTCTGAACCACGTCCCACCATGCTGGTGTTTGACTCCGGCGTCGGTGGGTTGTCGGTCTATGACGAGATCCGGCAtc                                                                                 <  2:21157/140‑1 (MQ=255)                                                           aCAGTCTATGGCTACCAAACTGCAGGACGGGAATACACCTTGTCTGGCAGCTACACCTTCTGAACCACGTCCCACCATGCTGGTGTTTGACTCCGGCGTCGGTGGGTTGTCGGTCTATGACGAGATCCGGCATCTCTTAc                                                                           <  1:686941/140‑1 (MQ=255)                                                             aGTCTATGGCTACCAAACTGCAGGACGGGAATACACCTTGTCTGGCAGCTACACCTTCTGAACCACGTCCCACCATGCTGGTGTTTGACTCCGGCGTCGGTGGGTTGTCGGTCTATGACGAGATCCGGCATCTCTTACCg                                                                         >  2:149523/1‑140 (MQ=255)                                                              gTCTATGGCTACCAAACTGCAGGACGGGAATACACCTTGTCTGGCAGCTACACCTTCTGAACCACGTCCCACCATGCTGGTGTTTGACTCCGGCGTCGGTGGGTTGTCGGTCTATGACGAGATCCGGCATCTCTTACCgg                                                                        <  2:496167/140‑1 (MQ=255)                                                               tCTATGGCTACCAAACTGCAGGACGGGAATACACCTTGTCTGGCAGCTACACCTTCTGAACCACGTCCCACCATGCTGGTGTTTGACTCCGGCGTCGGTGGGTTGTCGGTCTATGACGAGATCCGGCAtctc                                                                               <  1:179182/132‑1 (MQ=255)                                                               tCTATGGCTACCAAACTGCAGGACGGGAATACACCTTGTCTGGCAGCTACACCTTCTGAACCACGTCCCACCATGCTGGTGTTTGACTCCGGCGTCGGTGGGTTGTCGGTCTATGACGAGATCCGGCAtctc                                                                               >  2:179182/1‑132 (MQ=255)                                                               tCTATGGCTACCAAACTGCAGGACGGGAATACACCTTGTCTGGCAGCTACACCTTCTGAACCACGTCCCACCATGCTGGTGTTTGACTCCGGCGTCGGTGGGTTGTCGGTCTATGACGAGATCCGGCATCTCTTACCGGa                                                                       >  1:171180/1‑140 (MQ=255)                                                               tCTATGGCTACCAAACTGCAGGACGGGAATACACCTTGTCTGGCAGCTACACCTTCTGAACCACGTCCCACCATGCTGGTGTTTGACTCCGGCGTCGGTGGGCTGTCGGTGTATGACGAGATCCCGCATCTCTTACCGGa                                                                       >  1:210572/1‑140 (MQ=255)                                                               tCTATGGCTACCAAACTGCAGGACGGGAATACACCTTGTCTGGCAGCTACACCTTCTGAACCACGTCCCACCATGCTGGTGTTTGACTCCGACGTCGGTGCGTTATCGGTCTATGACGAGATGCGGCATCTCTTACCGGa                                                                       >  2:222847/1‑140 (MQ=255)                                                                cTATGGCTACCAAACTGCAGGACGGGAATACACCTTGTCTGGCAGCTACACCTTCTGAACCACGTCCCCCCATGCTGGTGTTTGACTCCGGCGTCGGTGGGTTGTCGGTCTATGACGAGATCCGGCATCTCTTACCGGAt                                                                      <  1:720828/140‑1 (MQ=255)                                                                cTATGGCTACCAAACTGCAGGACGGGAATACACCTTGTCTGGCAGCTACACCTTCTGAACCACGTCCCACCATGCTGGTGTTTGACTCCGGCGTCGGTGGGTTGTCGGTCTATGACGAGATCCGGCATCTCTTACCGGAt                                                                      <  2:561687/140‑1 (MQ=255)                                                                cTATGGCTACCAAACTGAAGGACGGGAATACACATAGTCTGGCAGCTAGTCCTTCTGAACCACGTCCCACCATGCTGGTGTTTGACTCCGGCGTCGGTGGGTTGTCGGTCTATAACGAGACCCGGCATCTCTTACCGGAt                                                                      <  2:353599/140‑1 (MQ=255)                                                                   tGGCTACCAAACTGCAGGACGGGAATACACCTTGTCTGGCAGCTACACCTTCTGAACCACGTCCCACCATGCTGGTGTTTGACTCCGGCGTCGGTGGGTTGTAGGTCTATTACGAGATCCGGCAACTCTTACCGGAtctc                                                                   >  1:324938/1‑140 (MQ=255)                                                                     gCTACCAAACTGCAGGACGGGAATACACCTTGTCTGGCAGCTACACCTTCTGAACCACGTCCCACCATGCTGGTGTTTGACTCCGGCGTCGGTGGGTTGTCGGTCTATGACGAGATCCGGCATCTCTTACCGGATCTCCa                                                                 <  1:381044/140‑1 (MQ=255)                                                                     gCTACCAAACTGCAGGACGGGAATACACCTTGTCTGGCAGCTACACCTTCTGAACCACGTCCCACCATGCTGGTGTTTGACTCCGGCGTCGGTGGGTTGTCGGTCTATGACGAGATCCGGCATCTCTTACCGGATCTCCa                                                                 >  1:11388/1‑140 (MQ=255)                                                                      cTACCAAACTGCAGGACGGGAATACACCTTGTCTGGCAGCTACACCTTCTGAACCACGTCCCACCATGCTGGTGTTTGACTCCGGCGTCGGTGGGTTGTCGGTCTATGACGAGATCCGGCATCTCTTACCGGATCTCCAt                                                                <  2:171180/140‑1 (MQ=255)                                                                      cTACCAAACTGCAGGACGGGAAGACCCCTTGTCTGGCAGCTACACCTTCTGAACCACGTCCCACCATGCTGGTGTTTGACTCCGGCGTCGGTGGGTTGTCGGTCTATGACGAGATCCGGCATCTCTTACCGGATCTCCAt                                                                <  1:323166/140‑1 (MQ=255)                                                                      cTACCAAACTGCAGGACGGGAAAACACCTTGCCTGCCAGATACACCTTCTGAACCACGTCCCACCATGCTGGTGTTTGACTCCGGCGTCGGTGGGTTGTCGGTCTATGACGAGATCCGGCATCTCTTACCGGATCTCCAt                                                                <  2:210572/140‑1 (MQ=255)                                                                         ccAAACTGCAGGACGGGAATACACCTTGTCTGGCAGCTACACCTTCTGAACCACGTCCCACCATGCTGGTGTTTGACTCCGGCGTCGGTGGGTTGTCGGTCTTTGACGAGATCCGGCGTCTCTTGCCGGATCTCCATTAc                                                             >  1:291034/1‑140 (MQ=255)                                                                         ccAAACTGCAGGACGGGAATACACCTTGTCTGGCAGCTACACCTTCTGAACCACGTCCCACCATGCTGGTGTTTGACTCCGGCGTCGGTGGGTTGTCGGTCTATGACGAGATCCGGCATCTCTTACCGGAt                                                                      >  1:610003/1‑131 (MQ=255)                                                                         ccAAACTGCAGGACGGGAATACACCTTGTCTGGCAGCTACACCTTCTGAACCACGTCCCACCATGCTGGTGTTTGACTCCGGCGTCGGTGGGTTGTCGGTCTATGACGAGATCCGGCATCTCTTACCGGAt                                                                      <  2:610003/131‑1 (MQ=255)                                                                         ccAAACTGCAGGACGGGAATACACCTTGTCTGGCAGCTACACCTTCTGAACCACGTCCCACCATGCTGGTGTTTGACTCCGGCGTCGGTGGGTTGTCGGTCTATGACGAGATCCGGCATCTCTTACCGGATCTCCAt                                                                >  2:19236/1‑137 (MQ=255)                                                                         ccAAACTGCAGGACGGGAATACACCTTGTCTGGCAGCTACACCTTCTGAACCACGTCCCACCATGCTGGTGTTTGACTCCGGCGTCGGTGGGTTGTCGGTCTATGACGAGATCCGGCATCTCTTACCGGATCTCCAt                                                                >  2:603009/1‑137 (MQ=255)                                                                         ccAAACTGCAGGACGGGAATACACCTTGTCTGGCAGCTACACCTTCTGAACCACGTCCCACCATGCTGGTGTTTGACTCCGGCGTCGGTGGGTTGTCGGTCTATGACGAGATCCGGCATCTCTTACCGGATCTCCAt                                                                <  1:19236/137‑1 (MQ=255)                                                                         ccAAACTGCAGGACGGGAATACACCTTGTCTGGCAGCTACACCTTCTGAACCACGTCCCACCATGCTGGTGTTTGACTCCGGCGTCGGTGGGTTGTCGGTCTATGACGAGATCCGGCATCTCTTACCGGATCTCCAt                                                                <  1:603009/137‑1 (MQ=255)                                                                         ccAAACTGCAGGACGGGAATACACCTTGTCTGGCAGCTACACCTTCTGAACCACGTCCCACCATGCTGGTGTTTGACTCCGGCGTCGGTGGGTTGTCGGTCTATGACGAGATCCGGCATCTCTTACCGGATCTCCATTAc                                                             >  2:168895/1‑140 (MQ=255)                                                                          cAAACTGCAGGACGGGAATACACCTTGTCTGGCAGCTACACCTTCTGAACCACGTCCCACCATGCTGGTGTTTGACTCCGGCGTCGGTGGGTTGTCGGTCTATGACGAGATCCGGCATCTCTTACCGGATCTCCATTACa                                                            <  1:655518/140‑1 (MQ=255)                                                                                  aGGACGGGAATACACCTTGTCTGGCAGCTACACCTTCTGAACCACGTCCCACCATGCTGGTGTTTGACTCCGGCGTCGGTGGGTTGTCGGTCTATGACGAGATCCGGCATCTCTTACCGGATCTCCATTACATTTATGCt                                                    <  1:5309/140‑1 (MQ=255)                                                                                    gACGGGAATACACCTTGTCTGGCAGCTACACCTTCTGAACCACGTCCCACCATGCTGGTGTTTGACTCCGGCGTCGGTGGGTTGTCGGTCTATGACGAGATCCGGCATCTCTTACCGGATCTCCATTACATTTATGCttt                                                  >  2:391648/1‑140 (MQ=255)                                                                                                 cTTGTCTGGCAGCTACACCTTCTGAACCACGTCCCACCATGCTGGTGTTTGACTCCGGCGTCGGTGGGTTGTCGGTCTATGACGAGATCCGGCATCTCTTACCGGATCTCCATTACATTTATGCTTTCGATAACGTCGCt                                     >  1:673851/1‑140 (MQ=255)                                                                                                   tGTCTGGCAGCTACACCTTCTGAACCACGTCCCACCATGCTGGTGTTTGACTCCGGCGTCGGTGGGTTGTCGGTCTATGACGAGATCCGGCCTCTCTTACGGGAACTCGATTACATTGATGCATTCGATAACGTCGCttt                                   >  2:205613/1‑140 (MQ=255)                                                                                                   tGTCTGGCAGCTACACCTTCTGAACCACGTCCCACCATGCTGGTGTTTGACTCCGGCGTCGGTGGGTTGTCGGTCTATGACGAGATCCGGCATCTCTTACCGGATCTCCATTACATTTATGCTTTCGATAACGTCGCttt                                   >  2:588934/1‑140 (MQ=255)                                                                                                   tGTCTGGCAGCTACACCTTCTGAACCACGTCCCACCATGCTGGTGTTTGACTCCGGCGTCGGTGGGTTGTCGGTCTATGACGAGATCCGGCATCTCTTACCGGATCTCCATTACATTTATGCTTTCGATAACGTCGCttt                                   >  1:395808/1‑140 (MQ=255)                                                                                                   tGTCTGGCAGCTACACCTTCTGAACCACGTCCCACCATGCTGGTGTTTGACTCCGGCGTCGGTGGGTTGTCGGTCTATGACGAGATCCGGCATCTCTTACCGGATCTCCATTACATTTATGCTTTCGATAACGTCGCttt                                   >  1:553486/1‑140 (MQ=255)                                                                                                    gTCTGGCAGCTACACCTTCTGAACCACGTCCCACCATGCTGGTGTTTGACTCCGGCGTCGGTGGGTTGTCGGTCTATGACGAGATCCGGCATCTCTTACCGGATCTCCATTACATTTATGCTTTCGATAACGTCGCTTTc                                  >  2:301926/1‑140 (MQ=255)                                                                                                     tCTGGCAGCTACACCTTCTGAACCACGTCCCACCATGCTGGTGTTTGACTCCGGCGTCGGTGGGTTGTCGGTCTATGACGAGATCCGGCATCTCTTACCGGATCTCCATTACATTTATGCTTTCGATAACGTCGCTTTcc                                 >  2:26464/1‑140 (MQ=255)                                                                                                        ggCAGCTACACCTTCTGAACCACGTCCCACCATGCTGGTGTTTGACTCCGGCGTCGGTGGGTTGTCGGTCTATGACGAGATCCGGCATCTCTTACCGGATCTCCATTACATTTATGCTTTCGATAACGTCGCTTTCCCGt                              >  1:335735/1‑140 (MQ=255)                                                                                                           aGCTACACCTTCTGAACCACGTCCCACCATGCTGGTGTTTGACTCCGGCGTCGGTGGGTTGTCGGTCTATGACGAGATCCGGCATCTCTTACCGGATCTCCATTACATTTATGCTTTCGATAACGTCGCTTTCCCGTATg                           <  1:698031/140‑1 (MQ=255)                                                                                                             cTACACCTTCTGAACCACGTCCCACCATGCTGGTGTTTGACTCCGGCGTCGGTGGGTTGTCGGTCTATGACGAGATCCGGCATCTCTTACCGGATCTCCATTACATTTATGCTTTCGATAACGTCGCTTTCCCGTATGGc                         >  2:719807/1‑140 (MQ=255)                                                                                                                   cTTCTGAACCACGTCCCACCATGCTGGTGTTTGACTCCGGCGTCGGTGGGTTGTCGGTCTATGACGAGATCCGGCATCTCTTACCGGATCTCCATTACATTTATGCTTTCGATAACGTCGCTTTCCCGTATGGCGaaaaa                   >  2:80941/1‑140 (MQ=255)                                                                                                                      cTGAACCACGTCCCACCATGCTGGTGTTTGACTCCGGCGTCGGTGGGTTGTCGGTCTATGACGAGATCCGGCATCTCTTACCGGATCTCCATTACATTTATGCTTTAGATAACGTCGATTTCCCGTATGGCTAAAAaagc                >  2:333996/1‑140 (MQ=255)                                                                                                                       tGAACCACGTCCCACCATGCTGGTGTTTGACTCCGGCGTCGGTGGGTTGTCGGTCTATGACGAGATCCGGCATCTCTTACCGGATCTCCATTACATTTATGCTTTCGATAACGTCGCTTTCCCGTATGGCGAAAAaagcg               <  2:11388/140‑1 (MQ=255)                                                                                                                       tGAACCACGTCCCACCATGCTGGTGTTTGACTCCGGCGTCGGTGGGTTGTCGGTCTATGACGAGATCCGGCATCTCTTACCGGATCTCCATTACATTTATGCTTTCGATAACGTCGCTTTCCCGTATGGCGAAAAaagcg               <  1:80941/140‑1 (MQ=255)                                                                                                                        gAACCACGTCCCACCATGCTGGTGTTTGACTCCGGCGTCGGTGGGTTGTCGGTCTATGACGAGATCCGGCATCTCTTACCGGATCTCCATTACATTTATGCTTTCGATAACGTCGCTTTCCCGTATGGCGAAAAaagcga              <  2:434934/140‑1 (MQ=255)                                                                                                                         aaCCACGTCCCACCATGCTGGTGTTTGACTCCGGCGTCGGTGGGTTGTCGGTCTATGACGAGATCCGGCATCTCTTACCGGATCTCCATTACATTTATGCTTTCGATAACGTCGCTTTCCCGTATGGCGAAAAaagcgaa             <  1:26464/140‑1 (MQ=255)                                                                                                                         aaCCACGTCCCACCATGCTGGTGTTTGACTCCGGCGTCGGTGGGTTGTCGGTCTATGACGAGATCCGGCATCTCTTACCGGATCTCCATTACATTTATGCTTTCGATAACGTCGCTTTCCCGTATGGCGAAAAaagcgaa             <  1:183435/140‑1 (MQ=255)                                                                                                                           ccACGTCCCACCATGCTGGTGTTTGACTCCGGCGTCGGTGGGTTGTCGGTCTATGACGAGATCCGGCATCTCTTACCGGATCTCCATTACATTTATGCTTTCGATAACGTCGCTTTCCCGTATGGCGAAAAaagcgaagc           <  1:719807/140‑1 (MQ=255)                                                                                                                                 cccACCATGCTGGTGTTTGACTCCGGCGTCGGTGGGTTGTCGGTCTATGACGAGATCCGGCATCTCTTACCGGATCTCCATTACATTTATGCTTTCGATAACGTCGCTTTCCCGTATGGCGAAAAaagcga              <  2:444596/131‑1 (MQ=255)                                                                                                                                 cccACCATGCTGGTGTTTGACTCCGGCGTCGGTGGGTTGTCGGTCTATGACGAGATCCGGCATCTCTTACCGGATCTCCATTACATTTATGCTTTCGATAACGTCGCTTTCCCGTATGGCGAAAAaagcga              >  1:444596/1‑131 (MQ=255)                                                                                                                                    aCCATGCTGGTGTTTGACTCCGGCGTCGGTGGGTTGTCGGTCTATGACGAGATCCGGCATCTCTTACCGGATCTCCATTACATTTATGCTTTCGATAACGTCGCTTTCCCGTATGGCGAAAAAAGCGAAGCGTTTAttgt  <  2:449721/140‑1 (MQ=255)                                                                                                                                       |                                                                                                                                         CACCTCTCACCTGACAGTTCGTGGTAAAATAGCCAACCTGTTCGACAAAGATTATGAGACAGTCTATGGCTACCAAACTGCAGGACGGGAATACACCTTGTCTGGCAGCTACACCTTCTGAACCACGTCCCACCGTGCTGGTGTTTGACTCCGGCGTCGGTGGGTTGTCGGTCTATGACGAGATCCGGCATCTCTTACCGGATCTCCATTACATTTATGCTTTCGATAACGTCGCTTTCCCGTATGGCGAAAAAAGCGAAGCGTTTATTGT  >  NC_000913/4165363‑4165633

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