Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,621,307	A→G	E37E (GAA→GAG)	purM →	phosphoribosylaminoimidazole synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,621,307	0	A	G	100.0%	251.0 / NA	70	E37E (GAA→GAG)	purM	phosphoribosylaminoimidazole synthetase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (37/33);  total (37/33)

TAACGCGTGGGGACCCAAGCAGTGACCGATAAAACCTCTCTTAGCTACAAAGATGCCGGTGTTGATATTGACGCGGGTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAAGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCA‑TTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACTGACGGCGTAGGTACCAAGCTGCGTCTGGCAATGGACTTAAAACGTCACGAC  >  NC_000913/2621176‑2621439                                                                                                                                    |                                                                                                                                      tAACGCGTGGGGACCCAAGCAGTGACCGATAAAACCTCTCTTAGCTACAAAGATGCCGGTGTTGATATTGACGCGGGTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATgg                                                                                                                                <  1:1878360/139‑1 (MQ=255)    cgcgTGGGGACCCAAGCAGTGACCGATAAAACCTCTCTTAGCTACAAAGATGCCGGTGTTGATATTGACGCGGGTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCg                                                                                                                             >  2:636980/1‑139 (MQ=255)        tGGGGACCCAAGCAGTGACCGATAAAACCTCTCTTAGCTACAAAGATGCCGGTGTTGATATTGACGCGGGTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCt                                                                                                                         >  1:661270/1‑139 (MQ=255)         ggggACCCAAGCAGTGACCGATAAAACCTCTCTTAGCTACAAAGATGCCGGTGTTGATATTGACGCGGGTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTg                                                                                                                        >  1:1543733/1‑139 (MQ=255)              cccAAGCAGTGACCGATAAAACCTCTCTTAGCTACAAAGATGCCGGTGTTGATATTGACGCGGGTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGggcgg                                                                                                                   <  2:661270/139‑1 (MQ=255)                 aaGCAGTGACCGATAAAACCTCTCTTAGCTACAAAGATGCCGGTGTTGATATTGACGCGGGTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCgg                                                                                                                            <  1:64936/127‑1 (MQ=255)                 aaGCAGTGACCGATAAAACCTCTCTTAGCTACAAAGATGCCGGTGTTGATATTGACGCGGGTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCgg                                                                                                                            >  2:64936/1‑127 (MQ=255)                 aaGCAGTGACCGATAAAACCTCTCTTAGCTACAAAGATGCCGGTGTTGATATTGACGCGGGTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCtt                                                                                                                >  2:422710/1‑139 (MQ=255)                     aGTGACCGATAAAACCTCTCTTAGCTACAAAGATGCCGGTGTTGATATTGACGCGGGTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGt                                                                                                            <  1:636980/139‑1 (MQ=255)                         aCCGATAAAACCTCTCTTAGCTACAAAGATGCCGGTGTTGATATTGACGCGGGTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTgcgc                                                                                                        >  1:480218/1‑139 (MQ=255)                         aCCGATAAAACCTCTCTTAGCTACAAAGATGCCGGTGTTGATATTGACGCGGGTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTgcgc                                                                                                        >  2:795606/1‑139 (MQ=255)                         aCCGATAAAACCTCTCTTAGCTACAAAGATGCCGGTGTTGATATTGACGCGGGTAATGATCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCAGGAGGTGATGGGCGGTCTGGGGGGATTGGGTgcgc                                                                                                        >  1:20928/1‑139 (MQ=255)                          ccGATAAAACCTCTCTTAGCTACAAAGATGCCGGTGTTGATATTGACGCGGGTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCt                                                                                                       <  2:1543733/139‑1 (MQ=255)                            gATAAAACCTCTCTTAGCTACAAAGATGCCGGTGTTGATATTGACGCGGGTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGt                                                                                                            >  1:1982393/1‑132 (MQ=255)                            gATAAAACCTCTCTTAGCTACAAAGATGCCGGTGTTGATATTGACGCGGGTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGt                                                                                                            <  2:1982393/132‑1 (MQ=255)                                     tctctTAGCTACAAAGATGCCGGTGTTGATATTGACGCGGGTAATTCTCTGGCTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTc                                                                                                               >  2:308614/1‑120 (MQ=255)                                     tctctTAGCTACAAAGATGCCGGTGTTGATATTGACGCGGGTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTc                                                                                                               <  1:308614/120‑1 (MQ=255)                                           aGCTACAAAGATGCCGGTGTTGATATTGACGCGGGTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGa                                                                                                                                   <  1:90272/94‑1 (MQ=255)                                           aGCTACAAAGATGCCGGTGTTGATATTGACGCGGGTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGa                                                                                                                                   >  2:90272/1‑94 (MQ=255)                                            gCTACAAAGATGCCGGTGTTGATATTGACGCGGGTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTgcg                                                                                                         >  1:947810/1‑119 (MQ=255)                                            gCTACAAAGATGCCGGTGTTGATATTGACGCGGGTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTgcg                                                                                                         <  2:947810/119‑1 (MQ=255)                                             cTACAAAGATGCCGGTGTTGATATTGACGCGGGTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTgcg                                                                                                         <  1:432216/118‑1 (MQ=255)                                             cTACAAAGATGCCGGTGTTGATATTGACGCGGGTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTgcg                                                                                                         >  2:432216/1‑118 (MQ=255)                                             cTACAAAGATGCCGGTGTTGATATTGACGCGGGTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCa                                                                                                <  2:1434448/127‑1 (MQ=255)                                             cTACAAAGATGCCGGTGTTGATATTGACGCGGGTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGGGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCa                                                                                                >  1:1434448/1‑127 (MQ=255)                                              tACAAAGATGCCGGTGTTGATATTGACGCGGGTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATggg                                                                                                                               >  2:48014/1‑95 (MQ=255)                                              tACAAAGATGCCGGTGTTGATATTGACGCGGGTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATggg                                                                                                                               <  1:48014/95‑1 (MQ=255)                                                   aGATGCCGGTGTTGATATTGACGCGGGTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCA‑TTGCCGCAAAAATATCGt                                                                             >  2:44829/1‑139 (MQ=255)                                                     aTGCCGGTGTTGATATTGACGCGGGTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCA‑TTGCCGCAAAAATATCGTGa                                                                           <  1:44829/139‑1 (MQ=255)                                                       gCCGGTGTTGATATTGACGCGGGTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCA‑TTGCCGCAAAAATATCGTGAAc                                                                         >  1:257933/1‑139 (MQ=255)                                                        ccGGTGTTGATATTGACGCGGGTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCA‑TTGCCGCAAAAATATCg                                                                              >  1:1540473/1‑133 (MQ=255)                                                        ccGGTGTTGATATTGACGCGGGTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCA‑TTGCCGCAAAAATATCg                                                                              <  2:1540473/133‑1 (MQ=255)                                                           gtgtTGATATTGACGCGGGTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCA‑TTGCCGCAAAAATATCGTGAACCCGt                                                                     >  1:197081/1‑139 (MQ=255)                                                                        cgcgGGTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCA‑TTGCCGCAAAAATATCGTGAACCCGTg                                                                    >  1:79125/1‑127 (MQ=255)                                                                        cgcgGGTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCA‑TTGCCGCAAAAATATCGTGAACCCGTg                                                                    <  2:79125/127‑1 (MQ=255)                                                                            ggTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCATTTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACt                                                     >  2:1213106/1‑139 (MQ=255)                                                                              tAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCA‑TTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACTGAc                                                  <  2:480218/139‑1 (MQ=255)                                                                              tAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCA‑TTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACTGAc                                                  <  2:257933/139‑1 (MQ=255)                                                                                   ctctGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCA‑TTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACTGACGGCGt                                             <  1:642490/139‑1 (MQ=255)                                                                                     ctGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCA‑TTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACTGACGGCGTAg                                           >  2:1866593/1‑139 (MQ=255)                                                                                       ggTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCA‑TTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACTGACg                                                 >  1:998138/1‑131 (MQ=255)                                                                                       ggTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCA‑TTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACTGACg                                                 <  2:998138/131‑1 (MQ=255)                                                                                       ggTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCA‑TTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACTGACGGCGTAGGt                                         >  1:1356498/1‑139 (MQ=255)                                                                                        gTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCA‑TTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACTGACGGCGTAGGTa                                        >  2:185484/1‑139 (MQ=255)                                                                                                  tCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCA‑TTGCCGCa                                                                                       <  2:1647384/82‑1 (MQ=255)                                                                                                  tCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCA‑TTGCCGCa                                                                                       >  1:1647384/1‑82 (MQ=255)                                                                                                  tCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCA‑TTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACTGACGGCGTAGGTACCAAGCTGCg                              >  2:919650/1‑139 (MQ=255)                                                                                                  tCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCA‑TTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACTGACGGCGTAGGTACCAAGCTGCg                              >  1:1061514/1‑139 (MQ=255)                                                                                                  tCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCA‑TTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACTGACGGCGTAGGTACCAAGCTGCg                              >  2:364453/1‑139 (MQ=255)                                                                                                   cAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATggg                                                                                                                               >  2:317563/1‑42 (MQ=255)                                                                                                   cAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATggg                                                                                                                               <  1:317563/42‑1 (MQ=255)                                                                                                    aaaGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCA‑TTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACTGACGGCGTAGGTACCAAGCTGCGTc                            <  1:795606/139‑1 (MQ=255)                                                                                                         cGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCA‑TTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACTGACGGCGTAGGTACCAAGCTGCGTCTGGCa                       <  1:422710/139‑1 (MQ=255)                                                                                                           tAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCA‑TTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACTGACGGCGTAGGTACCAAGCTGCGt                             >  1:1293177/1‑131 (MQ=255)                                                                                                           tAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCA‑TTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACTGACGGCGTAGGTACCAAGCTGCGt                             <  2:1293177/131‑1 (MQ=255)                                                                                                           tAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCA‑TTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACTGACGGCGTAGGTACCAAGCTGCGTCTGGCaa                      >  2:165363/1‑138 (MQ=255)                                                                                                           tAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCA‑TTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACTGACGGCGTAGGTACCAAGCTGCGTCTGGCaa                      <  1:165363/138‑1 (MQ=255)                                                                                                            aGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCA‑TTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACTGACGGCGTAGGTACCAAGCTGCGTCTGGCAATg                    >  2:563197/1‑139 (MQ=255)                                                                                                                aagaaAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCA‑TTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACTGACGGCGTAGGTACCAAGCTGCGTCTGGCAATGGACt                >  2:1969332/1‑139 (MQ=255)                                                                                                                  gaaAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCA‑TTGCCGCaaaaa                                                                                   <  1:2041767/70‑1 (MQ=255)                                                                                                                  gaaAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCA‑TTGCCGCaaaaa                                                                                   >  2:2041767/1‑70 (MQ=255)                                                                                                                     aaCGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCA‑TTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCa                                                       <  2:652609/95‑1 (MQ=255)                                                                                                                     aaCGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCA‑TTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCa                                                       >  1:652609/1‑95 (MQ=255)                                                                                                                     aaCGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCA‑TTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACTGACGGCGTAGGTACCAAGCTGCGTCTGGCAATGGACTTaaaa           <  1:364453/139‑1 (MQ=255)                                                                                                                        gcgTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCA‑TTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACTGACGGCGTAGGTACCAAGCTGCGTCTGGCAATGGACTTAAAACGt        <  1:1866593/139‑1 (MQ=255)                                                                                                                          gtcgtcCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCA‑TTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACTGACGGCGTAGGTACCAAGCTGCGTCTGGCaa                      >  2:664776/1‑123 (MQ=255)                                                                                                                          gtcgtcCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCA‑TTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACTGACGGCGTAGGTACCAAGCTGCGTCTGGCaa                      <  1:664776/123‑1 (MQ=255)                                                                                                                             gtcCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCA‑TTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACTGACGGCGTAGGTACCAAGCTGCGTCTGGCAATGGACTTAAAACGTCACGa   <  1:563197/139‑1 (MQ=255)                                                                                                                             gtcCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCA‑TTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACTGACGGCGTAGGTACCAAGCTGCGTCTGGCAATGGACTTAAAACGTCACGa   <  2:1356498/139‑1 (MQ=255)                                                                                                                               ccGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCATTTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACTGACGGCGTAGGTACCAAGCTGCGTCTGGCAATGGACTTAAAACGTCACGac  <  1:1213106/139‑1 (MQ=255)                                                                                                                                    |                                                                                                                                      TAACGCGTGGGGACCCAAGCAGTGACCGATAAAACCTCTCTTAGCTACAAAGATGCCGGTGTTGATATTGACGCGGGTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAAGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCA‑TTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACTGACGGCGTAGGTACCAAGCTGCGTCTGGCAATGGACTTAAAACGTCACGAC  >  NC_000913/2621176‑2621439

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