Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	66,266	A→G	H95H (CAT→CAC)	araD ←	L‑ribulose‑5‑phosphate 4‑epimerase AraD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	66,266	0	A	G	100.0%	120.2 / NA	37	H95H (CAT→CAC)	araD	L‑ribulose‑5‑phosphate 4‑epimerase AraD
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (17/20);  total (17/20)

CCCCATGCAAACGGGCCGTGGGAATGGACCAGAACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCAAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCGCCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCG‑TGGGTGGTGCCGGTTGCTGGAATCGACTGACCCGCCTGCGCCCAGATGGTGGCGTGGCGCGAGTGCGTATGCACAATGCCGCCAATGGAGGGGAATGCCTGATAGAGCAGCCGGTGAGTTGGCGTGTCGGAGGAGGGCTTTTTCGTACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAG  >  NC_000913/66021‑66505                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                 gcccATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCTTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAAGGGTGCCGTAGAAATAGTCGGCG‑TGGGTGGTGCCGGTGGCTGGAATCGTCTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGca                                                                                                                                                                                                                                              <  2:334298/249‑1 (MQ=255)        aaaatgggCCGGGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATTCCATGTTTTGCGACGGTTTCTACGATGACGTTACCGGTTTCCCACTCAAATTCTCCTGTGATTTCTGCGTCGGTCATTTCGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCG‑TGGGTGGTGCCGTTTGCTGGAACCGACGGGTCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATgccg                                                                                                                                                                                                                                      <  1:71902/246‑1 (MQ=255)                ccGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCG‑TGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGGGTGTGTGCACAATGCCGCCCATgg                                                                                                                                                                                                                               >  1:260598/1‑250 (MQ=255)                      gAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATGCCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCG‑TGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCCCAATGCCGCCAATGGACGGGa                                                                                                                                                                                                                         >  1:418595/1‑250 (MQ=255)                          ggACCAGGACGCCGGGCATTTGCGCTTCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTGCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCG‑TGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGcc                                                                                                                                                                                                                    <  1:335307/251‑1 (MQ=255)                                     ccGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCG‑TGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAAAGCCGCCAATGGACGGGAATGCCTGATAGagc                                                                                                                                                                                                           >  1:193923/1‑249 (MQ=255)                                                  cTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGGTGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGCCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCG‑TGGGTGGTGCCGGTTGCTGGCATCGACTGGCCCGCCTGCGCCCAGATGGTGGCGAGGCGGGAGTGTGTGTGCACAATGCCGCCAATGTACGGGCATTCCTGATAGAGCAGTCGGTGAGTTggg                                                                                                                                                                                            >  1:259780/1‑250 (MQ=255)                                                      aTCGATTCCCTGTTTTTCGAAGGTTTCTACGATGACTTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCG‑TGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCgtat                                                                                                                                                                                        <  1:128284/251‑3 (MQ=255)                                                      aTCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGTTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCG‑TGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGGGCACAATGCCCCCAATGGACGGGAATGCCTGATAGAGGCGGCGGTGAGTTggggtat                                                                                                                                                                                        >  2:198964/1‑246 (MQ=255)                                                         gATACCCTTTTTGTCGAAGTTTTCTACGATGCCTGTACCGTTTTCACACTCAGCTGCACCGTTGATTGCTGCGTCGGTTATTGTGCGGGTGCAGGGAATGGTGCCGTAGAAAGAGTCGGCG‑TGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGGGGCGTGGCGCGAGTGTGTGGGCACAATGCCGTCAATGGACTGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCgt                                                                                                                                                                                     <  2:36664/251‑2 (MQ=255)                                                                  tttttCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCG‑TGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGGGGCGGGCt                                                                                                                                                                            >  1:429037/1‑251 (MQ=255)                                                                         aaGGTTTCTTCGATGACGTTTCCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGTTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCT‑TGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCgc                                                                                                                                                                     <  1:324286/251‑2 (MQ=255)                                                                                  aCGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCG‑TGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAAc                                                                                                                                                            >  1:385792/1‑251 (MQ=255)                                                                                                                 caccGTTGATTGCTGCGTCGTTCATTTGGCGGGGGCAGGGAATGGTTCCGTAGAAATAGTCGGCGGTGGG‑GGTGCCGGTTGCTGGAATCGCCTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGCCGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGAc                                                                                                                             <  2:170752/248‑1 (MQ=255)                                                                                                                  accGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCG‑TGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGGGCACAATGCCGCCAATGGACGGGAATGCCTGGTAGAGCAGTCGGTGAGGTGGCGTATCGGGGGAGGGCTTTTTAGCACCTTCAACCCCTT‑CACCGGTTTCGGTGCTAaccaccaccc                                                                                                                           >  2:90375/2‑250 (MQ=255)                                                                                                                  accGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCG‑TGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGGGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATCGAGCAGGCGGTGAGTTGGCGGATCGGGGGGGGGCTTTTTCGCACCTTCAACCACTT‑CCCCGGTTTCGATGCTAACCACGAcc                                                                                                                            >  2:420630/2‑250 (MQ=255)                                                                                                                   ccGTTGGTTTCTGCGTCAGTCATTTTGCGGGTGCAGGGGATGATGCCGTAGAAATAGTCGGCG‑TGGGTGGTGCCGGTTGCTGGAATCGACTGGTCCGCCTGCGCCCAGATAGAGGCGTGGCGCGCGGGTGTGTGGACCATGCCGCCAATGGACGGGAATGCCTGATAGAGGAGTCGGTGAGGTGGCCTATCGGAGGAGGGGTTTTTCGCACCTTCAACCACTTTCACGGGTTT‑GATGcgaaacccggccc                                                                                                                           >  2:80916/1‑238 (MQ=255)                                                                                                                   ccGTTGATTTCTGGGTCGGTCATTTTGCGGGGGCAGGGAATGGTGCCTTAGAAAGAGTCGGCG‑TGGGTGGTGCCGGTTGCTGGAAGCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGACCat                                                                                                                          <  1:289008/251‑1 (MQ=255)                                                                                                                        gATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCG‑TGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAggg                                                                                                                                                                              >  1:443060/1‑195 (MQ=255)                                                                                                                                        aTTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCG‑TGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGGATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGGGGAGGGCTTTTTCGCACCTTCAACCACTT‑CACCGGCTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCt                                                                                                     >  1:157203/1‑251 (MQ=255)                                                                                                                                                    cAGGGAATGGTGCCGTAGAAATAGTCGGCG‑TGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGcc                                                                                         <  1:124492/251‑1 (MQ=255)                                                                                                                                                               gCCGTAGAAATAGTCGGCG‑TGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGa                                                                              <  1:413474/251‑1 (MQ=255)                                                                                                                                                               gCCGTAGAAATAGTCGGCG‑TGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGa                                                                              <  1:198964/251‑1 (MQ=255)                                                                                                                                                                     gAAATAGTCGGCG‑TGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATTGACGCCGGAAGGTTTGATcaca                                                                         >  1:133066/1‑250 (MQ=255)                                                                                                                                                                      aaaTAGTCGGCT‑TGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACaaa                                                                       <  1:379824/251‑1 (MQ=255)                                                                                                                                                                                 g‑tgggtggTGCCGGTTGTTGGATGCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGTCGCGAGGGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCCCCTTCAACCACTT‑CACCGGTTTCTATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCt                                                            <  1:423751/251‑1 (MQ=255)                                                                                                                                                                                          tgCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGAAGCCGCGCTCGCGATc                                                     >  1:293713/1‑250 (MQ=255)                                                                                                                                                                                                     gaaaccgcttggcccGCTTGCGCCCGGATAGTGGCGTGGCGTGAGTGTGTGTGCACAATGCCGCCAATGGACGGTAATGCCTGATAGAGCAGTCGGTGAGGTGGCGTATCGGAGGAGGGCTTTTTCGCAGTTTCAACCACTT‑CACCGTGTTCGATGCTAACCACGTCCATAGCGTCATCGCTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGAc                                         <  1:427298/239‑1 (MQ=255)                                                                                                                                                                                                             ctggcccGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGcccc                                 >  2:160563/5‑251 (MQ=255)                                                                                                                                                                                                                    gCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAg                          <  1:41979/251‑1 (MQ=255)                                                                                                                                                                                                                         cgcCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGa                     <  1:160563/251‑1 (MQ=255)                                                                                                                                                                                                                             cAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGACCATATAGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGCTCAACGGCGCTGACGTTTCCCCATGTGAGCGTGAcc                   >  1:88820/1‑249 (MQ=255)                                                                                                                                                                                                                                atagtggcgtggcGCGAGTGTTGGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGtt              <  1:438626/248‑1 (MQ=255)                                                                                                                                                                                                                                        gtggcgCGAGGGGGTGTGCACAAGGCCGACAATGGACGGGAATGACTGATAGAGCAGTCGGTGAGTTGGCGTATAGGAGGAGGGCTTTTGCGCACCTTCAACCACTT‑CACCGGTTGCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGTCCCATGTGAGCGTGACCAGGTTATGTTTTg      <  1:185370/251‑1 (MQ=255)                                                                                                                                                                                                                                          ggcgCGATTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTgg     <  1:114107/250‑1 (MQ=255)                                                                                                                                                                                                                                            cgcgAGTGTGTGTGCCCAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTTAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAg  <  1:2956/251‑1 (MQ=255)                                                                                                                                                                                                                                            cgagAGTGTGTGTGCACAATGCCGCCAATGGAAGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAGAGACGCAGCGCTCGCGTTCAACGGCGCTGACGTTGCCCCATGGGAGCGTGACAAGGTTATGTTTTGGc    >  2:180700/4‑249 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                 CCCCATGCAAACGGGCCGTGGGAATGGACCAGAACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCAAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCGCCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCG‑TGGGTGGTGCCGGTTGCTGGAATCGACTGACCCGCCTGCGCCCAGATGGTGGCGTGGCGCGAGTGCGTATGCACAATGCCGCCAATGGAGGGGAATGCCTGATAGAGCAGCCGGTGAGTTGGCGTGTCGGAGGAGGGCTTTTTCGTACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAG  >  NC_000913/66021‑66505

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