Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,881,694	C→T	R39H (CGC→CAC)	dmlR ←	transcriptional activator of dmlA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,881,694	0	C	T	100.0%	220.6 / NA	62	R39H (CGC→CAC)	dmlR	transcriptional activator of dmlA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (33/29);  total (33/29)

AAAGTTCGTCCATCATCTGACCGACATCCTGCAAAATCCGCTGCGCCCATTCATAAATCCGCTCTCCTTCTTCGGTAATGGTGACGCGGCGGGTGGTGCGGTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCGATGCGCTTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGCCACGGCGGCAAAACCAGCCCGGCGAGCCACCAGCATAAAGACGCGCAAATCATTCAGCAGCGGTAAATTATTCATGATTCGTGTTTTATGTTTCACCAGTTACGG  >  NC_000913/1881551‑1881839                                                                                                                                                |                                                                                                                                                  aaaGTTCGTCCATCATCTGACCGACATCCTGCAAAATCCGCTGCGCCCATTCATAAATCCGCTCTCCTTCTTCGGTAATGGTGACGCGGCGGGTGGTGCGGTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCGATGtgcttgct                                                                                                                                            >  1:1047877/1‑151 (MQ=255) aaaGTTCGTCCATCATCTGACCGACATCCTGCAAAATCCGCTGCGCCCATTCATAAATCCGCTCTCCTTCTTCGGTAATGGTGACGCGGCGGGTGGTGCGGTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCGATGtgcttgct                                                                                                                                            <  2:861439/151‑1 (MQ=255)  aaGTTCGTCCATCATCTGACCGACATCCTGCAAAATCCGCTGCGCCCATTCATAAATCCGCTCTCCTTCTTCGGTAATGGTGACGCGGCGGGTGGTGCGGTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCGATGTGCTTGCTg                                                                                                                                           <  2:160186/151‑1 (MQ=255)     ttCGTCCATCATCTGACCGACATCCTGCAAAATCCGCTGCGCCCATTCATAAATCCGCTCTCCTTCTTCGGTAATGGTGACGCGGCGGGTGGTGCGGTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCGATGTGCTTGCTGACg                                                                                                                                        >  1:900614/1‑151 (MQ=255)               atcTGACCGACATCCTGCAAAATCCGCTGCGCCCATTCATAAATCCGCTCTCCTTCTTCGGTAATGGTGACGCGGCGGGTGGTGCGGTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCGATGTGCTTGCTGACGAACGCCGGTg                                                                                                                              >  1:952237/1‑151 (MQ=255)                   gACCGACATCCTGCAAAATCCGCTGCGCCCATTCATTAATCCGCTCTCCTTCTTCGGTAATGGTGACGCGGCGGGAGGTTCGGTGTAACAACACCCCGTTTAGCGTTTGCTCCAGCAAGGCGATGTGCTTGCTGACGAACGCCGGTGAAAc                                                                                                                          >  2:156712/1‑151 (MQ=255)                   gACCGACATCCTGCAAAATCCGCTGCGCCCATTCATAAATCCGCTCTCCTTCTTCGGTAATGGTGACGCGGCGGGTGGTGCGGTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCGATGTGCTTGCTGACGAACGCCGGTGaa                                                                                                                            <  1:606239/149‑1 (MQ=255)                      cGACATCCTGCAAAATCCGCTGCGCCCATTCATAAATCCGCTCTCCTTCTTCGGTAATGGTGACGCGGCGGGTGGTGCGGTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCGATGTGCTTGCTGACGAACGCCGGTGAAACGcc                                                                                                                       >  2:782926/1‑151 (MQ=255)                       gACATCCTGCAAAATCCGCTGCGCCCATTCATAAATCCGCTCTCCTTCTTCGGTAATGGTGACGCGGCGGGTGGTGCGGTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCGATGTGCTTGCTGACGAACGCCGGTGAAACGccc                                                                                                                      <  1:62986/151‑1 (MQ=255)                            ccTGCAAAATCCGCTGCGCCCATTCATAAATCCGCTCTCCTTCTTCGGTAATGGTGACGCGGCGGGTGGTGCGGTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCGATGTGCTTGCTGACGAACGCCGGTGAAACGCCCAGttc                                                                                                                 <  2:490507/151‑1 (MQ=255)                              tGCAAAATCCGTTGCGCCCATTCATAAATCCGCTCTCCTTCTTCGGTAATGGTGACGCGGCGGGTGGTGCGGTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCGATGTGCTTGCTGACGAACGCCGGTGAAACGCCCAGttctt                                                                                                               >  1:19663/1‑151 (MQ=255)                               gCAAAATCCGCTGCGCCCATTCATAAATCCGCTCTCCTTCTTCGGTAATGGTGACGCGGCGGGTGGTGCGGTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCGATGTGCTTGCTGACGAACGCCGGTGAAACGCCCAGttcttc                                                                                                              <  1:260010/151‑1 (MQ=255)                               gCAAAATCCGCTGCGCCCATTCATAAATCCGCTCTCCTTCTTCGGTAATGGTGACGCGGCGGGTGGTGCGGTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCGATGTGCTTGCTGACGAACGCCGGTGAAACGCCCAGttcttc                                                                                                              <  1:514402/151‑1 (MQ=255)                                  aaaTCCGCTGCGCCCATTCATAAATCCGCTCTCCTTCTTCGGTAATGGTGACGCGGCGGGTGGTGCGGTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCGATGTGCTTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGc                                                                                                           >  1:670996/1‑151 (MQ=255)                                   aaTCCGCTGCGCCCATTCATAAATCCGCTCTCCTTCTTCGGTAATGGTGACGCGGCGGGTGGTGCGGTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCGATGTGCTTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGcc                                                                                                          >  1:813839/1‑151 (MQ=255)                                   aaTCCGCTGCGCCCATTCATAAATCCGCTCTCCTTCTTCGGTAATGGTGACGCGGCGGGTGGTGCGGTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCGATGTGCTTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGcc                                                                                                          >  1:1050388/1‑151 (MQ=255)                                     tCCGCTGCGCCCATTCATAAATCCGCTCTCCTTCTTCGGTAATGGTGACGCGGCGGGTGGTGCGGTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCGATGTGCTTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGCCAc                                                                                                        >  1:792742/1‑151 (MQ=255)                                       cGCTGCGCCCATTCATAAATCCGCTCTCCTTCTTCGGTAATGGTGACGCGGCGGGTGGTGCGGTGTAACTACACCACGTTTAGCGTTTGCTCCAGCAAGGCGATGTGCTTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGCCAcgg                                                                                                      <  2:21061/151‑1 (MQ=255)                                       cGCTGCGCCCATTCATAAATCCGCTCTCCTTCTTCGGTAATGGTGACGCGGCGGGTGGTGCGGTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCGATGTGCTTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGCCAcgg                                                                                                      <  1:419948/151‑1 (MQ=255)                                           gcgcCCATTCATAAATCCGCTCTCCTGCTTCGGTAATGGTGACGCGGCGGGTGGTGCGGTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCGATGTGCTTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGCCACggcggc                                                                                                  <  2:585972/151‑1 (MQ=255)                                            cgcCCATTCATAAATCCGCTCTCCTTCTTCGGTAATGGTGACGCGGCGGGTGGTGCGGTGTAACAACACCCCGTTTAGCGTTTGCTCCAGCAAGGCGATGTGCTTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGCCACggcggc                                                                                                  >  2:225808/1‑150 (MQ=255)                                               ccATTCATAAATCCGCTCTCCTTCTTCGGTAATGGTGACGCGGCGGGTGGTGCGGGGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCGATGTGCTTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGCCACGGCGGCaaa                                                                                               >  2:154227/1‑150 (MQ=255)                                                cATTCATAAATCCGCTCTCCTTCTTCGGTAATGGTGACGCGGCGGGTGGTGCGGTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCGATGTGCTTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGCCACGGCGGCAAAAc                                                                                             >  2:1034071/1‑151 (MQ=255)                                                       aaaTCCGCTCTCCTTCTTCGGTAATGGTGACGCGGCGGGTGGTGCGGTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCGATGTGCTTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGCCACGGCGGCAAAACCAGCCCg                                                                                      <  1:3796/151‑1 (MQ=255)                                                         aTCCGCTCTCCTTCTTCGGTAATGGTGACGCGGCGGGTGGTGCGGTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCGATGTGCTTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGCCACGGCGGCAAAACCAGCCCgg                                                                                     >  1:1891/1‑150 (MQ=255)                                                         aTCCGCTCTCCTTCTTCGGTAATGGTGACGCGGCGGGTGGTGCGGTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCGATGTGCTTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGCCACGGCGGCAAAACCAGCCCGGc                                                                                    >  1:554496/1‑151 (MQ=255)                                                            cGCTCTCCTTCTTCGGTAATGGTGACGCGGCGGGTGGTGCGGTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCGATGTGCTTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGCCACGGCGGCAAAACCAGCCCGGCGAg                                                                                 <  2:310321/151‑1 (MQ=255)                                                                  ccTTCTTCGGTAATGGTGACGCGGCGGGTGGTGCGGTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCGATGTGCTTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGCCACGGCGGCAAAACCAGCCCGGCGAGccacca                                                                           <  2:338010/151‑1 (MQ=255)                                                                      cttcGGTAATGGTGACGCGGCGGGTGGTGCGGTGTAACAACACCACGTTTATCGTTTGCTCCAGCAAGGCGTTGTGCTTGATGACGAAAGACGGCGAAACGCCCAGTTCTTCCGCCACGGCGGCAAAACCAGCCCGGCGAACCACCCGCa                                                                        >  2:633915/1‑150 (MQ=255)                                                                          ggTAATGGTGACGCGGCGGGTGGTGCGGTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCGATGTGCTTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGCCACGGCGGCAAAACCAGCCCGGCGAGCCACCAGCATaa                                                                     >  1:588100/1‑149 (MQ=255)                                                                            tAATGGTGACGCGGCGGGTGGTGCGGTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCGATGTGCTTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGCCACGGCGGCAAAACCAGCCCGGCGAGCCACCAGCATAAAGAc                                                                 >  1:616478/1‑151 (MQ=255)                                                                              aTGGTGACGCGGCTGGTGGGGCGGTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCGATGTGCTTGCTGACGAACGCCGGTGAAACGCCCAGTGCTTCCGCCACGGCGGCAAAACCAGCCCGGCGAGCCACCAGCATAAAg                                                                   <  2:540765/147‑1 (MQ=255)                                                                              aTGGTGACGCGGCGGGTGGTGCGGTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCGATGTGCTTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGCCACGGCGGCAAAACCAGCCCGGCGAGCCACCAGCATAAAGAcgc                                                               <  2:695401/151‑1 (MQ=255)                                                                                  tGACGCGGCGGGTGGTGCGGTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCGATGTGCTTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGCCACGGCGGCAAAACCAGCCCGGCGAGCCACCAGCATAAAGACgcgc                                                             >  1:817795/1‑149 (MQ=255)                                                                                  tGACGCGGCGGGTGGTGCGGTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCGATGTGCTTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGCCACGGCGGCAAAACCAGCCCGGCGAGCCACCAGCATAAAGACGCGCa                                                            >  1:817162/1‑150 (MQ=255)                                                                                             gtggtgCGGTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCGATGTGCTTGCTGACTAACGCCGGTGAAACGCCCAGTTCTTCCGCCACGGCGGCAAAACCAGCCCGGCGAGCCACCAGCATAAAGACGCGCAAATCATgcagca                                                <  2:791046/151‑1 (MQ=255)                                                                                             gtggtgCGGTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCGATGTGCTTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGCCACGGCGGCAAAACCAGCCCGGCGAGCCACCAGCATAAAGACGCGCAAATCATTcagca                                                <  1:53643/151‑1 (MQ=255)                                                                                              tggtgCGGTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCGATGTGCTTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGCCACGGCGGCAAAACCAGCCCGGCGAGCCACCAGCATAAAGACGCGCAAATCATTcagc                                                 >  1:1042865/1‑149 (MQ=255)                                                                                                 tgCGGTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCGATGTGCTTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGCCACGGCGGCAAAACCAGCCCGGCGAGCCACCAGCATAAAGACGCGCAAATCATTCAGCAGCgg                                            >  1:24985/1‑151 (MQ=255)                                                                                                  gCGGTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCGATGTGCTTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGCCACGGCGGCAAAACCAGCCCGGCGAGCCACCAGCATAAAGACGCGCAAATCATTCAGCAGCGGt                                           <  1:600390/151‑1 (MQ=255)                                                                                                    ggTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCGATGTGCTTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGCCACGGCGGCAAAACCAGCCCGGCGAGCCACCAGCATAAAGACGCGCAAATCATTCAGCAGCGGTaa                                         <  1:796580/151‑1 (MQ=255)                                                                                                    ggTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCGATGTGCTTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGCCACGGCGGCAAAACCAGCCCGGCGAGCCACCAGCATAAAGACGCGCAAATCATTCAGCAGCGGTaa                                         <  1:885790/151‑1 (MQ=255)                                                                                                     gtgtAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCGATGTGCTTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGCCACGGCGGCAAAACCAGCCCGGCGAGCCACCAGCATAAAGACGCGCAAATCATTCAGCAGCGGTaa                                         >  2:943110/1‑150 (MQ=255)                                                                                                          acaacaCCACGCTTAGCGTTTGCTCCAGCAAGGCGATGTGCTTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGCCACGGCGCCAAAACCAGCACGCCGAGCCCCCAGCATAAAGACGCGCACATCACTCAGCAGCGGTAAattatt                                   >  2:680072/1‑151 (MQ=255)                                                                                                           caacaCCACGTTTAGCGTTTGCTCCAGCAAGGCGATGTGCTTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGCCACGGCGGCAAAACCAGCCCGGCGAGCCACCAGCATAAAGACGCGCAAATCATTCAGCAGCGGTAAATTATTc                                  >  1:654954/1‑151 (MQ=255)                                                                                                           caacaCCACGTTTAGCGTGTGCTCAAGCAAGGCGATGTGCTGGCTGACGAAAGCCGGGGAAACGCTCAGCTCTTCCGCCCCGGCGGCAAAACCAGCCCGGCGAGCCACCAGCATAAACACGCGGAAATCATTCTGCAGCGGTCAatgatt                                   >  2:268410/1‑150 (MQ=255)                                                                                                                 cacGTTTAGCGTTTGCTCCAGCAAGGCGATGTGCTTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGCCACGGCGGCAAAACCAGCCCGGCGAGCCACCAGCATAAAGACGCGCAAATCATTCAGCAGCGGTAAATTATTCATGAt                             <  2:326839/150‑1 (MQ=255)                                                                                                                   cGTTTAGCGTTTGCTCCAGCAAGGCGATGTGCTTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGCCACGGCGGCAAAACCAGCCCGGCGAGCCACCAGCATAAAGACGCGCAAATCATTCAGCAGCGGTAAATTATTCATGATTCg                          <  1:310791/151‑1 (MQ=255)                                                                                                                    gTTTAGCGTTTGCTCCAGCAAGGCGATGTGCTTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGCCACGGCGGCAAAACCAGCCCGGCGAGCCACCAGCATAAAGACGCGCAAATCATTCAGCAGCGGTAAATTATTCATGATTc                           <  2:920610/149‑1 (MQ=255)                                                                                                                      ttAGCGTTTGCTCCAGCAAGGCGATGTGCTTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGCCACGGCGGCAAAACCAGCCCGGCGAGCCACCAGCATAAAGACGCGCAAATCATTCAGCAGCGGTAAATTATTCATGATTCgtgt                       >  2:458861/1‑151 (MQ=255)                                                                                                                      ttAGCGTTTGCTCCAGCAAGGCGATGTGCTTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGCCACGGCGGCAAAACCAGCCCGGCGAGCCACCAGCATAAAGACGCGCAAATCATTCAGCAGCGGTAAATTATTCATGATTCgtgt                       >  2:181482/1‑151 (MQ=255)                                                                                                                       tAGCGTTTGCTCCAGCAAGGCGATGTGCTTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGCCACGGCGTCAAAACCAGCCCGGCGAGCCACCAGCATAAAGACGCGCAAATCATTCAGCAGCGGTAAATTATTCATGATTCgtgt                       <  2:550502/150‑1 (MQ=255)                                                                                                                        aGCGTTTGCTCCAGCAAGGCGATGTGCTTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGCCACGGCGGCAAAACCAGCCCGGCGAGCCACCAGCATAAAGACGCGCAAATCATTCAGCAGCGGTAAATTATTCATGATTCGTGttt                     <  2:952188/151‑1 (MQ=255)                                                                                                                         gCGTTTGCTCCAGCAAGGCGATGTGCTTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGCCACGGCGGCAAAACCAGCCCGGCGAGCCACCAGCATAAAGACGCGCAAATCATTCAGCAGCGGTAAATTATTCATGATTCGTGtttt                    >  2:397365/1‑151 (MQ=255)                                                                                                                          cGTTTGCTCCAGCAAGGCGATGTGCTTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGCCACGGCGGCAAAACCAGCCCGGCGAGCCACCAGCATAAAGACGCGCAAATCATTCAGCAGCGGTAAATTATTCATGATTCGTGTTTTa                   <  2:464087/151‑1 (MQ=255)                                                                                                                           gTTTGCTCCAGCAAGGCGATGTGCTTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGCCACGGCGGCAAAACCAGCCCGGCGAGCCACCAGCATAAAGACGCGCAAATCATTCAGCAGCGGTAAATTATTCATGATTCGTGTTTTAt                  >  1:785698/1‑151 (MQ=255)                                                                                                                            tttGCTCCAGCAAGGCGATGTGCTTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGCCACGGCGGCAAAACCAGCCCGGCGAGCCACCAGCATAAAGACGCGCAAATCATTCAGCAGCGGTAAATTATTCATGATTCGTGTTTTATg                 >  1:785737/1‑151 (MQ=255)                                                                                                                               gCTCCAACAAGTCTATGTGCTTGCTGACGAAAGCCGGTGAAACGCCCAATTCTTCCGCCACGGCGGCAAAACCAGCCCGGCGAGCCACCAGCATAAACACGCGCAAATCATTCAGCCGCGGTAAATTATTCATGATTCGTGTTTTATGtt               <  2:800510/150‑1 (MQ=255)                                                                                                                                  ccAGCAAGGCGATGTGCTTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGCCACGGCGGCAAAACCAGCCCGGCGAGCCACCAGCATAAAGACGCGCAAATCATTCAGCAGCGGTAAATTATTCATGATTCGTGTTTTATGTTTCAc           <  2:813793/151‑1 (MQ=255)                                                                                                                                     gCAAGGCGATGTGCTTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGCCACGGCGGCAAAACCAGCCCGGCGAGCCACCAGCATAAAGACGCGCAAATCATTCAGCAGCGGTAAATTATTCATGATTCGTGTTTTATGTTTCACCAg        <  2:540469/151‑1 (MQ=255)                                                                                                                                      cAAGGCGATGTGCTTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGCCACGGCGGCAAAACCAGCCCGGCGAGCCACCAGCATAAAGACGCGCAAATCATTCAGCAGCGGTAAATTATTCATGATTCGTGTTTTATGTTTCACCAGt       >  2:988926/1‑151 (MQ=255)                                                                                                                                      cAAGGCGATGTGCTTGCTGACGAACGCCGGTGAAACGCCCAGTTCCTCCGCCACGGCGGCAAAACCAGCCCGGTGAGCCACCAGCATAAAGACGCGCAAATCATTCAGCAGCGGTAAATTATTCATGATTCGTGTTTGATGTTTCACCa         <  2:1006739/149‑1 (MQ=255)                                                                                                                                           cGATGTGCTTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGCCACGGCGGCAAAACCAGCCCGGCGCGCCACCAGCATAAAGACGCGCAAATCATTCAGCAGCGGTAAATTATTCATGATTCGTGTTTTATGTTTCACCAGTTACgg  >  2:749289/1‑151 (MQ=255)                                                                                                                                            gATGTGCTTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGCCACGGCGGCAAAACCAGCCCGGCGAGCCACCAGCATAAAGACGCGCAAATCATTCAGCAGCGGTAAATTATTCATGATTCGTGTTTTATGTTTCACCAGTTACgg  >  2:649356/1‑150 (MQ=255)                                                                                                                                                |                                                                                                                                                  AAAGTTCGTCCATCATCTGACCGACATCCTGCAAAATCCGCTGCGCCCATTCATAAATCCGCTCTCCTTCTTCGGTAATGGTGACGCGGCGGGTGGTGCGGTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCGATGCGCTTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGCCACGGCGGCAAAACCAGCCCGGCGAGCCACCAGCATAAAGACGCGCAAATCATTCAGCAGCGGTAAATTATTCATGATTCGTGTTTTATGTTTCACCAGTTACGG  >  NC_000913/1881551‑1881839

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