Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,026,812	(C)6→7	coding (204/459 nt)	mgsA ←	methylglyoxal synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,026,806	1	.	C	98.4%	218.0 / ‑3.5	62	coding (210/459 nt)	mgsA	methylglyoxal synthase
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (32/29);  total (32/30)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

ACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCA‑‑CCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGTTGCATACAGTACGTGTTGTTCCAGTAACGGTTGATGCCGTTCCACCCAG  >  NC_000913/1026669‑1026932                                                                                                                                           |                                                                                                                                aCCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCA‑‑cc                                                                                                                              >  1:485783/1‑140 (MQ=255)  ccGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCA‑‑ccc                                                                                                                             >  2:770489/1‑140 (MQ=255)            ccATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCAC‑CCCCCCATTggg                                                                                                                    <  2:167154/140‑1 (MQ=255)               tACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCAC‑CCCCCCATTGGGCCa                                                                                                                 <  1:1022233/140‑1 (MQ=255)                aCCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCAC‑CCCCCCATTGGGCCAc                                                                                                                >  1:797319/1‑140 (MQ=255)                    tCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCCCATTGGGCCACTCa                                                                                                             >  1:943321/1‑140 (MQ=255)                    tCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCAC‑CCCCCCATTGGGCCACTCaa                                                                                                            >  1:961183/1‑140 (MQ=255)                    tCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCAC‑CCCCCCATTGGGCCACTCaa                                                                                                            >  2:952508/1‑140 (MQ=255)                     cGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCAC‑CCCCCCATTGGGCCACTCAAc                                                                                                           >  2:948324/1‑140 (MQ=255)                            cGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCA‑‑CCCCCCATTGGGCCACTCAACATCGCGtt                                                                                                   <  2:485783/140‑1 (MQ=255)                               agcaAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCCCATTGGGCCACTCAACATCGCGTTg                                                                                                  <  2:943321/140‑1 (MQ=255)                                 caAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCAC‑CCCCCCATTGGGCCACTCAACATCGCGTTGACg                                                                                               >  2:176846/1‑140 (MQ=255)                                     gCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCAC‑CCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCa                                                                                           >  2:468284/1‑140 (MQ=255)                                      cTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCAC‑CCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCAt                                                                                          >  1:15149/1‑140 (MQ=255)                                       tttCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCAC‑CCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATg                                                                                         >  2:472267/1‑140 (MQ=255)                                         tCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCAC‑CCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGcc                                                                                       >  2:526099/1‑140 (MQ=255)                                                aGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCAC‑CCCCCCATTGGGCCACACAACATCGCGTTGACGTTCATGCCGGTcgcg                                                                                <  1:181151/140‑1 (MQ=255)                                                 ggATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCAC‑CCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTcgcgc                                                                               <  2:961183/140‑1 (MQ=255)                                                       tGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCAC‑CCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGaaa                                                                         >  1:1048155/1‑140 (MQ=255)                                                       tGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCAC‑CCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGaaa                                                                         >  1:753011/1‑140 (MQ=255)                                                       tGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCAC‑CCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGaaa                                                                         >  2:726660/1‑140 (MQ=255)                                                       tGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCAC‑CCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGaaa                                                                         >  1:168863/1‑140 (MQ=255)                                                       tGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCAC‑CCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGaaa                                                                         >  1:1268/1‑140 (MQ=255)                                                            cacggcaTTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCAC‑CCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAg                                                                    >  1:588329/1‑140 (MQ=255)                                                              cggcaTTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCAC‑CCCCCCATTGGGCCACTCAACATCCCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGtt                                                                  <  2:927289/140‑1 (MQ=255)                                                                gcaTTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCAC‑CCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTAc                                                                <  1:468284/140‑1 (MQ=255)                                                                 caTTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCAC‑CCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTAcc                                                               >  1:13216/1‑140 (MQ=255)                                                                   tttAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCAC‑CCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCgg                                                             >  1:217879/1‑140 (MQ=255)                                                                       gTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCAC‑CCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGt                                                         <  1:968848/140‑1 (MQ=255)                                                                         ggATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCAC‑CCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGc                                                       <  2:768668/140‑1 (MQ=255)                                                                         ggATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCAC‑CCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGc                                                       <  2:797319/140‑1 (MQ=255)                                                                         ggATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCAC‑CCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGc                                                       <  2:753011/140‑1 (MQ=255)                                                                         ggATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCAC‑CCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGc                                                       <  2:15149/140‑1 (MQ=255)                                                                         ggATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCAC‑CCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGc                                                       <  1:176846/140‑1 (MQ=255)                                                                          gATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCAC‑CCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGcc                                                      >  1:538267/1‑140 (MQ=255)                                                                           aTCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCAC‑CCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTAc                                                                >  1:187333/1‑129 (MQ=255)                                                                           aTCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCAC‑CCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTAc                                                                <  2:187333/129‑1 (MQ=255)                                                                           aTCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCAC‑CCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCt                                                     <  2:168863/140‑1 (MQ=255)                                                                            tCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCAC‑CCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTg                                                    <  2:1268/140‑1 (MQ=255)                                                                              ccAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCAC‑CCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGtt                                                  <  2:588329/140‑1 (MQ=255)                                                                                agaagaAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCAC‑CCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGTTGc                                                >  2:287195/1‑140 (MQ=255)                                                                                  aagaaAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCAC‑CCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCt                                                     >  2:862941/1‑133 (MQ=255)                                                                                  aagaaAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCAC‑CCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCt                                                     <  1:862941/133‑1 (MQ=255)                                                                                  aagaaAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCAC‑CCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGTTGCAt                                              <  1:726660/140‑1 (MQ=255)                                                                                     aaaaTCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCAC‑CCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGTTGCATACa                                           <  1:19505/140‑1 (MQ=255)                                                                                         tCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCAC‑CCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGTTGCATACAGTAc                                       >  1:851612/1‑140 (MQ=255)                                                                                            aTACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTAAC‑CCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGTTGCATACAGTACgtg                                    >  2:603809/1‑140 (MQ=255)                                                                                             tACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCAC‑CCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGTTGCATACAGTACgtgt                                   <  1:948324/140‑1 (MQ=255)                                                                                                aTCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCAC‑CCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGTTGCATACAGTACGtgttgt                                <  1:770489/140‑1 (MQ=255)                                                                                                       ttCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCAC‑CCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGTTGCATACAGTACGTGTTGTTCCAGTa                         <  1:472267/140‑1 (MQ=255)                                                                                                        tCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCAC‑CCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGTTGCATACAGTACGTGTTGTTCCAGTaa                        <  2:1048155/140‑1 (MQ=255)                                                                                                         cccTTCTGAGATCAATGCGCCAACCTGCTGGTCAC‑CCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGTTGCATACAGTACGTGTTGTTCCAGTAAc                       <  2:96839/140‑1 (MQ=255)                                                                                                          ccTTCTGAGATCAATGCGCCAACCTGCTGGTCAC‑CCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGTTGCATACAGTACGTGTTGTTCCAGTAACg                      >  1:670103/1‑140 (MQ=255)                                                                                                               tGAGATCAATGCGCCAACCTGCTGGTCAC‑CCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGTTGCATACAGTACGTGTTGTTCCAGTAACGGTTGa                 >  2:797032/1‑140 (MQ=255)                                                                                                               tGAGATCAATGCGCCAACCTGCTGGTCAC‑CCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGTTGCATACAGTACGTGTTGTTCCAGTAACGGTTGa                 <  1:952508/140‑1 (MQ=255)                                                                                                                gagaTCAATGCGCCAACCTGCTGGTCAC‑CCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGTTGCATACAGTACGTGTTGTTCCAGTAACGGTTGAt                <  2:538267/140‑1 (MQ=255)                                                                                                                    tCAATGCGCCAACCTGCTGGTCAC‑CCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGTTGCATACAGTACGTGTTGTTCCAGTAACGGTTGATGCCg            >  2:559207/1‑140 (MQ=255)                                                                                                                      aaTGCGCCAACCTGCTGGTCAC‑CCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGTTGCATACAGTACGTGTTGTTCCAGTAACGGTTGATGCCGtt          <  2:217879/140‑1 (MQ=255)                                                                                                                        tGCGCCAACCTGCTGGTCAC‑CCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGTTGCATACAGTACGTGTTGTTCCAGTAACGGTTGATGCCGTTcc        >  1:573864/1‑140 (MQ=255)                                                                                                                        tGCGCCAACCTGCTGGTCAC‑CCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGTTGCATACAGTACGTGTTGTTCCAGTAACGGTTGATGCCGTTcc        >  1:666492/1‑140 (MQ=255)                                                                                                                        tGCGCCAACCTGCTGGTCAC‑CCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGTTGCATACAGTACGTGTTGTTCCAGTAACGGTTGATGCCGTTcc        >  2:329830/1‑140 (MQ=255)                                                                                                                         gcgcCAACCTGCTGGTCAC‑CCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGTTGCATACAGTACGTGTTGTTCCAGTAACGGTTGATGCCGTTCCa       <  1:329830/140‑1 (MQ=255)                                                                                                                         gcgcCAACCTGCTGGTCAC‑CCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGTTGCATACAGTACGTGTTGTTCCAGTAACGGTTGATGCCGTTCCa       >  2:385237/1‑140 (MQ=255)                                                                                                                             cAACCTGCTGGTCAC‑CCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGTTGCATACAGTACGTGTTGTTCCAGTAACGGTTGATGCCGTTCCACCCa   >  1:169046/1‑140 (MQ=255)                                                                                                                              aaCCTGCTGGTCAC‑CCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGTTGCATACAGTACGTGTTGTTCCAGTAACGGTTGATGCCGTTCCACCCAg  <  1:385237/140‑1 (MQ=255)                                                                                                                                           |                                                                                                                                ACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCA‑‑CCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGTTGCATACAGTACGTGTTGTTCCAGTAACGGTTGATGCCGTTCCACCCAG  >  NC_000913/1026669‑1026932

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