Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,688,510	A→G	Y222H (TAC→CAC)	bcsC ←	cellulose biosynthesis protein BcsC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,688,510	0	A	G	100.0%	173.5 / NA	48	Y222H (TAC→CAC)	bcsC	cellulose biosynthesis protein BcsC
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (25/23);  total (25/23)

CAGCTGTTTTTGCTGTTCTGCCAGTTGCGATTGCGCAGCCGCCACGCTATCGCCATCACTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGCTACTAAACAGCAATAGCGCCAGATTGTTTTGCAGGCCCGTATTGCCCGGTGC  >  NC_000913/3688379‑3688645                                                                                                                                    |                                                                                                                                        cAGCTGTTTTTGCTGTTCTGCCAGTTGCGATTGCGCAGCCGCCACGCTATCGCCATCACTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTGCCAGAtt                                                                                                                                  >  2:298564/1‑139 (MQ=255) cAGCTGTTTTTGCTGTTCTGCCAGTTGCGATTGCGCAGCCGCCACGCTATCGCCATCACTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTGCCAGAtt                                                                                                                                  <  1:298564/139‑1 (MQ=255)    cTGTTTTTGCTGTTCTGCCAGTTGCGATTGCGCAGCCGCCACGCTATCGCCATCACTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTGCCAGATTTTag                                                                                                                              >  2:222065/1‑140 (MQ=255)         tttGCTGTTCTGCCAGTTGCGATTGCGCAGCCGCCACGCTATCGCCATCACTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTGCCAGATTTTAGAGGcc                                                                                                                         <  1:60394/140‑1 (MQ=255)                   tGCCAGTTGCGATTGCGCAGCCGCCACGCTATCGCCATCACTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTGCCAGATTTTAGAGGCCCCTTCGCGcc                                                                                                               <  2:239062/140‑1 (MQ=255)                    gCCAGTTGCGATTGCGCAGCCGCCACGCTATCGCCATCACTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTGCCAGATTTTAGAGGCCCCTTCGCGccc                                                                                                              <  1:467934/140‑1 (MQ=255)                          tGCGATTGCGCAGCCGCCACGCTATCGCCATCACTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTGCCAGATTTTAGAGGcc                                                                                                                         <  2:235470/123‑1 (MQ=255)                          tGCGATTGCGCAGCCGCCACGCTATCGCCATCACTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTGCCAGATTTTAGAGGcc                                                                                                                         >  1:235470/1‑123 (MQ=255)                          tGCGATTGCGCAGCCGCCACGCTATCGCCATCACTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTGCCAGATTTTAGAGGCCCCTTCGCGCCCGGCGtt                                                                                                        >  2:231807/1‑140 (MQ=255)                            cGATTGCGCAGCCGCCACGCTATCGCCATCACTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTGCCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCg                                                                                                      <  1:222065/140‑1 (MQ=255)                                      gccgccACGCTATCGCCATCACTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTGCCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCAt                                                                                            <  1:149955/140‑1 (MQ=255)                                      gccgccACGCTATCGCCATCACTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTGCCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCAt                                                                                            <  1:239838/140‑1 (MQ=255)                                      gccgccACGCTATCGCCATCACTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTGCCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCAt                                                                                            <  1:328977/140‑1 (MQ=255)                                                 aTCGCCATCACTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTGCCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGa                                                                                 >  1:836105/1‑140 (MQ=255)                                                 aTCGCCATCACTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTGCCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGa                                                                                 >  2:605548/1‑140 (MQ=255)                                                       aTCACTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTGCCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCaa                                                                           <  2:836105/140‑1 (MQ=255)                                                        tCACTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTGCCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCa                                                                            >  1:448973/1‑138 (MQ=255)                                                        tCACTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTGCCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCa                                                                            <  2:448973/138‑1 (MQ=255)                                                              aaGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTGCCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTc                                                                    >  1:458238/1‑140 (MQ=255)                                                              aaGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTGCCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTc                                                                    >  2:188440/1‑140 (MQ=255)                                                                    gagagaTATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTGCCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACg                                                              >  2:660195/1‑140 (MQ=255)                                                                    gagagaTATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTGCCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTT‑CCATCTGTTCCAGGACGGCAAAACCTTCGTCACgg                                                             >  2:826222/1‑140 (MQ=255)                                                                         atatTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTGCCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGAt                                                         <  1:188440/140‑1 (MQ=255)                                                                             tttttCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTGCCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCa                                                                <  1:710121/129‑1 (MQ=255)                                                                             tttttCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTGCCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCa                                                                >  2:710121/1‑129 (MQ=255)                                                                                     cgcCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTGCCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGCTACTAAACa                                             >  2:326064/1‑140 (MQ=255)                                                                                       ccGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTGCCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGAt                                                         <  2:23033/126‑1 (MQ=255)                                                                                       ccGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTGCCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGAt                                                         >  1:23033/1‑126 (MQ=255)                                                                                              acTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTGCCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGCTACTAAACAGCAATAgcg                                    <  2:667198/140‑1 (MQ=255)                                                                                              acTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTGCCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGCTACTAAACAGCAATAgcg                                    <  1:605548/140‑1 (MQ=255)                                                                                                tGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTGCCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGCTACTAAACAGCAATAGCGcc                                  <  2:458238/140‑1 (MQ=255)                                                                                                  gCATCACTGACGGGCATGTCTTTAATCTGCCCGTGCCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGCTACTAAACAGCAATAGCGCCAg                                >  2:466677/1‑140 (MQ=255)                                                                                                     tCACTGACGGGCATGTCTTTAATCTGCCCGTGCCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGCTACTAAACAGCAATAGCGCCa                                 <  2:230958/136‑1 (MQ=255)                                                                                                     tCACTGACGGGCATGTCTTTAATCTGCCCGTGCCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGCTACTAAACAGCAATAGCGCCa                                 >  1:230958/1‑136 (MQ=255)                                                                                                      cACTGACGGGCATGTCTTTAATCTGCCCGTGCCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGCTACTAAACAGCAATAGCGCCAGATTg                            <  1:231807/140‑1 (MQ=255)                                                                                                            cGGGCATGTCTTTAATCTGCCCGTGCCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGCTACTAAACAGCAATAGCGCCAGATTGTTTTGc                      >  1:32438/1‑140 (MQ=255)                                                                                                            cGGGCATGTCTTTAATCTGCCCGTGCCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGCTACTAAACAGCAATAGCGCCAGATTGTTTTGc                      >  1:636264/1‑140 (MQ=255)                                                                                                            cGGGCATGTCTTTAATCTGCCCGTGCCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGCTACTAAACAGCAATAGCGCCAGATTGTTTTGc                      >  1:418034/1‑140 (MQ=255)                                                                                                              ggCATGTCTTTAATCTGCCCGTGCCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGCTACTAAACAGCAATAGCGCCAGATTGTTTTGCAg                    >  1:37832/1‑140 (MQ=255)                                                                                                                cATGTCTTTAATCTGCCCGTGCCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGCTACTAAACAGCAATAGCGCCAGATTGTTTTGCAGGc                  <  2:37832/140‑1 (MQ=255)                                                                                                                cATGTCTTTAATCTGCCCGTGCCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGCTACTAAACAGCAATAGCGCCAGATTGTTTTGCAGGc                  >  2:862345/1‑140 (MQ=255)                                                                                                                  tGTCTTTAATCTGCCCGTGCCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGCTACTAAACAGCAATAGCGCCAGATTGTTTTGCAGGccc                >  1:157711/1‑140 (MQ=255)                                                                                                                         aaTCTGCCCGTGCCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGCTACTAAACAGCAATAGCGCCAGATTGTTTTGCAGGCCCGTATTGc         <  2:217214/140‑1 (MQ=255)                                                                                                                         aaTCTGCCCGTGCCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGCTACTAAACAGCAATAGCGCCAGATTGTTTTGCAGGCCCGTATTGc         <  2:109542/140‑1 (MQ=255)                                                                                                                          aTCTGCCCGTGCCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGCTACTAAACAGCAATAGCGCCAGATTGTTTTGCAGGCCCGTATTGcc        <  1:466677/140‑1 (MQ=255)                                                                                                                              gCCCGTGCCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGCTACTAAACAGCAATAGCGCCAGATTGTTTTGCAGGCCCGTATTGCCCGGt    >  1:106852/1‑140 (MQ=255)                                                                                                                               cccGTGCCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGCTACTAAACAGCAATAGCGCCAGATTGTTTTGCAGGCCCGTATTGCCCGGTg   >  1:739066/1‑140 (MQ=255)                                                                                                                                ccGTGCCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGCTACTAAACAGCAATAGCGCCAGATTGTTTTGCAGGCCCGTATTGCCCGGTGc  >  2:490081/1‑140 (MQ=255)                                                                                                                                    |                                                                                                                                        CAGCTGTTTTTGCTGTTCTGCCAGTTGCGATTGCGCAGCCGCCACGCTATCGCCATCACTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGCTACTAAACAGCAATAGCGCCAGATTGTTTTGCAGGCCCGTATTGCCCGGTGC  >  NC_000913/3688379‑3688645

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