Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,761,850	C→T	A36T (GCG→ACG)	yibF ←	glutathione S‑transferase homolog

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,761,850	0	C	T	100.0%	299.9 / NA	92	A36T (GCG→ACG)	yibF	glutathione S‑transferase homolog
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (48/44);  total (48/44)

AACATCGCCGGAGCGACATTCATTAATTCAATATATTCAGCGATGATCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGCGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATAT  >  NC_000913/3761717‑3761983                                                                                                                                      |                                                                                                                                      aaCATCGCCGGAGCGACATTCATTAATTCAATATATTCAGCGATGATCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTtat                                                                                                                                  >  1:156519/1‑139 (MQ=255)     tCGCCGGAGCGACATTCATTAATTCAATATATTCAGCGATGATCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAggg                                                                                                                              <  2:364882/139‑1 (MQ=255)     tCGCCGGAGCGACATTCATTAATTCAATATATTCAGCGATGATCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAggg                                                                                                                              <  2:1210965/139‑1 (MQ=255)      cGCCGGAGCGACATTCATTAATTCAATATATTCAGCGATGATCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGc                                                                                                                             <  1:350521/139‑1 (MQ=255)      cGCCGGAGCGACATTCATTAATTCAATATATTCAGCGATGATCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGc                                                                                                                             <  2:156519/139‑1 (MQ=255)                aCATTCATTAATTCAATATATTCAGCGATGATCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAggg                                                                                                                              <  1:772069/128‑1 (MQ=255)                aCATTCATTAATTCAATATATTCAGCGATGATCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAggg                                                                                                                              >  2:772069/1‑128 (MQ=255)                      attAATTCAATATATTCAGCGATGATCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAAtt                                                                                                             >  1:1541552/1‑139 (MQ=255)                           ttCAATATATTCAGCGATGATCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGa                                                                                                          <  1:1631682/137‑1 (MQ=255)                           ttCAATATATTCAGCGATGATCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGa                                                                                                          >  2:1631682/1‑137 (MQ=255)                              aaTATATTCAGCGATGATCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGtt                                                                                                     <  2:1566654/139‑1 (MQ=255)                                tatatTCAGCGATGATCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTAt                                                                                                   >  2:680352/1‑139 (MQ=255)                                tatatTCAGCGATGATCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTAt                                                                                                   >  1:857929/1‑139 (MQ=255)                                tatatTCAGCGATGATCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTAt                                                                                                   >  2:169704/1‑139 (MQ=255)                                tatatTCAGCGATGATCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTAt                                                                                                   >  1:919352/1‑139 (MQ=255)                                tatatTCAGCGATGATCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTAt                                                                                                   >  2:985420/1‑139 (MQ=255)                                  tcttCAGAGATGATCCGCGAATCAAACCACCATTCGCCCTCTTCGGTCACCTGCACCGGCACTTTTCCTAACGGTTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGc                                                                                                 <  2:919352/137‑1 (MQ=255)                                  tatTCAGCGATGATCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGc                                                                                                 <  1:680352/139‑1 (MQ=255)                                    ttCAGCGATGATCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGccc                                                                                               >  2:919015/1‑139 (MQ=255)                                     tCAGCGATGATCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAggg                                                                                                                              >  2:512535/1‑107 (MQ=255)                                     tCAGCGATGATCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAggg                                                                                                                              <  1:512535/107‑1 (MQ=255)                                     tCAGCGATGATCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCt                                                                                              <  1:985420/139‑1 (MQ=255)                                       aGCGATGATCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCttt                                                                                            >  2:1259179/1‑139 (MQ=255)                                        gCGATGATCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCtttt                                                                                           >  2:449608/1‑139 (MQ=255)                                         cGATGATCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTta                                                                                                                                   <  1:1357760/98‑1 (MQ=255)                                         cGATGATCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTta                                                                                                                                   >  2:1357760/1‑98 (MQ=255)                                         cGATGATCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTc                                                                                          >  1:877507/1‑139 (MQ=255)                                          gatgatCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCt                                                                                         >  2:123959/1‑139 (MQ=255)                                          gatgatCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCt                                                                                         >  2:495158/1‑139 (MQ=255)                                          gatgatCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCt                                                                                         >  1:153088/1‑139 (MQ=255)                                          gatgatCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCt                                                                                         <  1:1259179/139‑1 (MQ=255)                                          gaggatCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCTCTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCt                                                                                         >  2:960006/4‑139 (MQ=255)                                            tgatCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTaa                                                                                       <  2:474823/139‑1 (MQ=255)                                              atCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGa                                                                                                          <  1:975951/118‑1 (MQ=255)                                              atCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGa                                                                                                          >  2:975951/1‑118 (MQ=255)                                              atCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTaaca                                                                                     >  2:1250935/1‑139 (MQ=255)                                              atCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTaaca                                                                                     >  2:738811/1‑139 (MQ=255)                                              atCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAACACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGa                                                                                                          >  1:1250251/1‑118 (MQ=255)                                              atCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAACACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGa                                                                                                          <  2:1250251/118‑1 (MQ=255)                                                cGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTaacaac                                                                                   <  1:1038753/139‑1 (MQ=255)                                                cGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTaacaac                                                                                   <  2:943640/139‑1 (MQ=255)                                                cGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTaacaac                                                                                   <  2:1541552/139‑1 (MQ=255)                                                cGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTaacaac                                                                                   <  1:1250935/139‑1 (MQ=255)                                                cGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTaacaac                                                                                   <  1:449608/139‑1 (MQ=255)                                                   cGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGa                                                                                                          >  1:1371970/1‑113 (MQ=255)                                                   cGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGa                                                                                                          <  2:1371970/113‑1 (MQ=255)                                                     aaTCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAAt                                                                              >  1:1221944/1‑139 (MQ=255)                                                        cAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGa                                                                           >  2:1001328/1‑139 (MQ=255)                                                                cATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGc                                                                   >  2:490134/1‑139 (MQ=255)                                                                      cccTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACaa                                                             <  2:1480116/139‑1 (MQ=255)                                                                          ctTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAg                                                                            <  1:1505909/120‑1 (MQ=255)                                                                          ctTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAg                                                                            >  2:1505909/1‑120 (MQ=255)                                                                          ctTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACgg                                                         <  1:490134/139‑1 (MQ=255)                                                                          ctTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACgg                                                         <  2:524825/139‑1 (MQ=255)                                                                            tCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGc                                                       <  1:236899/139‑1 (MQ=255)                                                                             cGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCt                                                      <  2:1221944/139‑1 (MQ=255)                                                                                   ccAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTa                                                >  2:219177/1‑139 (MQ=255)                                                                                     aGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTaacaac                                                                                   >  2:51503/1‑102 (MQ=255)                                                                                     aGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTaacaac                                                                                   <  1:51503/102‑1 (MQ=255)                                                                                     aGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGc                                              >  1:1389691/1‑139 (MQ=255)                                                                                         ccGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGCGTAGCTAcc                                          >  2:975727/1‑139 (MQ=255)                                                                                          cGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACcg                                         >  1:1298005/1‑139 (MQ=255)                                                                                            gCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAggg                                                                                                                              <  1:822101/52‑1 (MQ=255)                                                                                            gCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAggg                                                                                                                              >  2:822101/1‑52 (MQ=255)                                                                                            gCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATGAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACcgac                                       >  2:1186985/1‑139 (MQ=255)                                                                                            gCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACcgac                                       <  2:877507/139‑1 (MQ=255)                                                                                                  ttCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACaaa                                                            >  1:623499/1‑112 (MQ=255)                                                                                                  ttCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACaaa                                                            <  2:623499/112‑1 (MQ=255)                                                                                                       aaCGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAACAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGc                            >  1:1117098/1‑139 (MQ=255)                                                                                                        aCGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGcc                           >  2:676457/1‑139 (MQ=255)                                                                                                        aCGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTACAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGcc                           <  1:1001328/139‑1 (MQ=255)                                                                                                         cGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCg                          <  1:975727/139‑1 (MQ=255)                                                                                                            gTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTGATAAATTCGAAAGTTATGCCCTTTTCTaacaac                                                                                   >  2:907281/1‑79 (MQ=255)                                                                                                            gTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTGATAAATTCGAAAGTTATGCCCTTTTCTaacaac                                                                                   <  1:907281/79‑1 (MQ=255)                                                                                                            gTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGa                                                                                                          >  2:1468108/1‑56 (MQ=255)                                                                                                            gTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGa                                                                                                          <  1:1468108/56‑1 (MQ=255)                                                                                                            gTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCt                                                                                              >  2:1623931/1‑68 (MQ=255)                                                                                                            gTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCt                                                                                              <  1:1623931/68‑1 (MQ=255)                                                                                                            gTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGt                       >  2:803570/1‑139 (MQ=255)                                                                                                               aaaTTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGa                                                                                                          <  1:928545/53‑1 (MQ=255)                                                                                                               aaaTTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGa                                                                                                          >  2:928545/1‑53 (MQ=255)                                                                                                               aaaTTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCt                    <  1:891168/139‑1 (MQ=255)                                                                                                                     cgcCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATg                                                                                                  >  2:503585/1‑55 (MQ=255)                                                                                                                     cgcCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATg                                                                                                  <  1:503585/55‑1 (MQ=255)                                                                                                                     cgcCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGaa                                                                          <  2:769402/79‑1 (MQ=255)                                                                                                                     cgcCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGaa                                                                          >  1:769402/1‑79 (MQ=255)                                                                                                                     cgcCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGTGc              >  2:1338922/1‑139 (MQ=255)                                                                                                                      gcCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACgg                                                         >  2:1220353/1‑95 (MQ=255)                                                                                                                      gcCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACgg                                                         <  1:1220353/95‑1 (MQ=255)                                                                                                                            ccGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGTGCGAGGaaa       <  2:1298005/139‑1 (MQ=255)                                                                                                                               ttGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAat    >  1:729189/1‑139 (MQ=255)                                                                                                                                 gTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAatat  <  1:919015/139‑1 (MQ=255)                                                                                                                                      |                                                                                                                                      AACATCGCCGGAGCGACATTCATTAATTCAATATATTCAGCGATGATCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGCGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATAT  >  NC_000913/3761717‑3761983

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