Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,246,914	T→C	intergenic (‑145/+149)	cirA ← / ← lysP	colicin IA outer membrane receptor and translocator; ferric iron‑catecholate transporter/lysine transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,246,914	0	T	C	100.0%	222.5 / NA	63	intergenic (‑145/+149)	cirA/lysP	colicin IA outer membrane receptor and translocator; ferric iron‑catecholate transporter/lysine transporter
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (31/32);  total (31/32)

TCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGA  >  NC_000913/2246770‑2247059                                                                                                                                                 |                                                                                                                                                  tCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCt                                                                                                                                             >  2:375053/1‑151 (MQ=255)  ccATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCGCGTCCCCGTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCtt                                                                                                                                            <  1:452394/151‑1 (MQ=255)   cATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCttt                                                                                                                                           >  2:268930/1‑151 (MQ=255)    atatCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCtttt                                                                                                                                          <  2:80078/151‑1 (MQ=255)         tcCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGCGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGAt                                                                                                                                     >  2:319699/1‑151 (MQ=255)            aTGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTgc                                                                                                                                  >  2:775890/1‑151 (MQ=255)              gAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTgcg                                                                                                                                 <  2:234109/150‑1 (MQ=255)                   aaCTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTa                                                                                                                            <  1:58890/150‑1 (MQ=255)                   aaCTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTa                                                                                                                            <  1:410549/150‑1 (MQ=255)                        cGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGAt                                                                                                                       >  2:481433/1‑150 (MQ=255)                        cGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGAt                                                                                                                       >  2:39129/1‑150 (MQ=255)                          aaaataaaatGGGTTATCGCTAACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAAt                                                                                                                    >  1:578899/1‑151 (MQ=255)                           aaataaaatGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATg                                                                                                                   <  2:447681/151‑1 (MQ=255)                           aaataaaatGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATg                                                                                                                   <  1:800158/151‑1 (MQ=255)                              taaaatGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCa                                                                                                                 >  1:686644/1‑150 (MQ=255)                               aaaatGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCaa                                                                                                                >  2:419475/1‑150 (MQ=255)                                 aatGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACg                                                                                                             <  2:227843/151‑1 (MQ=255)                                 aatGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACg                                                                                                             <  1:264038/151‑1 (MQ=255)                                 aatGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACg                                                                                                             <  2:7869/151‑1 (MQ=255)                                      gTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAAc                                                                                                        <  1:481397/151‑1 (MQ=255)                                        tATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACa                                                                                                       >  1:507610/1‑150 (MQ=255)                                        tATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACa                                                                                                       <  2:643202/150‑1 (MQ=255)                                         aTCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAAt                                                                                                     >  1:589137/1‑151 (MQ=255)                                          tCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAAtt                                                                                                    >  1:164936/1‑151 (MQ=255)                                           cGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAAtt                                                                                                    >  2:576411/1‑150 (MQ=255)                                           cGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAAtt                                                                                                    <  1:73868/150‑1 (MQ=255)                                               cacaTCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATc                                                                                                 <  1:514538/149‑1 (MQ=255)                                                acaTCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAAt                                                                                              <  2:310472/151‑1 (MQ=255)                                                 caTCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGTGGCGTAACGATAATGCAAACGATAACAATTATCAATc                                                                                             >  2:355677/1‑151 (MQ=255)                                                 caTCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAAt                                                                                              <  2:326512/150‑1 (MQ=255)                                                 caTCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATc                                                                                             >  2:8179/1‑151 (MQ=255)                                                 caTCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATc                                                                                             >  2:121459/1‑151 (MQ=255)                                                     ttcttcACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCaa                                                                                          >  2:613125/1‑150 (MQ=255)                                                     ttcttcACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCaa                                                                                          >  2:254767/1‑150 (MQ=255)                                                      tcttcACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCa                                                                                           <  1:806010/148‑1 (MQ=255)                                                        ttcACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGt                                                                                      >  1:283331/1‑151 (MQ=255)                                                         tcACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATg                                                                                       <  2:132138/149‑1 (MQ=255)                                                             gTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAAt                                                                                 >  2:752011/1‑151 (MQ=255)                                                               ccccTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAAttt                                                                               >  1:298862/1‑151 (MQ=255)                                                                cccTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAAtttt                                                                              <  1:375028/151‑1 (MQ=255)                                                                 ccTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGACGAAACAGCCTCTGATTTGCCAGCCTTTTTTGAGTGCGGCGTAACGATAAGGCAAGCGATAACAATTAGCAATCCAAATGTTAAAtttt                                                                              <  1:65414/150‑1 (MQ=255)                                                                      cGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATcc                                                                        <  2:49266/151‑1 (MQ=255)                                                                          cGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCAt                                                                    <  1:771227/151‑1 (MQ=255)                                                                                caTCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTg                                                              >  2:277263/1‑151 (MQ=255)                                                                                 aTCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGt                                                             <  1:609060/151‑1 (MQ=255)                                                                                   cgcgGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTc                                                            <  1:158573/150‑1 (MQ=255)                                                                                   cgcgGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTc                                                            >  2:149525/1‑150 (MQ=255)                                                                                       gTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGa                                                       <  1:218659/151‑1 (MQ=255)                                                                                        tGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGa                                                       >  2:356436/1‑150 (MQ=255)                                                                                          gACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAAt                                                    <  1:854252/151‑1 (MQ=255)                                                                                            cTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATgg                                                  >  1:20283/1‑151 (MQ=255)                                                                                                 tttttATGCAGGTGATCATCCCAGAAATCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATTGAAtt                                              <  1:813955/150‑1 (MQ=255)                                                                                                 tttttATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTg                                             >  1:18504/1‑151 (MQ=255)                                                                                                    ttATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGaa                                           <  2:406882/150‑1 (MQ=255)                                                                                                    ttATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGaa                                           >  2:506245/1‑150 (MQ=255)                                                                                                      aTGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGaaaa                                         >  2:611583/1‑150 (MQ=255)                                                                                                                    cccAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGTGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTga                          <  1:180657/151‑1 (MQ=255)                                                                                                                           aGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGGAACGATAATGCAAACGATAAGAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGCGAGGGCt                   >  1:487001/1‑151 (MQ=255)                                                                                                                            gCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCtt                  <  1:226400/151‑1 (MQ=255)                                                                                                                                 aaCAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCa              >  2:658433/1‑150 (MQ=255)                                                                                                                                 aaCAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATACAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCaa             <  1:172502/151‑1 (MQ=255)                                                                                                                                   cAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCaagg           <  2:182672/151‑1 (MQ=255)                                                                                                                                             tttGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGggagga  >  1:715631/1‑150 (MQ=255)                                                                                                                                                 |                                                                                                                                                  TCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGA  >  NC_000913/2246770‑2247059

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