Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	669,398	Δ1 bp	coding (510/612 nt)	cobC ←	putative adenosylcobalamin phosphatase/alpha‑ribazole phosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	669,395	0	A	.	100.0%	277.8 / NA	66	coding (513/612 nt)	cobC	putative adenosylcobalamin phosphatase/alpha‑ribazole phosphatase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (34/32);  total (34/32)

CAATGGTCGTTGCGGCTTACCCGTAAAAAAACAGAAAAGTCATGCATTCTCGACCCCGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCGAAAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTAAAATATTCTGATAGTGCTGAAATTCACTAAGCCTTGCGATAAAGCGTTCCA  >  NC_000913/669257‑669525                                                                                                                                           |                                                                                                                                   cAATGGTCGTTGCGGCTTACCCGTAAAAAAACAGAAAAGTCATGCATTCTCGACCCCGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCGaa                                                                                                                                   <  1:607671/140‑1 (MQ=255)    tGGTCGTTGCGGCTTACCCGTAAAAAAACAGAAAAGTCATGCATTCTCGACCCCGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCGaaatg                                                                                                                                <  2:1122419/140‑3 (MQ=255)     ggTCGTTGCGGCTTACCCGTAAAAAAACAGAAAAGTCATGCATTCTCGACCCCGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCGaaatgc                                                                                                                               >  1:1254782/1‑137 (MQ=255)          ttGCGGCTTACCCGTAAAAAAACAGAAAAGTCATGCATTCTCGACCCCGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCG‑AAATGCCACAt                                                                                                                         >  2:135986/1‑140 (MQ=255)          ttGCGGCTTACCCGTAAAAAAACAGAAAAGTCATGCATTCTCGACCCCGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCG‑AAATGCCACAt                                                                                                                         >  1:748735/1‑140 (MQ=255)          ttGCGGCTTACCCGTAAAAAAACAGAAAAGTCATGCATTCTCGACCCCGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCG‑AAATGCCACAt                                                                                                                         >  2:616224/1‑140 (MQ=255)            gCGGCTTACCCGTAAAAAAACAGAAAAGTCATGCATTCTCGACCCCGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCGaaa                                                                                                                                  >  2:1314134/1‑130 (MQ=255)            gCGGCTTACCCGTAAAAAAACAGAAAAGTCATGCATTCTCGACCCCGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCGaaa                                                                                                                                  <  1:1314134/130‑1 (MQ=255)                cTTACCCGTAAAAAAACAGAAAAGTCATGCATTCTCGACCCCGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCG‑AAATGCCACATAGCTTc                                                                                                                   >  2:1216665/1‑140 (MQ=255)                     ccGTAAAAAAACAGAAAAGTCATGCATTCTCGACCCCGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCG‑AAATGCCACATAGCTTCggcag                                                                                                              >  1:123340/1‑140 (MQ=255)                         aaaaaaaCAGAAAAGTCATGCATTCTCGACCCCGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCG‑AAATGCCACATAGCTTCGGCAGGCAt                                                                                                          >  2:764616/1‑140 (MQ=255)                               aCAGAAAAGTCATGCATTCTCGACCCCGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCG‑AAATGCCACATAGCTTCGGCAGGCATGCCAAt                                                                                                    >  1:549325/1‑140 (MQ=255)                               aCAGAAAAGTCATGCATTCTCGACCCCGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCG‑AAATGCCACATAGCTTCGGCAGGCATGCCAAt                                                                                                    >  1:1357366/1‑140 (MQ=255)                                      aGTCATGCATTCTCGACCCCGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCG‑AAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACgg                                                                                             <  1:135986/140‑1 (MQ=255)                                          aTGCATTCTCGACCCCGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCG‑AAATGCCACATAGCTTCGGCAGGCATGCCAATTAAAc                                                                                               >  1:1237935/1‑134 (MQ=255)                                          aTGCATTCTCGACCCCGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCG‑AAATGCCACATAGCTTCGGCAGGCATGCCAATTAAAc                                                                                               <  2:1237935/134‑1 (MQ=255)                                             cATTCTCGACCCCGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCG‑AAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTa                                                                                      >  1:892388/1‑140 (MQ=255)                                             cATTCTCGACCCCGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCG‑AAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTa                                                                                      <  2:1357366/140‑1 (MQ=255)                                               ttCTCGACCCCGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCG‑AAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAAc                                                                                    <  2:570990/140‑1 (MQ=255)                                               ttCTCGACCCCGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCG‑AAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAAc                                                                                    <  1:1216665/140‑1 (MQ=255)                                               ttCTCGACCCCGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCG‑AAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAAc                                                                                    <  1:616224/140‑1 (MQ=255)                                                tctcGACCCCGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCG‑AAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACa                                                                                   >  2:602380/1‑140 (MQ=255)                                                 ctcGACCCCGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCG‑AAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAg                                                                                  >  1:575847/1‑140 (MQ=255)                                                  tcGACCCCGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCG‑AAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGa                                                                                 <  1:764616/140‑1 (MQ=255)                                                  tcGACCCCGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCG‑AAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGa                                                                                 >  1:589505/1‑140 (MQ=255)                                                   cGACCCCGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCG‑AAATGCCACATAGCTTCGGCAGGCATGCCATTTAAACGGGCGATTAACAGAc                                                                                <  2:1427886/140‑1 (MQ=255)                                                   cGACCCCGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCG‑AAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGAc                                                                                <  2:1332630/140‑1 (MQ=255)                                                   cGACCCCGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCG‑AAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGAc                                                                                <  1:221277/140‑1 (MQ=255)                                                      ccccGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCG‑AAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCa                                                                             >  1:723218/1‑140 (MQ=255)                                                         cGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCG‑AGATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGAc                                                                                >  1:1443937/1‑134 (MQ=255)                                                         cGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCG‑AAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGAc                                                                                <  2:1443937/134‑1 (MQ=255)                                                         cGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCG‑AAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAGTa                                                                          <  2:1182152/140‑1 (MQ=255)                                                         cGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCG‑AAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAGTa                                                                          <  1:602380/140‑1 (MQ=255)                                                         cGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCG‑AAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAGTa                                                                          <  2:748735/140‑1 (MQ=255)                                                          gATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCG‑AAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAGTac                                                                         <  2:123340/140‑1 (MQ=255)                                                          gATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCG‑AAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAGTac                                                                         <  2:777713/140‑1 (MQ=255)                                                          gATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCG‑AAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAGTac                                                                         >  2:896409/1‑140 (MQ=255)                                                                 cGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCG‑AAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAGTACACCCTGa                                                                  >  2:1006006/1‑140 (MQ=255)                                                                      aTTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCG‑AAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCt                                                             <  1:896409/140‑1 (MQ=255)                                                                        tGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCG‑AAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTga                                                           <  2:448834/140‑1 (MQ=255)                                                                        tGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCG‑AAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTga                                                           >  1:1352537/1‑140 (MQ=255)                                                                           ggACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCG‑AAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTgacga                                                        <  2:1254782/140‑1 (MQ=255)                                                                           ggACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCG‑AAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTgacga                                                        <  2:1352537/140‑1 (MQ=255)                                                                           ggACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCG‑AAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTgacga                                                        <  2:780361/140‑1 (MQ=255)                                                                            gACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCG‑AAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTgacgac                                                       >  1:13110/1‑140 (MQ=255)                                                                              cgcgTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCG‑AAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTa                                                     >  2:1436447/1‑140 (MQ=255)                                                                              cgcgTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCG‑AAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTa                                                     >  2:1153323/1‑140 (MQ=255)                                                                               gcgTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCG‑AAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTaa                                                    <  1:1436447/140‑1 (MQ=255)                                                                                 gTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCG‑AAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTaaaa                                                  >  2:139211/1‑140 (MQ=255)                                                                                 gTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCG‑AAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTaaaa                                                  <  2:892388/140‑1 (MQ=255)                                                                                  tAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCG‑AAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTgac                                                          <  1:80137/131‑1 (MQ=255)                                                                                  tAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCG‑AAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTgac                                                          >  2:80137/1‑131 (MQ=255)                                                                                      gTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCG‑AAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTg                                                            <  2:740559/125‑1 (MQ=255)                                                                                      gTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCG‑AAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTg                                                            >  1:740559/1‑125 (MQ=255)                                                                                            aaTTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCG‑AAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTAAAATATTCTGATAg                                       <  1:1153323/140‑1 (MQ=255)                                                                                                    gTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCG‑AAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTAAAATATTCTGATAGTGCTGaaa                               <  2:498516/140‑1 (MQ=255)                                                                                                         tatCAATGGCGCTCCAACCCCCCTGGTCGACGCG‑AAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTAAAATATTCTGa                                          <  2:225514/124‑1 (MQ=255)                                                                                                         tatCAATGGCGCTCCAACACCCCTGGTCGACGCG‑AAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTAAAATATTCTGa                                          >  1:225514/1‑124 (MQ=255)                                                                                                          atCAATGGCGCTCCAACACCCCTGGTCAACGCG‑AAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAGTACGCCCTGATGGCTGACGACTAAAATATTCTGATAGTGCTGAAATTCACt                         <  2:656582/140‑1 (MQ=255)                                                                                                               tGGCGCTCCAACACCCCTGGTCAACGCG‑AAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTAAAATATTCTGATAGTGCTGAAATTCACTAAGcc                    <  2:549325/140‑1 (MQ=255)                                                                                                                  cgcTCCAACACCCCTGGTCAACGCG‑AAATGCCACATAGCTTCGGCAGGCATGTCAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTAAAATATTCTGATAGTGCTGAAATTCACTAAGCCTTg                 >  2:296931/1‑140 (MQ=255)                                                                                                                  cgcTCCAACACCCCTGGTCAACGCG‑AAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTAAAATATTCTGATAGTGCTGAAATTCACTAAGCCTTg                 >  2:1295962/1‑140 (MQ=255)                                                                                                                    cTCCAACACCCCTGGTCAACGCG‑AAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTAAAATATTCTGATAGTGCTGAAATTCACTAAGCCTTGCg               >  2:782739/1‑140 (MQ=255)                                                                                                                    cTCCAACACCCCTGGTCAACGCG‑AAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTAAAATATTCTGATAGTGCTGAAATTCACTAAGCCTTGCg               >  2:1427396/1‑140 (MQ=255)                                                                                                                       cAACACCCCTGGTCAACGCG‑AAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTAAAATATTCTGATAGTGCTGAAATTCACTAAGCCTTGCGATa            >  1:1115721/1‑140 (MQ=255)                                                                                                                          cacCCCTGGTCAACGCG‑AAATGCCACATAGCTTCGGCAGGCATGTCAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTAAAATATTCTGATAGTGCTGAAATTCACTAAGCCTTGCGATAAAg         <  1:296931/140‑1 (MQ=255)                                                                                                                          cacCCCTGGTCAACGCG‑AAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTAAAATATTCTGATAGTGCTGAAATTCTCTAAGCCTTGCGATAAAg         >  2:20181/1‑140 (MQ=255)                                                                                                                          cacCCCTGGTCAACGCG‑AAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTAAAATATTCTGATAGTGCTGAAATTCACTAAGCCTTGCGATAAAg         <  1:1006006/140‑1 (MQ=255)                                                                                                                              ccTGGTCAACGCG‑AAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTAAAATATTCTGATAGTGCTGAAATTCACTAAGCCTTGCGATAAAGCGtt     >  1:440974/1‑140 (MQ=255)                                                                                                                               cTGGTCAACGCG‑AAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTAAAATATTCTGATAGTGCTGAAATTCACTAAGCCTTGCGATAAAGCGTTc    >  1:1261952/1‑140 (MQ=255)                                                                                                                                 ggTCAACGCG‑AAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTAAAATATTCTGATAGTGCTGAAATTCACTAAGCCTTGCGATAAAGCGTTCca  >  1:47489/1‑140 (MQ=255)                                                                                                                                           |                                                                                                                                   CAATGGTCGTTGCGGCTTACCCGTAAAAAAACAGAAAAGTCATGCATTCTCGACCCCGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCGAAAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTAAAATATTCTGATAGTGCTGAAATTCACTAAGCCTTGCGATAAAGCGTTCCA  >  NC_000913/669257‑669525

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