Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,937,168	G→A	*297* (TAG→TAA)	rbsB →	D‑ribose ABC transporter periplasmic binding protein, ribose chemotaxis receptor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,937,168	0	G	A	100.0%	256.5 / NA	78	*297* (TAG→TAA)	rbsB	D‑ribose ABC transporter periplasmic binding protein, ribose chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (34/44);  total (34/44)

GCAAACTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAA  >  NC_000913/3937032‑3937292                                                                                                                                         |                                                                                                                             gCAAACTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAAtt                                                                                                                            <  2:1882885/139‑1 (MQ=255)         gcagcGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAg                                                                                                                    <  1:1066900/139‑1 (MQ=255)           agcGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCa                                                                                                                     >  2:1958751/1‑136 (MQ=255)           agcGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCa                                                                                                                     <  1:1958751/136‑1 (MQ=255)           agcGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGt                                                                                                                  >  1:698732/1‑139 (MQ=255)               aCTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTa                                                                                                              <  2:460754/139‑1 (MQ=255)                cTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTAt                                                                                                             <  1:793193/139‑1 (MQ=255)                 tATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTaa                                                                                                                        >  2:1419547/1‑127 (MQ=255)                 tATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTaa                                                                                                                        <  1:1419547/127‑1 (MQ=255)                   tCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGAc                                                                                                          >  2:1489631/1‑139 (MQ=255)                         aGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATg                                                                                                    >  1:305088/1‑139 (MQ=255)                           cTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGt                                                                                                  >  2:616558/1‑139 (MQ=255)                           cTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGt                                                                                                  <  1:616558/139‑1 (MQ=255)                           cTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGt                                                                                                  <  1:765500/139‑1 (MQ=255)                             aCCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTa                                                                                                              >  1:1784454/1‑125 (MQ=255)                             aCCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTa                                                                                                              <  2:1784454/125‑1 (MQ=255)                              cccGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGAc                                                                                               >  1:739000/1‑139 (MQ=255)                              cccGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGAc                                                                                               <  2:811217/139‑1 (MQ=255)                               ccGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTg                                                                                                       >  2:418261/1‑130 (MQ=255)                               ccGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTg                                                                                                       <  1:418261/130‑1 (MQ=255)                                    cAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAAt                                                                                         <  2:305088/139‑1 (MQ=255)                                       aTTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACg                                                                                      <  1:1941619/139‑1 (MQ=255)                                         tGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTc                                                                                    >  2:802653/1‑139 (MQ=255)                                         tGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTc                                                                                    >  2:1635738/1‑139 (MQ=255)                                           gcgcGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCAt                                                                                  <  1:1283983/139‑1 (MQ=255)                                            cgcgAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTa                                                                                                              >  1:1759203/1‑110 (MQ=255)                                            cgcgAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTa                                                                                                              <  2:1759203/110‑1 (MQ=255)                                            cgcgAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATc                                                                                 >  1:880733/1‑139 (MQ=255)                                                 aaGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATg                                                                                                            <  2:399864/107‑1 (MQ=255)                                                 aaGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATg                                                                                                            >  1:399864/1‑107 (MQ=255)                                                 aaGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAg                                                                            >  2:813053/1‑139 (MQ=255)                                                  aGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACCTCATCgcccga                                                                           <  1:1194638/139‑6 (MQ=255)                                                   ggCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGAt                                                                          <  2:739000/139‑1 (MQ=255)                                                   ggCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGAt                                                                          <  1:1635738/139‑1 (MQ=255)                                                     cGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAATTATCCGGTTGATCTGAAATTGGTTGTAAAGCAGTAATTTTAATCAGGATGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGa                                                                        <  1:1614326/139‑1 (MQ=255)                                                     cGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTa                                                                                                              >  2:1569614/1‑101 (MQ=255)                                                     cGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTa                                                                                                              <  1:1569614/101‑1 (MQ=255)                                                     cGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGa                                                                        <  2:880733/139‑1 (MQ=255)                                                     cGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGa                                                                        <  1:802653/139‑1 (MQ=255)                                                     cGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGa                                                                        <  1:1489631/139‑1 (MQ=255)                                                     cGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGa                                                                        <  2:1953262/139‑1 (MQ=255)                                                       tCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAAc                                                                      >  2:1033247/1‑139 (MQ=255)                                                          aaaCCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACgtg                                                                   <  1:1033247/139‑1 (MQ=255)                                                          aaaCCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACgtg                                                                   <  2:728256/139‑1 (MQ=255)                                                            aCCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTa                                                                 <  2:698732/139‑1 (MQ=255)                                                               gCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAata                                                              <  1:813053/139‑1 (MQ=255)                                                                cAGATAAAGTGCTGAAAGGCGAGAAAGTTGAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCAt                                                                                  <  2:1303078/118‑1 (MQ=255)                                                                cAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTAAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCAt                                                                                  >  1:1303078/1‑118 (MQ=255)                                                                        gTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCa                                                                                                                     >  2:1880208/1‑75 (MQ=255)                                                                        gTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCa                                                                                                                     <  1:1880208/75‑1 (MQ=255)                                                                             gAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATg                                                                         >  2:134595/1‑114 (MQ=255)                                                                             gAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATg                                                                         <  1:134595/114‑1 (MQ=255)                                                                                aGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCa                                             <  2:629906/139‑1 (MQ=255)                                                                                                     gTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTgg                        >  2:1019936/1‑139 (MQ=255)                                                                                                       aTCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGgcg                      <  1:1019936/139‑1 (MQ=255)                                                                                                          cGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACgcgc                                        >  2:642432/1‑118 (MQ=255)                                                                                                          cGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACgcgc                                        <  1:642432/118‑1 (MQ=255)                                                                                                                aTCTGAAACTGGTTGTTAAGCAGTAATTTTAATCGGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCAttt                 >  2:649950/1‑135 (MQ=255)                                                                                                                aTCTGAAACTGGTTGTTAAGCAGTAATTTTAATCGGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCAttt                 <  1:649950/135‑1 (MQ=255)                                                                                                                aTCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATACATACGTCATCGACAGAt                                                                          <  1:1176223/78‑1 (MQ=255)                                                                                                                aTCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATc                                                                                 <  2:1511440/71‑1 (MQ=255)                                                                                                                aTCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATc                                                                                 >  1:1511440/1‑71 (MQ=255)                                                                                                                aTCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGAt                                                                          >  2:1176223/1‑78 (MQ=255)                                                                                                                aTCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGa                                                                        <  2:1786617/80‑1 (MQ=255)                                                                                                                aTCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGa                                                                        >  1:1786617/1‑80 (MQ=255)                                                                                                                aTCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGa                                                                        <  2:682971/80‑1 (MQ=255)                                                                                                                aTCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGa                                                                        >  1:682971/1‑80 (MQ=255)                                                                                                                          ggTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGa   >  2:1798123/1‑139 (MQ=255)                                                                                                                           gttgttAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACgtg                                                                   >  2:408100/1‑74 (MQ=255)                                                                                                                           gttgttAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACgtg                                                                   <  1:408100/74‑1 (MQ=255)                                                                                                                           gttgttAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCAcgcg                                         >  1:898507/1‑100 (MQ=255)                                                                                                                           gttgttAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCAcgcg                                         <  2:898507/100‑1 (MQ=255)                                                                                                                           gttgttAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGaa  >  2:477178/1‑139 (MQ=255)                                                                                                                           gttgttAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGaa  >  1:1540778/1‑139 (MQ=255)                                                                                                                           gttgttAAGCAGTAATTTTAATCACGTTGTATGACCTGATGGTGACATAAAt                                                                                         >  2:1186950/1‑52 (MQ=255)                                                                                                                           gttgttAAGCAGTAATTTTAATCACGTTGTATGACCTGATGGTGACATAAAt                                                                                         <  1:1186950/52‑1 (MQ=255)                                                                                                                            ttgttAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGgcgc                     <  2:55262/119‑1 (MQ=255)                                                                                                                            ttgttAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGgcgc                     >  1:55262/1‑119 (MQ=255)                                                                                                                                     agtaatTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGgc                       >  1:2039905/1‑108 (MQ=255)                                                                                                                                     agtaatTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGgc                       <  2:2039905/108‑1 (MQ=255)                                                                                                                                         |                                                                                                                             GCAAACTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAA  >  NC_000913/3937032‑3937292

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