Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,687,875	C→T	S433S (TCG→TCA)	bcsC ←	cellulose synthase outer membrane channel

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,687,875	0	C	T	100.0%	214.1 / NA	60	S433S (TCG→TCA)	bcsC	cellulose synthase outer membrane channel
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (27/33);  total (27/33)

GCTCCGGGTCGTTCGACTTCTGCTGCGCCAGATTGCGCATTAACGTATCGGCCTGGCTGCGTTGTCCGGCCTGCCAGAGATCCTGCGAAAGTCGGTAAGTAATCCATACGCTGCCGGGGTCCAGCGCCAGTCGTTGCCGCTGAAGTGCTGCCGCCTGCGCCCATTTGCCCTGGTTTTCCAGTGCCTCTGCCTGCTGTGCCAGACGGTCGTTTTGCAGGCTGCGTTCGATATCATCAATGCTACGCCGCTGACTGGCAGAGAGCGAGGCGATAAACGCTTCAGCTTTTTCTGGCGATTGCTGGCGGTAAATATTTGCCAGCCCGCGCACGGCGTTAGTGTTGCCGCTGTCCATACGCAAGGTCTGCTGATAATAACGTTCGGCGGCGGGATAATCTTTTCGCGCCATCGCCACATCGCCCAGCCCCAGCACTGCATAACTGTCGGTGTTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTCAGCAG  >  NC_000913/3687583‑3688141                                                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                           gCTCCGGGTCGTTCGACTTCTGCTGCGCCAGATTGCGCATTAACGTATCGGCCTGGCTGCGTTGTCCGGCCTGCCAGAGATCCTGCGAAAGTCGGTAAGTAATCCATACGCTGCCGGGGTCCAGCGCCAGTCGTTGCCGCTGAAGTGCTGCCGCCTGCGCCCATTTGCCCTGGTTTTCCAGTGCCTCTGCCTGCTGTGCCAGACGGTCGTTTTGCAGGCTGCGTTCGATATCATCAATGCTACGCCGCTGACTGGCAGAGAGCGAGGCGATAAACGCTTCAGCTTTTTCTGGTGATTGCt                                                                                                                                                                                                                                                                     >  2:192388/1‑300 (MQ=255) gCTCCGGGTCGTTCGACTTCTGCTGCGCCAGATTGCGCATTAACGTATCGGCCTGGCTGCGTTGTCCGGCCTGCCAGAGATCCTGCGAAAGTCGGTAAGTAATCCATACGCTGCCGGGGTCCAGCGCCAGTCGTTGCCGATGAAGTGCTGCCGCCTGCGCCCATTTGCCCTGGTTTTCCAGTGCCTCTGCCTGCTGTGCCAGACGGTCGTTTTGCAGGCTGCGTTCGATATCATCAATGCTACGCCGCTGACTGGCAGAGAGCGAGGCGATAAACGCTTCAGCTTTTTCTGGTGATTGCt                                                                                                                                                                                                                                                                     >  2:742280/1‑300 (MQ=255)      gggTCGTTCGACTTCTGCTGCGCCAGATTGCGCATTAACGTATCGGCCTGGCTGCGTTGTCCGGCCTGCCAGAGATCCTGCGAAAGTCGGTAAGTAATCCATACGCTGCCGGGGTCCAGCGCCAGTCGTTGCCGCTGAAGTGCTGCCGCCTGCGCCCATTTGCCCTGGTTTTCCAGTGCCTCTGCCTGCTNTGCCAGACGGTCGTTTTGCAGGCTGCGTTCGATATCATCAATGCTACGCCGCTGACTGGCAGAGAGCGAGGCGATAAACGCTTCAGCTTTTTCTGGTGATTGCTGGCgg                                                                                                                                                                                                                                                                >  1:642428/1‑300 (MQ=255)                  ttCTGCTGCGCCAGATTGCGCATTAACGTATCGGCCTGGCTGCGTTGTCCGGCCTGCCAGAGATCCTGCGAAAGTCGGTAAGTAATCCATACGCTGCCGGGGTCCAGCGCCAGTCGTTGCCGATGAAGTGCTGCCGCCTGCGCCCATTTGCCCTGGTTTTCCAGTGCCTCTGCCTGCTGTGCCAGACGGTCGTTTTGCAGGCTGCGTTCGATATCATCAATGCTACGCCGCTGACTGGCAGAGAGCGAGGCGATAAACGCTTCAGCTTTTTCTGGTGATTGCTGGCGGTAAATATTTGcc                                                                                                                                                                                                                                                    <  1:758716/300‑1 (MQ=255)                         gcgcCAGATTGCGCATTAACGTATCGGCCTGGCTGCGTTGTCCGGCCTGCCAGAGATCCTGCGAAAGTCGGTAAGTAATCCATACGCTGCCGGGGTCCAGCGCCAGTCGTTGCCGCTGAAGTGCTGCCGCCTGCGCCCATTTGCCCTGGTTTTCCAGTGCCTCTGCCTGCTGTGCCAGACGGTCGTTTTGNAGGCTGCGTTCGATATCATCAATGCTACGCCGCTGACTGGCAGAGAGCGAGGCGATAAACGCTTCAGCTTTTTCTGGTGATTGCTGGCGGTAAATATTTGCCAGCCcgc                                                                                                                                                                                                                                             >  1:707672/1‑300 (MQ=255)                            ccAGATTGCGCATTAACGTATCGGCCTGGCTGCTTTGTCCGGCCTGCCAGAGATCCTTCGAAAGTCGGTAAGTAATCCATACTCTGCCGGGGTCCAGCGCCAGTCGTTGCCGCTGAAGTGCTGCCGCCTGCGCCCATTTGCCCTGGTTTTCCAGTGCCTCTGCCTGCTGTGCCAGACGGTCGTTTTGCAGGCTGCGTTCGATATCATCAATGCTACGCCGCTGACTGGCAGAGAGCGAGGCGATAAACGCTTCAGCTTTTTCTGGTGATTGCTGGCGGTAAATATTTGCCAGCCCGCGCa                                                                                                                                                                                                                                          <  2:553262/300‑1 (MQ=255)                            ccAGATTGCGCATTAACGTATCGGCCTGGCTGCGTTGTCCGGCCTGCCAGAGATCCTGCGAAAGTCGGTAAGTAATCCATACGCTGCCGGGGTCCAGCGCCAGTCGTTGCCGCTGAAGTGCTGCCGCCTGCGCCCATTTGCCCTGGTTTTCCAGTGCCTCTGCCTGCTGTGCCAGACGGTCGTTTTGCAGGCTGCGTTCGATATCATCAATGCTACGCCGCTGACTGGCAGAGAGCGAGGCGATAAACGCTTCAGCTTTTTCTGGTGATTGCTGGCGGTAAATATTTGCCAGCCCGCGCa                                                                                                                                                                                                                                          <  1:750897/300‑1 (MQ=255)                             cAGATTGCGCATTAACGTATCGGCCTGGCTGCGTTGTCCGGCCTGCCAGAGATCCTGCGAAAGTCGGTAAGTAATCCATACGCTGCCGGGGTCCAGCGCCAGTCGTTGCCGCTGAAGTGCTGCCGCCTGCGCCCATTTGCCCTGGTTTTCCAGTGCCTCTGCCTGCTGTGCCAGACGGTCGTTTTGCAGGCTGCGTTCGATATCATCAATGCTACGCCGCTGACTGGCAGAGAGCGAGGCGATAAACGCTTCAGCTTTTTCTGGTGATTGCTGGCGGTAAATATTTGCCAGCCCGCGCAc                                                                                                                                                                                                                                         <  2:317146/300‑1 (MQ=255)                                     gcATTAACGTATCGGCCTGGCTGCGTTGTCCGGCCTGCCAGAGATCCTGCGAAAGTCGGTAAGTAATCCATACGCTGCCGGGGTCCAGCGCCAGTCGTTGCCGCTGAAGTGCTGCCGCCTGCGCCCATTTGCCCTGGTTTTCCAGTGCCTCTGCCTGCTGTGCCAGACGGTCGTTTTGCAGGCTGCGTTCGATATCATCAATGCTACGCCGCTGACTGGCAGAGAGCGAGGCGATAAACGCTTCAGCTTTTTCTGGTGATTGCTGGCGGTAAATATTTGCCAGCCCGCGCACGGCGTTAg                                                                                                                                                                                                                                 >  1:5117/1‑300 (MQ=255)                                      cATTAACGTATCGGCCTGGCTGCGTTGTCCGGCCTGCCAGAGATCCTGCGAAAGTCGGTAAGTAATCCATACGCTGCCGGGGTCCAGCGCCAGTCGTTGCCGCTGAAGTGCTGCCGCCTGCGCCCATTTGCCCTGGTTTTCCAGTGCCTCTGCCTGCTGTGCCAGACGGTCGTTTTGCAGGCTGCGTTCGATATCATCAATGCTACGCCGCTGACTGGCAGAGAGCGAGGCGATAAACGCTTCAGCTTTTTCTGGTGATTGCTGGCGGTAAATATTTGCCAGCCCGCGCACGGCGTTAgt                                                                                                                                                                                                                                >  1:903708/1‑300 (MQ=255)                                       aTTAACGTATCGTCCTGGCTGCGTTGTCCGGCCTGCCAGAGATCCTGCGAAAGTCGGTAAGTAATCCATACGCTGCCGGGGTCCAGCGCCAGTCGTTGCCGCTGAAGTGCTGCCGCCTGCGCCCATTTGCCCTGGTTTTCCAGTGCCTCTGCCTGCTGTGCCAGACGGTCGTTTTGCAGGCTGCGTTCGATATCATCAATGCTACGCCGCTGACTGGCAGAGAGCGAGGCGATAAACGCTTCAGCTTTTTCTGGTGATTGCTGGCGGTAAATATTTGCCAGCCCGCGCACGGCGTTAgtg                                                                                                                                                                                                                               <  2:335897/300‑1 (MQ=255)                                       aTTAACGTATCGGCCTGGCTGCGTTGTCCGGCCTGCCAGAGATCCTGCGAAAGTCGGTAAGTAATCCATACGCTGCCGGGGTCCAGCGCCAGTCGTTGCCGCTGAAGTGNTGCCGCCTGCGCCCATTTGCCCTGGTTTTCCAGTGCCTCTGCCTGCTGTGCCAGACGGTCGTTTTGCAGGCTGCGTTCGATATCATCAATGCTACGCCGCTGACTGGCAGAGAGCGAGGCGATAAACGCTTCAGCTTTTTCTGGTGATTGCTGGCGGTAAATATTTGCCAGCCCGCGCACGGCGTTAgtg                                                                                                                                                                                                                               <  1:640381/300‑1 (MQ=255)                                                   gCCTGGCTGCGTTGTCCGGCCTGCCAGAGATCCTGCGAAAGTCGGTAAGTAATCCATACGCTGCCGGGGTCCAGCGCCAGTCGTTGCCGCTGAAGTGCTGCCGCCTGCGCCCATTTGCCCTGGTTTTCCAGTGCCTCTGCCTGCTGTGCCAGACGGTCGTTTTGCAGGCTGCGTTCGATATCATCAATGCTACGCCGCTGACTGGCAGAGAGCGAGGCGATAAACGCTTCAGCTTTTTCTGGTGATTGCTGGCGGTAAATATTTGCCAGCCCGCGCACGGCGTTAGTGTTGCCGCTGTcc                                                                                                                                                                                                                   >  1:330425/1‑300 (MQ=255)                                                        gCTGCGTTGTCCGGCCTGCCAGAGATCCTGCGAAAGTCGGTAAGTAATCCATACGCTGCCGGGGTCCAGCGCCATTCGTTGCCGCTGAAGTGCTGCCGCCTGCGCCCATTTGCCCTGGTTTTCCAGTGCCTCTGCCTGCTGTGCCAGACGGTCGTTTTGCAGGCTGCGTTCGATATCATCAATGCTACGCCGCTGACTGGCAGAGAGCGAGGCGATAAACGCTTCAGCTTTTTCTGGTGATTGCTGGCGGTAAATATTTGCCAGCCCGCGCACGGCGTTAGTGTTGCCGCTGTCCAGACg                                                                                                                                                                                                              <  2:535510/300‑1 (MQ=255)                                                               tGTCCGGCCTGCCAGAGATCCTGCGAAAGTCGGTAAGTAATCCATACGCTGCCGGGGTCCAGCGCCAGTCGTTGCCGCTGAAGTGCTGCCGCCTGCGCCCATTTGCCCTNGTTTTCCAGTGCCTCTGCCTGCTGTGCCAGACGGTCGTTTTGCAGGCTGCGTTCGATATCATCAATGCTACGCCGCTGACTGGCAGAGAGCGAGGCGATAAACGCTTCAGCTTTTTCTGGTGATTGCTGGCGGTAAATATTTGCCAGCCCGCGCACGGCGTTAGTGTTGCCGCTGTCCATACGCAAGGTc                                                                                                                                                                                                       <  1:701367/300‑1 (MQ=255)                                                               tGTCCGGCCTGCCAGAGATCCTGCGAAAGTCGGTAAGTAATCCATACGCTGCCGGGGTCCAGCGCCAGTCGTTGCCGCTGAAGTGCTGCCGCCTGCGCCCATTTGCCCTGGTTTTCCAGTGCCTCTGCCTGCTGTGCCAGACGGTCGTTTTGCAGGCTGCGTTCGATATCATCAATGCTACGCCGCTGACTGGCAGAGAGCGAGGCGATAAACGCTTCAGCTTTTTCTGGTGATTGCTGGCGGTAAATATTTGCCAGCCCGCGCACGGCGTTAGTGTTGCCGCTGTCCATACGCAAGGTc                                                                                                                                                                                                       <  1:365665/300‑1 (MQ=255)                                                                      ccTGCCAGAGATCCTGCGAAAGTCGGTAAGTAATCCATACGCTGCCTGGGTCCAGCGCCAGTCGTTGCCGCTGAAGTGCTGCCGCCTGCGCCCATTTGCCCTGGTTTTCCAGTGCCTCTGCCTGCTGTGCCANACGGTCGTTTTGCAGGCTGCGTTCGATATCATCAATGCTACGCCGCTGACTGGCAGAGAGCGAGGCGATAAACGCTTCAGCTTTTTCTGGTGATTGCTGGCGGTAAATATTTGCCAGCCCGCGCACGGCGTTAGTGTTGCCGCTGTCCATACGCAAGGTCTGCTGat                                                                                                                                                                                                <  2:216988/300‑1 (MQ=255)                                                                       cTGCCAGAGATCCTGCGAAAGTCGGTAGGTAATCCATACGCTGCCGGGGTCCAGCGCCAGTCGTTGCCGCTGAAGTGCTGCCGCCTGCGCCCATTTGCCCTGGTTTTCCAGTGCCTCTGCCTGCTGTGCCAGACGGTCGTTTTGCAGGCTGCGTTCGATATCATCAATGCTACGCCGCTGACTGGCAGAGAGCGAGGCGATAAACGCTTCAGCTTTTTCTGGTGATTGCTGGCGGTAAATATTTGCCAGCCCGCGCACGGCGCTAGTGTTGCCGCTGCCCATACGCAAGGTCTGCTGata                                                                                                                                                                                               >  2:315786/1‑300 (MQ=255)                                                                         gCCAGAGATCCTGCGAAAGTCGGTAAGTAATCCATACGCTGCCGGGGTCCAGCGCCAGTCGTTGCCGCTGAAGTGCTGCCGCCTGCGCCCATTTGCCCTGGTTTTCCAGTGCCTCTGCCTGCTGTGCCAGACGGTCGTTTTGCAGGCTGCGTTCGATATCATCAATGCTACGCCGCTGACTGGCAGAGAGCGAGGCGATAAACGCTTCAGCTTTTTCTGGTGATTGCTGGCGGTAAATATTTGCCAGCCCGCGCACGGCGTTAGTGTTGCCGCTGTCCATACGCAAGGTCTGCTGataat                                                                                                                                                                                             >  1:184556/1‑300 (MQ=255)                                                                         gCCAGAGATCCTGCGAAAGTCGGTAAGTAATCCATACGCTGCCGGGGTCCAGCGCCAGTCGTTGCCGCTGAAGTGCTGCCGCCTGCGCCCATTTGCCCTGGTTTTCCAGNGCCTCTGCCTGCTGTGCCAGACGGTCGTTTTGCAGGCTGCGTTCGATATCATCAATGCTACGCCGCTGACTGGCAGAGAGCGAGGCGATAAACGCTTCAGCTTTTTCTGGTGATTGCTGGCGGTAAATATTTGCCAGCCCGCGCACGGCGTTAGTGTTGCCGCTGTCCATACGCAAGGTCTGCTGataat                                                                                                                                                                                             <  1:645257/300‑1 (MQ=255)                                                                           cAGAGATCCTGCGAAAGTCGGTAAGTAATCCATACGCTGCCGGGGTCCAGCGCCAGTCGTTGCCGCTGAAGTGCTGCCGCCTGCGCCCATTTGCCCTGGTTTTCCAGTGCCTCTGCCTGCTGTGCCAGACGGTCGTTTTGCAGGCTGCGTTCGATATCATCAATGCTACGCCGCTGACTGGCAGAGAGCGAGGCGATAAACGCTTCAGCTTTTTCTGGTGATTGCTGGCGGTAAATATTTGCCAGCCCGCGCACGGCGTTAGTGTTGCCGCTGTCCATACGCAAGGTCTGCTGAtaataa                                                                                                                                                                                           >  1:541948/1‑300 (MQ=255)                                                                            agagATCCTGCGAAAGTCGGTAAGTAATCCATACGCTGCCGGGGTCCAGCGCCAGTCGTTGCCGCTGAAGTGCTGCCGCCTGCGCCCATTTGCCCTGGTTTTCCAGTGCCTCTGCCTGCTGTGCCAGACGGTCGTTTTGCAGGCTGCGTTCGATATCATCAATGCTACGCCGCTGACTGGCAGAGAGCGAGGCGATAAACGCTTCAGCTTTTTCTGGTGATTGCTGGCGGTAAATATTTGCCAGCCCGCGCACGGCGTTAGTGTTGCCGCTGTCCATACGCAAGGTCTGCTGATAATAAc                                                                                                                                                                                          <  1:515026/300‑1 (MQ=255)                                                                                                     aaTCCATACGCTGCCGGGGTCCAGCGCCAGTCGTTGCCGCTGAAGTGCTGCCGCCTGCGCCCATTTGCCCTGGTTTTCCAGTGCCTCTGCCTGCTGTGCCAGACGGTCGTTTTGCAGGCTGCGTTCGATATCATCAATGCTACGCCGCTGACTGGCAGAGAGCGAGGCGATAAACGCTTCAGCTTTTTCTGGTGATTGCTGGCGGTAAATATTTGCCAGCCCGCGCACGGCGTTAGTGTTGCCGCTGTCCATACGCAAGGTCTGCTGATAATAACGTTCGGCGGCGGGATAATCTTTTcg                                                                                                                                                                 >  1:794532/1‑300 (MQ=255)                                                                                                              gCTGCCGGGGTCCAGCGCCAGTCGTTGCCGCTGAAGTGCTGCCGCCTGCGCCCATTTGCCCTGGTTTTCCAGTGCCTCTGCCTGCTGTGCCAGACGGTCGTTTTGCAGGCTGCGTTCGATATCATCAATGCTACGCCGCTGACTGGCAGAGAGCGAGGCGATAAACGCTTCAGCTTTTTCTGGTGATTGCTGGCGGTAAATATTTGCCAGCCCGCGCACGGCGTTAGTGTTGCCGCTGTCCATACGCAAGGTCTGCTGATAATAACGTTCGGCGGCGGGATAATCTTTTCGCGCCATCGc                                                                                                                                                        <  1:572155/300‑1 (MQ=255)                                                                                                                      ggTCCAGCGCCATTCGTTGCCGCTGAAGTGCTGCCGCCTGCGCCCATTTGCCCTGGTTTTCCAGTGCCTCTGCCTGCTGTGCCAGACGGTCGTTTTGCAGGCTGCGTTCGATATCATCAATGCTACGCCGCTGACTGGCAGAGAGCGAGGCGATAAACGCTTCAGCTTTTTCTGGTGATTGCTGGCGGTAAATATTTGCCAGCCCGCGCACGGCGTTAGTGTTGCCGCTGTCCATACGCAAGGTCTGCTGATAATAACGTTCGGCGGCGGGATAATCTTTTCGCGCCATCGCCACATCGc                                                                                                                                                <  2:911041/300‑1 (MQ=255)                                                                                                                       gTCCAGCGCCAGTCGTTGCCGCTGAAGTGCTGCCGCCTGCGCCCATTTGCCCTGGTCTTCCAGTGCCTCTGCCTGCTGTGCCAGACGGTCGTTTTGCAGGCTGCGTTCGATATCATCAATGCTACGCCGCTGACTGGCAGAGAGCGAGGCGATAAACGCTTAAGCTTTTTCTGGTGATTGCTGGCGGTAAATATTTGCCAGCCCGCGCACGGCGTTAGTGTTGCCGCTGTCCATACGCAAGGTCTGCTGATAATAACGTTCGGCGGCGGGATAATCTTTTCGCGCCATCGCCACATCgac                                                                                                                                               >  2:601616/1‑298 (MQ=255)                                                                                                                              gcCAGTCGTTGCCGCTTAAGTGCTGCCGCCTGCGCCCATTTGCCCTGGTTTTCCAGTGCCTCTGCCTGCTGTGCCAGACGGTCGTTTTGCAGGCTGCGTTCGATATCATCAATGCTACGCCGCTGACTGGCAGAGAGCGAGGCGATAAACGCTTCAGCTTTTTCTGGTGATTGCTGGCGGTAAATATTTGCCAGCCCGCGCACGGCGTTAGTGTTGCCGCTGTCCATACGCAAGGTCTGCTGATAATAACGTTCGGCGGCGGGATAATCTTTTCGCGCCATCGCCACATCGcccagcccc                                                                                                                                        <  1:618973/300‑1 (MQ=255)                                                                                                                                   tCGTTGCCGCTGAAGTGCTGCCGCCTGCGCCCATTTGCCCTGGTTTTCCAGTGCCTCTGCCTGCTGTGCCAGACGGTCGTTTTGCAGGCTGCGTTCGATATCATCNATGCTACGCCGCTGACTGGCAGAGAGCGAGGCGATAAACGCTTCAGCTTTTTCTGGTGATTGCTGGCGGTAAATATTTGCCAGCCCGCGCACGGCGTTAGTGTTGCCGCTGTCCATACGCAAGGTCTGCTGATAATAACGTTCGGCGGCGGGATAATCTTTTCGCGCCATCGCCACATCGCCCAGCCCCAGCAc                                                                                                                                   <  2:840004/300‑1 (MQ=255)                                                                                                                                   tCGTTGCCGCTGAAGTGCTGCCGCCTGCGCCCATTTGCCCTGGTTTTCCAGTGCCTCTGCCTGCTGTGCCAGACGGTCGTTTTGCAGGCTGCGTTCGATATCATCAATGCTACGCCGCTGACTGGCAGAGAGCGAGGCGATAAACGCTTCAGCTTTTTCTGGTGATTGCTGGCGGTAAATATTTGCCAGCCCGCGCACGGCGTTAGTGTTGCCGCTGTCCATACGCAAGGTCTGCTGATAATAACGTTCGGCGGCGGGATAATCTTTTCGCGCCATCGCCACATCGCCCAGCCCCAGCAc                                                                                                                                   <  1:821631/300‑1 (MQ=255)                                                                                                                                   tCGTTGCCGCTGAAGTGCTGCCGCCTGCGCCCATCTGCCCTGGTTTTCCAGTGCCTCTGCCTGCTGTGCCAGACGGTCGTTTTGCAGGCTGCGTTCGATATCATCAATGCTACGCCGCTGACTGGCAGAGAGCGAGGCGATAAACGCTTCAGCTTTTTCTGGTGATTGCTGGCGGTAAATATTTGCCAGCCCGCGCACGGCGTTAGTGTTGCCGCTGTCCATACGCAAGGTCTGCTGATAATAACGTTCGGCGGCGGGATAATCTTTTCGCGCCAGCGCCACATCGCCCAGCCCCAGCAc                                                                                                                                   <  2:354576/300‑1 (MQ=255)                                                                                                                                   tCGTTGCCGCTGAAGTGCGGCCGCCGGCGCCCATTTGCCCTGGTTGTCCAGTGCCTCTGCCTGCTGTGCCAGACGGTCGTTTTGCAGGCTGCGTTCGATATCATCAATGCTACGCCGCTGACTGGCAGAGAGCGAGGCGATAAACGCTTCAGCTTTTTCTGGTGATTGCTGGCGGTAAATATTTGCCAGCCCGCGCACGGCGTTAGTGTTGCCGCTGTCCATACGCAAGGTCTGCTGATAATAACGTTCGGCGGCGGGATAATCTTTTCGCGCCATCGCCACATCGCCCAGCCCCAGCAc                                                                                                                                   <  1:895113/300‑1 (MQ=255)                                                                                                                                         ccGCTGAAGTGCTGCCGCCTGCGCCCATTTGCCCTGGTTTTCCAGTGCCTCTGCCTGCTGTGCCAGACGGTCGTTTTGCAGGCTGCGTTCGATATCATCAATGCTACGCCGCTGACTGGCAGAGAGCGAGGCGATAAACGCTTCAGCTTTTTCTGGTGATTGCTGGCGGTAAATATTTGCCAGCCCGCGCACGGCGTTAGTGTTGCCGCTGTCCATACGCAAGGTCTGCTGATAATAACGTTCGGCGGCGGGATAATCTTTTCGCGCCATCGCCACACCGCCCAGCCCCAGCACTGCaca                                                                                                                             >  2:102723/1‑298 (MQ=255)                                                                                                                                           gcggAAGTGCTGCCGCCTGCGCCCATTTGCCCTGGTTTTCCAGTGCCTCTGCCTGCTGTGCCAGACGGTCGTTTTGCAGGCTGCGTTCGATATCATCAATGCTACGCCGCTGACTGGCAGAGAGCGAGGCGATAAACGCTTCAGCTTTTTCTGGTGATTGCTGGCGGTAAATATTTGCCAGCCCGCGCACGGCGTTAGTGTTGCCGCTGTCCATACGCAAGGTCTGCTGATAATAACGTTCGGCGGCGGGATAATCTTTTCGCGCCATCGCCACATCGCCCAGCCCCAGCACTGCATAAc                                                                                                                           <  1:932089/297‑1 (MQ=255)                                                                                                                                           gCTGAAGTGCTGCCGCCTGCGCCCATTTGCCCTGGTTTTCCAGTGCCTCTGCCTGCTGTGCCAGACGGTCGTTTTGCAGGCTGCGTTCGATATCATCAATGCTACGCCGCTGACTGGCAGAGAGCGAGGCGATAAACGCTTCAGCTTTTTCTGGTGATTGCTGGCGGTAAATATTTGCCAGCCCGCGCACGGCGTTAGTGTTGCCGCTGTCCATACGCAAGGTCTGCTGATAATAACGTTCGGCGGCGGGATAATCTTTTCGCGCCATCGCCACATCGCCCAGCCCCAGCACTGCATAAc                                                                                                                           <  1:271785/300‑1 (MQ=255)                                                                                                                                           gCTGAAGTGCTGCCGCCTGCGCCCATTTGCCCTGGTTTTCCAGTGCCTCTGCCTGCTGTGCCAGACGGTCGTTTTGCAGGCTGCGTTCGATATCATCAATGCTACGCCGCTGACTGGCAGAGAGCGAGGCGATAAACGCTTCAGCTTTTTCTGGTGATTGCTGGCGGTAAATATTTGCCAGCCCGCGCACGGCGTTAGTGTTGCCGCTGTCCATACGCAAGGTCTGCTGATAATAACGTTCGGCGGCGGGATAATCTTTTCGCGCCATCGCCACATCGCCCAGCCCCAGCACTGCATAAc                                                                                                                           <  2:871804/300‑1 (MQ=255)                                                                                                                                             tGAAGTGCTGCCGCCTGCGCCCATTTGCCCTGGTTTTCCAGTGCCTCTGCCTGCTGTGCCAGACGGTCGTTTTGCAGGCTGCGTTCGATATCATCAATGCTACGCNGCTGACTGGCAGAGAGCGAGGCGATAAACGCTTCAGCTTTTTCTGGTGATTGCTGGCGGTAAATATTTGCCAGCCCGCGCACGGCGTTAGTGTTGCCGCTGTCCATACGCAAGGTCTGCTGATAATAACGTTCGGCGGCGGGATAATCTTTTCGCGCCATCGCCACATCGCCCAGCCCCAGCACTGCATAACTg                                                                                                                         <  2:827044/300‑1 (MQ=255)                                                                                                                                                 gTGCTGCCGCCTGCGCCCATTTGCCCTGGTTTTCCAGTGCCTCTGCCTGCTGTGCCAGACGGTCGTTTTGCAGGCTGCGTTCGATATCATCAATGCTACGCCGCTGACTGGCAGAGAGCGAGGCGATAAACGCTTCAGCTTTTTCTGGTGATTGCTGGCGGTAAATATTTGCCAGCCCGCGCACGGCGTTAGTGTTGCCGCTGTCCATACGCAAGGTCTGCTGATAATAACGTTCGGCGGCGGGATAATCTTTTCGCGCCATCGCCACATCGCCCAGCCCCAGCACTGCATAACTGTCgg                                                                                                                     >  1:620600/1‑300 (MQ=255)                                                                                                                                                         gccTGCGCCCATTTGCCCTGGTTTTCCAGTGCCTCTGCCTGCTGTGCCAGACGGTCGTTTTGCAGGCTGCGTTCGATATCATCAATGCTACGCCGCTGACTGGCAGAGAGCGAGGCGATAAACGCTTCAGCTTTTTCTGGTGATTGCTGGCGGTAAATATTTGCCAGCCCGCGCACGGCGTTAGTGTTGCCGCTGTCCATACGCAAGGTCTGCTGATAATAACGTTCGGCGGCGGGATAATCTTTTCGCGCCATCGCCACATCGCCCAGCCCCAGCACTGCATAACTGTCGGTGTTATCg                                                                                                             <  1:281072/300‑1 (MQ=255)                                                                                                                                                                                 tCCAGTGCCTCTGCCTGCTGTGCCAGACGGTCGTTTTGCAGGCTGCGTTCGATATCATCAATGCTACGCCGCTGACTGGCAGAGAGCGAGGCGATAAACGCTTCAGCTTTTTCTGGTGATTGCTGGCGGTAAATATTTGCCAGCCCGCGCACGGCGTTAGTGTTGCCGCTGTCCATACGCAAGGTCTGCTGATAATAACGTTCGGCGGCGGGATAATCTTTTCGCGCCATCGCCACATCGCCCAGCCCCAGCACTGCATAACTGTCGGTGTTATCGACATTACGCGCCTGCTGGAACAgg                                                                                     >  1:933607/1‑300 (MQ=255)                                                                                                                                                                                  ctaGTGCCTCTGCCTGCTGTGCCAGACGGTCGTTTTGCATGTTGCGTTCGATATCATCAATGCTACGCCGCTGACTTGCAGAGAGCGAGGCGATAAACGCTTCAGCTTTTTCTGGTGATTGCTGGCGGTAAATATTTGCCAGCCCGCGCACGGCGTTAGTGTTGCCGCTGTCCATACGCAAGGTCTGCTGATAATAACGTTCGGCGGCGGGATAATCTTTTCGCGCCATCGCCACATCGCCCAGCCCCAGCACTGCATAACTGTCGGTGTTATCGACATTACGCGCCTGCTGGAATAGGc                                                                                    <  2:296508/298‑1 (MQ=255)                                                                                                                                                                                   cAGTGCCTCTGCCTGCTGTGCCAGACGGTCGTTTTGCAGGCTGCGTTCGATATCATCAATGCTACGCCGCTGACTGGCAGAGAGCGAGGCGATAAACGCTTCAGCTTTTNCTGGTGATTGCTGGCGGTAAATATTTGCCAGCCCGCGCACGGCGTTAGTGTTGCCGCTGTCCATACGCAAGGTCTGCTGATAATAACGTTCGGCGGCGGGATAATCTTTTCGCGCCATCGCCACATCGCCCAGCCCCAGCACTGCATAACTGTCGGTGTTATCGACATTACGCGCCTGCTGGAACAGGCg                                                                                   <  1:677299/300‑1 (MQ=255)                                                                                                                                                                                       gCCTCTGCCTGCTGTGCCAGACGGTCGTTTTGCAGGCTGCGTTCGATATCATCAATGCTACGCCGCTGAATGGCAGAGAGCGAGGCGATAAACGCTTCAGCTTTTTCTGGTGATTGCTGGCGGTAAATATTTGCCAGCCCGCGCACGGCGTTAGTGTTGTCGCTGTCCCTACGCCAGGCTTGCTGATCCTAACGCTCGTCCGCAGGATAGCTTTTTCGCGCCATCGCCACATCGCACAGCACCATTACTGAATAACTGTCGGTGTTATCGAAATTAAGAGCCTGCTGGAAAAGGAGTTCt                                                                               >  2:375843/1‑300 (MQ=255)                                                                                                                                                                                          tctGCCTGCTGTGCCAGACGGTCGTTTTGCAGGCTGCGTTCGATATCATCAATGCTACGCCGCTGACTGGCAGAGAGCGAGGCGATAAACGCTTCAGCTTTTTCTGGTGATTGCTGGCGGTAAATATTTGCCAGCCCGCGCACGGCGTTAGTGTTGCCGCTGTCCATACGCAAGGTCTGCTGATAATAACGTTCGGCGGCGGGATAATCTTTTCGCGCCATCGCCACATCGCCCAGCCCCAGCACTGCATAACTGTCGGTGTTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGcc                                                                            >  1:719845/1‑300 (MQ=255)                                                                                                                                                                                                    gtgCCAGACGGTCGTTTTGCAGGCTGCGTTCGATATCATCAATGCTACGCCGCTGACTGGCAGAGAGCGAGGCGATAAACGCTTCAGCTTTTTCTGGTGATTGCTGGCGGTAAATATTTGCCAGNCCGCGCACGGCGTTAGTGTTGCCGCTNTCCATACGCAAGGTCTGCTGATAATAACGTTCGGCGGCGGGATAATCTTTTCGCGCCATCGCCACATCGCCCAGCCCCAGCACTGCATAACTGTCGGTGTTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGat                                                                  >  2:651060/1‑300 (MQ=255)                                                                                                                                                                                                    gtgCCAGACGGTCGTTTTGCAGGCTGCGTTCGATATCATCAATGCTACGCCGCTGACTGGCAGAGAGCGAGGCGATAAACGCTTCAGCTTTTTCTGGTGATTGCTGGCGGTAAATATTTGCCAGCCCGCGCACGGCGTTAGTGTTGCCGCTGTCCATACGCAAGGTCTGCTGATAATAACGTTCGGCGGCGGGATAATCTTTTCGCGCCATCGCCACATCGCCCAGCCCCAGCACTGCATAACTGTCGGTGTTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGat                                                                  >  1:769717/1‑300 (MQ=255)                                                                                                                                                                                                     tgCCAGACGGTCGTTTTGCAGGCTGCGTTCGATATCATCAATGCTACGCCGCTGACTGGCAGAGAGCGAGGCGATAAACGCTTCAGCTTTTTCTGGTGATTGCTGGCGGTAAATATTTGCCAGCCCGCGCACGGCGTTAGTGTTGCCGCTGTCCATACGCAAGGTCTGCTGATAATAACGTTCGGCGGCGGGATAATCTTTTCGCGCCATCGCCACATCGCCCAGCCCCAGCACTGCATAACTGTCGGTGTTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGatt                                                                 <  1:410504/300‑1 (MQ=255)                                                                                                                                                                                                         agacattCGTTTTGCAGGCTGCGTTCGATATCATCAATGCTACGCCGCTGACTGGCAGAGAGCGAGTCGATAAACGCTTCAGCTTTTTCTGGTGATTGCTGGCGGTAAATATTTGCCAGCCCGCGCACTGCGTTAGTGTTGCCGCTGTCCATACGCAAGGTCTGCTAATAATAACGTTCGGCGGCGGGATAATCTTTTCGCGCCATCGCCACATCGCCCAGCCCCAGCACTGCATAACTGTCGGTGTTATCGACATTACGCGCCTGTTGGAACAGGCGTTCTGCCCGGTCAGGATTATTg                                                             <  2:139604/294‑1 (MQ=255)                                                                                                                                                                                                         aGACGGTCGTTTTGCAGGCTGCGTTCGATATCATCAATGCTACGCCGCTGACTGGCAGAGAGCGAGGCGATAAACGCTTCAGCTTTTTCTGGTGATTGCTGTCGGTAAATATTTGCCAGCCCGCGCACGGCGTTAGTGTTGCCGCTGTCCATACGCAAGGTCTGCTGATAATAACGTTCGGCGGCGGGATAATCTTTTCGCGCCATCGCCACATCGCCCAGCCCCAGCACTGCATAACTGTCGGTGTTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGattatgg                                                             <  1:796955/300‑3 (MQ=255)                                                                                                                                                                                                         aGACGGTCGTTTTGCAGGCTGCGTTCGATATCATCAATGCTACGCCGCTGACTGGCAGAGAGCGAGGCGATAAACGCTTCAGCTTTTTCTGGTGATTGCTGGCGGTAAATATTTGCCAGCCCGCGCACGGCGTTAGTGTTGCCGCTGTCCATACGCAAGGTCTGCTGATAATAACGTTCGGCGGCGGGATAATCTTTTCGCGCCATCGCCACATCGCCCAGCCCCAGCACTGCATAACTGTCGGTGTTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGATTATTg                                                             <  2:336100/300‑1 (MQ=255)                                                                                                                                                                                                                             gCGTTCGATATCATCAATGCTACGCCGCTGACTGGCAGAGAGCGAGGCGATAAATGCTTCAGCTTTTTCTGGTGATTGCTGGCGGTAAATATTTGCCAGCCCGCGCACGGCGTTAGTGTTGCCGCTGTCCATACGCAAGGTCTGCTGATAATAACGTTCGGCGGCGGGATAATCTTTTCGCGCCATCGCCACATCGCCCAGCCCCAGCACTGCATAACTGTCGGTGTTATCGACATTACGCGCATGCTGGAACAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGcc                                         >  2:252941/1‑300 (MQ=255)                                                                                                                                                                                                                             gCGTTCGATATCATCAATGCTACGCCGCTGACTGGCAGAGAGCGAGGCGATAAACGCTTCAGCTTTTTCTGGTGATTGCTGGCGGTAAATATTTGCCAGCCCGCGCACGGCGTTAGTGTTGCCGCTGTCCATACGCAAGGTCTGCTGATAATAACGTTCGGCGGCGGGATAATCTTTTCGCGCCATCGCCACATCGCCCAGCCCCAGCACTGCATAACTGTCGGTGTTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGcc                                         >  1:851017/1‑300 (MQ=255)                                                                                                                                                                                                                              cGTTCGATATCATCAATGCTACGCCGCTGACTGGCAGAGAGCGAGGCGATAAACGCTTCAGCTTTTTCTGGTGATTGCTGGCGGTAAATATTTGCCAGCCCGCGCACGGCGTTAGTGTTGCCGCTGTCCATACGCAAGGTCTGCTGATAATAACGTTCGGCGGCGGGATAATCTTTTCGCGCCATCGCCACATCGCCCAGCCCCAGCACTGCATAACTGTCGGTGTTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGccc                                        <  1:460830/300‑1 (MQ=255)                                                                                                                                                                                                                               gTTCGATATCATCAATGCTACGCCGCTGACTGGCAGAGAGCGAGGCGATAAACGCTTCAGCTTTTTCTGGTGATTGCTGGCGGTAAATATTTGCCAGCCCGCGCACGGCGTTAGTGTTGCCGCTGTCCATACGCAAGGTCTGCTGATAATAACGTTCGGCGGCGGGATAATCTTTTCGCGCCATCGCCACATCGCCCAGCCCCAGCACTGCATAACTGTCGGTGTTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGCCct                                       >  2:736137/1‑300 (MQ=255)                                                                                                                                                                                                                                      atcatcAATGCTACGCCGCTGACTGGCAGAGAGCGAGGCGATAAACGCTTCAGCTTTTTCTGGTGATTGCTGGCGGTAAATATTTGCCAGCCCGCGCACGGCGTTAGTGTTGCCGCTGTCCATACGCAAGGTCTGCTGATAATAACGTTCGGCGGCGGGATAATCTTTTCGCGCCATCGCCACATCGCCCAGCCCCAGCACTGCATAACTGTCGGTGTTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGAt                                >  2:620471/1‑300 (MQ=255)                                                                                                                                                                                                                                        catcaATGCTACGCCGCTGACTGGCAGAGAGCGAGGCGATAAACGCTTCAGCTTTTTCTGGTGATTGCTGGCGGTAAATATTTGCCAGCCCGCGCACGGCGTTAGTGTTTCCGCTGTCCATACGCACGGTCTGCTGATAATAACGTTCGGCGGCGAGATAATCTTTTCGCGCCATCGCCACATCGCCCAGCCCCAGCACTGCATAACTGTAGGTGTTATCGACATTACGCGCATGCTGGAAAAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGATCg                              >  1:726453/1‑300 (MQ=255)                                                                                                                                                                                                                                        catcaATGCTACGCCGCTGACTGGCAGAGAGCGAGGCGATAAACGCTTCAGCTTTTTCTGGTGATTGCTGGCGGTAAATATTTGCCAGCCCGCGCACGGCGTTAGTGTTGCCGCTGTCCATACGCAAGGTCTGCTGATAATAACGTTCGGCGGCGGGATAATCTTTTCGCGCCATCGCCACATCGCCCAGCCCCAGCACTGCATAACTGTCGGTGTTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGATCg                              >  2:529151/1‑300 (MQ=255)                                                                                                                                                                                                                                                      cgcTGACTGGCAGAGAGCGAGGCGATAAACGCTTCAGCTTTTTCTGGTGATTGCTGGCGGTAAATATTTGCCAGCCCGCGCACGGCGTTAGTGTTGCCGCTGTCCATACGCAAGGTCTGCTGATAATAACGTTCGGCGGCGGGATAATCTTTTCGCGCCATCGCCACATCGCCCAGCCCCAGCACTGCATAACTGTCGGTGTTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGATCGCCAGCCAGTAGCgg                <  1:418581/300‑1 (MQ=255)                                                                                                                                                                                                                                                          gACTGGCAGAGAGCGAGGCGATAAACGCTTCAGCTTTTTCTGGTGATTGCTGGCGGTAAATATTTGCCAGCCCGCGCACGGCGTTAGTGTTGCCGCTGTCCATACGCAAGGTCTGCTGATAATAACGTTCGGCGGCGGGATAATCTTTTCGCGCCATCGCCACATCGCCCAGCCCCAGCACTGCATAACTGTCGGTGTTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTa            >  2:904231/1‑300 (MQ=255)                                                                                                                                                                                                                                                                  gagagCGAGGCGATAAACGCTTCAGCTTTTTCTGGTGATTGCTGGCGGTAAATATTTGCCAGCCCGCGCACGGCGTTAGTGTTGCCGCTGTCCATACGCAAGGTCTGCTGATAATAACGTTCGGCGGCGGGATAATCTTTTCGCGCCATCGCCACATCGCCCAGCCCCAGCACTGCATAACTGTCGGTGTTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTcagc    <  1:893212/300‑1 (MQ=255)                                                                                                                                                                                                                                                                    gagCGAGGCGATAAACGCTTCAGCTTTTTCTGGTGATTGCTGGCGGTAAATATTTGCCAGCCCGCGCACGGCGTTAGTGTTGCCGCTGTCCATACGCAAGGTCTGCTGATAATAACGTTCGGCGGCGTGATAATCTTTTCGCGCCATCGCCACATCGCCCAGCCCCAGCACTGCATAACTGTCGGTGTTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTcagcag  >  1:454326/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                           GCTCCGGGTCGTTCGACTTCTGCTGCGCCAGATTGCGCATTAACGTATCGGCCTGGCTGCGTTGTCCGGCCTGCCAGAGATCCTGCGAAAGTCGGTAAGTAATCCATACGCTGCCGGGGTCCAGCGCCAGTCGTTGCCGCTGAAGTGCTGCCGCCTGCGCCCATTTGCCCTGGTTTTCCAGTGCCTCTGCCTGCTGTGCCAGACGGTCGTTTTGCAGGCTGCGTTCGATATCATCAATGCTACGCCGCTGACTGGCAGAGAGCGAGGCGATAAACGCTTCAGCTTTTTCTGGCGATTGCTGGCGGTAAATATTTGCCAGCCCGCGCACGGCGTTAGTGTTGCCGCTGTCCATACGCAAGGTCTGCTGATAATAACGTTCGGCGGCGGGATAATCTTTTCGCGCCATCGCCACATCGCCCAGCCCCAGCACTGCATAACTGTCGGTGTTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTCAGCAG  >  NC_000913/3687583‑3688141

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