Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,559,680	G→A	G77G (GGC→GGT)	ddpC ←	D,D‑dipeptide ABC transporter permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,559,680	0	G	A	100.0%	220.9 / NA	68	G77G (GGC→GGT)	ddpC	D,D‑dipeptide ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (31/37);  total (31/37)

AAGCACACCGGATAGACATCCGAGTAGCGAACCAATCATACCCGCAATGGCAACCACGACTAATCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGTGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGC  >  NC_000913/1559548‑1559804                                                                                                                                     |                                                                                                                             aaGCACACCGGATAGACATCCGAGTAGCGAACCAATCATACCCGCAATGGCAACCACGACTAATCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAAc                                                                                                                        >  1:539792/1‑139 (MQ=255)        ccGGATAGACATCCGAGTAGCGAACCAATCATACCCGCAATGGCAACCACGACTAATCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTgcg                                                                                                                 <  1:1355817/139‑1 (MQ=255)           gATAGACATCCGAGTAGCGAACCAATCATACCCGCAATGGCAACCACGACTAATCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACTAGTGcgccg                                                                                                              <  2:762529/139‑1 (MQ=255)                 cATCCGAGTAGCGAACCAATCATACCCGCAATGGCAACCACGACTAATCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAg                                                                                                        <  2:1541546/139‑1 (MQ=255)                     cGAGTAGCGAACCAATCATACCCGCAATGGCAACCACGACTAATCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAggcgg                                                                                                    >  2:1303645/1‑139 (MQ=255)                              aaCCAATCATACCCGCAATGGCAACCACGACTAATCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAgg                                                                                                       >  1:747361/1‑127 (MQ=255)                              aaCCAATCATACCCGCAATGGCAACCACGACTAATCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAgg                                                                                                       <  2:747361/127‑1 (MQ=255)                              aaCCAATCATACCCGCAATGGCAACCACGACTAATCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAggcg                                                                                           >  2:781684/1‑139 (MQ=255)                               aCCAATCATACCCGCAATGGCAACCACGACTAATCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAggcgg                                                                                          <  1:1303645/139‑1 (MQ=255)                                ccAATCATACCCGCAATGGCAACCACGACTAATCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAggcggg                                                                                         <  2:299214/139‑1 (MQ=255)                                   aTCATACCCGCAATGGCAACCACGACTAATCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAGAggcg                                                                                           <  1:376309/134‑1 (MQ=255)                                   aTCATACCCGCAATGGCAACCACGACTAATCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAggcg                                                                                           >  2:376309/1‑134 (MQ=255)                                            gCAATGGCAACCACGACTAATCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCaa                                                                             <  2:109890/139‑1 (MQ=255)                                            gCAATGGCAACCACGACTAATCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCaa                                                                             <  1:506921/139‑1 (MQ=255)                                             cAATGGCAACCACGACTAATCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAAt                                                                            >  2:273128/1‑139 (MQ=255)                                             cAATGGCAACCACGACTAATCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAAt                                                                            >  2:776806/1‑139 (MQ=255)                                             cAATGGCAACCACGACTAATCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAAt                                                                            >  1:714355/1‑139 (MQ=255)                                              aaTGGCAACCACGACTAATCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTa                                                                                   >  1:123251/1‑131 (MQ=255)                                              aaTGGCAACCACGACTAATCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTa                                                                                   <  2:123251/131‑1 (MQ=255)                                              aaTGGCAACCACGACTAATCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATg                                                                           >  2:270802/1‑139 (MQ=255)                                              aaTGGCAACCACGACTAATCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATg                                                                           >  1:1118758/1‑139 (MQ=255)                                              aaTGGCAACCACGACTAATCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGC‑AAAGGCGGGCGGTTAAATCAATgg                                                                          >  2:1287761/1‑139 (MQ=255)                                                  gCAACCACGACTAATCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGt                                                                       >  1:1588092/1‑139 (MQ=255)                                                   cAACCACGACTAATCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCaa                                                                             >  2:1019487/1‑132 (MQ=255)                                                   cAACCACGACTAATCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCaa                                                                             <  1:1019487/132‑1 (MQ=255)                                                    aaCCACGACTAATCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTg                                                                     >  1:468818/1‑139 (MQ=255)                                                           aCTAATCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAAt                                                                            <  1:1149417/125‑1 (MQ=255)                                                           aCTAATCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAAt                                                                            >  2:1149417/1‑125 (MQ=255)                                                           aCTAATCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCg                                                                        >  1:1262759/1‑129 (MQ=255)                                                           aCTAATCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCg                                                                        <  2:1262759/129‑1 (MQ=255)                                                            cTAATCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCaa                                                                             >  2:105862/1‑123 (MQ=255)                                                            cTAATCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCaa                                                                             <  1:105862/123‑1 (MQ=255)                                                            cTAATCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATg                                                             <  2:714355/139‑1 (MQ=255)                                                                tCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTc                                                         <  2:610118/139‑1 (MQ=255)                                                                    gCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAggcgg                                                                                          <  2:1587166/102‑1 (MQ=255)                                                                    gCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAggcgg                                                                                          >  1:1587166/1‑102 (MQ=255)                                                                          attgattgCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACg                                               <  1:302605/139‑1 (MQ=255)                                                                            tgattgCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGt                                             <  2:1588092/139‑1 (MQ=255)                                                                              attgCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGTGa                                           >  1:833962/1‑139 (MQ=255)                                                                              attgCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCC‑GCCACGGTGaa                                          >  2:1245800/1‑139 (MQ=255)                                                                                tgCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGTGaaa                                         <  2:468818/139‑1 (MQ=255)                                                                                  cTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGTGaaaaa                                       >  2:1071063/1‑139 (MQ=255)                                                                                             ccAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCC‑GCCACGGTGAAAAAATCATCATCAAt                           <  1:1245800/139‑1 (MQ=255)                                                                                              cAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGTGAAAAAATCATCATCAAt                           <  2:182992/139‑1 (MQ=255)                                                                                              cAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGTGAAAAAATCATCATCAAt                           <  1:958164/139‑1 (MQ=255)                                                                                                gTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGTGAAAAAATCATCATCAATAg                         <  1:855108/139‑1 (MQ=255)                                                                                                gTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGTGAAAAAATCATCATCAATAg                         <  1:781684/139‑1 (MQ=255)                                                                                                gTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGTGAAAAAATCATCATCAATAg                         <  2:1118758/139‑1 (MQ=255)                                                                                                gTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGTGAAAAAATCATCATCAATAg                         <  1:270802/139‑1 (MQ=255)                                                                                                gTACGCGGCTAAACAGATAGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGTGAAAAAATCATCATCAATAg                         >  1:895427/1‑139 (MQ=255)                                                                                                     cgGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGTGAAAAAATCATCATCAATAGCATTa                    <  2:833962/139‑1 (MQ=255)                                                                                                      ggCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATg                                                             <  2:864151/97‑1 (MQ=255)                                                                                                      ggCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATg                                                             >  1:864151/1‑97 (MQ=255)                                                                                                      ggCTAAACAGATCGCGTCCCACTTCATCGGTACC‑AACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGTGAAAAAATCATCATCAATAGCATTAAt                  <  1:1371762/139‑1 (MQ=255)                                                                                                       gCTAAACTGATGGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGTGAAAAAATCATCATCAATAGCATTAAt                  <  1:866474/139‑1 (MQ=255)                                                                                                       gCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGTGAAAAAATCATCATCAATAGCATTAAt                  <  2:1135393/139‑1 (MQ=255)                                                                                                                aTCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGTGAAAAAATCATCATCAATAGCATTAATACAatgatg         >  2:108066/1‑139 (MQ=255)                                                                                                                  cgcgTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAggcgg                                                                                          <  1:1448815/56‑1 (MQ=255)                                                                                                                  cgcgTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAggcgg                                                                                          >  2:1448815/1‑56 (MQ=255)                                                                                                                  cgcgTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGTGAAAAAATCATCATCAAt                           >  1:425106/1‑119 (MQ=255)                                                                                                                  cgcgTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGTGAAAAAATCATCATCAAt                           <  2:425106/119‑1 (MQ=255)                                                                                                                  cgcgTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGTGAAAAAATCATCATCAATAGCATTAATACAATGATGAc       >  2:576471/1‑139 (MQ=255)                                                                                                                     gTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGTGAAAAAATCATCATCAATAGCATTAGTACAATGATGACTgc    <  1:273128/139‑1 (MQ=255)                                                                                                                     gTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGTGAAAAAATCATCATCAATAGCATTAATACAATGATGACTgc    >  1:1390837/1‑139 (MQ=255)                                                                                                                     gTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGTGAAAAAATCATCATCAATAGCATTAATACAATGATGACTgc    <  1:108066/139‑1 (MQ=255)                                                                                                                      tCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGTGAAAAAATCATCATCAATAGCATTAATACAATGATGACTgcg   <  2:213999/139‑1 (MQ=255)                                                                                                                       cccACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGTGAAAAAATCATCATCAATAGCATTAATACAATGATGACTgcgc  >  2:1379161/1‑139 (MQ=255)                                                                                                                                 cGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGt                                                                       >  2:733560/1‑60 (MQ=255)                                                                                                                                 cGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGt                                                                       <  1:733560/60‑1 (MQ=255)                                                                                                                                     |                                                                                                                             AAGCACACCGGATAGACATCCGAGTAGCGAACCAATCATACCCGCAATGGCAACCACGACTAATCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGTGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGC  >  NC_000913/1559548‑1559804

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