Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	496,382	C→G	S421S (TCC→TCG)	htpG →	chaperone protein HtpG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	496,382	0	C	G	92.1%	108.6 / 3.2	38	S421S (TCC→TCG)	htpG	chaperone protein HtpG
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base G (20/15);  total (21/17)
Fisher's exact test for biased strand distribution p-value = 5.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.13e-01

AGCACGGTAACGCGTAACCTGCGCAATGCGCTGACCAAGCGTGTGCTGCAAATGCTGGAAAAACTGGCGAAAGACGACGCGGAAAAATACCAGACCTTCTGGCAACAGTTTGGCCTGGTACTGAAAGAAGGTCCGGCGGAAGATTTCGCTAACCAGGAAGCGATCGCCAAACTGCTGCGTTTTGCTTCTACCCATACCGATTCTTCTGCGCAGACCGTATCTCTGGAAGACTACGTTTCCCGCATGAAAGAAGGGCAGGAGAAAATCTACTACATCACCGCAGACAGCTATGCGGCAGCGAAGAGCAGCCCGCACCTGGAACTGCTGCGTAAGAAAGGCATCGAAGTTCTGCTGCTTTCCGACCGCATCGATGAGTGGATGATGAACTATCTGACTGAGTTCGACGGTAAACCGTTCCAGTCGGTGTCTAAAGTTGACGAGTCGCTTGAAAAACTGGCTGACGAAGTTGATGAGAGCGCGAAAG  >  NC_000913/496143‑496626                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                     aGCACGGTAACGCGTAACCTGCGCAATGCGCTGACCAAGCGTGTGCTGCAAATGCTGGAAAAACTGGCGAAAGACGACGCGGAAAAATACCAGACCTTCTGGCAACAGTTTGGCCTGGTACTGAAAGAAGGTCCGGCGGAAGATTTCGCTAACCAGGAAGCGATCGCCAAACTGCTGCGTTTTGCTTCTACCCATACCGATTCTTCTGCGCAGACCGTATCTCTGGAAGACTACGTTTCGCGCATGAaaga                                                                                                                                                                                                                                           <  1:135808/251‑1 (MQ=255)                     gcgcAATGCGCTGACCAAGCGTGTGCTGCAAATGCTGGAAAAACTGGCGAAAGACGACGCGGAAAAATACCAGACCTTCTGGCAACAGTTTGGCCTGGTACTGAAAGAAGGTCCGGCGGAAGATTTCGCTAACCAGGAAGCGATCGCCAAACTGCTGCGTTTTGCTTCTACCCATACCGATTCTTCTGCGCAGACCGTATCTCTGGAAGACTACGTTTCCCGCATGAAAGAAGGGCAGGAGAAAATctac                                                                                                                                                                                                                        <  1:113778/250‑1 (MQ=255)                      cgcAATGCGCTGACCAAGCGTGTGCGGCAAATGCTGGAAAAACTGGCGAAAGACGACGCGGAAAAATACCAGACCTTCTGGCAACAGTTTGGCCTGGTACTGAAAGAAGGTCCGGCGGAAGATGTCGCTAACCAGGAAGCGATCGCCAAACTGCTGCGTTTTGCTTCTACCCATACCGATTCTTCTGCGCAGACCGTATCTCTGGAAGACTACGTTTCGCGCATGAAAGAAGGGCAGGAGAAAATctacta                                                                                                                                                                                                                      <  1:293729/251‑1 (MQ=255)                      cgcAATGAGATGACAAAGCGTGAGCTGCAACTGCTGTCCAAACTGGAGAAAGCCGACGCGGAAAAATACCAGACCTTCTGGCAACAGGTTTGCCTGGTACTGAAAGAAGTTCCGACGGTAGATTTCGCTAACCCGGAAGCGCTCTCCACACTGCTGGGTTTTGCTTCTCCCCATCACGATTCTTTTTCGCAGACACTCTCTCTGGTATACTACGTTTCGTGCATGAATAACGGGCAGGAGAATATccactc                                                                                                                                                                                                                      >  2:280006/1‑250 (MQ=255)                       gcAATGCGCTGACCAAGCGTGGGCTGCAAATGCTGGAAAAACTGGCGAAAGACGACGCGGAAAAATACCAGACCTTCTGGCAACAGTTTGGCCTGGTACTGAAAGAAGGTCCGGCGGAAGATTTCGCTAACCAGGAAGCGATCGCCAAACTGCTGCGTTTTGCTTCTACCCATACCGATTCTTCTGCGCAGACCGTATCTCTGGAAGACTACGTTTCGCGCATGAAAGAAGGGCAGGAGAAAATCtactac                                                                                                                                                                                                                     <  1:181879/251‑1 (MQ=255)                          aTGCGCTGACCAAGCGTGTGCTGCAAATGCTGGAAAAACTGGCGAAAGACGACGCGGAAAAATACCAGACCTTCTGGCAACAGTTTGGCCTGGTACTGAAAGAAGGTCCGGCGGAAGATTTCGCTAACCAGGAAGCGATCGCCAAACTGCTGCGTTTTGCTTCTACCCATACCGATTCTTCTGCGCAGACCGTATCTCTGGAAGACTACGTTTCGCGCATGAAAGAAGGGCAGGAGAAAATCTACTACATc                                                                                                                                                                                                                  >  1:102645/1‑251 (MQ=255)                                tGACAAAGCGTGTGATGCAAATGATGGAAAAACTGGAGAAAGACGAAGCGGAAAAATACAAGAACTTCTGGTAAAAGTTTGGACTGGTACTGAAAGAAGGTCCGGAGGAAGATTTCGATAATCAGGAAGTGATAGCAAAAGTGCTGCGTTTTGCTTCTAACTATATCGATTCTTCTGCGCAGACCGTATCTCTGGAAGTTTAGGTTTTGAGCATGAAAGATGGGAAGGAGAAAATctattaaatctacgc                                                                                                                                                                                                             >  2:134015/1‑238 (MQ=255)                                               tgcAAATGCTGGAAAAACTGGCGAAAGACGACGCGGAAAAATACCAGACCTTCTGGCAACAGTTTGGCATGGTACTGAAAGAAGGTCCGGCGGAAGATTTCGCTAACCAGGAAGCGATCGCCAAACTGCTGCGTTTTGCTTCTACCCATACCGATTCTTCTGCGCAGACCGTATCTCTGGAAGACTACGTTTCGCGCATGAAAGAAGGGCAGGAGAAAATCTACTACATCACCGCAGACAGCTATGCGGc                                                                                                                                                                                              >  1:304/1‑250 (MQ=255)                                                        tGGAAAAACTGGCGAAAGACGACGCGGAAAAATACCAGACCTTCTGGCAACAGTTTGGCCTGGTACTGAAAGAAGGTCCGGCGGAAGATTTCGCTAACCAGGAAGCGATCGCCAAACTGCTGCGTTTTGCTTCTACCCATACCGATTCTTCTGCGCAGACCGTATCTCTGGAAGACTACGTTTCGCGCATGAAAGAAGGGCAGGAGAAAATCTACTACATCACCGCAGACAGCTATGCGGCAGCGAAGagc                                                                                                                                                                                    >  1:104509/1‑251 (MQ=255)                                                          gAAAAACTGGCGAAAGACGACGCGGAAAAATACCAGACCTTCTGGCAACAGTTTGGCCTGGTACTGAAAGAAGGTCCGGCGGAAGATTTCGCTAACCAGGAAGCGATCGCCAAACTGCTGCGTTTTGCTTCTACCCATACCGATTCTTCTGCGCAGACCGTATCTCTGGAAGACTACGTTTCGCGCATGAAAGAAGGGCAGGAGAAAATCTACTACATCACCGCAGACAGCTATGCGGCAGCGAAGagcag                                                                                                                                                                                  <  2:102645/251‑1 (MQ=255)                                                                      aaaGACGAGTCGGAAAAATACCAGAGCTTCTGGAAACAGTTTGGCGTGGTACGGAAAGAAGGTCCGGCGTAATATTTCGCTAACCAGTAAGCGATCGCCAAACTGCTGCGTTTTGCTTCTACCCATACCGAGTCTGCTGCGCAGACCGTATCTCTGGAAGACTCCGTTTCGCGCATGAAAGAAGGGCAGGAGAAAAGCTACTACATCACCTCAGACAGCTAGGCGGCAGCGAAGAGCAGCCCGCACCTggt                                                                                                                                                                      <  2:64990/251‑2 (MQ=255)                                                                      aaaGACGACGCGGAAAAATACCAGACCTTCTGGCAACAGTTTGGCCTGGTACTGAAAGAAGGTCCGGCGGAAGATTTCGCTAACCAGGAAGCGATCGCCAAACTGCTGCGTTTTGCTTCTACCCATACCGATTCTTCTGCGCAGACCGTATCTCTGGAAGACTACGTTTCGCGCATGAAAGAAGGGCAGGAGAAAATCTACTACATCACCGCAGACAGCTATGCGGCAGCGAAGAGCAGCCCGCACCTgg                                                                                                                                                                       >  1:61486/1‑250 (MQ=255)                                                                       aaGACGACGCGGAAAAATACCAGCCCTTCTGGCAACAGTTTGGACTGGTACTGAAAGAAGGTACGGCGGACGATTTAGCTAACCAGGAAGCGATCGCCAAACTGCTGCGTTTTGCTTCTGCCCATGCCGATTCTTCTGCGCAGACCGTATCTCTGGAAGACTACGTTTAGCGCATGCAAGAAGGGCAGGTGAAAATCTACTAAATCACCGCCGACAGCACTGCGGCAGCGAAGAGCAGCCCGccacagga                                                                                                                                                                      >  2:169005/1‑243 (MQ=255)                                                                                 ggAAAAATACCAGACCTTCTGGCAACAGTTTGGCCTGGTACTGAAAGAAGGTCCGGCGGAAGATTTCGCTAACCAGGAAGCGATCGCCAAACTGCTGCGTTTTGCTTCTACCCATACCGATTCTTCTGCGCAGACCGTATCTCTGGAAGACTACGTTTCGCGCATGAAAGAAGGGCAGGAGAAAATCTACTACATCACCGCAGACAGCTATGCGGCAGCGAAGAGCAGCCCGCACCTGGAACTGCTGCGt                                                                                                                                                            >  1:84798/1‑250 (MQ=255)                                                                                   aaaaaTACCAGACCTTCTGGCAACAGTTTGGCCTGGTACTGAAAGAAGGTCCGGCGGAAGATTTCGCTAACCAGGAAGCGATCGCCAAACTGCTGCGTTTTGCTTCTACCCATACCGATTCTTCTGCGCAGACCGTATCTCTGGAAGACTACGTTTCGCGCATGAAAGAAGGGCAGGAGAAAATCTACTACATCACCGCAGACAGCTATGCGGCAGCGAAGAGCAGCCCGCACCTGGAACTGCTGCGTAAg                                                                                                                                                         <  1:84949/251‑1 (MQ=255)                                                                                                     ggCAACAGTTTGGCCTGGTACTGAAAGAAGGTCCGGCGGAAGATTTCGCTAACCAGGAAGCGATCGCCAAACTGCTGCGTTTTGCTTCTACCCATACCGATTCTTCTGCGCAGACCGTATCTCTGGAAGACTACGTTTCGCGCATGAAAGAAGGGCAGGAGAAAATCTACTACATCACCGCAGACAGCTATGCGGCAGCGAAGAGCAGCCCGCACCTGGAACTGCTGCGTAAGAAAGGCATCGAAGTTctg                                                                                                                                       >  1:152505/1‑251 (MQ=255)                                                                                                              ttGGCCTGGTACTGAAAGAAGGTCCGGCGGAAGATTTCGCTAACCAGGAAGCGATCGCCAAACTGCTGCGTTTTGCTTCTACCCATACCGATTCTTCTGCGCAGACCGTATCTCTGGAAGACTACGTTTCGCGCATGAAAGAAGGGCAGGAGAAAATCTACTACATCACCGCAGACAGCTATGCGGCAGCGAAGAGCAGCCCGCACCTGGAACTGCTGCGTAAGAAAGGCATCGAAGTTCTGCTGCTTTcc                                                                                                                              >  1:150110/1‑251 (MQ=255)                                                                                                                   cTGGTACTCAAAGAAGGTCCGGTGGATGATTTCGCTAACCAGGAAACGATCGCCAAACTGCTGCGTTTTGCTTCTACCCATACCGATTCTTCTGCGCTGCGCGTATCTCTGCAAGACTACGTTTCGCGCATGCAAGAAGGGCAGGAGAAAATCTACTACATCACCGCAGACAGCTAGACAGCAGCGAAGATCAGCCCTCACATGCAACTGATGCGTCAGAAATGCTTCGACTTTCTGATGCTTTCCGAc                                                                                                                           >  2:251342/1‑249 (MQ=255)                                                                                                                        aCTGAAAGAAGGTTCGGCGGAAGATTTCGCTAACCAGGAAGCGATCGCCAAACTGCTGCGTTTTGCTTCTACCCATACCGATTCTTCTGCGCAGACCGTATCTCTGGAAGACTACGTTTCGCGCATGAAAGAAGGGCAGGAGAAAATCTACTACATCACCGCAGACAGCTATGCGGCAGCGAAGAGCAGCCCGCACCTGGAACTGCTGCGTAAGAAAGGCATCGAAGTTCTGCTGCTTTCCGACCGCATc                                                                                                                     <  1:16448/250‑1 (MQ=255)                                                                                                                                 agGTCCGGCGGAAGATTTCGCTAACCAGGAAGCGATCGCCAAACTGCTGCGTTTTGCTTCTACCCATACCGATTCTTCTGCGCAGACCGTATCTCTGGAAGACTACGTTTCGCGCATGAAAGAAGGGCAGGAGAAAATCTACTACATCACCGCAGACAGCTATGCGGCAGCGAAGAGCAGCCCGCACCTGGAACTGCTGCGTAAGAAAGGCATCGAAGTTCTGCTGCTTTCCGACCGAATCGATGAGTggg                                                                                                           >  2:213691/1‑250 (MQ=255)                                                                                                                                         cggAAGATTTCGCTAACCAGGAAGCGATCGCCAAACTGCTGCGTTTTGCTTCTACCCATACCGATTCTTCTGCGCAGACCGTATCTCTGGAAGACTACGTTTCGCGCATGAAAGAAGGGCAGGAGAAAATCTACTACATCACCGCAGACAGCTATGCGGCAGCGAAGAGCAGCCCGCACCTGGAACTGCTGCGTAAGAAAGGCATCGAAGTTCTGCTGCTTTCCGACCGCATCGATGAGTGGATGATGAAc                                                                                                   <  1:6190/251‑1 (MQ=255)                                                                                                                                                 tgcGCTAACCAGGAAGCGATCGCCAAACTTCTGCGTTTTGCTTCTACCCATACCGATTCTTCTGCGCAGACCGTATCTCTGGAAGACTACGTTTCGCGCATGAAAGAAGGGCAGGAGAAAATCTACTACATCACCGCAGACAGCTATGCGGCAGCGAAGAGCAGCCCGCACCTGGAACTGCTGCGTAAGAAAGGCATCGAAGTTCTGCTGCTTTCCGACCGCATCGATGAGTGGATGATGAACTATctga                                                                                            <  2:152505/248‑1 (MQ=255)                                                                                                                                                  tCGCTAACCAGGAAGCGATCGCCAAACTGCTGCGTTTTGCTTCTACCCATACCGATTCTTCTGCGCAGACCGTATCTCTGGAAGACTACGTTTCGCGCATGAAAGAAGGGCAGGAGAAAATCTACTACATCACCGCAGACAGCTATGCGGCAGCGAAGAGCAGCCCGCACCTGGAACTGCTGCGTAAGAAAGGCATCGAAGTTCTGCTGCTTTCCGACCGCATCGATGAGTGGATGATGAACTATctgact                                                                                          >  1:196705/1‑251 (MQ=255)                                                                                                                                                        aCCAGGAAGCGATCGCCAAACTGCTGCGTTTTGCTTCTACCCATACCGATTCTTCTGCGCAGACCGTATCTCTGGAAGACTACGTTTCGCGCATGAAAGAAGGGCAGGAGAAAATCTACTACATCACCGCAGACAGCTATGCGGCAGCGAAGAGCAGCCCGCACCTGGAACTGCTGCGTAAGAAAGGCATCGAAGTTCTGCTGCTTTCCGACCGCATCGATGAGTGGATGATGAACTATCTGACTGAGTTc                                                                                    <  1:134015/251‑1 (MQ=255)                                                                                                                                                                 cGATCGCCAAACTGCTGCGTTTTGCTTCTACCCATACCGATTCTTCTGCGCAGACCGTATCTCTGGAAGACTACGTTTCCCGCATGAAAGAAGGGCAGGAGAAAATCTACTACATCACCGCAGACAGCTATGCGGCAGCGAAGAGCAGCCCGCACCTGGAACTGCTGCGTAAGAAAGGCATCGAAGTTCTGCTGCTTTCCGACCGCATCGATGAGTGGATGATGAACTATCTGACTGAGTTCGACGGTaaa                                                                           <  1:59332/251‑1 (MQ=255)                                                                                                                                                                          aaCTGCTGCGTTTTGCTTCTACCCATACCGATTCTTCTGCGCAGACCGTATCTCTGGAAGACTACGTTTCGCGCATGAAAGAAGGGCAGGAGAAAATCTACTACATCACCGCAGACAGCTATGCGGCAGCGAAGAGCAGCCCGCACCTGGAACTGCTGCGTAAGAAAGGCATCGAAGTTCTGCTGCTTTCCGACCGCATCGATGAGTGGATGATGAACTATCTGACTGAGTTCGACGGTAAACCGTTCCAg                                                                  <  2:196705/251‑1 (MQ=255)                                                                                                                                                                             tgctgcGTTTTGCTTCTACCCATACCGATTCTTCTGCGCAGACCGTATCTCTGGAAGACTACGTTTCGCGCATGAAAGAAGGGCAGGAGAAAATCTACTACATCACCGCAGACAGCTATGCGGCAGCGAAGAGCAGCCCGCACCTGGAACTGCTGCGTAAGAAAGGCATCGAAGTTCTGCTGCTTTCCGACCGCATCGATGAGTGGATGATGAACTATCTGACTGAGTTCGACGGTAAACCGTTCCAGTCg                                                               >  1:156887/1‑251 (MQ=255)                                                                                                                                                                              gctgcGTTTTGCTTCTACCCATACCGATTCTTCTGCGCAGACCGTATCTCTGGAAGACTACGTTTCCCGCATGAAAGAAGGGCAGGAGAAAATCTACTACATCACCGCAGACAGCTATGCGGCAGCGAAGAGCAGCCCGCACCTGGAACTGCTGCGTAAGAAAGGCATCGAAGTTCTGCTGCTTTCCGACCGCATCGATGAGTGGATGATGAACTATCTGACTGAGTTCGACGGTAAACCGTTCCAGTCgg                                                              >  2:259162/1‑251 (MQ=255)                                                                                                                                                                               ctgcGTTTTGCTTCTACCCATACCGATTCTTCTGCGCAGACCGTATCTCTGGAAGACTACGATTCGCGCATGACAGAAGGGCAGGAGAAAATCTACTACATCACCGCAGACAGCTATGCAGAAGCGAAGAGCAGCCCGCACCTGGAACTGCTGCGTAAGAAAGGCATCGAAAATCTGCTGCTTTCCGACCGCATCGATGAGTGGCTGATGAACTATCTGACTGAGTTCGACGGTGAACCGTTCCAGTCGgt                                                             >  2:242402/1‑251 (MQ=255)                                                                                                                                                                                     tttGCTTCTACCCATACCGATTCTTCTGCGCAGACCGTATCTCTGGAAGACTACGTTTCGCGCATGAAAGAAGGGCAGGAGAAAATCTACTACATCACCGCAGACAGCTATGCGGCAGCGAAGAGCAGCCCGCACCTGGAACTGCTGCGTAAGAAAGGCATCGAAGTTCTGCTGCTTTCCGACCGCATCGATGAGTGGATAATGAACTATCTGACTGCGTTCGACGGTCAACCGTTACAGTCGGTGTCTa                                                        >  2:60046/1‑250 (MQ=255)                                                                                                                                                                                              aCCCATACCGATTCTTCTGCGCAGACCGTATCTCTGGAAGACTACGTTTCGCGCATGAAAGAAGGGCAGGAGAAAATCTACTACATCACCGCAGACAGCTATGCGGCAGCGAAGAGCAGCCCGCACCTGGAACTGCTGCGTAAGAAAGGCATCGAAGTTCTGCTGCTTTCCGACCGCATCGATGAGTGGATGATGAACTATCTGACTGAGTTCGACGGTAAACCGTTCCAGTCGGTGTCTAAAGTTGACg                                               >  1:111793/1‑250 (MQ=255)                                                                                                                                                                                                  aTACCGATTCGTCTGCGCAGACCGTATCTCTGGAAGACTACGTTTCGCGCATGAAAGAAGGGCAGGAGAAAATCTACTACATCACCGCAGACAGCTATGCGGCAGCGAAGAGCAGCCCGCACCTGGAACTGCTGCGTAAGAAAGGCATCGAAGTTCTGCTGCTTTCCGACCGCATCGATGAGTGGATGATGAACTATCTGACTGAGTTCGACGGTAAACCGTTCTAGTCGGTGTCTAAAGTTGACGAGTCg                                          <  1:106975/251‑1 (MQ=255)                                                                                                                                                                                                            ttctGCGCAGACCGTATCTCTGGAAGACTACGTTTCGCGCATGAAAGAAGGGCAGGAGAAAATCTACTACATCACCGCAGACAGCTATGCGGCAGCGAAGAGCAGCCCGCACCTGGAACTGCTGCGTAAGAAAGGCATCGAAGTTCTGCTGCTTTCCGACCGCATCGATGAGTGGATGATGAACTATCTGACTGAGTTCGACGGTAAACCGTTCCAGTCGGTGTCTAAAGTTGACGAGTc                                           <  2:107054/240‑1 (MQ=255)                                                                                                                                                                                                            ttctGCGCAGACCGTATCTCTGGAAGACTACGTTTCGCGCATGAAAGAAGGGCAGGAGAAAATCTACTACATCACCGCAGACAGCTATGCGGCAGCGAAGAGCAGCCCGCACCTGGAACTGCTGCGTAAGAAAGGCATCGAAGTTCTGCTGCTTTCCGACCGCATCGATGAGTGGATGATGAACTATCTGACTGAGTTCGACGGTAAACCGTTCCAGTCGGTGTCTAAAGTTGACGAGTc                                           >  1:107054/1‑240 (MQ=255)                                                                                                                                                                                                                               ctGGAAGACTACGTTTCGCGCATGAAAGAAGGGCAGGAGAAAATCTACTACATCACCGCAGACAGCTATGCGGCAGCGAAGAGCAGCCCGCACCTGGAACTGCTGCGTAAGAAAGGCATCGAAGTTCTGCTGCTTTCCGACCGCATCGATGAGTGGATGATGAACTATCTGACTGAGTTCGACGGTAAACCGTTCCAGTCGGTGTCTAAAGTTGACGAGTCGCTTGAAAAACTGGCTGACGAAGTtga                <  1:243833/248‑1 (MQ=255)                                                                                                                                                                                                                                  gAAGACTACGTTTCGCGCATGAAAGAAGGGCAGGAGAAACTCTACTACATCACCGCAGACAGCTATGCGGCAGCGAAGAGCAGCCCGCACCTGGAACTGCTGCGTAAGAAAGGCATCGAAGTTCTGCTGCTTTCCGACCGCATCGATGAGTGGCTGATGAACTATCTGACTGAGTTCGACGGTAAACCGTTCCAGTCGGTGTCTAAAGTTGACGAGTCGCTTGAAAAACTGGCTGACGAAGTTGTTgag            >  1:175402/1‑249 (MQ=255)                                                                                                                                                                                                                                     gACTACGTTTCGCGCATGAAAGAAGGGCAGGAGAAAATCTACTACATCACCGCAGACAGCTATGCGGCAGCGAAGAGCAGCCCGCACCTGGAACTGCTGCGTAAGAAAGGCATCGAAGTTCTGCTGCTTTCCGACCGCATCGATGAGTGGATGATGAACTATCTGACTGAGTTCGACGGTAAACCGTTCCAGTCGGTGTCTAAAGTTGACGAGTCGCTTGAAAAACTGGCTGTCGAAGTTGATGAGAgcgc       >  1:20056/1‑251 (MQ=255)                                                                                                                                                                                                                                          cGTTTCGCGCATGAAAGAAGGTCAGGAGAAAATCTACTACATAACCGAAGAAAGCTATTCGGCAGCGAAGAGCAGCCCGCGCCTGGAACTGATGCGTAAGAAAGGCATTGAAGGTCTGCTTCTTTCCGACCGCATCGATGAGTGGATGATGAACTAGCTGACGTAGTTCGACGGTAAACCGTTCCAGTCGGGTTCGAAAGTTGATGAGTCGCTTGAAAAACTGGCTTACGAAGTTGATGAGAGCGTGAaag  <  2:78145/251‑1 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                     AGCACGGTAACGCGTAACCTGCGCAATGCGCTGACCAAGCGTGTGCTGCAAATGCTGGAAAAACTGGCGAAAGACGACGCGGAAAAATACCAGACCTTCTGGCAACAGTTTGGCCTGGTACTGAAAGAAGGTCCGGCGGAAGATTTCGCTAACCAGGAAGCGATCGCCAAACTGCTGCGTTTTGCTTCTACCCATACCGATTCTTCTGCGCAGACCGTATCTCTGGAAGACTACGTTTCCCGCATGAAAGAAGGGCAGGAGAAAATCTACTACATCACCGCAGACAGCTATGCGGCAGCGAAGAGCAGCCCGCACCTGGAACTGCTGCGTAAGAAAGGCATCGAAGTTCTGCTGCTTTCCGACCGCATCGATGAGTGGATGATGAACTATCTGACTGAGTTCGACGGTAAACCGTTCCAGTCGGTGTCTAAAGTTGACGAGTCGCTTGAAAAACTGGCTGACGAAGTTGATGAGAGCGCGAAAG  >  NC_000913/496143‑496626

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