Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,391,321	G→C	6.4%	V236V (GTC→GTG)	tldD ←	putative peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,391,321	0	G	C	6.4%	207.1 / 4.1	62	V236V (GTC→GTG)	tldD	putative peptidase
Reads supporting (aligned to +/- strand):  ref base G (31/27);  new base C (2/2);  total (33/29)
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.80e-01

CCTGTCCACTAAATACTGAAGTGCCACGGCGGTTGAAGTCGCCTTCCAGACCGTGACCAACCGCTTCATGCAACAGCACGCCCGGCCAACCTGCGCCAAGTACTACCGGCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGAACG  >  NC_000913/3391176‑3391465                                                                                                                                                  |                                                                                                                                                 ccTGTCCACTAAATACTGAAGTGCCACGGCGGTTGAAGTCGCCTTCCAGACCGTGACCAACCGCTTCATGCAACAGCACGCCCGGCCAACCTGCGCCAAGTACTACCGGCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGAcc                                                                                                                                               <  8:66095/149‑1 (MQ=255)  cTGTCCACTAAATCCTGAAGTGCCACGGCGGTTGAAGTCGCCTTCCAGCCCGTGACCAACCGCTTCATGCAACAGCACGCCCGGCCAACCTGCGCCAAGTACTACCGGCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCa                                                                                                                                              <  5:317866/149‑1 (MQ=255)  cTGTCCACTAAATACTGAAGTGCCACGGCGGTTGAAGTCGCCTTCCAGACCGTGACCAACCGCTTCATGCAACAGCACGCCCGGCCAACCTGCGCCAAGTACTACCGGCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCa                                                                                                                                              <  2:243988/149‑1 (MQ=255)   tGTCCACTAAATAATGAAGTGCCACGGCGGTTGAAGTCGCCTTCCAGACCGTGACAAACCGCTTCATGCAACAGCACGCCCGGCCAACCTGCGCCAAGTACTACCGGCATGGTGCCCGCTGGTGAACCAACGGCAGAAAGATTGACCac                                                                                                                                             >  3:285752/1‑148 (MQ=255)     tCCACTAAATACTGAAGTGCCACGGCGGTTGAAGTCGCCTTCCAGACCGTGACCAACCGCTTCATGCAACAGCACGCCCGGCCAACCTGCGCCAAGTACTACCGGCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAgcg                                                                                                                                           >  3:77357/1‑149 (MQ=255)       cACTAAATACTGAAGTGCCACGGCGGTTGAAGTCGCCTTCCAGACCGTGACCAACCGCTTCATGCAACAGCACGCCCGGCCAACCTGCGCCAAGTACTACCGGCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAgcgc                                                                                                                                          >  1:365491/1‑148 (MQ=255)        aCTAAATACTGAAGTGCCACGGCGGTTGAAGTCGCCTTCCAGACCGTGACCAACCGCTTCATGCAACAGCACGCCCGGCCAACCTGCGCCAAGTACTACCGGCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCGCCa                                                                                                                                        >  7:115994/1‑149 (MQ=255)        aCTAAATACTGAAGTGCCACGGCGGTTGAAGTCGCCTTCCAGACCGTGACCAACCGCTTCATACAACAGCACGCCCGGCCAACCTGCGCCAAGTACTACCGGCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCGCCa                                                                                                                                        <  8:273914/149‑1 (MQ=255)          tAAATACTGAAGTGCCACGGCGGTTGAAGTCGCCTTCCAGACCGTGACCAACCGCTTCATGCAACAGCACGCCCGGCCAACCTGCGCCAAGTACTACCGGCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCGCCATa                                                                                                                                      >  3:147761/1‑149 (MQ=255)           aaaTACTGAAGTGCCACGGCGGTTGAAGTCGCCTTCCAGACCGTGACCAACCGCTTCATGCAACAGCACGCCCGGCCAACCTGCGCCAAGTACTACCGGCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTCACCAGCGCCATAc                                                                                                                                     >  7:112992/1‑149 (MQ=255)           aaaTACTGAAGTGCCACGGCGGTTGAAGTCGCCTTCCAGACCGTGACCAACCGCTTCATGCAACAGCACGCCCGGCCAACCTGCGCCAAGTACTACCGGCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTCACCAGCGCCATAc                                                                                                                                     >  7:107314/1‑149 (MQ=255)                          aCGGCGGTTGAAGTCGCCTTCCAGACCGTGACCAACCGCTTCATGCAACAGCACGCCCGGCCGACCTGCGCCAAGTACTACCGGCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTg                                                                                                                      <  7:97227/149‑1 (MQ=255)                              cggTTGAAGTCGCCTTCCAGACCGTGACCAACCGCTTCATGCAACAGCACGCCCGGCCAACCTGCGCCAAGTACTACCGGCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCa                                                                                                                  <  7:135657/149‑1 (MQ=255)                              cggTTGAAGTCGCCTTCCAGACCGTGACCAACCGCTTCATGCAACAGCACGCCCGGCCAACCTGCGCCAAGTACTACCGGCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCa                                                                                                                  >  7:231671/1‑149 (MQ=255)                               ggTTGAAGTCGCCTTCCAGACCGTGACCAACCGCTTCATGCGACAGCACGCCCGGCCAACCTGCGCCAAGTACTACCGGCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCAt                                                                                                                 <  4:13958/149‑1 (MQ=255)                                 ttGAAGTCGCCTTCCAGACCGTGACCAACCGCTTCATGCAACAGCACGCCCGGCCAACCTGCGCCAAGTACTACCGGCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGc                                                                                                               <  5:29596/149‑1 (MQ=255)                                    aaGTCGCCTTCCAGCCCGTGACCAACCGCGTCATGCATCAGCACGCCCGGCCCACCTGCGCACAGTACTACCGGCAAGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCAt                                                                                                             <  5:97882/148‑1 (MQ=255)                                       tCGCCTTCCAGACCGTGACCAACCGCTTCATGCAACAGCACGCCCGGCCAACCTGCGCCAAGTACTACCGGCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCAt                                                                                                             >  4:328283/1‑145 (MQ=255)                                       tCGCCTTCCAGACCGTGACCAACCGCTTCATGCAACAGCACGCCCGGCCAACCTGCGCCAAGTACTACCGGCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCAt                                                                                                             <  3:328283/145‑1 (MQ=255)                                          ccTTCCAGACCGTGACCAACCGCTTCATGCAACAGCACGCCCGGCCAACCTGCGCCAAGTACTACCGGCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACg                                                                                                      <  6:50200/149‑1 (MQ=255)                                           cTTCCAGACAGTGACCAACCGCTTCATGCAACAGCACGCCCGGCCAACCTGCGCCAAGTACTACCGGCATCGTGACCGCTGCAGCAGCAACTGCAGAAAGATTGCCCAGCGCCATACGCACTGCTTCATTTGCCCATGCATCCTCACg                                                                                                      >  4:300462/1‑148 (MQ=255)                                              ccAGACCGTGACCAACCGCTTCATGCAACAGCACGCCCGGCCAACCTGCGCCAAGTACTACCGGCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACt                                                                                                  >  3:282744/1‑149 (MQ=255)                                              ccAGACCGTGACCAACCGCTTCATGCAACAGCACGCCCGGCCAACCTGCGCCAAGTACTACCGGCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCCCCATACGCACAGCGTCTTTTGCCCATGCATCCGCACGCACt                                                                                                  >  3:330620/1‑149 (MQ=255)                                                aGACCGTGACCAACCGCTTCCTGCAACAGCACGCCCGGCCAACCTGCGCCAAGTACTACCGGCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTCACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTc                                                                                                <  8:112992/149‑1 (MQ=255)                                                aGACCGTGACCAACCGCTTCATGCAACAGCACGCCCGGCCAACCTGCGCCAAGTACTACCGGCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTCACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTc                                                                                                <  8:107314/149‑1 (MQ=255)                                                      tGACCAACCGCTTCATGCAACAGCACGCCCGGCCAACCTGCGCCAAGTACTACCGGCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTc                                                                                          <  2:15210/149‑1 (MQ=255)                                                      tGACCAACCGCTTCATGCAACAGCACGCCCGGCCAACCTGCGCCAAGTACTACCGGCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTc                                                                                          <  8:115994/149‑1 (MQ=255)                                                           aaCCGCTTCATGCAACAGCACGCCCGGCCAACCTGCGCCAAGTACTACCGGCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGAt                                                                                     <  4:77357/149‑1 (MQ=255)                                                           aaCCGCTTCATGCAACAGCACGCCCGGCCAACCTGCGCCAAGTACTACCGGCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGAt                                                                                     <  4:86767/149‑1 (MQ=255)                                                                   cATGCAACAGCACGCCCGGCCAACCTGCGCCAAGTACTACCGGCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAg                                                                              <  8:231671/148‑1 (MQ=255)                                                                        aaCAGCACGCCCGGCCAACCTGCGCCAAGTACTACCGGCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGgaagaa                                                                        >  4:11085/1‑149 (MQ=255)                                                                         aCAGCACGCCCGGCCAACCTGCGCCAAGTACTACCGGCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAAt                                                                       <  2:321931/149‑1 (MQ=255)                                                                         aCAGCACGCCCGGCCAACCTGCGCCAAGTACTACCGGCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAAt                                                                       <  3:300462/149‑1 (MQ=255)                                                                              aCGACCGGCCAACCTGCGCCAAGTACTACCGGCATGGTGCCAGCTGGTTCAGCAACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACTGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATa                                                                  <  4:285752/149‑1 (MQ=255)                                                                                gCCCGGCCAACCTGCGCCAAGTACTACCGGCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCCGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATCAc                                                                >  3:2943/1‑149 (MQ=255)                                                                                  ccGGCCAACCTGCGCCAAGTACTACCGGCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAAcc                                                               >  4:55876/1‑148 (MQ=255)                                                                                          ccTGCGCCAAGTACTACCGGCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGAcc                                                      <  2:365491/149‑1 (MQ=255)                                                                                          ccTGCGCCAAGTACTACCGGCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGAcc                                                      >  5:338351/1‑149 (MQ=255)                                                                                            tGCGCCAAGTACTACCGGCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGAccgc                                                    >  1:364461/1‑149 (MQ=255)                                                                                               gcCAAGTACTACCGGCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGAccgccgc                                                 >  1:128719/1‑149 (MQ=255)                                                                                                          ccGGCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCa                                      >  6:28680/1‑149 (MQ=255)                                                                                                          ccGGCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGCACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCa                                      >  6:70579/1‑149 (MQ=255)                                                                                                               aTGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACcgcg                                 >  5:239862/1‑149 (MQ=255)                                                                                                                tGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCaa                                                             >  8:49475/1‑120 (MQ=255)                                                                                                                tGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCaa                                                             <  7:49475/120‑1 (MQ=255)                                                                                                                tGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACcgcgt                                >  7:293522/1‑149 (MQ=255)                                                                                                                tGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACcgcgt                                >  4:151506/1‑149 (MQ=255)                                                                                                                  gTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACcgcgttc                              <  4:282744/149‑1 (MQ=255)                                                                                                                   tGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGATATCATAACCAAAACGACCGCCGCCGCCACTGGCACcgcgttcg                             >  2:102607/1‑149 (MQ=255)                                                                                                                   tGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCa                                            <  6:34564/134‑1 (MQ=255)                                                                                                                   tGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCa                                            >  5:34564/1‑134 (MQ=255)                                                                                                                       cGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACcgcgttcgcgtt                         >  2:343202/1‑149 (MQ=255)                                                                                                                         cTGGTGCAGCAACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCTAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTg                       <  3:11085/149‑1 (MQ=255)                                                                                                                               cagcaACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCAtct                 >  6:177436/1‑149 (MQ=255)                                                                                                                                agcaACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCAtctt                <  1:82495/149‑1 (MQ=255)                                                                                                                                 gcaACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCAtcttc               >  4:108289/1‑149 (MQ=255)                                                                                                                                  caACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCAtcttct              >  4:166372/1‑149 (MQ=255)                                                                                                                                  caACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCAtcttc               >  6:69748/1‑148 (MQ=255)                                                                                                                                   aaCGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATcttctt             <  7:331120/149‑1 (MQ=255)                                                                                                                                         aGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTcgacga       <  5:70579/149‑1 (MQ=255)                                                                                                                                             aGATTGACCAGCGCCATACGCTCTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGAAc   >  7:141416/1‑149 (MQ=255)                                                                                                                                              gATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGAACg  >  6:228896/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                 CCTGTCCACTAAATACTGAAGTGCCACGGCGGTTGAAGTCGCCTTCCAGACCGTGACCAACCGCTTCATGCAACAGCACGCCCGGCCAACCTGCGCCAAGTACTACCGGCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGAATTCATAACCAAAACGACCGCCGCCGCCACTGGCACCGCGTTCGCGTTTGCCATCTTCTTCGACGAGAACG  >  NC_000913/3391176‑3391465

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: ‑