Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_bme_pgi	4,122,984	C→T	R119H (CGT→CAT)	ftsN ←	essential cell division protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_bme_pgi	4,122,984	0	C	T	100.0%	200.5 / NA	63	R119H (CGT→CAT)	ftsN	essential cell division protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (26/37);  total (26/37)

CTGTTCCGCTAATTGCTGCGCCTGACGTTGGCGCTGTAGCGTTTGCTGACGCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGACGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTGGTGGTAGTCCGTTTCCGGT  >  NC_000913_3_bme_pgi/4122845‑4123129                                                                                                                                            |                                                                                                                                                  cTGTTCCCCTCATTGCTGCACCTGACGTTGGCCCTGTAGCGTTTGCTGACGCTATTCCGCCGTCTGTTCATTCCACGCAACTTCAACCACCTGCATTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTg                                                                                                                                          >  2:232765/1‑149 (MQ=255)   gTTCCGCTAATTGCTGCGCCTGACGTTGGCGCTGTAGCGTTTGCTGACGCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGgt                                                                                                                                        <  1:952604/149‑1 (MQ=255)   gTTCCGCTAATTGCTGCGCCTGACGTTGGCGCTGTAGCGTTTGCTGACGCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGgt                                                                                                                                        <  2:197924/149‑1 (MQ=255)        gCTAATTTCTGCGCCGGAGGTTGGGGCTGTAGCGTTTGCTCACCCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCGTTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAg                                                                                                                                   <  2:216909/149‑1 (MQ=255)        gCTAATTGCTGCGCCTGACGTTGGCGCTGTAGCGTTTGCTGACGCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAg                                                                                                                                   <  1:129974/149‑1 (MQ=255)                ctgcGCCTGACGTTGGCGCTGTAGCGTTTGCTGAGGCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCCGTTGCTcct                                                                                                                           >  2:64589/1‑148 (MQ=255)                  gcgcCTGACGTTGGCGCTGTAGCGTTTGCTGACGCTGTTCCGGCGTCTGTTCATTCCACGGCACTTCAACCAGCTGCGTTGTCCGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGc                                                                                                                         >  1:237230/1‑149 (MQ=255)                    gcCTGACGTTGGCGCTGTAGCGGTTGCTGACGCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGt                                                                                                                       <  1:490080/149‑1 (MQ=255)                       tGACGTTGGCGCTGTAGCGTTTGCTGACGCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGtttt                                                                                                                    >  2:125667/1‑149 (MQ=255)                       tGACGTTGGCGCTGTAGCGTTTGCTGACGCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTCGTTCTGGTGTCAGTTGCTCCGGCGtttt                                                                                                                    >  2:64550/1‑149 (MQ=255)                         aCGTTGGCGCTGTAGCGTTTGCTGACGCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTttca                                                                                                                  >  1:767509/1‑149 (MQ=255)                           gTTTGCGCTGTAGCGTTTGCTGACGCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCCGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTTGTGTCAGTTGCTCCGGCGTTttcact                                                                                                                <  1:771793/149‑1 (MQ=255)                            ttGGGGGTGTAGCGTTTGCTGACGTTGTTCGGGGGTCTGTTCATTTCACGGAACTTCAACGAGCGGCGTTGGCTGTTCTCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCACGGGCTCCGGCGATttctact                                                                                                               <  2:809973/149‑5 (MQ=255)                             tGGCGCTGTAGCGTTTGCTGACGCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTttcacttc                                                                                                              <  1:285122/149‑1 (MQ=255)                               gcgcTGTAGCGTTTGCTGACGCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTTATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTttcacttcac                                                                                                            <  1:151790/149‑1 (MQ=255)                               gcgcTGTAGCGGTTGCTGACGCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTACTGGCGCATATCAGCCTGCATTTGTTCCAGCTGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTttcacttcac                                                                                                            <  2:5599/149‑1 (MQ=255)                                  cTGTAGCGTTTGCTGACGCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTcaccac                                                                                                         >  2:465731/1‑149 (MQ=255)                                    gTAGCGGTTGCGGACGCTGGTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGGTGCTGGCGCAGATCAGCCTGCATTTGTTCAAGAAGCCGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCg                                                                                                       <  2:197965/149‑1 (MQ=255)                                         gTTTGCTGACGCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGa                                                                                                  >  1:407944/1‑149 (MQ=255)                                                gACGCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCt                                                                                           <  2:947113/149‑1 (MQ=255)                                                gACGCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCt                                                                                           <  1:172435/149‑1 (MQ=255)                                                  cGCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCtgt                                                                                         <  1:316463/149‑1 (MQ=255)                                                   gCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTGTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCtgtg                                                                                        >  1:533751/1‑149 (MQ=255)                                                    gtGTTCCGGCGTCTGTTCATTCCCCGGAACTTCAACCAGCTGCGTTGGCTGCTAGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTgg                                                                                       <  2:766903/148‑1 (MQ=255)                                                    cTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGACTCCATTTGTGCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAGGTTTCTGTgg                                                                                       >  2:416432/1‑149 (MQ=255)                                                         ccGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGcgca                                                                                  >  2:219466/1‑149 (MQ=255)                                                         ccGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGcgca                                                                                  >  1:842854/1‑149 (MQ=255)                                                                 tGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTcc                                                                          <  2:237230/149‑1 (MQ=255)                                                                  gTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGGTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTccc                                                                         >  2:812278/1‑149 (MQ=255)                                                                      aTTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCt                                                                     <  1:373071/149‑1 (MQ=255)                                                                       ttCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTg                                                                    <  1:235760/149‑1 (MQ=255)                                                                         ccACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTCGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGc                                                                  <  1:64550/149‑1 (MQ=255)                                                                           actgAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGa                                                                <  1:424512/146‑1 (MQ=255)                                                                                     aaCCAGCTGCGTTGGCTGCTGGCGAATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGct                                                      >  1:654433/1‑149 (MQ=255)                                                                                      aCCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGctc                                                     >  1:198792/1‑149 (MQ=255)                                                                                       ccAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGctct                                                    <  1:315279/149‑1 (MQ=255)                                                                                       ccAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGctct                                                    >  1:303083/1‑149 (MQ=255)                                                                                             gCGTTGGCTGCTGGCGCATATCGGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATg                                              <  1:982888/149‑1 (MQ=255)                                                                                               gTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTTTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTa                                            >  1:450180/1‑149 (MQ=255)                                                                                                  ggCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGAGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTTgcg                                         >  2:419806/1‑149 (MQ=255)                                                                                                    ctgctgTCGCATATCCGGCTGCATTTGTAAAAGAATCTGGTGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGACGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGcc                                       <  2:381677/149‑1 (MQ=255)                                                                                                     tgctgGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCa                                      >  2:649680/1‑149 (MQ=255)                                                                                                       ctgGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGc                                    <  1:381247/149‑1 (MQ=255)                                                                                                       ctgGCGCATATCAGCCGGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGc                                    <  1:823802/149‑1 (MQ=255)                                                                                                          gcgcATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGtt                                 >  1:124704/1‑149 (MQ=255)                                                                                                            gcATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGttct                               >  1:455977/1‑149 (MQ=255)                                                                                                             cATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGttctt                              <  2:842854/149‑1 (MQ=255)                                                                                                             cATATCAGCCGGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCTGCGttctg                              <  2:10742/149‑2 (MQ=255)                                                                                                               tatCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTtcttct                            <  2:366424/149‑1 (MQ=255)                                                                                                                atCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTg                           <  1:465731/149‑1 (MQ=255)                                                                                                                atCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTg                           >  2:726299/1‑149 (MQ=255)                                                                                                                 tCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTgg                          <  1:939858/149‑1 (MQ=255)                                                                                                                  cAGCCTGCATTTGTTCAGGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGt                         <  2:654433/149‑1 (MQ=255)                                                                                                                      cTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTtg                     <  1:705018/149‑1 (MQ=255)                                                                                                                       tGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTtgg                    >  1:367541/1‑149 (MQ=255)                                                                                                                        gCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTtggt                   <  1:649680/149‑1 (MQ=255)                                                                                                                         taTTTGTTGAAGAAGCGGATGTTGTTCTGGTGTCAGGTGCTCCGGCGTTTTCACGTCACCACCGGCAGAAGGGTCTGTGGGCGCACGCGCTCCCGGCGGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTtggtg                  <  2:474721/148‑1 (MQ=255)                                                                                                                          aTTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTtggtgg                 >  2:331945/1‑149 (MQ=255)                                                                                                                          aTTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTtggtgg                 <  1:219466/149‑1 (MQ=255)                                                                                                                          aTTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTtggtgg                 <  2:997037/149‑1 (MQ=255)                                                                                                                               ttCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTTGTTTTGGTGGTAGTc            <  2:373349/149‑1 (MQ=255)                                                                                                                               ttCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGAACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTGTTCTGGTTTTGGTGGTAGTc            >  1:742012/1‑149 (MQ=255)                                                                                                                                         tGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTGGTGGTAGTCCGTTTCCGGt  >  1:727965/1‑149 (MQ=255)                                                                                                                                            |                                                                                                                                                  CTGTTCCGCTAATTGCTGCGCCTGACGTTGGCGCTGTAGCGTTTGCTGACGCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGACGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTGGTGGTAGTCCGTTTCCGGT  >  NC_000913_3_bme_pgi/4122845‑4123129

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