Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,486,216	G→T	6.8%	Q33H (CAG→CAT)	crp →	cAMP‑activated global transcription factor, mediator of catabolite repression

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,486,216	0	G	T	6.8%	195.1 / 2.5	63	Q33H (CAG→CAT)	crp	cAMP‑activated global transcription factor, mediator of catabolite repression
Reads supporting (aligned to +/- strand):  ref base G (29/29);  new base T (3/2);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

TTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAA‑AGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAG  >  NC_000913/3486075‑3486360                                                                                                                                              |                                                                                                                                                  ttATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGaaa                                                                                                                                            >  1:716080/1‑149 (MQ=255)     cTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTTTTCACCCTGGTGAATAAGc                                                                                                                                        <  1:470788/149‑1 (MQ=255)     cTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGc                                                                                                                                        <  1:749838/149‑1 (MQ=255)     cTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGc                                                                                                                                        <  1:483381/149‑1 (MQ=255)     cTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGc                                                                                                                                        <  1:199683/149‑1 (MQ=255)      tGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCg                                                                                                                                       <  2:649013/149‑1 (MQ=255)        gCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTTTTCACCCTAGTGAATAAGCgtg                                                                                                                                     <  2:448504/149‑3 (MQ=255)        gCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGa                                                                                                                                     <  2:59658/149‑1 (MQ=255)                 aaaGCTTATAACCGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGta                                                                                                                            <  1:106691/149‑1 (MQ=255)                        ataACAGAGGATAATCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGGAACGCTGTACTACATc                                                                                                                     <  1:491011/149‑1 (MQ=255)                          aaCAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTCCCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGt                                                                                                                   <  1:222927/149‑1 (MQ=255)                          aaCAGAGGATAACCGCGCATGGCGCTTGGCAACCCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTTCCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGt                                                                                                                   <  1:222794/149‑1 (MQ=255)                           aCAGCGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGGCCCGACTCTCGGATGGTTCTTGTCTCATTGCCACAATCATAAGTGCCCATCCAAGAGCACGCTTATTCACCAGGGAGAAAAGGCGGAAACGCTGTACTACATCGta                                                                                                                  >  1:86961/1‑148 (MQ=255)                              gagGATAACCGCGCATGGTGCTTGGCAAGCCGCAAACAGCCCCGACCATCGAATGGTTCTTGTCTCATTGCCACATTCATAGTTACGCATCCAAGAGCACGCTTATTCACCAGGGGGAAAAAGCGGAAACGCTGTACTACATCGTTaaa                                                                                                               <  2:133129/149‑1 (MQ=255)                                     aCCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCGGTGCTACATCGTTAAAGGCTCtg                                                                                                        >  2:549615/1‑149 (MQ=255)                                      ccGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCtgt                                                                                                       <  1:353834/149‑1 (MQ=255)                                        gcgcATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTgg                                                                                                     <  2:396688/149‑1 (MQ=255)                                              tgggCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGc                                                                                               <  1:206247/146‑1 (MQ=255)                                              ggTGCGTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGc                                                                                               <  1:348426/149‑1 (MQ=255)                                                   ttGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATc                                                                                          >  1:250179/1‑149 (MQ=255)                                                       cAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCGGATCAA‑Ag                                                                                     >  2:603469/1‑149 (MQ=255)                                                               aaaCAGACCCGACTGTCGAATGGTTATTGTCTCGTTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAGAGCGGAATCGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAA‑AGACGAagag                                                                             <  1:477811/149‑1 (MQ=255)                                                               aaaCAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTCCATCGTAAAAGGCTCTGTGGCAGTGCTGACCAA‑AGCCGGagag                                                                             >  2:56876/1‑149 (MQ=255)                                                               aaaCAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAA‑AGACGAagag                                                                             <  2:14033/149‑1 (MQ=255)                                                                          aCTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAA‑AGACGAAGAGGGTAAAGAAAt                                                                  <  1:117408/149‑1 (MQ=255)                                                                          aCTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAA‑AGACGAAGAGGGTAAAGAAAt                                                                  <  1:760945/149‑1 (MQ=255)                                                                          aCTCTCGAATGGTTCTTGTCTCATTGCAACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAA‑AGACGAAGAGGGTAAAGAAAt                                                                  <  2:571795/149‑1 (MQ=255)                                                                              tcGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAA‑AGACGAAGAGGGTAAAGAAATGATc                                                              <  1:580406/149‑1 (MQ=255)                                                                                gAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAA‑AGACAAAGAGGGTAAAGAAATGATccg                                                            >  1:600486/1‑148 (MQ=255)                                                                                 aaTGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAA‑AGACGAAGAGGGTAAAGAAATGATCctc                                                           >  1:298076/1‑149 (MQ=255)                                                                                   tggaccTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACAGCGTTAAAGGCTCTGTGGCAGTGCTGATCAA‑AGACGAAGAGGGTAAAGAAATGATCCtctc                                                         <  1:97884/144‑1 (MQ=255)                                                                                           gTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAA‑AGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGa                                                 >  1:633900/1‑149 (MQ=255)                                                                                            tctcATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAA‑AGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGaa                                                >  1:630321/1‑149 (MQ=255)                                                                                            tctcATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAA‑AGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGaa                                                >  1:749990/1‑149 (MQ=255)                                                                                            tctcATTGCCACATTCATAAGTACCCAGCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAA‑AGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGaa                                                <  1:763247/149‑1 (MQ=255)                                                                                             ctcATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAA‑AGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAAt                                               >  1:122186/1‑149 (MQ=255)                                                                                              tcATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTTTTTATTTTGAAAAAGCGGAAACGCCGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAA‑AGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATc                                              >  2:777541/1‑149 (MQ=255)                                                                                              tcATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAA‑AGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATc                                              >  2:692569/1‑149 (MQ=255)                                                                                              tcATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAA‑AGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATc                                              >  2:143257/1‑149 (MQ=255)                                                                                                  tGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAA‑AGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAggg                                          >  1:739434/1‑149 (MQ=255)                                                                                                    ccACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAA‑AGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTg                                        >  1:540347/1‑149 (MQ=255)                                                                                                      acaTTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAA‑AGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGAt                                      >  1:630801/1‑149 (MQ=255)                                                                                                      acaTTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAA‑AGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGAt                                      >  1:91470/1‑149 (MQ=255)                                                                                                      acaTTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAA‑AGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGAt                                      >  1:128150/1‑149 (MQ=255)                                                                                                       cagtCATAAGTACCCACCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCACTGCCGTCGTGCTG‑TCAATACAATTAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGAtt                                     <  1:186238/146‑1 (MQ=255)                                                                                                          tCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAA‑AGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTa                                  >  1:559759/1‑149 (MQ=255)                                                                                                            ctAAGCACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAA‑AGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCCGAATCAGGGTGATTTTAtt                                <  2:534009/148‑1 (MQ=255)                                                                                                            aTAAGTACCCATCCAAGAGCACGCTTCTTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAA‑CGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTAtt                                <  2:515730/149‑1 (MQ=255)                                                                                                            aTAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAA‑AGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTAtt                                <  2:413604/149‑1 (MQ=255)                                                                                                                    ccATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAA‑AGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGGGAACt                        >  1:387027/1‑149 (MQ=255)                                                                                                                      aTCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAA‑AGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTgg                      >  1:84372/1‑149 (MQ=255)                                                                                                                      aTCCAAGAGCACGCTTATTCACCAGGGTGAAAAAACGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAA‑AGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTgg                      >  2:436503/1‑149 (MQ=255)                                                                                                                        ccAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAA‑AGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGc                    >  1:444869/1‑149 (MQ=255)                                                                                                                              gCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAA‑AGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGttt              <  1:551107/149‑1 (MQ=255)                                                                                                                                 cgcgTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAA‑AGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGaa           <  2:569712/149‑1 (MQ=255)                                                                                                                                 cacTTATACACCAGGGGGAAAAAGCGGAAACGCTGTACTACATCGTGCAAGGCTCTGTGGCAGTGCTGATCAA‑AGACGAAGTGTGTAAAGAAATGATCCTCTCCTATCTGAATCAGGTTGATTTTATTGGCGAACTGGGCCTGTTTGaa           <  2:356022/147‑1 (MQ=255)                                                                                                                                  gCTTTTTCACCCTGGTGAATAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCCGTGCTGATCAA‑AGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGAGTTTATTGGCGAACTGGGCCTGTTTGCag          >  1:147789/1‑149 (MQ=255)                                                                                                                                  gCTTTTTCACCCTGGTGAATAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAA‑AGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTCGATCAGGGTGATGTTATTGGCGAACTGGGCCTGTTTGAag          >  1:18893/1‑149 (MQ=255)                                                                                                                                    ttATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAA‑AGACGAAGAGGGTAAAGAAATGATCCTCTGCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGCagcg        >  1:103980/1‑147 (MQ=255)                                                                                                                                    ttATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAA‑AGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAagag        >  1:621897/1‑149 (MQ=255)                                                                                                                                          aCCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAA‑AGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAg  >  1:74128/1‑149 (MQ=255)                                                                                                                                          aCCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAA‑AGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAg  >  1:365787/1‑149 (MQ=255)                                                                                                                                          aCCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAA‑AGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAg  >  1:197817/1‑149 (MQ=255)                                                                                                                                              |                                                                                                                                                  TTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAA‑AGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAG  >  NC_000913/3486075‑3486360

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