Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,723,515	T→G	6.5%	K17Q (AAG→CAG)	htpG ←	chaperone protein HtpG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,723,515	0	T	G	6.5%	191.7 / 3.6	61	K17Q (AAG→CAG)	htpG	chaperone protein HtpG
Reads supporting (aligned to +/- strand):  ref base T (26/31);  new base G (2/2);  total (28/33)
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.91e-01

GGTCAGCATCGCCCTCGAACAGCTCTGGCTTGGCCAGGGCTTCGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTCGAAACTCCGATCTGATGGCGTTGGCGCCTGGCGGCCGGTTCGGCCAGCGTTAACGGCGGATGTCATGCAGATGGGGGCTGCCCGGATGATT  >  NC_002947/4723368‑4723661                                                                                                                                                    |                                                                                                                                                   ggTCAGCATCGCCCTCGAACAGCTCTGGCTTGGCCAGGGCTTCGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTc                                                                                                                                                   <  1:56795/149‑1 (MQ=255)  gTCAGCATCGCCCTCGAACAGCTCTGGCTTGGCCAGGGCTTCGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCa                                                                                                                                                  >  8:166038/1‑149 (MQ=255)       cATCGCCCTCGAACAGCTCTGGCTTGGCCAGGGCTTCGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGAAGCAACTGCTGCACCTca                                                                                                                                             >  1:219086/1‑148 (MQ=255)         tCGCCCTCGAACAGCTCTGGCTTGGCCAGGGCTTCGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGAAGCAACTGCTGCACCTcatg                                                                                                                                           <  2:219086/149‑4 (MQ=255)           gCCCTCGAACAGCTCTGGCTTGGCCAGGGCTTCGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTc                                                                                                                                          <  2:178597/148‑1 (MQ=255)           gCCCTCGAACAGCTCTGGCTTGGCCAGGGCTTCGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCt                                                                                                                                         <  6:176484/149‑1 (MQ=255)                 gAACAGCTCTGGCTTGGCCAGGGCTTCGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAg                                                                                                                                    <  2:18959/148‑1 (MQ=255)                 gAACAGCTCTGGCTTGGCCAGGGCTTCGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGc                                                                                                                                   >  7:28349/1‑149 (MQ=255)                     aGCTCTGGCTTGGCCAGGGCTTCGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAg                                                                                                                               >  4:166848/1‑149 (MQ=255)                     aGCTCTGGCTTGGCCAGGGCTTCGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAg                                                                                                                               <  8:229954/149‑1 (MQ=255)                          tGGCTTGGCCAGGGCTTCGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGttt                                                                                                                          <  7:166038/149‑1 (MQ=255)                                 gCCAGGGCTTCGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCtttt                                                                                                                    <  5:298489/148‑1 (MQ=255)                                 gCCAGGGCTTCGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTtg                                                                                                                   <  5:243016/149‑1 (MQ=255)                                 gCCAGGGCTTCGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTtg                                                                                                                   >  6:243016/1‑149 (MQ=255)                                  ccAGGGCTTCGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTtgt                                                                                                                  <  6:255368/149‑1 (MQ=255)                                  ccAGGGCTTCGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTtgt                                                                                                                  <  3:166848/149‑1 (MQ=255)                                    aGGGCTTCGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTgtgt                                                                                                                >  5:49938/1‑149 (MQ=255)                                     gggCTTCGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGtt                                                                                                               <  8:165564/149‑1 (MQ=255)                                     gggCTTCGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGtt                                                                                                               >  4:101983/1‑149 (MQ=255)                                      ggCTTCGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGttt                                                                                                              <  2:43514/149‑1 (MQ=255)                                        cTTCGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGTTTc                                                                                                             >  5:284090/1‑148 (MQ=255)                                        cTTCGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGTTTc                                                                                                             >  1:116072/1‑148 (MQ=255)                                           cGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCtttt                                                                                                                    <  2:167263/138‑1 (MQ=255)                                           cGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCtttt                                                                                                                    >  1:167263/1‑138 (MQ=255)                                               aCGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGTTTCCACAGTCa                                                                                                     >  2:244332/1‑149 (MQ=255)                                                 gCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGTTTCCACAGTCATg                                                                                                   <  3:101983/149‑1 (MQ=255)                                                  cAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCt                                                                                                                       >  6:58699/1‑128 (MQ=255)                                                           ggcggcGTCGGAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCg                                                                                                                                             <  2:110849/97‑1 (MQ=255)                                                           ggcggcGTCGGAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCg                                                                                                                                             >  1:110849/1‑97 (MQ=255)                                                            gcggcgTCGGAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTCGAAAc                                                                                        >  1:116169/1‑149 (MQ=255)                                                               gcgTCGGAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTCGAAACTcc                                                                                     >  3:263356/1‑149 (MQ=255)                                                                 gTCGGAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCGTTTTGTGTTTCCACAGTCATGGTCTTCGAAACTCCGa                                                                                   <  8:28349/149‑1 (MQ=255)                                                                    ggAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTtgt                                                                                                                  <  8:281075/115‑1 (MQ=255)                                                                    ggAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTtgt                                                                                                                  >  7:281075/1‑115 (MQ=255)                                                                        gCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTCGAAACTCCGATCTGATg                                                                            >  5:69435/1‑149 (MQ=255)                                                                           ttGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGTTTCCACAGTCATGCTCTTCGAAACTCCGATCTGATGGCg                                                                         >  1:66555/1‑149 (MQ=255)                                                                              gAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTCGAAACTCCGATCTGATGGCGTTg                                                                      <  4:84297/149‑1 (MQ=255)                                                                                aaTCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTCGAAACTCCGATCTGATGGCGTTGgc                                                                    <  6:69435/149‑1 (MQ=255)                                                                                  tCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTCGAAACTCCGATCTGATGGCGTTGgcgc                                                                  <  2:116072/149‑1 (MQ=255)                                                                                    aGCTCCCGCAGGAAGATTTCCTTGTTCGAGTATAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTCGAAACTCCGATCTGATGGCGTTGGCGCCt                                                                <  7:213771/149‑1 (MQ=255)                                                                                           gCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTCGAAACTCCGATCTGATGGCGTTGGCGCCTggcggc                                                          >  1:159319/1‑148 (MQ=255)                                                                                            cAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTCGAAACTCCGATCTGATGGCGTTGGCGCCTGGCGGCCg                                                        <  6:284090/149‑1 (MQ=255)                                                                                             aGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGTTTCCACAGTCATGCTCTTCGAAACTCCGATCTGATGGCGTTGGCGCCTGGCGGCCgg                                                       <  2:66555/149‑1 (MQ=255)                                                                                               gAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTCGAAACTCCGATCTGATGGCGTTGGCGCCTGGCGGCCGGtt                                                     <  1:244332/149‑1 (MQ=255)                                                                                                 aGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTCGAAACTCCGATCTGATGGCGTTGGCGCCTGGCGGCCGGTTCg                                                   <  2:116169/149‑1 (MQ=255)                                                                                                    tttCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTCGAAACTCCGATCTGATGGCGTTGGCGCCTGGCGGCCGGTTCGGc                                                 <  8:157720/148‑1 (MQ=255)                                                                                                        cTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTCGAAACTCCGATCTGATGGCGTTGGCGCCTGGCGGCCGGTTCGGCCAGCg                                            >  5:256788/1‑149 (MQ=255)                                                                                                             tCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTCGAAACTCCGATCTGATGGCGTTGGCGCCTGGCGGCCGGTTCGGCCAGCGTTaa                                        <  6:256788/148‑1 (MQ=255)                                                                                                               gAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTCGAAACTCCGATCTGATGGCGTTGGCGCCTggc                                                             >  7:113921/1‑125 (MQ=255)                                                                                                               gAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTCGAAACTCCGATCTGATGGCGTTGGCGCCTggc                                                             <  8:113921/125‑1 (MQ=255)                                                                                                                  tACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTCGAAACTCCGATCTGATGGCGTTGGCGCCTGGCGGCCGGTTCGGCCAGCGTTAAcggcg                                   >  1:288050/1‑148 (MQ=255)                                                                                                                     aGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTCGAAACTCCGATCTGATGGCGTTGGCGCCTGGCGGCCGGTTCGGCCAGCGTTAACGGCGGATg                               >  3:15328/1‑149 (MQ=255)                                                                                                                               gaccgAGGTGAAGCAACTGCTGCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTCGAAACTCCGATCTGATGGCGTTGGCGCCTGGCGGCCg                                                        >  7:151596/5‑114 (MQ=39)                                                                                                                               gaccgAGGTGAAGCAACTGCTGCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTCGAAACTCCGATCTGATGGCGTTGGCGCCTGGCGGCCg                                                        <  8:151596/110‑1 (MQ=39)                                                                                                                                cATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTCGAAACTCCGATCTGATGGCGTTGGCGCCTGGCGGCCGGTTCGGCCAGCGTTAACGGCGGATGTCATGCAGATg                    >  4:148849/1‑149 (MQ=255)                                                                                                                                 aTGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTCGAAACTCCGATCTGATGGCGTTGGCGCCTGGCGGCCGGTTCGGCCAGCGTTAACGGCGGATGTCATGCAGATgg                   <  4:15328/149‑1 (MQ=255)                                                                                                                                  tGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTCGAAACTCCGATCTGATGGCGTTg                                                                      <  7:108441/97‑1 (MQ=255)                                                                                                                                  tGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTCGAAACTCCGATCTGATGGCGTTg                                                                      >  8:108441/1‑97 (MQ=255)                                                                                                                                      gTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTCGAAACTCCGATCTGATGGCGTTGGCGCCTGGCGGCCGGTTCGGCCAGCGTTAACGGCGGATGTCATGCAGATGGGGGCt              >  4:244383/1‑149 (MQ=255)                                                                                                                                                 gattCACCTCGGTCTGTAAGCCCAGCGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTCGAAACTCCGATCTCATGGCGTTGGCGCCTGCCGGCCGGTTCCCCCAGCGTAAACGGCGGAAGTCAAGCAGATGGGGGATGCCCGgatgat   >  6:23349/3‑149 (MQ=255)                                                                                                                                                  cTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTCGAAACTCCGATCTGATGGCGTTGGCGCCTGGCGGCCGGTTCGGCCAGCGTTAACGGCGGATGTCATGCAGATGGGGGCTGCCCGGATGAtt  <  3:244383/149‑1 (MQ=255)                                                                                                                                                  cTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTCGAAACTCCGATCTGATGGCGTTGGCGCCTGGCGGCCGGTTCGGCCAGCGTTAACGGCGGATGTCATGCAGATGGGGACTGCCCGGATGAtt  >  4:180053/1‑149 (MQ=255)                                                                                                                                                   ttCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTCGAAACTCCGATCTGATGGCGTTGGCGCCTGGCGGCCGGTTCGGCCAGCGTTAACGGCGGATGTCATGCAGATGGGGGCTGCCCGGATGAtt  >  3:238855/1‑148 (MQ=255)                                                                                                                                                    |                                                                                                                                                   GGTCAGCATCGCCCTCGAACAGCTCTGGCTTGGCCAGGGCTTCGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCAGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCTGGAAGCCCAGCGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTCGAAACTCCGATCTGATGGCGTTGGCGCCTGGCGGCCGGTTCGGCCAGCGTTAACGGCGGATGTCATGCAGATGGGGGCTGCCCGGATGATT  >  NC_002947/4723368‑4723661

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