Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,686,705	G→A	F823F (TTC→TTT)	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,686,705	0	G	A	97.9%	138.1 / ‑3.7	47	F823F (TTC→TTT)	bcsC	cellulose synthase outer membrane channel
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (26/20);  total (26/21)
Fisher's exact test for biased strand distribution p-value = 4.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.16e-01

GCCCGATATCATCGCTGTAACTG‑ATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGTCATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTCATATTGACGAAATCACTGCGAAAGAACA‑TCCGCCCGTCAGAATACGGCGCATCCACCTGCAACATGGTAGTGTGCGCTTTCAGATCGGAGTAACCACCGGTGCCGCTCGAACCCCAGTAATCGTGCTCAAGGGTGACGTTAAGATCCTGCTGGCGATAGAGGTCCGCCGCATCGCTGCGCACGCCACGTTTCAGCCAGTCATCTTTCTCGTCGTTACGGGTCAGTCGGGTAAAGGTGT‑CGTTATCCTGCGGACGCGTCGT  >  NC_000913/3686458‑3686941                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                               gCCCGATATCACCGCTGTAA‑TGCAGGCCGCCGACCCCATCCGCCCCTTGGAACCCCATCGTCGTGGTACCGATATCCCAGCTCATGACGTCATTTCGCCATCCAACCGCCACGCTGGCTCCGGACTTTGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGTCATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTCATATTGACGAAATCACTGCGAAAAAAc                                                                                                                                                                                                                                             <  2:5576/251‑1 (MQ=255)   ccGATATCATCGCTGTAACTG‑ATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACATGTGCCCCAGTTGTCATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTCATATTGACGAAATCCCTGCGAAAAAACA‑t                                                                                                                                                                                                                                          >  1:340620/1‑251 (MQ=255)            ccGCTGTAATGG‑ATGCCGCCGACCCCATCCACCACGTTGAAGCCCACCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGTCATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTCATATTGACGAAATCACTGCGAAAAAACA‑TCCGCCCGt                                                                                                                                                                                                                                  <  1:259153/249‑1 (MQ=255)                   aaCTG‑ATGCCGCTGACCGCATCCACCACGTTTAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGTCATCCCATTTGCCAGCGGCATTAGTGGAGAAACTGCCGACGTTCATATTGACGAAATCACTGCGAAAAAACA‑TCCGCCCGTCAGAATa                                                                                                                                                                                                                           <  1:334406/250‑1 (MQ=255)                    aCTG‑ATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGTCATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTCATATTGACGAAATCACTGCGAAAAAACA‑TCCGCCCGTCAGAATACg                                                                                                                                                                                                                         >  1:133096/1‑251 (MQ=255)                                  cccccaTCCCCCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGTCATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTCATATTGACGAAATCACTGCGAAAAAACA‑TCCGCCCGTCAGAATACGGCGCATCCACCTg                                                                                                                                                                                                            <  1:220880/247‑1 (MQ=255)                                      caTCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGTCATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTCATATTGACGAAATCACTGCGAAAAAACA‑TCCGCCCGTCAGAATACGGCGCATCCACCTGCAAc                                                                                                                                                                                                        >  1:267472/1‑251 (MQ=255)                                          caccacGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGTCATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTCATATTGACGAAATCACTGCGAAAAAACA‑TCCGCCCGGCAGAATACGGCGCATCCACCTGCAACATgg                                                                                                                                                                                                    >  1:262541/1‑251 (MQ=255)                                                gTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGTCATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTCATATTGACGAAATCACTGCGAAAAAACA‑TCCGCCCGTCAGAATACGGCGACTCCACCTGCAACATgtgagtgg                                                                                                                                                                                              >  2:103084/1‑244 (MQ=255)                                                             ggCGTGGTACCGATATCCAGGCTCCATACTTCATTCCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGTCATCCCATTTGCCATCGGCATTAGTGTAGAAACTGCCGACGTTCATATTGACGAAATCACTGCGAAAAAACA‑TCCGCCCGTCAGAATACGGCGCATCCACCTGCAACATGGTAGTGTGCGCTTTCAGATc                                                                                                                                                                                 <  1:16189/251‑1 (MQ=255)                                                                    tACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGTCATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTCATATTGACGAAATCACTGCGAAAAAACA‑TCCGCCCGTCAGAATACGGCGCATCCACCTGCAACATGGTAGTGTGCGCTTTCAGATCGGAGTa                                                                                                                                                                           >  1:71409/1‑250 (MQ=255)                                                                            tCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGTCATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTCATATTGACGAAATCACTGCGAAAAAACA‑TCCGCCCGTCAGAATACGGCGCATCCACCTGCACCATGGTAGTGTGCGCTTTCAGATCGGAGTAACCACCggg                                                                                                                                                                  >  1:46924/1‑250 (MQ=255)                                                                            ccccAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGTCATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTCATATTGACGAAATCACTGCGAAAAAACA‑TCCGCCCGTCAGAATACGGCGCATCCACCTGCAACATGGTAGTGTGCGCTTTCAGATCGGAGTAACCACCGGt                                                                                                                                                                  <  1:214033/250‑1 (MQ=255)                                                                               cagcccatgccgttaTTTCGCCTGCCCACCGCCACGTTTGCACATTAATCCGACGTGCTGCGGGTGCCGCTACTGTCCTGTAATGTACAGGTGCCCCAGTTGTTATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTACATATTGACGAAATCACTGCGAAAGAACA‑CCCGGCCGTCAGAATACGGCGATTCCGCCTGCAACATGGTAGTGTGCGCTTTCAGATCTGAGTAACCACCGGTGcc                                                                                                                                                               <  2:294781/237‑1 (MQ=255)                                                                                  cacccgacttcagttgccaggCCGACCGCCACGCTGGCTCCGGAAGCCGACTTGCTGCTGTTGCCGCTACAGATCTTTAATCGAAAGGTTCCCCAGGTGTCATCCCATTTGCCATCGGCAGTAGTGGAGAAACTGCCGACGTTCATATTGACTAAACCACTGCGAAAACACA‑TCCGCCCGTCAGAATACGGCGCATCCCCCTGCAACATGGTAGTGTGCGCTGTCAGATCGTAGTAACCTCCGGTGCCGCt                                                                                                                                                            <  2:299982/231‑1 (MQ=255)                                                                                  cTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGTCATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTCATATTGACGAAATCACTGCGAAAAAACA‑TCCGCCCGTCAGAATACGGCGCATCCACCTGCAACATGGTAGTGTGCGCTTTCAGATCGGAGTAACCACCGGTGCCGCt                                                                                                                                                            >  1:83342/1‑251 (MQ=255)                                                                                   cccAGACGTCATTTCGCCAGCCAACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGTCATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTCATATTGACGAAATCACTGCGAAAAAACA‑TCCGCCCGTCAGAATACGGCGCATCCACCTGCAACATGGTAGTGTGCGCTTTCAGATCGGAGTAACCACCGGTGCCGCTc                                                                                                                                                           <  1:252853/250‑1 (MQ=255)                                                                                    ccAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGTCATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTCATATTGACGAAATCACTGCGAAAAAACA‑TCCGCCCGTCAGAATACGGCGCATCCACCTGCAACATGGTAGTGTGCGCTTTCAGATCGGAGTAACCACCGGTGCCGCTCg                                                                                                                                                          >  1:287388/1‑251 (MQ=255)                                                                                                   ccAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGTCATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTCATATTGACGAAATCACTGCGAAAAAACA‑TCCGCCCGTCAGAATACGGCGCATCCACCTGCAACATGGTAGTGTGCGCTTTCAGATCGGAGTAACCACCGGTGCCGCTCGAACCCCAGTAATCGt                                                                                                                                           >  1:164902/1‑251 (MQ=255)                                                                                                         gACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGTCATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTCATATTGACGAAATCACTGCGAAAAAACA‑TCCGCCCGTCAGAATACGGCGCATCCACCTGCAACATGGTAGTGTGCGCTTTCAGATCGGAGTAACCACCGGTGCAGCTCGAACCCCAGTAATCGTGCTCa                                                                                                                                      >  1:178445/1‑250 (MQ=255)                                                                                                                     ggCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTGGTCATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTCATATTGACGAAATCACTGCGAAAAAACA‑TCCGCCCGTCAGAATACGGCGCATCCACCTGCAACATGGTAGTGTGCGCTTTCAGATCGGAGTAACCACCGGTGCCGCGCGATCCCCAGTAATCGTGCTCAAGGGTGACGTTaa                                                                                                                         <  2:178445/251‑1 (MQ=255)                                                                                                                      gCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGTCATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTCATATTGACGAAATCACTGCGAAAAAACA‑TCCGCCCGTCAGAATACGGCGCATCCACCTGCAACATGGTAGTGTGCGCTTTCAGATCGGAGTAACCACCGGTGCCGCTCGAACCCCAGTAATCGTGCTCAAGGGTGACGTTa                                                                                                                          >  1:328782/1‑249 (MQ=255)                                                                                                                      gCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGTCATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTCATATTGACGAAATCACTGCGAAAAAACA‑TCCGCCCGTCAGAATACGGCGCATCCACCTGCAACATGGTAGTGTGCGCTTTCAGATCGGAGTAACCACCGGTGCCGCTCGAACCCCAGTAATCGTGCTCAAGGGTGACGTTa                                                                                                                          >  2:320684/1‑249 (MQ=255)                                                                                                                       cACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGTCATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTCATATTGACGAAATCACTGCGAAAAAACA‑TCCGCCCGTCAGAATACGGCGCATCCACCTGCAACATGGTAGTGTGCGCTTTCAGATCGGAGTAACCACCGGTGCCGCTCGAACCCCAGTAATCGTGCTCAAGGGTGACGTTAAGa                                                                                                                       <  1:154306/251‑1 (MQ=255)                                                                                                                           ggAATCCTACAGGTTGCTGTTGCCGCTACAGTTCGGTAATGTTCAGGTGCCCCAGTTGTCATTCCTGGTGCCATCGGCAATAGTGGTGAAAATGCCGACGTTCATATTGAAGAATTCACTGCGAAAAAACA‑TCCGCACTTCAGAAGACGGCGCGTCCACCTGCAACATGGTAGTGTGCGCTTTCAGATCGGAGTAACCACCGGTGCCGCTCGTACCCCAGTAATCTTGCTCAAGGGTGACGTTCAGAtgct                                                                                                                   <  2:273286/251‑4 (MQ=255)                                                                                                                           ggAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGTCATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTCATATTGACGAAATCACTGCGAAAAAACA‑TCCGCCCGTCAGAATACGGCGCATCCACCTGCAACATGGTAGTGTGCGCTTTCAGATCGGAGTAACCACCGGTGCCGCTCGAACCCCAGTAATCGTGCTCAAGGGTGACGTTAAGATCct                                                                                                                   <  1:118581/251‑1 (MQ=255)                                                                                                                                  gACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGTCATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTCATATTGACGAAATCACTGCGAAAAAACA‑TCCGCCCGTCAGAATACGGCGCATCCACCTGCAACATGGTAGTGTGCGCTTTCAGATCGGAGTAACCACCGGTGCCGCTCGAACCCCAGTAATCGTGCTCAAGGGTGACGTTAAGATCCTGCTGGCg                                                                                                            >  1:301035/1‑251 (MQ=255)                                                                                                                                           cGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGTCATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTCATATTGACGAAATCACTGCGAAAAAACA‑TCCGCCCGTCAGAATACGGCGCATCCACCTGCAACATGGTAGTGTGCGCTTTCAGATCGGAGTAACCACCGGTGCCGCTCGAACCCCCGTAATCGTGCTCAAGGGTGACGTTAAGATCCTGCTGGCGATAGAGGTc                                                                                                   >  1:4677/1‑251 (MQ=255)                                                                                                                                                               tAATGTACAGGTGCACCAGTTGTCATACCATTTGCCATCGGCATTAGTGGAGAAAATGCCGCCGTTCATATTGACGAAATCACTGCGAAAAAAAA‑TCCGACCGTCAGAATACGGCGCATCCACCTGCAACATGGTAGAGGCCGATTTCAGATCGGAGTAAACACCGGTGCCGCTCGAACCCCACTAATCGCGCTCCAGGGTGTCGTTAAGCTCCGGCTGTCGACAGAGGTCAACCGCATCGCAGCGCTCGc                                                                               >  2:308152/1‑251 (MQ=255)                                                                                                                                                                aaTGTACAGGCGCCCCAGTTGTCATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTCATATTGACGAAATCACTGCGAAAAAACA‑TCCGCCCGTCAGAATACGGCGCATCCACCTGCAACATGGTAGTGTGCGCTTTCAGATCGGAGTAACCACCGGTGCCGCTCGCACCCCAGTAGTCGTGCTCAAGGGTGCCGTTAAGATCCTGCTGGCGATAGAGGTACGCCGCATCGcgggtcacgcc                                                                              >  2:58330/1‑241 (MQ=255)                                                                                                                                                                        ggtgtcccatttcttaTCCCATTTGACACCGGCACTAGTGTAGAAACTGCAGTCCTTTAAAATGACAATATCCCTACGAAAAAACA‑TCCGTCCGTCAGAATACGGCGGTTCCCCCTGCAACATGGTTGTGTGCGCTTTCAGATCTGAGTAACCCCCGGGGCCGCCCGAACCCCAGTAATCGTGCTCGAGGGTGACGTTAAGATTCTGTTGGAGATATAGGTCCGCCGCATCGCTGCGCACGCCACGTTTCa                                                                      <  2:97454/236‑1 (MQ=255)                                                                                                                                                                           gCCCCAGTTGTCATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTCATATTGACGAAATCACTGCGAAAAAACA‑TCCGCCCGTCAGAATACGGCGCATCCACCTGCAACATGGTAGTGTGCGCTTTCAGATCGGAGTAACCACCGGTGCCGCTCGAACCCCAGTAATCGTGCTCAAGGGTGACGTTAAGATCCTGCTGGCGCTAGAGGTCCGCCGCATc                                                                                          >  2:213902/1‑228 (MQ=255)                                                                                                                                                                           gCCCCAGTTGTCATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTCATATTGACGAAATCACTGCGAAAAAACA‑TCCGCCCGTCAGAATACGGCGCATCCACCTGCAACATGGTAGTGTGCGCTTTCAGATCGGAGTAACCACCGGTGCCGCTCGAACCCCAGTAATCGTGCTCAAGGGTGACGTTAAGATCCTGCTGGCGATAGAGGTCCGCCGCATc                                                                                          <  1:213902/228‑1 (MQ=255)                                                                                                                                                                             cccAGTTGTCATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTCATATTGACGAAATCACTGCGAAAAAACA‑TCCGCCCGTCAGAATACGGCGCATCCACCTGCAACATGGTAGTGTGCGCTTTCAGATCGGAGTAACCACCGGTGCCGCTCGAACCCCAGTAATCGTGCTCAAGGGTGACGTTAAGATCCTGCTGGCGATAGAGGTCCGCCGCATCGCTGCGCACGCCACGTTTCAGcc                                                                   >  1:259796/1‑249 (MQ=255)                                                                                                                                                                                   tGTCATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTCATATTGACGAAATCACTGCGAAAAAACA‑TCCGCCCGTCAGAATACGGCGCATCCACCTGCAACATGGTAGTGTGCGCTTTCAGATCGGAGTAACCACCGGTGCCGCTCGAACCCCAGTAATCGTGCTCAAGGGTGACGTTAAGATCCTGCTGGCGATAGAGGTCCGCCGCATCGCTGCGCACGCCACGTTTCAGCCAGTCATc                                                            <  1:213190/250‑1 (MQ=255)                                                                                                                                                                                           cATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTCATATTGACGAAATCACTGCGAAAAAACA‑TCCGCCCGTCAGAATACGGCGCATCCACCTGCAACATGGTAGTGTGCGCTTTCAGATCGGAGTAACCACCGGTGCCGCTCGAACCCCAGTAATCGTGCTCAAGGGTGACGTTAAGATCCTGCTGGCGATAGAGGTCCGCCGCATCGCTGCGCACGCCACGTTTCAGCCAGTCATCTTTCtcgtc                                                   <  1:221975/251‑1 (MQ=255)                                                                                                                                                                                                         aTTAGTGGAGAAACTGCCGACGTTCATATTGACGAAATCACTGCGAAAAAACA‑TCCGCCCGTCAGAATACGGCGCATCCACCTGCAACATGGTAGTGTGCGCTTTCAGATCGGAGTAACCACCGGTGCCGCTCGAACCCCAGTAATCGTGCTCAAGGGTGACGTTAAGATCCTGCTGGCGATAGAGGTCCGCCGCATCGCTGCGCACGCCACGTTTCAGCCAGTAATCTTTCTCGTCGTTACGGGTCAGtg                                     >  2:286539/1‑250 (MQ=255)                                                                                                                                                                                                                   aaaatGCCGACGTTCATATAGACGAAATCACTGCGAAAAAACA‑TCCGCACGTCAGAATACGGCGCATACACCTACAACATGGTAGTGAGCGCTTTCAGATCGGAGTAACCACCGGTGACGCTCGAACCCCAGTAATCGTGCTCAAGGGGGCCGCTAAGATCATGCTGGCGAAAGAGGTCAGCCGCATCGCTGCGCACGCCACGTTTAGGCCAGTCATCTTTCTCGTCGTTCCGGGACCGTCGGGCAAAgg                            >  2:286559/5‑250 (MQ=255)                                                                                                                                                                                                                        gCCGACGTTCATATTGACGAAATCACTGCGAAAAAACA‑TCCGCCCGTCAGAATACGGCGCATCCACCTGCAACATGGTAGTGTGCGCTTTCAGATCGGAGTAACCACCGGTGCCGCTCGAACCCCAGTAATCGTGCTCAAGGGTGACGTTAAGATCCTGCTGGCGATAGAGGTCCGCCGCATCGCTGCGCACGCCACGTTTCAGCCAGTCATCTTTCTCGTCGTTACGGGTCAGTCGGGTAAAGGTGT‑CGt                     <  1:103084/251‑1 (MQ=255)                                                                                                                                                                                                                         ccGACGTTCATATTGACGAAATCACTGCGAAAAAACA‑TCCGCCCGTCAGAATACGGCGCATCCACCTGCAACATGGTAGTGTGCGCTTTCAGATCGGCGTAACCACCGGTGCCGCTCGAACCCCCGTAATCGTGCTCAAGGGTGACGTTAAGATCCTGCGGGCGCTAGAGGTCCGCCGCATCGCTGCGCACGCCACGTTTCAGCCAGTCATCTTTCTCGTCGTTACGGGGCAGTCGGGTAAAGGTGT‑CGtt                    >  2:246766/1‑251 (MQ=255)                                                                                                                                                                                                                         ccGACGTTCATATTGACGAAATCACTGCGAAAAAACA‑TCCGCCCGTCAGAATACGGCGCATCCACCTGCAACATGGTAGTGTGCGCTTTCAGATCGGAGTAACCACCGGTGCCGCTCGAACCCCAGTAATCGTGCTCAAGGGTGACGTTAAGATCCTGCTGGCGATAGAGGTCCGCCGCATCGCTGCGCACGCCACGTTTCAGCCAGTCATCTTTCTCGTCGTTACGGGTCAGTCGGGTAAAGGTGT‑CGtt                    <  1:286539/251‑1 (MQ=255)                                                                                                                                                                                                                            cccctaatatTAACGAAATCCCTGCGCAAAAACACTCC‑CCCGTCATAATACGGCGCACCCACCCGCAACATGGTAGTGGGCTCTTTCAGATCGGAGTACCCACCGGTGCCGCTCGAACCCCATTAAGCGTTCTCAAGGGTGCGGTTCAGATCCTGCTGGCGATAGAGGTCCGCCGCATCGCTGCGCCCGCCACGGTTCAGCCAGTCATCTTTCTCTTCGTTACGGGTCAGTCGGGTAAAGGTGT‑CGttctc                 <  2:340620/245‑4 (MQ=255)                                                                                                                                                                                                                             cgtgcATATTGACGAAATCACTGCGAAAAAACA‑TCCGCCCGTCAGAATACGGCGCATCCACCTGCAACATGGTAGTGTGCGCTTTCAGAGCGGAGTAACCCCCGGTGCCGCTCGAACCCCAGTAATc                                                                                                                                             <  2:199918/123‑1 (MQ=255)                                                                                                                                                                                                                             cGTTCATATTGACGAAATCACTGCGAAAAAACA‑TCCGCCCGTCAGAATACGGCGCATCCACCTGCAACATGGTAGTGTGCGCTTTCAGATCGGAGTAACCACCGGTGCCGCTCGAACCCCAGTAATc                                                                                                                                             >  1:199918/1‑127 (MQ=255)                                                                                                                                                                                                                                 cATATTGACGAAATCACTGCGAAAAAACC‑TCCACCCGTCAGAATACGGCGCATCCACCTGCAACATGGTAGTGTGCGCTTTCAGATCGGAGTAACCACCGGTGCCGCTCGAACCCCAGTAATCGTGCTCAAGGGCGACGTTAAGAGCCTACTGGCGATAGAGGTCCGCCGCATCGCTGCGCACGCCACGTTTCAGCCAGTCATCTTTTTCGTCGTTACGGGTCAGTCGGGAAAAGG‑GTGCTTTATCCTGCg            >  2:261347/1‑251 (MQ=255)                                                                                                                                                                                                                                 cATATTGACGAAATCACTGCGAAAAAACA‑TCCGCCCGTCAGAATACGGCGCATCCACCTGCAACAAGGTAGTGTGCGCTTTCAGATAGGAGTAACCACCGGCGCCGATCGAACCCCAATAATCGTGCTCAAGGGTGACGTTAAGATCCTGCTGGCGATAGAGGTCCGCCGAATCGCTGCGCACGCCACGTTTCAGCCAGACATCTTTCTCGTGATTGCGGGTCAGTCGGGTAACGGGGT‑CATGATCCTggg            >  1:153343/1‑249 (MQ=255)                                                                                                                                                                                                                                           aaaTCACTGCGAAAAAACA‑TCCGCCCGTCAGAATACGGCGCATCCACCTGCAACATGGTAGTGTGCGCTTTCAGATCGGAGTAACCACCGGTGCCGCTCGAACCCCAGTAATCGTGCTCAAGGGTGACGTTAAGATCCTGCTGGCGATAGAGGTCCGCCGCATCGCTGCGCACGCCACGTTTCAGCCAGTCATCTTTCTCGTCGTTACGGGTCAGTCGGGTAAAGGTGT‑CGTTATCCTGCGGACGcgtcgt  >  1:39528/1‑251 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                               GCCCGATATCATCGCTGTAACTG‑ATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGTCATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTCATATTGACGAAATCACTGCGAAAGAACA‑TCCGCCCGTCAGAATACGGCGCATCCACCTGCAACATGGTAGTGTGCGCTTTCAGATCGGAGTAACCACCGGTGCCGCTCGAACCCCAGTAATCGTGCTCAAGGGTGACGTTAAGATCCTGCTGGCGATAGAGGTCCGCCGCATCGCTGCGCACGCCACGTTTCAGCCAGTCATCTTTCTCGTCGTTACGGGTCAGTCGGGTAAAGGTGT‑CGTTATCCTGCGGACGCGTCGT  >  NC_000913/3686458‑3686941

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