Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,650,674	A→C	7.8%	H173Q (CAT→CAG) ‡	PP_4959 ←	diguanylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,650,674	0	A	C	7.8%	190.3 / 6.9	64	H173Q (CAT→CAG) ‡	PP_4959	diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base A (31/28);  new base C (3/2);  total (34/30)
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.03e-01

CCCCTGGTAGGCCTTGTAGAAGTCCGACTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGG‑AAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTGGACCTTGGGCTTGAGGA  >  NC_002947/5650527‑5650822                                                                                                                                                     |                                                                                                                                                     ccccTGGTAGGCCTTGTAGAAGTCCGACTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGG‑AAATGGCCGGATTGGCAAAACGGATGATGCCAt                                                                                                                                                     <  1:293246/149‑1 (MQ=255)  cccTGGTAGGCCTTGTAGAAGTCCGACTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGG‑AAATGGCCGGATTGGCAAAACGGATGATGCCATg                                                                                                                                                    >  4:240887/1‑149 (MQ=255)    cTGGTATGCCTTGTATAAGTCCGTCTCGCGCCACATGAGGGTACTGGTCATTTGCACAACGTCGCTCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAAATGGCCGGATTGGCAAAACGGATGATGCCATgg                                                                                                                                                   <  1:48903/149‑1 (MQ=255)    cTGGTAGGGGTTGTATCAGTCCGACTGGCGGCACAGGCTGGTACTGGTCAGTTGCACAACGTGGTGCAATTGGGTGCCCTGCAGGTGTTGGACCGGTGCCGCCAGCAAACGGGTTAATGGCCGGATTGGCAAAACGGATGATGCCATgg                                                                                                                                                   <  1:44249/149‑1 (MQ=255)      ggTAGGCCTTGTAGAAGTCCGACTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGG‑AAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCa                                                                                                                                                >  8:74187/1‑149 (MQ=255)           gCCTTGTAGAAGTCCGACTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGG‑AAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAac                                                                                                                                           <  6:295892/149‑1 (MQ=255)              ttGTAGAAGTCCGACTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGG‑AAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAaccac                                                                                                                                        >  3:12873/1‑149 (MQ=255)                ttAGAAGTCCTACTAGCGCAACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGG‑AAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAaccac                                                                                                                                        <  6:227396/146‑1 (MQ=255)                gTAGAAGTCCGACTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGG‑AAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAaccac                                                                                                                                        <  7:74187/147‑1 (MQ=255)                gTAGAAGTCCGACTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGG‑AAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAAccacca                                                                                                                                      >  7:2880/1‑149 (MQ=255)                    aaGTCCGACTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGG‑AAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAAt                                                                                                                                  <  8:35261/149‑1 (MQ=255)                    aaGTCCGACTAGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCAAGCAAACGGG‑AAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAAt                                                                                                                                  >  4:83684/1‑149 (MQ=255)                      gTCCTACTCGCGACACAGGCTGGTACTTGTCAGTTGAACAACGTCGAGCAATTGGGTGACCTGCAGGTGTTGAACAGGTGCCGCCAGCAAACGGG‑AAATGGCCGGATTGGCAAAACGTATGATGCCATGGGCATCAACCACCACAATGc                                                                                                                                >  4:233291/1‑149 (MQ=255)                      gTCCGACTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGG‑AAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGc                                                                                                                                <  3:245084/149‑1 (MQ=255)                         cGACTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGG‑AAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCtt                                                                                                                             >  3:189871/1‑149 (MQ=255)                            cTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGG‑AAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTcgg                                                                                                                          >  2:53812/1‑149 (MQ=255)                                   cacaGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGG‑AAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCAtt                                                                                                                   >  4:66809/1‑149 (MQ=255)                                        gCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGG‑AAATGGCCGGTTTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCa                                                                                                              <  1:93248/149‑1 (MQ=255)                                        gCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGG‑AAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTc                                                                                                                >  2:226858/1‑147 (MQ=255)                                         cTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGG‑AAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAg                                                                                                             <  3:83684/149‑1 (MQ=255)                                           ggTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGG‑AAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGa                                                                                                           <  3:240887/149‑1 (MQ=255)                                               cttgTCAGTTGCACAACGTCTAGCAATTGGGTTCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGG‑AAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATgg                                                                                                        <  8:2880/145‑1 (MQ=255)                                                  gTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGG‑AAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGc                                                                                                    <  3:66809/149‑1 (MQ=255)                                                  gTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAAAGGG‑AAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGgacggacg                                                                                                     >  4:285393/1‑148 (MQ=255)                                                      gTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGG‑AAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCAtt                                                                                                 >  1:55434/1‑148 (MQ=255)                                                           acaACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGG‑AAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGc                                                                                           >  6:108648/1‑149 (MQ=255)                                                            caACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGG‑AAATGGCCGGATTGGCAAAACGGATGATGCCCATGGCATCAtccgtt                                                                                                                                       <  7:151693/104‑4 (MQ=255)                                                            caACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGG‑AAATGGCCGGATTGGCAAAACGGATGATGCCCATGGCATCAtccgtt                                                                                                                                       >  8:151693/1‑101 (MQ=255)                                                                   gAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGG‑AAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGcc                                                                                    <  3:156764/148‑1 (MQ=255)                                                                     gCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGG‑AAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGtc                                                                                 >  1:181043/1‑149 (MQ=255)                                                                      cAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGG‑AAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGtct                                                                                >  1:273764/1‑149 (MQ=255)                                                                          tGGGTGCCCTGCAGGTTTTGGACCGGTGCCGCCAGCAAACGTG‑AAATTGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCa                                                                             <  4:189871/148‑1 (MQ=255)                                                                           gggTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGG‑AAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAg                                                                            <  6:139390/148‑1 (MQ=255)                                                                           gggTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGG‑AAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAg                                                                            <  2:305468/148‑1 (MQ=255)                                                                                cccTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGG‑AAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCAc                                                                       <  1:53812/148‑1 (MQ=255)                                                                                         tgtTGCACCGGTGCCGCCAGCAAACGGG‑AAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAgctg                                                             >  3:289637/1‑149 (MQ=255)                                                                                             gCACCGGTGCCGCCAGCAAACGGG‑AAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGctgctg                                                          >  1:202873/1‑148 (MQ=255)                                                                                              cACCGGTGCCGCCAGCAAACGGG‑AAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGctgctg                                                          >  4:294744/1‑147 (MQ=255)                                                                                              cACCGGTGCCGCCAGCAAACGGG‑AAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGctgctg                                                          >  4:290018/1‑147 (MQ=255)                                                                                                 cGGTGCCGCCAGCAAACGGG‑AAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTgcagca                                                     >  2:268046/1‑149 (MQ=255)                                                                                                  ggTGCCGCCAGCAAACGGG‑AAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCACACGACTCAGCTGCTGCAGCAt                                                    >  1:242881/1‑149 (MQ=255)                                                                                                            gCAAACGGG‑AAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTcc                                                                              >  5:288260/1‑113 (MQ=255)                                                                                                            gCAAACGGG‑AAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTcc                                                                              <  6:288260/113‑1 (MQ=255)                                                                                                            gCAAACGGG‑AAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCATCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTAAGATGCTGCAGCATGCGCCGGATg                                          >  4:117249/1‑149 (MQ=255)                                                                                                                  ggg‑AAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCgctgc                                    <  3:294744/149‑1 (MQ=255)                                                                                                                  ggg‑AAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCgctgc                                    <  3:290018/149‑1 (MQ=255)                                                                                                                  ggg‑AAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCgct                                      >  4:170843/1‑147 (MQ=255)                                                                                                                  ggg‑AAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCgct                                      <  3:170843/147‑1 (MQ=255)                                                                                                                    g‑AAATGGCCGGATTGTCAAACCGGATGATGCCATGGGCATCAACCACCAGAATGCCTTAGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGctgctg                                  <  8:248404/149‑1 (MQ=255)                                                                                                                    g‑AAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGctgctg                                  <  5:108648/149‑1 (MQ=255)                                                                                                                           ggcGGATTGGCAAAACGTATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGTTGGACGCATTGAAGGCCCTGTCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGag                            <  3:233291/147‑1 (MQ=255)                                                                                                                                    gCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTCGACCAGCAgcgcaac                   <  2:93014/149‑4 (MQ=255)                                                                                                                                      aaaaCGGATTATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCgg                                             <  8:155712/121‑1 (MQ=255)                                                                                                                                      aaaaCGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTgca                                                        <  1:295556/110‑1 (MQ=255)                                                                                                                                      aaaaCGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTgca                                                        >  2:295556/1‑110 (MQ=255)                                                                                                                                      aaaaCGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCgg                                             >  7:155712/1‑121 (MQ=255)                                                                                                                                      aaaaCGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTGGa                 >  6:319636/1‑149 (MQ=255)                                                                                                                                       aaaCGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTGGAc                >  3:109030/1‑149 (MQ=255)                                                                                                                                         gtggttgatgCCCATGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCATGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTGGACCt              >  6:111773/6‑149 (MQ=255)                                                                                                                                         aCGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTGGACCt              >  4:244894/1‑149 (MQ=255)                                                                                                                                          cGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTGGACCtt             >  8:9058/1‑149 (MQ=255)                                                                                                                                             ttgatgCCCATGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTgg                  >  8:285419/2‑141 (MQ=255)                                                                                                                                             ttgatgCCCATGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTCCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTgg                  <  7:285419/140‑1 (MQ=255)                                                                                                                                               gatgCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTGGACCTTGGGCt        >  4:77771/1‑149 (MQ=255)                                                                                                                                                     aTGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGATCACGACTAAGCTGCTGCAGCATGCGCCGGTTGCGCTCCTGGTCGAGCAGCGCCTGGACCTTGGGCTTGAgca  >  7:119990/1‑147 (MQ=255)                                                                                                                                                     |                                                                                                                                                     CCCCTGGTAGGCCTTGTAGAAGTCCGACTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGG‑AAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTGGACCTTGGGCTTGAGGA  >  NC_002947/5650527‑5650822

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