Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,692,114	C→T	S165N (AGC→AAC)	bcsB ←	cellulose synthase periplasmic subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,692,114	0	C	T	100.0%	131.6 / NA	40	S165N (AGC→AAC)	bcsB	cellulose synthase periplasmic subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (22/18);  total (22/18)

CATCCGGCGCACCCGCAAAGACCATCGGCAAGGTGTTGGTGCGGTTATCGCTCGGGTCAAAGAATGGCACCGGGAAGTGTGACAGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGCTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCATGCTGCCCGGCG  >  NC_000913/3691887‑3692360                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       cATCCAGCGAACCCGCAAAGAACATCGGCAAGGTGTTCGTGCGGATATCGCTCGAGTCAAAAAATAGCACCGGGAAGTGTGACAGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGCTCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTAGCCGGTTTTTCGCACACGGCCTGATAATGGCCGACAAACTCCAGCCGTACTCGGTTGAAGTCGTTAATaaaacgagagttactggcc                                                                                                                                                                                                                                 >  2:171212/1‑235 (MQ=255)      ggCGCACCCGCAAAGACCATCGGCAAGGTGTTGGTGCGGTTATCGCTCGGGTCAAAGAATGGCACCGGGAAGTGTGACAGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTg                                                                                                                                                                                                                            <  1:149427/251‑1 (MQ=255)                   aGACCATCGGCAAGGTGTTGGTGCGGTTATCGCTCGGGTCAAAGAATGGCACCGGGAAGTGTGACAGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTCCACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGttttt                                                                                                                                                                                                               >  2:5220/1‑251 (MQ=255)                       cATCGGCAAGGTGTTGGTGCGGTTATCGCTCGGGTCAAAGAATGGCACCGGGAAGTGTGACCGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTAc                                                                                                                                                                                                           >  2:176180/1‑251 (MQ=255)                                     tGGTGCGGTTATCGCTCGGGTCAAAGAATGGCACCGGGAAGTGTGACAGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCttt                                                                                                                                                                                              <  2:21980/250‑1 (MQ=255)                                                   atcggggcaatgcacggCACCGGGAAGTGTGACAGGCCCTGCGTCACAGTCAGGGTCTGATAGGTCAGATCGCGTCCACTGCTCCGCTTAACATCCAGCCAAAGCTCGGTGCTGGCCTGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCAGTTGCTCCATCGTTTTTTTTCCCAACTGTCCTTTGGTCACTGGAAGCAc                                                                                                                                                                               <  2:110699/236‑1 (MQ=255)                                                      gggTCAAAGAATGGCACCGGGAAGTGTGACAGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGc                                                                                                                                                                                                                     >  1:163198/1‑210 (MQ=255)                                                      gggTCAAAGAATGGCACCGGGAAGTGGGACAGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGGTAATGGGCATTTGCGCCAGc                                                                                                                                                                                                                     <  2:163198/210‑1 (MQ=255)                                                        gTCAAAGAATGGCACCGGGAAGTGTGACAGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCa                                                                                                                                                                          <  2:40119/251‑1 (MQ=255)                                                          cAAAGAATGGCACCGGGAAGTGTGACAGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGGCCTGATAATGGCCGACAAACGCCAGCCGTGCACGGTTGAAGGCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGgc                                                                                                                                                                                    <  1:272467/239‑1 (MQ=255)                                                          cAAAGAATGGCACCGGGAAGTGTGACAGGTCATTCTTAACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGgc                                                                                                                                                                                    >  2:272467/1‑239 (MQ=255)                                                          cAAAGAATGGCACCGGGAAGTGTGACAGGTCATTATTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATc                                                                                                                                                                        >  1:242070/1‑251 (MQ=255)                                                           aaaGAATGGCACCGGGAAGTGTGACAGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATc                                                                                                                                                                        >  2:74328/1‑250 (MQ=255)                                                            aaGATTGGCACCGGGAAGTGTGACAGGTCATTCTTCACATCCAGGGCCTGATAGGTCGGACTCAGTCCACCGCTCCGCCCAACACCCAGCCAAAGCGAGGTGCTGGCCGGTTTTTCGCACACATACAGATAATGGACAACAAACTCCAGCCGTACACGAATGAAGTCGTTGATAAGCAGTGGGGTAATGGGCTGTCGCGCCAGCGTTTTTTTTCACAACAGTTCTTTGGTCAGTGGCAGCACGCCCATCaa                                                                                                                                                                      <  2:65148/251‑2 (MQ=255)                                                                      ccGGGAAGTGTGACAGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTg                                                                                                                                                            >  1:45057/1‑251 (MQ=255)                                                                                ggACAGGTCAGTCGTCTCATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCGCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTGTCGCACACGGCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTGAATAAACAGTGGGTTAAGGGGGAGTTGCGCCAGCGTGGTTTTACCCAACTGTTCTTTGGTCTCTGGCAGCACGCCCATCAGTTGATCATTGAGGTAAACCt                                                                                                                                                  <  2:271561/250‑1 (MQ=255)                                                                                           tcttcACATTCAGGGTATGATAGGTAAGTTCAAGTCCACTGCTCAGCCCAACATCCAGACAAAGCGTGGTACTGGCCGGTTTTTGGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGCACACGGTTGAAGACGTTAATAAACAGGGTGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCACATCAGTTCATCATTGAGATAACCCTTTAACTGCGAc                                                                                                                                       >  2:88699/1‑251 (MQ=255)                                                                                              tcACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTTGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGAc                                                                                                                                       <  1:80103/248‑1 (MQ=255)                                                                                                  aTTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGAc                                                                                                                                  >  1:141653/1‑249 (MQ=255)                                                                                                       gggTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATGGAGATAAACCTTTAACTGCGACTGGACAGGCAGt                                                                                                                           <  1:244261/251‑1 (MQ=255)                                                                                                          tCTGATAGGTCAGCTCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATACACCTTTAACTGCGACTGGACAGGCAGTAAc                                                                                                                        >  2:145284/1‑251 (MQ=255)                                                                                                                    cAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAAATGCGACTGGACAGGCAGTAAcgatggcggtg                                                                                                              >  2:1048/1‑248 (MQ=255)                                                                                                                     atatccattacACTGCTCCGCCAACCCTCCAGACAAAGCGTAGAGCCGGCGGGTTTTTAGCACACTTCCCGATAATGGACGCAAAACTCCAGACGTCCACGGTTGACGTCCTTAATAAACAGTGGGTTAATGGGCATTTGCACCAGCGTTTTTTTACCCAACTGTTCATTGGTCACTGGCAGAACGCCCATCAGATCATCATTGAAATACACATTTAACAGCGACTGGACAGGCAGTAAcgactgcgatgg                                                                                                             >  2:264761/11‑251 (MQ=255)                                                                                                                      gATCCCGTCCCCTGCTCCGCCCAACATCCAGCCAAAGCGTGGAGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCAATGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACGGGCCAGGCATTAACGCTGGCGTTggg                                                                                                            >  2:110735/1‑251 (MQ=255)                                                                                                                                  tGCTCCGCCCAACCTCACGCCCAAGCGTGGTCCTGCCCGGTTTTTAGCACACGTCCTGATACTGGCCGACACACTCCAGCCGTACACGGTTGACGTCGTTACTAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACCGTTCTTTGGTCACTGGCAGCACGCCCCTCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGAAGGCGGTGGGTCGTATTCGAgg                                                                                                >  2:134157/1‑251 (MQ=255)                                                                                                                                   gCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCATCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAgg                                                                                                <  1:178528/250‑1 (MQ=255)                                                                                                                                                        aaGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAGCTCCAGCCGTACACGGTTGAAGTCGTTCATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTAACCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAAcgatgg                                                                                                                   >  1:129313/1‑210 (MQ=255)                                                                                                                                                                     gCCGGTGTTTCGCCCACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGGTGAAGGCGTTAATAAACAGTGGGTTAATGGGCATTTGCTCCAGCGTGTTGTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACTAGGGCGATGGGGTGTATGCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTgcg                                                             <  2:122618/251‑1 (MQ=255)                                                                                                                                                                           ttttCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATAcc                                                       >  2:86092/1‑251 (MQ=255)                                                                                                                                                                               cGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCaaa                                                    <  2:9872/250‑1 (MQ=255)                                                                                                                                                                                     cGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTACCGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAAtg                                             >  2:118546/1‑251 (MQ=255)                                                                                                                                                                                     cGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAAtg                                             >  2:29682/1‑251 (MQ=255)                                                                                                                                                                                                  gCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGt                                >  2:43910/1‑251 (MQ=255)                                                                                                                                                                                                                  cGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCAc                >  2:130357/1‑251 (MQ=255)                                                                                                                                                                                                                     acacGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGAGAAACCTTTAACTGCGACTGCACAGGCAGTATCGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCCCTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCGCGTAGCACCa              <  1:124728/250‑1 (MQ=255)                                                                                                                                                                                                                       acGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCAtgc           <  1:19206/251‑1 (MQ=255)                                                                                                                                                                                                                            tGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCATGCTGcc       <  1:86092/250‑1 (MQ=255)                                                                                                                                                                                                                            tGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCATGCTGcc       <  1:111666/250‑1 (MQ=255)                                                                                                                                                                                                                               aGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCGGTTCGGGTTAATGCCACGTAGCACCATGCTGCCcggc   <  1:176180/251‑1 (MQ=255)                                                                                                                                                                                                                                gTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCATGCTGCCcggcg  >  1:194930/1‑251 (MQ=255)                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       CATCCGGCGCACCCGCAAAGACCATCGGCAAGGTGTTGGTGCGGTTATCGCTCGGGTCAAAGAATGGCACCGGGAAGTGTGACAGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGCTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCATGCTGCCCGGCG  >  NC_000913/3691887‑3692360

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