Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,164,123	C→G	A162G (GCA→GGA)	btuB →	cobalamin outer membrane transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,164,123	0	C	G	100.0%	135.8 / NA	43	A162G (GCA→GGA)	btuB	cobalamin outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (22/21);  total (22/21)

CAGGTCAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTAC‑GCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAG‑TTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCTTTATGGCGCGCTGGAGCATAACTTTACTGATGCCTGGAGCGGCTTTGTGCGCGGCT  >  NC_000913/4163877‑4164359                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                              cAGGTCAGCTCTCAGCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGGTGGTGTTAAGTGAGGGCGTAC‑GCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTGGAATATACCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAggg                                                                                                                                                                                                                                           <  2:708709/251‑1 (MQ=255)                  aTTTTTATTCGCGGTACAAATGCCAGTCATGTTTTGGTGTTAATTGATGGCGTAC‑GCCTAAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATTAACCCCGAACGGAAATTTCAGGAGGGTGGGGAATCAATa                                                                                                                                                                                                                              <  1:190686/247‑1 (MQ=255)                  aTTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTAC‑GCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATa                                                                                                                                                                                                                              <  1:13458/247‑1 (MQ=255)                       tATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTAC‑GCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAAt                                                                                                                                                                                                                               >  1:160431/1‑241 (MQ=255)                       tATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTAC‑GCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAAt                                                                                                                                                                                                                               <  2:160431/241‑1 (MQ=255)                             cgGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTAC‑GCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAG‑TTATCAGAACTatg                                                                                                                                                                                                              >  1:486112/1‑251 (MQ=255)                             cgGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTAC‑GCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGGGGTGAATATCATTACGGAGCGCCGTGAACCCCGAACGGGAGTTTCAGGAGGGGGGGGAAGGAAAAG‑TTATCAGaatcatg                                                                                                                                                                                                              >  2:36778/1‑246 (MQ=255)                                       gcagttcagtgtttGGTGTTAATTGAGGGC‑TACCGCCTTAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCTCTTGTCCAGCGTGTTGAATATACCCGTGGGTCGCGCTCCGCTGTTTATGGGTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAG‑TTACCAGATCTATGATGTCTCTAc                                                                                                                                                                                                    <  2:210399/239‑1 (MQ=255)                                            tcctgtgTTGGTGTTAATTGATGGCGTAC‑GCCTGACTCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTCTGGTTCCGATGCAATAGGCGGGGGGGAGAATATACTCACGCCGCGCGCGGAACCCGGAACGGGAATTTCCGGAGGGTGGGGACGCAACAGGTTTTC‑GAACTATGATGTCTCTAccaagc                                                                                                                                                                                               >  2:285474/4‑246 (MQ=255)                                               tgtgTTGGTGTTAATTGATGGCGTAC‑GCCTGAATATGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGCTGAACCCGGAACGGAAATTTGAGGAGGGGGGGGAAGCCATAG‑GTATCAGAACTATGAAGTCTCTACGCAGcaac                                                                                                                                                                                            >  1:537072/1‑251 (MQ=255)                                                gtgtTGGTGTTAATTGATGGCGTAC‑GCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCAGCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAG‑TTATCAGAACTATGATGTCTCTACGCAGcaaca                                                                                                                                                                                           <  1:458518/251‑1 (MQ=255)                                                       tgtTAATTGATGGCGTAC‑GCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGTATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAG‑TTATCAGAACTATGATGTCTCTACGCAGCAACAACTgggg                                                                                                                                                                                    >  1:298664/1‑251 (MQ=255)                                                       tgtTAATTGATGGCGTAC‑GCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGGAGCAATAG‑TTATCAGAACTATGATGTCTCTACGCAGCAACAACTgggg                                                                                                                                                                                    >  1:462059/1‑251 (MQ=255)                                                           aaTTGATGGCGTAC‑GCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGGGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAG‑TTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGAt                                                                                                                                                                                 >  1:185665/1‑250 (MQ=255)                                                                aTGGCGTAC‑GCCTGAATCTGGCGGGGGTGAGTGGTTATGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGGGGTGAATATCATCATGAGGCGCGATGAACCCGGAACGGAACATTCAGGAGGGTGGGGAAGCAATAG‑TTATAAGAACTATGATGGCTATACGCAGCAAAGACGGGGGGATAAGac                                                                                                                                                                            >  1:527626/1‑250 (MQ=255)                                                                     gTAC‑GTCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAG‑TTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGt                                                                                                                                                                      <  1:552827/251‑1 (MQ=255)                                                                     gTAC‑GCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGGGGGGAAGCAATAG‑TTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGt                                                                                                                                                                      >  1:134465/1‑251 (MQ=255)                                                                                        gtggTGAGTGTTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTGGGCCAGCGTGTTGAATATATCCGGGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCGCGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAG‑TTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTa                                                                                                                                                    <  2:559341/249‑1 (MQ=255)                                                                                                        ccGACCTTAGCAAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGATCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGGGGTAAATATCATCACGACGCGCGATTAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAG‑TTATCAGAACTATGCTGTCTCTACGCAGCTACAACTGGGGGATAAGACACGGGCAACGCTGTTGGGAGCTTATGCCCATACTCAAGGt                                                                                                                                    >  2:400140/1‑251 (MQ=255)                                                                                                         cgtccTTAGCCAGTTCCCTATTGCGCTGGTCCAGCGTGTTTAATAGATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACTCGGAACGGAAATTTCAGGAGGGTGGGGAAGCTATAG‑TTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGtt                                                                                                                                   <  1:49579/248‑1 (MQ=255)                                                                                                          gACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAG‑TTATCAGAACTATGATGTCTCTTCGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTTTTGGGCGATTATGCCCATACTCATGGttt                                                                                                                                  >  1:565824/1‑250 (MQ=255)                                                                                                                  ccAGTTCCCTATTGCGCGTGTCCAGCGTGTTGTATAAATCCGTGGGCCGCGCGCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAG‑TTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGAtgttg                                                                                                                          <  2:298664/251‑1 (MQ=255)                                                                                                                              tGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAG‑TTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGt                                                                                                               <  1:477102/250‑1 (MQ=255)                                                                                                                                 gcTTGTCCAGCGTGTTGACTATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGAGAATATCATCACGACGCGCGATAAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAAAAG‑TTATCAGGACTCTGATGTCTCTAAGCAGCAACAACTGGGGGGAAAAACACGGGGAACGATGTTGGGGGATTATGCCCATACTCATGGGTATGGTGTTGTTGCCTATGGTaaa                                                                                                            >  2:653501/1‑249 (MQ=255)                                                                                                                                   ttGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGGAATTTCAGGAGGGTGGGGAAGCAATAG‑TTATCAGAACTATGATGTCTATACGCCGCAACAACTGGGGGAAAAGACACGGGTAACGCTGTTGGGGGATTATGCACATACTCATGGTTATGATGTTGTTGCATATGGTacaacc                                                                                                         >  2:343404/1‑246 (MQ=255)                                                                                                                                     gTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAG‑TTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCg                                                                                                        <  1:227051/250‑1 (MQ=255)                                                                                                                                                               gCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAG‑TTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATg                                                                             <  1:468684/251‑1 (MQ=255)                                                                                                                                                                      tCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAG‑TTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGtttttt                                                                      >  1:159310/1‑251 (MQ=255)                                                                                                                                                                            gTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAG‑TTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCCAGCGCAGACAGATAACGATGGTTTTTTAAGTa                                                                 >  1:484802/1‑250 (MQ=255)                                                                                                                                                                                aTGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAG‑TTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAAc                                                             <  2:134465/250‑1 (MQ=255)                                                                                                                                                                                     tCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAG‑TTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTGGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCttt                                                        <  1:285474/250‑1 (MQ=255)                                                                                                                                                                                                    ggggTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAG‑TTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCTTTATGGCGCGCTGGAGct                                        >  1:65305/1‑250 (MQ=255)                                                                                                                                                                                                     gggTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAG‑TTATCAGAACTATGATGTCTCTACGCAGCAACAACTGTGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCTTTATGGCGCGCTGGAGCAt                                       <  1:400140/251‑1 (MQ=255)                                                                                                                                                                                                         ggttaaTATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAG‑TTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCTTTATGGCGCGCTGGAGCATAACt                                   <  2:241054/247‑1 (MQ=255)                                                                                                                                                                                                                      acgacgCGGGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAG‑TTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCTTTATGGCGCGCTGGAGCATAACTTTACTGATGCCTg                      <  2:90425/251‑1 (MQ=255)                                                                                                                                                                                                                        gacgCGCGATGAACCCGGAACGGAATTGTCAGGCGGGTGGGGACGCGATCG‑TTATCAGAACTATGATGTCTCTACGCAGCAACCACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCACATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCTTTATGGCGCGCTGGAGCATAACTTTACTgtgacctgg                     >  2:102191/1‑242 (MQ=255)                                                                                                                                                                                                                                aTGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAG‑TTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGCCAGATAACGATGGTTTTTTAAGTAAAACGCTTTATGGCGCGCTGGAGCATAACTTTACTGATGccgttcggaagagc            >  2:175923/1‑239 (MQ=255)                                                                                                                                                                                                                                aTGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAG‑TTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCTTTATGGCGCGCTGGAGCATAACTTTACTGATGc                         <  1:175923/238‑1 (MQ=255)                                                                                                                                                                                                                                    actcggaacggaaATTTCAGGAGGGTGGGGAAGCAATAG‑TTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCTTTATGGCGCGCTGGAGCATAACTTTACTGATGCCTGGAGCGGCTTTGTgc        <  1:169589/248‑1 (MQ=255)                                                                                                                                                                                                                                     cccGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAG‑TTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCTTTATGGCGCGCTGGAGCATAACTTTACTGATGCCTGGAGCGGCTTTGTgcg       >  1:258062/1‑251 (MQ=255)                                                                                                                                                                                                                                       cggaacggaaATTTCAGGAGGGTGGGGAAGCAATAG‑TTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCTTTATGGCGCGCTgg                                            >  1:182523/1‑212 (MQ=255)                                                                                                                                                                                                                                       cggaacggaaATTTCAGGAGGGTGGGGAAGCAATAG‑TTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCTTTATGGCGCGCTgg                                            <  2:182523/212‑1 (MQ=255)                                                                                                                                                                                                                                          aacggtaATTTCAGGAGGGTGGGGAAGCAATAG‑TTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCTTTATGGCGCGCTGGAGCATAACTTTACTGATGCCTGGAGCGGCTTTGTGCGCggct  >  2:592945/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                              CAGGTCAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTAC‑GCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAG‑TTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCTTTATGGCGCGCTGGAGCATAACTTTACTGATGCCTGGAGCGGCTTTGTGCGCGGCT  >  NC_000913/4163877‑4164359

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