Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,650,676	G→T	8.2%	H173N (CAT→AAT) ‡	PP_4959 ←	diguanylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,650,676	0	G	T	8.2%	183.2 / 6.3	61	H173N (CAT→AAT) ‡	PP_4959	diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base G (31/25);  new base T (3/2);  total (34/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.27e-01

CCCTGGTAGGCCTTGTAGAAGTCCGACTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGG‑AAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTGGACCTTGGGCTTGAGGA  >  NC_002947/5650528‑5650822                                                                                                                                                      |                                                                                                                                                   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)                                                                                                                                                      |                                                                                                                                                   CCCTGGTAGGCCTTGTAGAAGTCCGACTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGG‑AAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTGGACCTTGGGCTTGAGGA  >  NC_002947/5650528‑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: ‑