Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	139,295	T→C	6.2%	Q644R (CAG→CGG)	gcd ←	glucose dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	139,295	0	T	C	6.2%	205.1 / 2.1	64	Q644R (CAG→CGG)	gcd	glucose dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (30/30);  new base C (2/2);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAATCGCCACTTCACGATTTGGATCAACGG  >  NC_000913/139156‑139439                                                                                                                                            |                                                                                                                                                 cgTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGa                                                                                                                                          <  3:229950/148‑1 (MQ=255) cgTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGa                                                                                                                                          <  3:119673/148‑1 (MQ=255)     ccAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGc                                                                                                                                      <  2:290606/148‑1 (MQ=255)       aaTACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCg                                                                                                                                    >  8:218370/1‑148 (MQ=255)        aTACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCgg                                                                                                                                   <  8:16383/148‑1 (MQ=255)         tACGTTTCTTCCACACCACTTCATTAGTTTTCCGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGa                                                                                                                                  >  5:271169/1‑148 (MQ=255)                   ccACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTAccc                                                                                                                        >  3:3480/1‑148 (MQ=255)                     acacCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGCATGCCGGATTCCGTACCCGt                                                                                                                      >  6:342/1‑148 (MQ=255)                            ttCATTAGTTTTCAGCTCCCGCGCCGTGATAAAACCCCATTCTGGCTGTTTACATGGCAGACCAATTGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTg                                                                                                               <  5:8056/148‑1 (MQ=255)                                  aGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATc                                                                                                          <  7:84468/147‑1 (MQ=255)                                   gTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCtt                                                                                                        <  4:99573/148‑1 (MQ=255)                                    ttttCAGAGCCAGCGCCGAGATATAACCCCATCCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCttt                                                                                                       <  6:18248/148‑1 (MQ=255)                                     tttCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTc                                                                                                      >  8:33958/1‑148 (MQ=255)                                          gATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGCCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGgctg                                                                                                 >  2:35310/1‑148 (MQ=255)                                           aTCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGgctg                                                                                                 >  4:67406/1‑147 (MQ=255)                                             ccAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATC‑TTCGGCTGCTc                                                                                              >  4:113266/1‑147 (MQ=255)                                                           taACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCa                                                                                 <  6:27497/147‑1 (MQ=255)                                                           taACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAg                                                                                >  4:19891/1‑148 (MQ=255)                                                             aCCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCCTCGTTCGGCTGCTCCATCGGGTTGCCAGGa                                                                              >  2:121938/1‑148 (MQ=255)                                                                ccATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCa                                                                           >  6:78126/1‑148 (MQ=255)                                                                    gctggctgTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATTCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCAcgcg                                                                       >  2:297875/1‑148 (MQ=255)                                                                    gctggctgTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCCTCTTTCTGCTGCTCCATCGGGTTGCCAGGACACcgcg                                                                       >  4:193368/1‑148 (MQ=255)                                                                    gctggctgTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGGCTGCTCCATCGGGt                                                                                      <  5:84119/133‑1 (MQ=255)                                                                    gctggctgTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGGCTGCTCCATCGGGt                                                                                      >  6:84119/1‑133 (MQ=255)                                                                            tttACATGGCAGACCAAATGGTGAGAGCAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGt                                                               >  6:315331/1‑148 (MQ=255)                                                                               aCATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGttt                                                             <  6:271169/147‑1 (MQ=255)                                                                                 aTGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGa                                                          <  4:3480/148‑1 (MQ=255)                                                                                 aTGCCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGa                                                          <  1:35310/148‑1 (MQ=255)                                                                                   ggCAGACCAAATGGTGAGAGGAAAGGGTTGAGCGTGACACCAAACGGTACACAGAACTGTGGCTGAATGCAGGATTCCGTACCCGTGCCTATGGCAAAATTCGGATGCTCCATCGCGTTGCCAGGACCACGCGGGATCAGTTTCTaaa                                                        >  8:5715/1‑148 (MQ=255)                                                                                        aCCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGaaacaaac                                                   <  3:19891/148‑1 (MQ=255)                                                                                         ccgaaTGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACCGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACg                                                  >  5:19557/4‑148 (MQ=255)                                                                                         ccAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACg                                                  >  7:86115/1‑148 (MQ=255)                                                                                         ccAAATGGTGAGAGGAACGGGTTGACCGTGACACCATACGGTACACCGTACCGTGGCTGAATGCCGGATTCCCTACCCGTGCCTTTGGAATCTTTCGGCAGCTCCATCGCCTTGCCAGGCCCACCCGGGATCAGTTTAGAAACAAACg                                                  >  5:19558/1‑148 (MQ=255)                                                                                           aaaTGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACgg                                                 >  1:122979/1‑147 (MQ=255)                                                                                           aaaTGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGc                                                <  4:284293/148‑1 (MQ=255)                                                                                           aaaTGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGAGCAGTTTCGAAACAAACGGc                                                <  3:67406/148‑1 (MQ=255)                                                                                            aaTGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACgg                                                 >  8:182488/1‑146 (MQ=255)                                                                                            aaTGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCa                                               <  8:86115/148‑1 (MQ=255)                                                                                            aaTGGTGAGAGCAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCa                                               <  5:315331/148‑1 (MQ=255)                                                                                                 tGAGAGCAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTTGCAGAATGCCGGATTGCGTACGCGTGCCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGAGCACGCGGGATCAGTTTCGAAACAAACGGCAGTGcc                                          <  7:5715/148‑1 (MQ=255)                                                                                                   agagGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCAt                                        >  2:127397/1‑148 (MQ=255)                                                                                                      ggAAAGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGACGGATTCCGTACACGTGCCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGACAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTgg                                     >  8:289995/1‑148 (MQ=255)                                                                                                       gAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACg                                                  <  7:263830/134‑1 (MQ=255)                                                                                                        aaCGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGt                                   >  1:307018/1‑148 (MQ=255)                                                                                                        aaCGGATTGAGCGTGACACCATACGGTACACAGTACTCTGGCTGAATGCCGGATTCAGTACCCGTGCCTTTGGCATCTTTCGGCTGCTCCATCGGTTTGCCAGCACCACGCGGGATCAGTTTCGCAACAAACgg                                                 >  8:263830/1‑134 (MQ=255)                                                                                                          cGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTg                                 <  7:33958/148‑1 (MQ=255)                                                                                                          cGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTg                                 >  6:17430/1‑148 (MQ=255)                                                                                                             gTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGc                               >  7:116596/1‑147 (MQ=255)                                                                                                             gTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATC‑TTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCaa                             <  3:113266/148‑1 (MQ=255)                                                                                                              ttGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCaa                             >  7:271650/1‑148 (MQ=255)                                                                                                                   cGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAATCGcc                        <  8:116596/148‑1 (MQ=255)                                                                                                                       acagCATACGTTTCACAGTACCGTGGCTGAATGCCGGATTCCGTACCCGTGGCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCTCGCGGGTTCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAATCGCCACtt                    <  6:19558/148‑1 (MQ=255)                                                                                                                       acacCATACGGTACACCGTACCGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAATCGCCACtt                    <  6:19557/148‑1 (MQ=255)                                                                                                                          ccATACGGTACACCGTACTGGGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAATCGCCACTTCAc                 <  7:218370/148‑1 (MQ=255)                                                                                                                             tACGGTACACCGTACTGTAGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCATGCCCACGCGGGATCAGTTTCGAAACAAACGTCAGTGCCATTGGGTTGGCAATCGCCCCTTAACGAt              >  6:103863/1‑148 (MQ=255)                                                                                                                                 gTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAATCGCCACTTCACGATTTgg          <  2:307018/148‑1 (MQ=255)                                                                                                                                  tACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAATCGCCACTTCACGATTTGGa         >  3:237068/1‑148 (MQ=255)                                                                                                                                      ccGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGtt                                  <  7:102973/119‑1 (MQ=255)                                                                                                                                      ccGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGtt                                  >  8:102973/1‑119 (MQ=255)                                                                                                                                      ccGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAATCGCCACTTCACGATTTGGATCaa     <  2:122979/148‑1 (MQ=255)                                                                                                                                      ccGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAATCGCCACTTCACGATTTGGATCaa     <  8:271650/148‑1 (MQ=255)                                                                                                                                      ccGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAATCGCCACTTCACGATTTGGATCaa     <  5:78126/148‑1 (MQ=255)                                                                                                                                       cGTTCTGTGGCTGAATGCCGGATTCCGTACCCGGGCCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAATCGCCACTTCACGATTTGGATCAAc    <  5:342/148‑1 (MQ=255)                                                                                                                                          aCTGTGGATGAATGCCGGATTCCGTACCAGTGCCTTTGGCATCTTTAGGCTGCTCCATCGGGTTGCCAGGACAACGCGGGATCAGTTTAGAAACAAACGGCAGTGACATTGGGTTGGCAATCGCCCCTTAACGATTTGGATCAACgg  >  7:140593/1‑147 (MQ=255)                                                                                                                                            |                                                                                                                                                 CGTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAATCGCCACTTCACGATTTGGATCAACGG  >  NC_000913/139156‑139439

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