Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,694,840	G→A	H133Y (CAT→TAT)	bcsA ←	cellulose synthase, catalytic subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,694,840	0	G	A	100.0%	273.7 / NA	84	H133Y (CAT→TAT)	bcsA	cellulose synthase, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (38/46);  total (38/46)

GCGCCAGCGGGTTAAACGGCTGAGTAACGCAGATTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATGCTGTGTTTTATGGCTAACGCGCTCAGGCAGCT‑CGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGAAAATGCCCGCCTGCGCGGTTTCGGCGTTTCTTTGCGCGATGCACCGATCA  >  NC_000913/3694711‑3694968                                                                                                                                  |                                                                                                                                  gcgcCAGCGGGTTAAACGGCTGAGTAACGCAGATTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGtttt                                                                                                                          >  1:1747453/1‑139 (MQ=255) gcgcCAGCGGGTTAAACGGCTGAGTAACGCAGATTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGtttt                                                                                                                          >  1:1419377/1‑139 (MQ=255) gcgcCAGCGGGTTAAACGGCTGAGTAACGCAGATTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGCGtttt                                                                                                                          >  2:1806093/1‑139 (MQ=255)   gcCAGCGGGTTAAACGGCTGAGTAACGCAGATTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTAt                                                                                                                        >  2:736229/1‑139 (MQ=255)    ccAGCGGGTTAAACGGCTGAGTAACGCAGATTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATg                                                                                                                       >  2:1076015/1‑139 (MQ=255)    ccAGCGGGTTAAACGGCTGAGTAACGCAGATTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATg                                                                                                                       >  1:248713/1‑139 (MQ=255)    ccAGCGGGTTAAACGGCTGAGTAACGCAGATTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATg                                                                                                                       >  1:535682/1‑139 (MQ=255)    ccAGCGGGTTAAACGGCTGAGTAACGCAGATTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATg                                                                                                                       >  2:769192/1‑139 (MQ=255)       gCGGGTTAAACGGCTGAGTAACGCAGATTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCt                                                                                                                    >  2:1232495/1‑139 (MQ=255)        cGGGTTAAACGGCTGAGTAACGCAGATTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTa                                                                                                                   <  1:769192/139‑1 (MQ=255)         gggTTAAACGGCTGAGTAACGCAGATTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTaa                                                                                                                  <  1:736229/139‑1 (MQ=255)         gggTTAAACGGCTGAGTAACGCAGATTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTaa                                                                                                                  <  1:1586480/139‑1 (MQ=255)           gTTAAACGGCTGAGTAACGCAGATTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAAcg                                                                                                                <  1:1232495/139‑1 (MQ=255)             tAAACGGCTGAGTAACGCAGATTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAAcgcg                                                                                                              <  1:385139/139‑1 (MQ=255)               aaCGGCTGAGTAACGCAGATTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCt                                                                                                            <  1:248296/139‑1 (MQ=255)                     tGAGTAACGCAGATTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCa                                                                                                      >  2:413940/1‑139 (MQ=255)                         tAACGCAGATTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCT‑c                                                                                                 >  2:139736/1‑139 (MQ=255)                         tAACGCAGATTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCT‑c                                                                                                 >  2:1680786/1‑139 (MQ=255)                           aCGCAGATTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCT‑CGt                                                                                               <  1:1035278/139‑1 (MQ=255)                              cAGATTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCT‑CGTTCa                                                                                            <  1:1680786/139‑1 (MQ=255)                               aGATTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCT‑CGTTCAt                                                                                           <  1:1768398/139‑1 (MQ=255)                               aGATTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCT‑CGTTCAt                                                                                           <  1:803660/139‑1 (MQ=255)                               aGATTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCT‑CGTTCAt                                                                                           <  1:1661895/139‑1 (MQ=255)                                gATTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCT‑CGTTCATc                                                                                          <  1:413940/139‑1 (MQ=255)                                 aTTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCT‑CGtt                                                                                              >  1:474076/1‑134 (MQ=255)                                 aTTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCT‑CGtt                                                                                              <  2:474076/134‑1 (MQ=255)                                 aTTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCT‑CGTTCATcc                                                                                         >  1:292404/1‑139 (MQ=255)                                 aTTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCT‑CGTTCATcc                                                                                         >  1:486277/1‑139 (MQ=255)                                   tAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGCTACTGTGTTTTATGGCTAACGCGCTCAGGCAGCT‑CGTTCATCCAt                                                                                       <  1:1806093/139‑1 (MQ=255)                                    aaCGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCT‑CGTTCATCCAtt                                                                                      <  1:1139940/139‑1 (MQ=255)                                     aCGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCT‑CGTTCATCCATTg                                                                                     <  1:53803/139‑1 (MQ=255)                                           gAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCT‑CGTTCATCCATTGATGGTa                                                                               <  1:1076015/139‑1 (MQ=255)                                                      aaGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCT‑CGTTCATCCATTGATGGTAACGTCCACGaa                                                                    >  2:1437287/1‑139 (MQ=255)                                                            aCGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCT‑CGTTCATCCATTGATGGTAACGTCCACGAATAtttt                                                              >  1:86135/1‑139 (MQ=255)                                                             cGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCT‑CGTTCATCCATTGATGGTAACGTCCACGAATATTTTg                                                             <  2:248713/139‑1 (MQ=255)                                                               atgatACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCT‑CGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCa                                                           <  1:139736/139‑1 (MQ=255)                                                                tgatACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGAGCGAGCTCTTTGTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCT‑CGTTCATCCATTGATCGTAACGTCCACGAATATTTTGCAg                                                          <  2:1316012/139‑1 (MQ=255)                                                                 gatACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCT‑CGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGa                                                         <  2:86135/139‑1 (MQ=255)                                                                          gATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGACAGCT‑CGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGAAAAt                                                >  1:1133272/1‑139 (MQ=255)                                                                              aaCCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTaa                                                                                                                  <  1:1094607/70‑1 (MQ=255)                                                                              aaCCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTaa                                                                                                                  >  2:1094607/1‑70 (MQ=255)                                                                               aCCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAAcgcg                                                                                                              <  1:1372535/73‑1 (MQ=255)                                                                               aCCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAAcgcg                                                                                                              >  2:1372535/1‑73 (MQ=255)                                                                               aCCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCT‑CGTTc                                                                                             <  2:398741/89‑1 (MQ=255)                                                                               aCCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCT‑CGTTc                                                                                             >  1:398741/1‑89 (MQ=255)                                                                               aCCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCT‑CGTTCATCCATTGATGGTAACGTCCACGa                                                                     >  1:1552000/1‑113 (MQ=255)                                                                               aCCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCT‑CGTTCATCCATTGATGGTAACGTCCACGa                                                                     <  2:1552000/113‑1 (MQ=255)                                                                               aCCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCT‑CGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGAAAATGccc                                            <  1:432893/138‑1 (MQ=255)                                                                               aCCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCT‑CGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGAAAATGccc                                            >  2:432893/1‑138 (MQ=255)                                                                               aCCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCT‑CGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGAAAATGCCCg                                           >  1:1233011/1‑139 (MQ=255)                                                                                ccgccgCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCT‑CGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGAAAATGccc                                            >  1:381943/1‑137 (MQ=255)                                                                                ccgccgCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCT‑CGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGAAAATGccc                                            <  2:381943/137‑1 (MQ=255)                                                                                   ccgcGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCT‑CGTTCATCCAt                                                                                       >  1:1685868/1‑91 (MQ=255)                                                                                   ccgcGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCT‑CGTTCATCCAt                                                                                       <  2:1685868/91‑1 (MQ=255)                                                                                     gcgcCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCT‑CGTTCATCCATTGATGGTAACg                                                                            <  1:218196/100‑1 (MQ=255)                                                                                     gcgcCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCT‑CGTTCATCCATTGATGGTAACg                                                                            >  2:218196/1‑100 (MQ=255)                                                                                     gcgcCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCT‑CGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGAAAATGCCCGCCTgcg                                     >  2:960770/1‑139 (MQ=255)                                                                                     gcgcCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCT‑CGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGAAAATGCCCGCCTgcg                                     >  2:63307/1‑139 (MQ=255)                                                                                      cgcCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAg                                                                                                     >  2:687884/1‑75 (MQ=255)                                                                                      cgcCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAg                                                                                                     <  1:687884/75‑1 (MQ=255)                                                                                       gcCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCT‑CGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGa                                                    <  2:1480330/122‑1 (MQ=255)                                                                                       gcCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCT‑CGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGa                                                    >  1:1480330/1‑122 (MQ=255)                                                                                             gCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCT‑CGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGAAAATGCCCGCCTGCGCGGttt                               <  1:1759629/137‑1 (MQ=255)                                                                                             gCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCT‑CGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGAAAATGCCCGCCTGCGCGGttt                               >  2:1759629/1‑137 (MQ=255)                                                                                                 tCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCTCCGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGAAAATGCCCGCCTGCGCGGTTTCGGCg                          <  2:1747453/139‑1 (MQ=255)                                                                                                       gACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCT‑CGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGAAAATGCCCGCCTGCGCGGTTTCGGCGTTTCttt                   <  1:1718540/139‑1 (MQ=255)                                                                                                        aCCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCT‑CGTTCAt                                                                                           <  2:434630/66‑1 (MQ=255)                                                                                                        aCCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCT‑CGTTCAt                                                                                           >  1:434630/1‑66 (MQ=255)                                                                                                          cGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCT‑CGTTCATCCATTGATGGTAACGTGCACGAATATTTTGCAGACCTGAAAATGCCCGCCTGCGCGGTTTCGGCGTTTCTTTgcg                <  2:292404/139‑1 (MQ=255)                                                                                                          cGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCT‑CGTTCATCCATTGATGGTAACGTCCACGaa                                                                    <  2:1009441/87‑1 (MQ=255)                                                                                                          cGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCT‑CGTTCATCCATTGATGGTAACGTCCACGaa                                                                    >  1:1009441/1‑87 (MQ=255)                                                                                                          cGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCT‑CGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCa                                                           <  2:214047/96‑1 (MQ=255)                                                                                                          cGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCT‑CGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCa                                                           >  1:214047/1‑96 (MQ=255)                                                                                                           gAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCT‑CGTTCATCCAt                                                                                       <  2:1252768/67‑1 (MQ=255)                                                                                                           gAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCT‑CGTTCATCCAt                                                                                       >  1:1252768/1‑67 (MQ=255)                                                                                                                   ttttCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCT‑CGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGAAAATGCCCGCCTGCGCGGTTTCGGCGTTTCTTTGCGCGATGCAcc       <  1:960770/139‑1 (MQ=255)                                                                                                                   ttttCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCT‑CGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGAAAATGCCCGCCTGCGCGGTTTCGGCGTTTCTTTGCGCGATGCAcc       <  2:535682/139‑1 (MQ=255)                                                                                                                        tcATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCT‑CGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGAAAATGCCCGCCTGCGCGGTTTCGGCGTTTCTTTGCGCGATGCACCGATCa  <  2:1320468/139‑1 (MQ=255)                                                                                                                        tcATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGACAGCT‑CGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGAAAATGCCCGCCTGCGCGGTTTCGGCGTTTCTTTGCGCGATGCACCGATCa  <  2:1133272/139‑1 (MQ=255)                                                                                                                        tcATCCAGATACTGTGTTGTATGGCTAACGCGCTCAGGCAGCT‑CGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGAAAATGCCCGCCTGCGCGGTTTCGGCGTTTCTTTGCGCGATGCACCGATCa  <  2:1233011/139‑1 (MQ=255)                                                                                                                          aTCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCT‑CGTTCATCCAt                                                                                       <  1:451174/52‑1 (MQ=255)                                                                                                                          aTCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCT‑CGTTCATCCAt                                                                                       >  2:451174/1‑52 (MQ=255)                                                                                                                          aTCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCT‑CGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGAAAATGCCCGCCTGCGCGGTTTCGGCGTTTCt                     <  2:1068185/118‑1 (MQ=255)                                                                                                                          aTCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCT‑CGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGAAAATGCCCGCCTGCGCGGTTTCGGCGTTTCt                     >  1:1068185/1‑118 (MQ=255)                                                                                                                                  |                                                                                                                                  GCGCCAGCGGGTTAAACGGCTGAGTAACGCAGATTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATGCTGTGTTTTATGGCTAACGCGCTCAGGCAGCT‑CGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGAAAATGCCCGCCTGCGCGGTTTCGGCGTTTCTTTGCGCGATGCACCGATCA  >  NC_000913/3694711‑3694968

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