Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,698,243	G→A	M10I (ATG→ATA)	bcsG →	cellulose phosphoethanolamine transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,698,243	0	G	A	97.4%	113.5 / NA	39	M10I (ATG→ATA)	bcsG	cellulose phosphoethanolamine transferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (27/11);  minor base C (1/0);  total (28/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAATTATTGTCGTTTGCGCACTGATATTTTTCCCGCTGGGCTATCTGGCGCGGCACTCTTTGCGACGCATTCGCGACACCTTACGTTTGTTCTTTGCTAAACCTCGTTATGTTAAACC‑GG‑CCGGGACGTTAC‑GCCGCACGGAAAAAGCCAGGGCAA‑CCAAAAAATGACTCAATTTACGCAAAAT‑ACCGCCATGCCTTCTTCCCTCTGGCAATACTGGCGCGGCCTTTCCGGCTGGAACTTCTATTTTCTGGTTAAGTTCGGCCTGTTGTGGGCGGGATATCTTAACTTCCATCCGCTCCTCAATTTGGTGTTTGCCGCGTTTCTGCTGATGCCCCTTCCGCGCTACAGCCTGCATCGCTTGCGCCACTGGATTGCCCTGCCGATCGGCTTTGCTTTGTTCTGGCATGACACCTGGTTGCCTGGCCCGGAAAGCATA  >  NC_000913/3697997‑3698486                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                    cTGTTAGATGATGCTGTGGAGCACTCATCATGATGACCATCAGCGATATCATTGAAATTATTGTCGTTTGAGCACTGATATTTTTCCCGCTGGGCTATCTGGCGCGGCACTCTTTACGCCGCATTCGCGACACCTTACGTTTGTTCTTTGCTAAACCTCGTTATGTTAAACC‑GG‑CCGGGACGTTAC‑GCCGAAAGGAAAAAGCCAGGGCAA‑CCAAAAAATGACTCAACTTACGCAAAAT‑ACCGCCATACctt                                                                                                                                                                                                                                                 >  1:42156/1‑251 (MQ=255)           atgCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAATTATTGTCGTTTGCGCACTGATATTTTTCCCGCTGGGCTATCTGGCGCGGCACTCTTTGCGACGCATTCGCGACACCTTACGTTTGTTCTTTGCTAAACCTCGTTATGTTAAACC‑GG‑CCGGGACGTTAC‑GCCGCACGGAAAAAGCCAGGGCAA‑CCAAAAAATGACTCAATTTACGCAAAAT‑ACCGCCATACCTTCTTCCCTCTg                                                                                                                                                                                                                                       >  2:202997/1‑251 (MQ=255)                gtgGAGCTCTAATCATTATAACAATCAGCGATATAATTGAATTTTTTGTCGTTTGCGCACTGATATTTTTCCCGCTGGGCTATGTGGCGCGGCACTCTTCGCGCCGCATGCGCGACCCCTTATGTTTGTTCTTTGCAAACCCGCGTTATTTAAAACC‑GG‑CCGGGAC‑TTACCGCCGCACGAAAAAACCCAGGGCAA‑CAAAAAAATGACTCATTTTCCGAAAAAT‑ACCGCCATACCTTCTTCCCTCTGGCAAt                                                                                                                                                                                                                                  <  2:81875/251‑1 (MQ=255)                            taatgatgACCATCAGCGATATCATTGAAATTATTGTCGTTTGCTCACTGATATTTTTCCCGCTGGGCTATCTGGCGCGGCACTCTTTGCGACGCATTCGCGACACCTTACGTTTGTTCTTTGCTAAACCTCGTTATGTTAAACA‑GG‑CCGGGACGTTAC‑GCCGCACGGAAAAAGCCAGGGCAA‑CCAAAAAATGACTCAATTTACGCACAAT‑ACCGCCATACCTTCTTCCCTCTGGCAATACTGGCGCGGcc                                                                                                                                                                                                                      >  2:38165/3‑251 (MQ=255)                                          aGCGATATAATTGAAATTATTGACGTTTGCGCACTGATATTTTTACCGCTTGGCTATCTGGCGCGGCACTCTTTGCGACGCATTCGCGACACCTTACGTTTGTTCTTTGCTAAACCTCGCTCTGTTAAACC‑GG‑CCGGGACGTTAC‑GCCGCACGGAAAAAGCCAGGGCAA‑CCAAAGAATGACTCAATTTACGCAAAAT‑CCCGCCATCCCTTCTTCCCTCTGGCAATACTGGCGAGGCCTTTCCGGCTGGAAc                                                                                                                                                                                                        >  1:370715/1‑251 (MQ=255)                                                              tGTCGTTTGCGCACTGATATTTTTCCCGCTGTGCTATCTGGCGCGGCACTCTTTGCGACGCATTCGCGCCACCTTACGTTTGTTCTTTGCTAAACCTCGTTATGTTAAACC‑GG‑CCGGGACGTTAC‑GCCGCACGGAAAAAGCCAGGGCAA‑CCAAAAAATGACTCAATTTACGCAAAAT‑ACCGCCATACCTTCTTCCCTCTGGCACTACTGGCGCGGCCTTTCCGGCTGGCCCTTCTTTTTTCTGGTTAAGtt                                                                                                                                                                                    >  2:44088/1‑251 (MQ=255)                                                                   tttGCCCACTGATTTTTTCCCCGCTGGGCTATCTGGCGCGGCCCTCTTTGCGCCGCATTCGCGCCCCCTTACTTTTTTTCTTTGCAAACCCTCGTTATGTTAACCC‑GG‑CCGGGACGTTAC‑GCCGCACGGAAAAAGCCAGGCAAA‑CAAAAAAAGGACTAAATTTACAAAAAAT‑ACCGCCATACCTTCTTCCCTCTGGCAATACTGGCGCGGCCTTTCCGGCTGGAACTTCTATTTTCTGGTTAAGTTCgg                                                                                                                                                                                 <  1:35808/249‑1 (MQ=255)                                                                        gcACTGATATTTTTCCCGCTGGGCTATCTGGCGCGGCACTCTTTGCGACGCATTCGCGACACCTTACGTTTGTTCTTTGCTAAACCTCGTTATGTTAAACC‑GG‑CCGGGACGTTAC‑GCCGCACGGAAAAAGCCAGGGCAA‑CCAAAAAATGACTCAATTTACGCAAAAT‑ACCGCCATACCTTCTTCCCTCTGGCAATACTGGCGCGGCCTTTCCGGCTGGAACTTCTATTTTCTGGTTAAGTTCGGCCtgttg                                                                                                                                                                          >  2:395738/1‑251 (MQ=255)                                                                          cctatttttttcCCCGCGGGGCTATCGGGCGCGCCTCTCTTTGCGCCCCATCCGCGACCCCTTCCGTTTGTTTTTTTCAAACCCTCTTTATGTAAACCC‑GG‑CGGGGACTTTAC‑GCCGCACGAAAAAAGCCGGGTCAA‑CAAAAAAATGACTAAATTTACGAAAAAA‑ACCGCCATACCTTCTTCCCTCTGGCAATACTGGCGCGGCCTTTCCGGCTGGAACTTCTATTTTCTGGTTAAGTTCGGCCTGTtgtg                                                                                                                                                                        <  1:395738/244‑1 (MQ=255)                                                                                   tttCCCGCTGGGCTATCTGGCGCGGCACTCTTTGCGACGCATTCGCGACACCTTACGTTTGTTCTTTGCTAAACCTCGTTATGTTAAACC‑GG‑CCGGGACGTTAC‑GCCGCACGGAAAAAGCCAGGGCAA‑CCAAAAAATGACTCAATTTACGCAAAAT‑ACCGCCATACCTTCTTCCCTCTGGCAATACTGGCGCGGCCTTTCCGGCTGGAACTTCTATTTTCTGGTTCAGTTCGGCCTGTTGTGGGCggggtc                                                                                                                                                               >  2:188253/1‑248 (MQ=255)                                                                                     tCCCGCTGGGCTATCTGGCGCGGCACTCTTTGCGACGCATTCGCGACACCTTACGTTTGTTCTTTGCTAAACCTCGTTATGTTAAACC‑GG‑CCGGGACGTTAC‑GCCGCACGGAAAAAGCCAGGGCAA‑CCAAAAAATGACTCAATTTACGCAAAAT‑ACCGCCATACCTTCTTCCCTCTGGCAATACTGGCGCGGCCTTTCCGGCTGGAACTTCTATTTTCTGGTTAAGTTCGGCCTGTTGTGGGCGGGATATc                                                                                                                                                             >  2:136724/1‑251 (MQ=255)                                                                                           tGGGCTATCTGGCGCGGCACTCCTTGCGACGCATTCGCGACACCTTGCGTTTGTTCTTTGCTAAACCTCGTTATGTTAAACC‑GG‑CCGGGACGTTAC‑GCCGCACGGAAAAAGCCAGGGCAA‑CCAAAAAATGACTCAATTTACGCACAAT‑ACCGCCATACCTTCTTCCCTCTGGCAATACTGGCGCGGCCTTTCCGGCTGGAACTTCTATTTTCTGGTTAAGTTCGGCCTGTTGTGGGCGGGATATCTTAACt                                                                                                                                                       >  1:379974/1‑251 (MQ=255)                                                                                                                    gcgcctcATTCGCGACTCTTTACTTTTGTTCTTGTCAAACCCTCGTATTGTGAAACC‑GG‑CGGGGACGTCAC‑CCCGTAAGAAAAAAGCCAGTGCAC‑CAAAAAAATGACACATTTTACAAAAAAC‑TCCGCCATACCTTCTTCCCTCTGGCAATACTGGCGCGGCCTTTCCGGCTGGAACTTCTATTTTCTGGTTAAGTTCGGCCTGTTGTGGGCGGGATATCTTAACTTCCATCCGCTCCTCAATTTGGTGtt                                                                                                                              <  1:219059/245‑1 (MQ=255)                                                                                                                                 gACACCTTACGTTTGTTCTTTGCTAAACCTCGTTATGTTAACCC‑GG‑CCGGGACGTTAC‑GCCGCACGGAAAAAGCCAGGGCAA‑CCAAAAAATGACTCAATTTACGCAAAAT‑ACCGCCATACCTTCTTCCCTCTGGCAATAATGGCGCGGCCTTTCCGGCTGGAACTTCTCTTTTCTGGTTAAGTTCGGCCTGTTGTGGGCGGGATATCTTAACTTCCATCCGATCCTCAATTTGGTGTTTGCCGCGTTTctt                                                                                                                 >  2:280585/1‑250 (MQ=255)                                                                                                                                   cacCTTACGTTTGTTCTTTGCTAAACCTCGTTATGTTAAACC‑GG‑CCGGGACGTTAC‑GCCGCACGGAAAAAGCCAGGGCAA‑CCAAAAAATGACTCAATTTACGCAAAAT‑ACCGCCATACCTTCTTCCCTCTGGCAATACTGGCGCGGCCTTTCCGGCTGGAACTTCTATTTTCTGGTTAAGTTCGGCCTGTTGTGGGCGGGATATCTTAACTTCCATCCGCTCCTCAATTTGGTGTTTGCCGCGTTTctgct                                                                                                               >  2:161728/1‑251 (MQ=255)                                                                                                                                    acCTTACGTTTGTTCTTTGCTAAACCTCGTTATGTTAAACC‑GG‑CCGGGACGTTAC‑GCCGCACGGAAAAAGCCAGGGCAA‑CCAAAAAATGACTCAATTTACGCAAAAT‑ACCGCCATACCTTCTTCCCTCTGGCAATACTGGCGCGGCCTTTCCGGCTGGAACTTCTATTTTCTGGTTAAGTTCGGCCTGTTGTGGGCGGGATATCTTAACTTCCATCCGCTCCTCAATTTGGTGTTTGCCGCGTTTctgctg                                                                                                              >  1:205579/1‑251 (MQ=255)                                                                                                                                     ccTTACGTTTGTTTTTTCCAAAACCTGGTTATGTTAACCC‑GG‑CCGGGCCGTTCC‑GCCGCACGGAAAAAGCCAGGGCAC‑CAAAAAAAGGACTAATTTTCCGCAAAAT‑CCCGCCATACCTTCTTCCCTCTGGCAATACTGGCGCGGCCTTTCCGGCTGGAACTTCTATTTTCTGGTTAAGTTCGGCCTGTTGTGGGCGGGATATCTTAACTTCCATCCGCTCCTCAATTTGGTGTTTGCCGCGTTTCTGCTGa                                                                                                             <  1:202997/251‑1 (MQ=255)                                                                                                                                       ttACGTTTGTTCTTTGCTAAACCTCGTTATGTTAAACC‑GG‑CCGGGACGTTAC‑GCCGCACGGAAAAAGCCAGGGCAA‑CCAAAAAATGACTCAATTTACGCAAAAT‑ACCGCCATACCTTCTTCCCTCTGGCAATACTGGCGCGGCCTTTCCGGCTGGCACTTCTATTTTCTGGTTAAGTTCGGCCTGTTGTGGGCGGGATATCTTAACTTCCATCCGCTCCTCTATTTGGTGTTTGCCGCGTTTCTGCTGATg                                                                                                           >  2:110074/1‑251 (MQ=255)                                                                                                                                        ttgttttgttttgtgcaaaccctgtttATGT‑AAACCAGG‑CGGGGACGTTGC‑CCCGCACGAAAAAAGCAAGGGAAA‑CAAAAAAATGACACAATTACCGAAAAAT‑ACCGCCATACCTTCTTCCCTCTGGCAATACTGGCGCGGCCTTTCCGGCTGGAACTTCTATTTTCTGGTTAAGTTCGGCCTGTTGTGGGCGGGATATCTTAACTTCCATCCGCTCCTCAATTTGGTGTTTGCCGCGTTTCTGCTGATg                                                                                                           <  2:108163/225‑1 (MQ=255)                                                                                                                                         gggggtgctctcccgtcccccTCGTT‑TGGTAAACC‑GGGCGGGGAGGTTAC‑GCCGAACGAAAAAAGCCAGGACAA‑CAAAAAAAGGACTAAATTACCGCAAAAG‑ACCGCCATACCTTCTTCCCTCTGGCAATACTGGCGCGGCCTTTCCGGCTGGAACTTCTATTTTCTGGTTAAGTTCGGCCTGTTGTGGGCGGGATATCTTAACTTCCATCCGCTCCTCAATTTGGTGTTTGCCGCGTTTCTGCTGATGcc                                                                                                         <  1:254899/232‑1 (MQ=255)                                                                                                                                                           aCCTCGTTATGTTAAACC‑GG‑CCGGGACGTTAC‑GCCGCACGGAAAAAGCCAGGGCAA‑CCAAAAAATGACTCAATTTACGCAAAAT‑ACCGCCATACCTTCTTCCCTCTGGCAATACTGGCGCGGCCTTTCCGGCTGGAACTTCTATTTTCTGGTTAAGTTCGGCCTGTTGTGGGCGGGATATCTTAACTTCCATCCGCTCCTCAATTTGGTGTTTGCCGCGTTTCTGCTGATGCCCCTTCCGCGCTACAGCCt                                                                                       >  2:414112/1‑251 (MQ=255)                                                                                                                                                            ccTCGTTATGTTAAACC‑GG‑CCGGGCCGTTAC‑GCCGCACGGAAAAAGCCAGGGCAA‑CCAAAAAATGCCTCAATTTACGCAAAAT‑ACCGCCATACCTTCTTCCCTCTGGCAATACTGGCGCGGCCTTTCCGGCTGGAACTTCTATTTTCTGGTTAAGTTCGGCCTGTTGTGGGCGGGATATCTTAACTTCCATCCGCTCCTCAATTTGGTGTTTGCCGCGTTTCTGCTGATGcc                                                                                                         >  2:254899/1‑232 (MQ=255)                                                                                                                                                               cGTTATGTTAAACC‑GG‑CCGGGACGTTAC‑GCCGCACGGAAAAAGCCAGGGCAA‑CCAAAAAATGACTCAATTTACGCAAAAT‑ACCGCCATACCTTCTTCCCTCTGGCAATACTGGCGCGGCCTTTCCGGCTGGAACTTCTATTTTCTGGTTAAGTTCGGCCTGTTGTGGGCGGGATATCTTAACTTCCATCCGCTCCTCAATTTGGTGTTTGCCGCGTTTCTGCTGATGCCCCTTCCGCGCTACAGCCTGCAt                                                                                   >  2:137409/1‑251 (MQ=255)                                                                                                                                                               cGTTATGTTAAACC‑GG‑CCGGGACGTTAC‑GCCGCACGGAAAAAGCCAGGGCAA‑CCAAAAAATGACTCAATTTACGCAAAAT‑ACCGCCATACCTTCTTCCCTCTGGCAATACTGGCGCGGCCTTTCCGGCTGGAACTTCTATTTTCTGGTTAAGTTCGGCCTGTTGTGGGCGGGATATCTTAACTTCCATCCGCTCCTCAATTTGGTGTTTGCCGCGTTTCTGCTGATGCCCCTTCCGCGCTACAGCCTGCAt                                                                                   >  1:316353/1‑251 (MQ=255)                                                                                                                                                                     gTTAAACC‑GG‑CCGGGACGTTAC‑GCCGCACGGAAAAAGCCAGGGCAA‑CCAAAAAATGACTCAATTTACGCAAAAT‑ACCGCCATACCTTCTTCCCTCTGGCAATACTGGCGCGGCCTTTCCGGCTGGAACt                                                                                                                                                                                                       >  2:263065/1‑129 (MQ=255)                                                                                                                                                                                     gACGTTAC‑GCCGCACGGAAAAAGCCAGGGCAA‑CCAAAAAATGACTCAATTTCCGCAAAAT‑ACCGCCATACCTTCTTCCCTCTGGCAATACTGGCGCGGCCTTTCCGGCTGGAACTTCTATTTTCTGGTTAAGTTCGGCCTGTTGTGGGCGGGATATCTTAACTTCCATCCGCTCCTCAATTTGGTGTTTGCCGCGTTTCTGCTGATGCCCCTTCCGCGCTACAGCATGCATCGCTTGCGCCACTGGATTGc                                                               >  2:263306/1‑251 (MQ=255)                                                                                                                                                                                          gtatgccgacaggAAAAAGCCCGGGCAAACAAAAAAATGACTAAATTT‑CCCAAAATTCCCGCCATACCTTCTTCCCTCTGGCAATACTGGCGCGGCCTTTCCGGCTGGAACTTCTATTTTCTGGTTAAGTTCGGCCTGTTGTGGGCGGGATATCTTAACTTCCATCCGCTCCTCAATTTGGTGTTTGCCGCGTTTCTGCTGATGCCCCTTCCGCGCTACAGCCTGCATCGCTTGCGCCACGGGATTGCCCt                                                            <  2:64399/240‑1 (MQ=255)                                                                                                                                                                                                 gcACGGAAAACGCCAGGGCAA‑CCAAAAAATGACTCAATTTACGCAACAT‑ACCGCCATACCTTCTTCCCTCTGGCAATACTGGCGCGGCCTTTCCGGCTGGAACTTCTATTTTCTGGTTAAGTTCGGCCTGTTGTGGGCGGGATATCTTAACTTCCATCCGCTCCTCAATTTGGTGTTTGCCGCGTTTCTGCTGATGCCCCTTCCGCGCTACAGCCTGCATCGCTTGCGCCACTGGATTGCCCTTCCGATCg                                                    >  2:408971/1‑251 (MQ=255)                                                                                                                                                                                                        aaaaGCCAGGGCAA‑CCAAAAAATGAATCAATTTACGCAAAAT‑ACCGCCATACCTTCTTCCCTATGGCAATACTGGCGCGGCCTTTCCGGCTGGAACTTCTATTTTCTGGTTAAGTTCGGCCTGTTGTGGGCGGGATATCTTAACTTCCATCCGCTCCTCAATTTGGTGTTTGCCGCGTTTCTGCTGATGCCCCTTCCGCGCTACAGCCTGCATCGCTTGCGCCACTGGCTTGCCCTGCCGCTCGgctttgc                                             >  2:404347/1‑251 (MQ=255)                                                                                                                                                                                                         aaaGCCAGGCAAA‑CAAAAAAATCACCCATTTCACGAAAAAT‑ACCGCCATACCTTCTTCCCTCTGGCAATACTGGCGCGGCCTTTCCGGCTGGAACTTCTATTTTCTGGTTAAGTTCGGCCTGTTGTGGGCGGGATATCTTAACTTCCATCCGCTCCTCAATTTGGTGTTTGCCGCGTTTCTGCTGATTCCCCTTCCGCGCTACAGCCTGCATCGCTTGCGCCACTGGATTGCCCTGCCGATCGgctttgct                                            <  2:168778/251‑1 (MQ=255)                                                                                                                                                                                                            gCCAGGGCAA‑CCAACGAATGACTCAATTTACGCAAAAT‑ACCGCCATACCTTCTTCCCTCTGGCAATACTGGCGCAGCCTTTACGGCTGGAACTTCTATTTTCTGGATCAGTTCGGCCTGTTGTGGGCGGGATATCTTAACTTCCATCCGCTCCTCAATTTGGTGTTTGCCGCGTTTCTGCTGATGCACCTTCCGCGCTTCAGCCTGCCTCGCTTGCGCCACTGGATTGCCCTGCCGATCGGctttgctttg                                         >  1:282898/1‑251 (MQ=255)                                                                                                                                                                                                            gCCAGGGCAA‑CCAAAAAATGACTCAATTTACGCAAAAT‑ACCGCCCTACCTTCTTCCCTCTGGCAATACTGGCGCGGGATTTCAGGCTGGCACTTCtctt                                                                                                                                                                                                 >  1:419887/1‑96 (MQ=255)                                                                                                                                                                                                                     a‑CCAAAAAATATCTCATTTTCCGAAAAAA‑ACCGCCATACCTTCTTCCCTCTGGCAATTCTGGCGCGGCCTTTCCGGCTGGAACTTCTATTTTCTGGTTAAGTTCGGCCTGTTGTGGGCGGGATATCTTAACTTCCATCCGCTCCTCAATTTGGTGTTTGCCGCGTTTCTGCTGATGCCCCTTCCGCGCTACAGCCTGCATCGCTTGCGCCACTGGATTGCCCTGCCGATCGGCTTTGCTTTGTTCTGGCAt                                <  1:302914/251‑1 (MQ=255)                                                                                                                                                                                                                       ccAAAAAATGACTCAATTTACGCAAAAT‑ACCGCCATACCTTCTTCCCTCTGGCAATACTGGCGCGGCCTTTCCGGCTGGAACTTCTATTTTCTGGTTAAGTTCGGCCTGTTGTGGGCGGGATATCTTAACTTCCATCCGCTCCTCAATTTGGTGTTTGCCGCGTTTCTGCTGATGCCCCTTCCGCGCTACAGCCTGCATCGCTTGCGCCACTGGATTGCCCTGCCGCTCGGCTTTGCTTTGTTCTGGCATg                               >  2:107636/1‑251 (MQ=255)                                                                                                                                                                                                                               tGACTCAATTTACGCAAAAT‑ACCGCCATACCTTCTTCCCTCTGGCAATACTGGCGCGGCCTTTCCGGCTGGAACTTCTATTTTCTGGTTAAGTTCGGCCTGTTGTGGGCGGGATATCTTAACTTCCATCCGCTCCTCAATTTGGTGTTTGCCGCGTTTCTGCTGATGCCCCTTCCGCGCTACAGCCTGCCTCGCTTGCGCCACTGGATTGCCCTGCCGATCGGCTTTGCTTTGTTCTGGCATGACACCTgg                       >  2:404767/1‑251 (MQ=255)                                                                                                                                                                                                                                       tttACGCACACT‑ACAGCCATACCTTCTTCCCTCTGGCAATACTGGCGCGGCCTTACCGGCCGGAACTTCTATTTTCTGGTTAAGTTCGGCCTGTTGTGGGCGGGATATCTTAACTTCCATCCGCTCCTCAATTTGGTGTTTGCCGCGTTTCTGCTGATGCCCCTTCCGCGCTCCAGCCTGCCTCGCTTGCGCCACTGGATTGCCCTGCCGATCGGCTTTGCTTTGTTCTGGCATGCCACCTGGTTGCCTgg               >  2:412837/1‑251 (MQ=255)                                                                                                                                                                                                                                              aaaaT‑ACCGCCATACCTTCTTCCCTCTGGCAATACTGGCGCGGCCTTTCCGGCTGGAACTTCTATTTTCTGGTTAAGTTCGGCCTGTTGTGGGCGGGATATCTTAACTTCCATCCGCTCCTCAATTTGGTGTTTGCCGCGTTTCTGCTGATGCCCCTTCCGCGCTACAGCCTGCATCGCTTGCGCCACTGGATTGCCCTGCCGATCGGCTTTGCTTTGTTCTGGCATGACACCTGGTTGCCTGGCCCGGa         >  2:330133/1‑250 (MQ=255)                                                                                                                                                                                                                                                 aattccGCCATACCTTCTTCCCTCTGGCAATACTGGCGCGGCCTTTCCGGCTGGAACTTCTATTTTCTGGTTAAGTTCGGCCTGTTGTGGGCGGGATATCTTAACTTCCATCCGCTCCTCAATTTGGTGTTTGCCGCGTTTCTGCTGATGCCCCTTCCGCGCTACAGCCTGCATCGCTTGCGCCACTGGATTGCCCTGCCGATCGGCTTTGCTTTGTTCTGGCATGACACCTGGTTGCCTGGCCCGGAAAg      <  1:73527/247‑1 (MQ=255)                                                                                                                                                                                                                                                      cGCCATACCTTCTTCCCTCTGGCACTACTGGCGCGGCCTTTCCGGCTGGAACTTCTCTTTTCTGGTTAAGTTCGGCCTGTTGTGGGCGGGATATCTTAACTTCCATCCGCTCCTCAATTTGGTGTTTGCCGCGTTTCTGCTGATGCCCCTTCCGCGCTACAGCCTGCATCGCTTGCGCCACTGGCTTGCCCTGCCGATCGGCTTTGCTTTGTTCTGGCATGACACCTGGTTGCCTGGCCCGGAAAGCATa  >  2:205226/1‑250 (MQ=255)                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                    CTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATCAGCGATATCATTGAAATTATTGTCGTTTGCGCACTGATATTTTTCCCGCTGGGCTATCTGGCGCGGCACTCTTTGCGACGCATTCGCGACACCTTACGTTTGTTCTTTGCTAAACCTCGTTATGTTAAACC‑GG‑CCGGGACGTTAC‑GCCGCACGGAAAAAGCCAGGGCAA‑CCAAAAAATGACTCAATTTACGCAAAAT‑ACCGCCATGCCTTCTTCCCTCTGGCAATACTGGCGCGGCCTTTCCGGCTGGAACTTCTATTTTCTGGTTAAGTTCGGCCTGTTGTGGGCGGGATATCTTAACTTCCATCCGCTCCTCAATTTGGTGTTTGCCGCGTTTCTGCTGATGCCCCTTCCGCGCTACAGCCTGCATCGCTTGCGCCACTGGATTGCCCTGCCGATCGGCTTTGCTTTGTTCTGGCATGACACCTGGTTGCCTGGCCCGGAAAGCATA  >  NC_000913/3697997‑3698486

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