Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,085,831	A→G	L411P (CTG→CCG)	plaP ←	putrescine importer, low affinity

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,085,831	0	A	G	100.0%	274.6 / NA	76	L411P (CTG→CCG)	plaP	putrescine importer, low affinity
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (38/38);  total (38/38)

GATGATTATGCTACGTCTTCGTACTGCGGAACCGGATTACGGAAGCTCTTGGTGACGCAAGCCAGGTAAATCAGACCGATAGCAGCCCAGATCAGACCCAGAACCATTGAGCTTTCTTCCAGGTTAACCCACAGCGCACCAACGGTCAGCGCACCACACATCGGCAGGAACAGATACTGGAAGTGATCTTTCAGCGTCTTGTTACGCTTCTCACGGATCCAGAACTGCGAGATGACCGACAGGTTAACGAAAGTGAACGCCACCA  >  NC_000913/2085699‑2085963                                                                                                                                     |                                                                                                                                     gatgatTATGCTACGTCTTCGTACTGCGGAACCGGATTACGGAAGCTCTTGGTGACGCAAGCCAGGTAAATCAGACCGATAGCAGCCCAGATCAGACCCAGAACCATTGAGCTTTCTTCCAGGTTAACCCAAGGCGCAc                                                                                                                                >  1:1233873/1‑139 (MQ=255)          gCTACGTCTTCGTACTGCGGAACCGGATTACGGAAGCTCTTGGTGACGCAAGCCAGGTAAATCAGACCGATAGCAGCCCAGATCAGACCCAGAACCATTGAGCTTTCTTCCAGGTTAACCCACGGCGCACCAACGGTCa                                                                                                                       <  1:77198/139‑1 (MQ=255)          gCTACGTCTTCGTACTGCGGAACCGGATTACGGAAGCTCTTGGTGACGCAAGCCAGGTAAATCAGACCGATAGCAGCCCAGATCAGACCCAGAACCATTGAGCTTTCTTCCAGGTTAACCCACGGCGCACCAACGGTCa                                                                                                                       <  2:192604/139‑1 (MQ=255)          gCTACGTCTTCGTACTGCGGAACCGGATTACGGAAGCTCTTGGTGACGCAAGCCAGGTAAATCAGACCGATAGCAGCCCAGATCAGACCCAGAACCATTGAGCTTTCTTCCAGGTTAACCCACGGCGCACCAACGGTCa                                                                                                                       <  2:1457664/139‑1 (MQ=255)              cGTCTTCGTACTGCGGAACCGGATTACGGAAGCTCTTGGTGACGCAAGCCAGGTAAATCAGACCGATAGCAGCCCAGATCAGACCCAGAACCATTGAGCTTTCTTCCAGGTTAACCCACGGCGCACCAACGGTCa                                                                                                                       <  1:151839/135‑1 (MQ=255)              cGTCTTCGTACTGCGGAACCGGATTACGGAAGCTCTTGGTGACGCAAGCCAGGTAAATCAGACCGATAGCAGCCCAGATCAGACCCAGAACCATTGAGCTTTCTTCCAGGTTAACCCACGGCGCACCAACGGTCa                                                                                                                       >  2:151839/1‑135 (MQ=255)                           cGGAACCGGATTACGGAAGCTCTTGGTGACGCAAGCCAGGTAAATCAGACCGATAGCAGCCCAGATCAGACCCAGAACCATTGAGCTTTCTTCCAGGTTAACCCACGGCGCACCAACGGTCAGCGCACCACACATCGGc                                                                                                      >  2:967974/1‑139 (MQ=255)                           cGGAACCGGATTACGGAAGCTCTTGGTGACGCAAGCCAGGTAAATCAGACCAATAGCAGCCCAGATCAGACCCAGAACCATTGAGCTTTCTTCCAGGTTAACCCACGGCGCACCa                                                                                                                              >  1:495430/1‑115 (MQ=255)                           cGGAACCGGATTACGGAAGCTCTTGGTGACGCAAGCCAGGTAAATCAGACCAATAGCAGCCCAGATCAGACCCAGAACCATTGAGCTTTCTTCCAGGTTAACCCACGGCGCACCa                                                                                                                              <  2:495430/115‑1 (MQ=255)                           cGGAACCGGATTACGGAAGCTCTTGGTGACGCAAGCCAGGTAAATCAGACCAATAGCAGCCCAGATCAGACCCAGAACCATTGAGCTTTCTTCCAGGTTAACCCACGGCGCACCa                                                                                                                              >  1:1156904/1‑115 (MQ=255)                           cGGAACCGGATTACGGAAGCTCTTGGTGACGCAAGCCAGGTAAATCAGACCAATAGCAGCCCAGATCAGACCCAGAACCATTGAGCTTTCTTCCAGGTTAACCCACGGCGCACCa                                                                                                                              <  2:1156904/115‑1 (MQ=255)                              aaCCGGATTACGGAAGCTCTTGGTGACGCAAGCCAGGTAAATCAGACCGATAGCAGCCCAGATCAGACCCAGAACCATTGAGCTTTCTTCCAGGTTAACCCACGGCGCACCAACGGTCAGCGCACCACACATCGGCAgg                                                                                                   >  2:1434832/1‑139 (MQ=255)                                ccGGATTACGGAAGCTCTTGGTGACGCAAGCCAGGTAAATCAGACCGATAGCAGCCCAGATCAGACCCAGAACCATTGAGCTTTCTTCCAGGTTAACCCACGGCGCACCAACGGTCAGCGCACCACACATCGGCAGGaa                                                                                                 <  1:484332/139‑1 (MQ=255)                                ccGGATTACGGAAGCTCTTGGTGACGCAAGCCAGGTAAATCAGACCGATAGCAGCCCAGATCAGACCCAGAACCATTGAGCTTTCTTCCAGGTTAACCCACGGCGCACCAACGGTCAGCGCACCACACATCGGCAGGaa                                                                                                 <  1:322157/139‑1 (MQ=255)                                             gCTCTTGGTGACGCAAGCCAGGTAAATCAGACCGATAGCAGCCCAGATCAGACCCAGAACCATTGAGCTTTCTTCCAGGTTAACCCACGGCGCa                                                                                                                                 <  1:1427957/94‑1 (MQ=255)                                             gCTCTTGGTGACGCAAGCCAGGTAAATCAGACCGATAGCAGCCCAGATCAGACCCAGAACCATTGAGCTTTCTTCCAGGTTAACCCACGGCGCa                                                                                                                                 >  2:1427957/1‑94 (MQ=255)                                             gCTCTTGGTGACGCAAGCCAGGTAAATCAGACCGATAGCAGCCCAGATCAGACCCAGAACCATTGAGCTTTCTTCCAGGTTAACCCACGGCGCACCAACGGTCAGCGCACCACACATCGGCAGGAACAGATACTGGAAg                                                                                    >  2:663001/1‑139 (MQ=255)                                              ctctTGGTGACGCAAGCCAGGTAAATCAGACCGATAGCAGCCCAGATCAGACCCAGAACCATTGAGCTTTCTTCCAGGTTAACCCACGGCGCACCAACGGTGAGCGCACCACACATCGGCAgg                                                                                                   >  1:289563/1‑123 (MQ=255)                                              ctctTGGTGACGCAAGCCAGGTAAATCAGACCGATAGCAGCCCAGATCAGACCCAGAACCATTGAGCTTTCTTCCAGGTTAACCCACGGCGCACCAACGGTGAGCGCACCACACATCGGCAgg                                                                                                   <  2:289563/123‑1 (MQ=255)                                              ctctTGGTGACGCAAGCCAGGTAAATCAGACCGATAGCAGCCCAGATCAGACCCAGAACCATTGAGCTTTCTTCCAGGTTAACCCACGGCGCACCAACGGTCAGCGCACCACACATCGGCAGGAACAGATACTGGAAGt                                                                                   >  2:188020/1‑139 (MQ=255)                                              ctctTGGTGACGCAAGCCAGGTAAATCAGACCGATAGCAGCCCAGATCAGACCCAGAACCATTGAGCTTTCTTCCAGGTTAACCCACGGCGCACCAACGGTCAGCGCACCACACATCGGCAGGAACAGATACTGGAAGt                                                                                   >  2:354535/1‑139 (MQ=255)                                              ctctTGGTGACGCAAGCCAGGTAAATCAGACCGATAGCAGCCCAGATCAGACCCAGAACCATTGAGCTTTCTTCCAGGTTAACCCACGGCGCACCAACGGTCAGCGCACCACACATCGGCAGGAACAGATACTGGAAGt                                                                                   >  2:897435/1‑139 (MQ=255)                                                ctTGGGGACGCAAGCCAGGTAAATCAGACCGATAGCAGCCCAGATCAGACCCAGAACCATTGAGCTTTCTTCCAGGTTAACCCACGGCGCACCAACGGTCAGCGCACCACACATCGGCAGGAACAGATACTGGAAGTGa                                                                                 >  2:1057158/1‑139 (MQ=255)                                                     tGACGCAAGCCAGGTAAATCAGACCGATAGCAGCCCAGATCAGACCCAGAACCATTGAGCTTTCTTCCAGGTTAACCCACGGCGCACCAACGGTCAGCGCACCACACATCGGCAGGAACAGATACTGGAAGTGATCttt                                                                            >  2:815248/1‑139 (MQ=255)                                                        cGCAAGCCAGGTAAATCAGACCGATAGCAGCCCAGATCAGACCCAGAACCATTGAGCTTTCTTCCAGGTTAACCCACGGCGCACCAACGGTCAGCGCACCACACATCGGCAGGAACAGATACTGGAAGTGATCTTTCAg                                                                         >  1:122560/1‑139 (MQ=255)                                                             gCCAGGTAAATCAGACCGATAGCAGCCCAGATCAGACCCAGAACCATTGAGCTTTCTTCCAGGTTAACCCACGGCGCACCAACGGTCAGCGCACCACACATCGGc                                                                                                      >  1:1491833/1‑105 (MQ=255)                                                             gCCAGGTAAATCAGACCGATAGCAGCCCAGATCAGACCCAGAACCATTGAGCTTTCTTCCAGGTTAACCCACGGCGCACCAACGGTCAGCGCACCACACATCGGc                                                                                                      <  2:1491833/105‑1 (MQ=255)                                                                    aaaTCAGACCGATAGCAGCCCAGATCAGACCCAGAACCATTGAGCTTTCTTCCAGGTTAACCCACGGCGCACCAACGGTCAGCGCACCACACATCGGCAGGAACAGATACTGGAAGTGATCTTTCAGCGTCTTGTTACg                                                             <  2:1233873/139‑1 (MQ=255)                                                                    aaaTCAGACCGATAGCAGCCCAGATCAGACCCAGAACCATTGAGCTTTCTTCCAGGTTAACCCACGGCGCACCAACGGTCAGCGCACCACACATCGGCAGGAACAGATACTGGAAGTGATCTTTCAGCGTCTTGTTACg                                                             <  1:897435/139‑1 (MQ=255)                                                                         aGACCGATAGCAGCCCAGATCAGACCCAGAACCATTGAGCTTTCTTCCAGGTTAACCCACGGCGCACCAACGGTCAGCGCACCACACATCGGCAGGAACAGATACTGGAAGTGATCTTTCAGCGTCTTGTTACGCTtct                                                        <  1:354535/139‑1 (MQ=255)                                                                            ccGATAGCAGCCCAGATCAGACCCAGAACCATTGAGCTTTCTTCCAGGTTAACCCACGGCGCACCAACGGTCAGCGCACCACACATCGGCAGGAACAGATACTGGAAGTGATCTTTCAGCGTCTTGTTACGCTTCTGAc                                                     <  1:967974/139‑1 (MQ=255)                                                                            ccGATAGCAGCCCAGATCAGACCCAGAACCATTGAGCTTTCTTCCAGGTTAACCCACGGCGCACCAACGGTCAGCGCACCACACATCGGCAGGAACAGATACTGGAAGTGATCTTTCAGCGTCTTGTTACGCTTCTCAc                                                     <  1:663001/139‑1 (MQ=255)                                                                            ccGATAGCAGCCCAGATCAGACCCAGAACCATTGAGCTTTCTTCCAGGTTAACCCACGGCGCACCAACGGTCAGCGCACCACACATCGGCAGGAACAGATACTGGAAGTGATCTTTCAGCGTCTTGTTACGCTTCTCAc                                                     <  1:1057158/139‑1 (MQ=255)                                                                                tAGCAGCCCAGATCAGACCCAGAACCATTGAGCTTTCTTCCAGGTTAACCCACGGCGCACCAACGGTCAGCGCACCACACATCGGCAGGa                                                                                                  >  1:1279720/1‑90 (MQ=255)                                                                                tAGCAGCCCAGATCAGACCCAGAACCATTGAGCTTTCTTCCAGGTTAACCCACGGCGCACCAACGGTCAGCGCACCACACATCGGCAGGa                                                                                                  <  2:1279720/90‑1 (MQ=255)                                                                                tAGCAGCCCAGATCAGACCCAGAACCATTGAGCTTTCTTCCAGGTTAACCCACGGCGCACCAACGGTCAGCGCACCACACATCGGCAGGAACAGATACTGGAAGTGATCTTTCAGCg                                                                       <  1:473196/117‑1 (MQ=255)                                                                                tAGCAGCCCAGATCAGACCCAGAACCATTGAGCTTTCTTCCAGGTTAACCCACGGCGCACCAACGGTCAGCGCACCACACATCGGCAGGAACAGATACTGGAAGTGATCTTTCAGCg                                                                       >  2:473196/1‑117 (MQ=255)                                                                                    agcCCAGATCAGACCCAGAACCATTGAGCTTTCTTCCAGGTTAACCCACGGCGCACCAACGGTCAGCGCACCACACATCGGCAGGAACAGATACTGGAAGTGATCTTTCAGCGTCTTGTTACGCTTCTCACGGATCCAg                                             <  2:122560/139‑1 (MQ=255)                                                                                     gcCCAGATCAGACCCAGAACCATTGAGCTTTCTTCCAGGTTAACCCACGGCGCACCAACGGTCAGCGCACCACACATCGGCAGGAACa                                                                                               >  1:403456/1‑88 (MQ=255)                                                                                     gcCCAGATCAGACCCAGAACCATTGAGCTTTCTTCCAGGTTAACCCACGGCGCACCAACGGTCAGCGCACCACACATCGGCAGGAACa                                                                                               <  2:403456/88‑1 (MQ=255)                                                                                     gcCCAGATCAGACCCAGAACCATTGAGCTTTCTTCCAGGTTAACCCACGGCGCACCAACGGTCAGCGCACCACACATCGGCAGGAACAGATACTGGAAGTGa                                                                                 >  1:1108619/1‑102 (MQ=255)                                                                                     gcCCAGATCAGACCCAGAACCATTGAGCTTTCTTCCAGGTTAACCCACGGCGCACCAACGGTCAGCGCACCACACATCGGCAGGAACAGATACTGGAAGTGa                                                                                 <  2:1108619/102‑1 (MQ=255)                                                                                          gATCAGACCCAGAACCATTGAGCTTTCTTCCAGGTTAACCCACGGCGCACCAACGGTCAGCGCACCACACATCGGCAGGAACAGATACTGGAAGTGATCTTTCAGCGTCTTGTTACGCTTCTCACGGATCCAGAACTGc                                       >  1:540896/1‑139 (MQ=255)                                                                                           aTCAGACCCAGAACCATTGAGCTTTCTTCCAGGTTAACCCACGGCGCACCAACGGTCAGCGCACCACACATCGGCAGGAACAGATACTGGAAGTGATCTTTCAGCGTCTTGTTACGCTTCTCACGGATCCAGAACTGCg                                      >  1:579987/1‑139 (MQ=255)                                                                                                aCCCAGAACCATTGAGCTTTCTTCCAGGTTAACCCACGgcg                                                                                                                                   >  2:1585574/1‑41 (MQ=255)                                                                                                aCCCAGAACCATTGAGCTTTCTTCCAGGTTAACCCACGgcg                                                                                                                                   <  1:1585574/41‑1 (MQ=255)                                                                                                aCCCAGAACCATTGAGCTTTCTTCCAGGTTAACCCACGGCGCACCAACGGTCAGCGCACCACACATCGGCAGGAACAGATACTGGAAGTGATCTTTCAGCGTCTTGTTACGCTTCTCACGGATCCAg                                             >  1:768825/1‑127 (MQ=255)                                                                                                aCCCAGAACCATTGAGCTTTCTTCCAGGTTAACCCACGGCGCACCAACGGTCAGCGCACCACACATCGGCAGGAACAGATACTGGAAGTGATCTTTCAGCGTCTTGTTACGCTTCTCACGGATCCAg                                             <  2:768825/127‑1 (MQ=255)                                                                                                aCCCAGAACCATTGAGCTTTCTTCCAGGTTAACCCACGGCGCACCAACGGTCAGCGCACCACACATCGGCAGGAACAGATACTGGAAGTGATCTTTCAGCGTCTTGTTACGCTTCTCACGGATCCAg                                             >  1:603821/1‑127 (MQ=255)                                                                                                aCCCAGAACCATTGAGCTTTCTTCCAGGTTAACCCACGGCGCACCAACGGTCAGCGCACCACACATCGGCAGGAACAGATACTGGAAGTGATCTTTCAGCGTCTTGTTACGCTTCTCACGGATCCAg                                             <  2:603821/127‑1 (MQ=255)                                                                                                aCCCAGAACCATTGAGCTTTCTTCCAGGTTAACCCACGGCGCACCAACGGTCAGCGCACCACACATCGGCAGGAACAGATACTGGAAGTGATCTTTCAGCGTCTTGTTACGCTTCTCACGGATCCAGAACTGCg                                      >  1:1142352/1‑134 (MQ=255)                                                                                                aCCCAGAACCATTGAGCTTTCTTCCAGGTTAACCCACGGCGCACCAACGGTCAGCGCACCACACATCGGCAGGAACAGATACTGGAAGTGATCTTTCAGCGTCTTGTTACGCTTCTCACGGATCCAGAACTGCg                                      <  2:1142352/134‑1 (MQ=255)                                                                                                aCCCAGAACCATTGAGCTTTCTTCCAGGTTAACCCACGGCGCACCAACGGTCAGCGCACCACACATCGGCAGGAACAGATACTGGAAGTGATCTTTCAGCGTCTTGTTACGCTTCTCACGGATCCAGAACTGCGAGATg                                 >  2:4201/1‑139 (MQ=255)                                                                                                     gAACCATTGAGCTTTCTTCCAGGTTAACCCACGGCGCACCAACGGTCAGCGCACCACACATCGGCAGGAACAGATACTGGAAGTGATCTTTCAGCGTCTTGTTACGCTTCTCACGGATCCAGAACTGCgaga                                   <  2:3817/132‑1 (MQ=255)                                                                                                     gAACCATTGAGCTTTCTTCCAGGTTAACCCACGGCGCACCAACGGTCAGCGCACCACACATCGGCAGGAACAGATACTGGAAGTGATCTTTCAGCGTCTTGTTACGCTTCTCACGGATCCAGAACTGCgaga                                   >  1:3817/1‑132 (MQ=255)                                                                                                     gAACCATTGAGCTTTCTTCCAGGTTAACCCACGGCGCACCAACGGTCAGCGCACCACACATCGGCAGGAACAGATACTGGAAGTGATCTTTCAGCGTCTTGTTACGCTTCTCACGGATCCAGAACTGCGAGATGACCGa                            <  2:579987/139‑1 (MQ=255)                                                                                                     gAACCATTGAGCTTTCTTCCAGGTTAACCCACGGCGCACCAACGGTCAGCGCACCACACATCGGCAGGAACAGATACTGGAAGTGATCTTTCAGCGTCTTGTTACGCTTCTCACGGATCCAGAACTGCGAGATGACCGa                            >  2:738940/1‑139 (MQ=255)                                                                                                     gAACCATTGAGCTTTCTTCCAGGTTAACCCACGGCGCACCAACGGTCAGCGCACCACACATCGGCAGGAACAGATACTGGAAGTGATCTTTCAGCGTCTTGTTACGCTTCTCACGGATCCAGAACTGCGAGATGACCGa                            <  1:738940/139‑1 (MQ=255)                                                                                                          aTTGAGCTTTCTTCCAGGTTAACCCACGGCGCACCAACGGTCAGCGCACCACACATCGGCAGGAACAGATACTGGAAGTGATCTTTCAGCGTCTTGTTACGCTTCTCAc                                                     <  2:21653/109‑1 (MQ=255)                                                                                                          aTTGAGCTTTCTTCCAGGTTAACCCACGGCGCACCAACGGTCAGCGCACCACACATCGGCAGGAACAGATACTGGAAGTGATCTTTCAGCGTCTTGTTACGCTTCTCAc                                                     >  1:21653/1‑109 (MQ=255)                                                                                                          aTTGAGCTTTCTTCCAGGTTAACCCACGGCGCACCAACGGTCAGCGCACCACACATCGGCAGGAACAGATACTGGAAGTGATCTTTCAGCGTCTTGTTACGCTTCTCACGGATCCAGAACTGCGAGATGACCGACAGGt                       <  1:1434832/139‑1 (MQ=255)                                                                                                          aTTGAGCTTTCTTCCAGGTTAACCCACGGCGCACCAACGGTCAGCGCACCACACATCGGCAGGAACAGATACTGGAAGTGATCTTTCAGCGTCTTGTTACGCTTCTCACGGATCCAGAACTGCGAGATGACCGACAGGt                       >  1:531893/1‑139 (MQ=255)                                                                                                           ttGAGCTTTCTTCCAGGTTAACCCACGGCGCACCAACGGTCAGCGCACCACACATCGGCAGGAACAGATACTg                                                                                        >  2:361705/1‑73 (MQ=255)                                                                                                           ttGAGCTTTCTTCCAGGTTAACCCACGGCGCACCAACGGTCAGCGCACCACACATCGGCAGGAACAGATACTg                                                                                        <  1:361705/73‑1 (MQ=255)                                                                                                           ttGAGCTTTCTTCCAGGTTAACCCACGGCGCACCAACGGTCAGCGCACCACACATCGGCAGGAACAGATACTGGAAGTGATCTTTCAGCGTCTTGTTACGCTTCTCAc                                                     <  1:1083410/108‑1 (MQ=255)                                                                                                           ttGAGCTTTCTTCCAGGTTAACCCACGGCGCACCAACGGTCAGCGCACCACACATCGGCAGGAACAGATACTGGAAGTGATCTTTCAGCGTCTTGTTACGCTTCTCAc                                                     >  2:1083410/1‑108 (MQ=255)                                                                                                                  ttcttcCAGGTTAACCCACGGCGCACCAACGGTCAGCGCACCACACATCGGCAGGAACAGATACTGGAAGTGATCTTTCAGCg                                                                       >  1:717250/1‑83 (MQ=255)                                                                                                                  ttcttcCAGGTTAACCCACGGCGCACCAACGGTCAGCGCACCACACATCGGCAGGAACAGATACTGGAAGTGATCTTTCAGCg                                                                       <  2:717250/83‑1 (MQ=255)                                                                                                                   tcttcCAGGTTAACCCACGGCGCACCAACGGTCAGCGCACCACACATCGGCAGGAACAGATACTGGAAGTGATCTTTCAGCGTCTTGTTACGCTTCTCACGGATCCAGAACTGCGAGATGACCGACAGGTTAACGAAAg              <  1:4201/139‑1 (MQ=255)                                                                                                                             tAACCCACGGCGCACCAACGGTCAGCGCACCACACATCgg                                                                                                       <  2:166236/40‑1 (MQ=255)                                                                                                                             tAACCCACGGCGCACCAACGGTCAGCGCACCACACATCgg                                                                                                       >  1:166236/1‑40 (MQ=255)                                                                                                                               aCCCACGGCGCACCAACGGTCAGCGCACCACACATCGGCAGGAACAGATACTGGAAGTGATCTTTCAGCGTCTTGTTGCGCTTCTCAc                                                     >  1:1057478/1‑88 (MQ=255)                                                                                                                               aCCCACGGCGCACCAACGGTCAGCGCACCACACATCGGCAGGAACAGATACTGGAAGTGATCTTTCAGCGTCTTGTTACGCTTCTCAc                                                     <  2:1057478/88‑1 (MQ=255)                                                                                                                               aCCCACGGCGCACCAACGGTCAGCGCACCACACATCGGCAGGAACAGATACTGGAAGTGATCTTTCAGCGTCTTGTTACGCTTCTCACGGATCCAGAACTGCGAGATGACCGACAGGTTaa                    <  2:172164/121‑1 (MQ=255)                                                                                                                               aCCCACGGCGCACCAACGGTCAGCGCACCACACATCGGCAGGAACAGATACTGGAAGTGATCTTTCAGCGTCTTGTTACGCTTCTCACGGATCCAGAACTGCGAGATGACCGACAGGTTaa                    >  1:172164/1‑121 (MQ=255)                                                                                                                               aCCCACGGCGCACCAACGGTCAGCGCACCACACATCGGCAGGAACAGATACTGGAAGTGATCTTTCAGCGTCTTGTTACGCTTCTCACGGATCCAGAACTGCGAGATGACCGACAGGTTAACGAAAGTGAACGccacca  >  2:1253444/1‑139 (MQ=255)                                                                                                                                     |                                                                                                                                     GATGATTATGCTACGTCTTCGTACTGCGGAACCGGATTACGGAAGCTCTTGGTGACGCAAGCCAGGTAAATCAGACCGATAGCAGCCCAGATCAGACCCAGAACCATTGAGCTTTCTTCCAGGTTAACCCACAGCGCACCAACGGTCAGCGCACCACACATCGGCAGGAACAGATACTGGAAGTGATCTTTCAGCGTCTTGTTACGCTTCTCACGGATCCAGAACTGCGAGATGACCGACAGGTTAACGAAAGTGAACGCCACCA  >  NC_000913/2085699‑2085963

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: ‑