Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,567,081	Δ1 bp	7.7%	coding (847/879 nt)	PP_2250 ←	LysR family transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,567,081	0	G	.	7.7%	204.1 / 10.4	65	coding (847/879 nt)	PP_2250	LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (29/31);  new base . (2/3);  total (31/34)
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.38e-01

GCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCG‑TGCA‑GCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCGCTGGCCAGCAGCGGGGCGACCCGCGACAACGGCAGATCTCCCCAGCC  >  NC_002947/2566937‑2567228                                                                                                                                                 |                                                                                                                                                      gcgcAGCGAGCGCATGACGCGGTGGGGGGCCTATCGACCGAGATATCGGCTTTAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCG‑TGCa                                                                                                                                                  <  4:282155/149‑1 (MQ=255)        gcgcTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCG‑TGCA‑GCctgct                                                                                                                                          >  8:275909/1‑149 (MQ=255)           cTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGC‑‑TGCACGCctgctgat                                                                                                                                       <  2:133673/149‑3 (MQ=255)            tCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCG‑TGCA‑GCCTGCTGCAg                                                                                                                                      >  3:172233/1‑149 (MQ=255)              aTGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGTAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCG‑TGCA‑GCCTGCTGCAGc                                                                                                                                     <  5:178551/148‑1 (MQ=255)               tGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCG‑TGCA‑GCCTGCTGCAGCgg                                                                                                                                   <  4:217884/149‑1 (MQ=255)                            ggCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCG‑TGCA‑GCCTGCTGCAGCGGCTGTTGTTTgcgc                                                                                                                      >  8:262914/1‑149 (MQ=255)                               cTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCG‑TGCA‑GCCTGCTGCAGCGGCTGTTGTTTGCGCca                                                                                                                    >  5:229571/1‑148 (MQ=255)                                  tCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCG‑TGCA‑GCCTGCTGCAGCGGCTGTTGTTTGCGCcaca                                                                                                                  <  7:325180/147‑1 (MQ=255)                                    gACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCG‑TGCA‑GCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGat                                                                                                              <  1:184163/149‑1 (MQ=255)                                     aCCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCG‑TGCA‑GCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGata                                                                                                             >  8:30556/1‑149 (MQ=255)                                     aCCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCG‑TGCA‑GCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGata                                                                                                             >  6:135882/1‑149 (MQ=255)                                         agaTATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTAACGGATCAGCCAGCTGCTCGATCAGCAGGCG‑TGCA‑GCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCa                                                                                                         >  5:366001/1‑149 (MQ=255)                                             atCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCG‑TGCA‑GCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACggg                                                                                                     >  1:63790/1‑149 (MQ=255)                                             atCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCG‑TGCA‑GCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACggg                                                                                                     >  1:23180/1‑149 (MQ=255)                                              tCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCG‑TGCA‑GCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACggg                                                                                                     >  4:307837/1‑148 (MQ=255)                                                ggCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCG‑TGCA‑GCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAg                                                                                                  >  5:262006/1‑149 (MQ=255)                                                ggCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCG‑TGCA‑GCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAg                                                                                                  <  3:307837/149‑1 (MQ=255)                                                     gAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGCCGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCG‑TGCA‑GCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGc                                                                                              <  7:275909/148‑1 (MQ=255)                                                        gggTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCG‑TGCA‑GCCTGCTGCAGCGGCTGTTGttt                                                                                                                          >  2:45666/1‑117 (MQ=255)                                                        gggTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCG‑TGCA‑GCCTGCTGCAGCGGCTGTTGttt                                                                                                                          <  1:45666/117‑1 (MQ=255)                                                        gggTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCG‑TGCA‑GCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAAtt                                                                                          <  5:135882/149‑1 (MQ=255)                                                          gTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCG‑TGCA‑GCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCg                                                                                        >  3:197044/1‑149 (MQ=255)                                                             aaaGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCG‑TGCA‑GCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGtt                                                                                      <  4:172233/148‑1 (MQ=255)                                                             aaaGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCG‑TGCA‑GCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGtt                                                                                      <  4:68472/148‑1 (MQ=255)                                                             aaaGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCG‑TGCA‑GCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTg                                                                                     >  3:342772/1‑149 (MQ=255)                                                             aaaGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCG‑TGCA‑GCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTg                                                                                     >  7:30075/1‑149 (MQ=255)                                                              aaGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCG‑TGCA‑GCCTTCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGTGCAGTTGCAATTCGTTGc                                                                                    <  1:341501/149‑1 (MQ=255)                                                                gATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCG‑TGCA‑GCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCgt                                                                                  <  7:30556/149‑1 (MQ=255)                                                                  ttGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCG‑TGCA‑GCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCgtg                                                                                 <  7:262914/148‑1 (MQ=255)                                                                      aGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCG‑TGCA‑GCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTa                                                                            <  4:342772/149‑1 (MQ=255)                                                                             aGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCG‑TGCA‑GCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAAtt                                                                     >  5:124387/1‑149 (MQ=255)                                                                             aGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCG‑TGCA‑GCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAAtt                                                                     >  5:180096/1‑149 (MQ=255)                                                                               ggAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCG‑TGCA‑GCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTcc                                                                   <  5:62816/149‑1 (MQ=255)                                                                                gAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCG‑TGCA‑GCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTATGGAATTc                                                                    >  4:130659/1‑147 (MQ=255)                                                                                gAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCG‑TGCA‑GCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTc                                                                    >  8:185572/1‑147 (MQ=255)                                                                                  aaaTTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCG‑TGCA‑GCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGc                                                                                              >  6:58463/1‑119 (MQ=255)                                                                                  aaaTTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCG‑TGCA‑GCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGc                                                                                              <  5:58463/119‑1 (MQ=255)                                                                                  aaaTTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCG‑TGCA‑GCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAg                                                                 <  6:137335/148‑1 (MQ=255)                                                                                          aCAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGC‑‑TGCACGCctgctgat                                                                                                                                       <  4:191973/70‑3 (MQ=37)                                                                                          aCAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGC‑‑TGCACGCctgctgat                                                                                                                                       >  3:191973/1‑68 (MQ=37)                                                                                            aGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCG‑TGCA‑GCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGTACGGGCAGTTGCAATTCGTTGAGTGTATTACGGAATTCAAGGCGCACCAGAc                                                      >  6:73526/1‑149 (MQ=255)                                                                                             gACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCG‑TGCA‑GCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACg                                                     <  6:229571/149‑1 (MQ=255)                                                                                                     gacggacgGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCG‑TGCA‑GCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGtt                                                                                      >  6:149735/1‑108 (MQ=255)                                                                                                     gacggacgGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCG‑TGCA‑GCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGtt                                                                                      <  5:149735/108‑1 (MQ=255)                                                                                                      acggacggCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCG‑TGCA‑GCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCGCTGGc                                            <  7:185572/149‑1 (MQ=255)                                                                                                            ggCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCG‑TGCA‑GCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGAATTACGGAATTCCAGGCGCACCAGACGCCCGCTGGCcagcag                                      >  5:372723/1‑149 (MQ=255)                                                                                                              cAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTTGCA‑GCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCGCTGGCCagcagc                                     >  1:246362/1‑149 (MQ=255)                                                                                                                 tCACGGATCAGCCAGCTGCTCGATCAGCAGGCG‑TGCA‑GCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCGCTGGc                                            >  1:206418/1‑138 (MQ=255)                                                                                                                 tCACGGATCAGCCAGCTGCTCGATCAGCAGGCG‑TGCA‑GCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCGCTGGc                                            <  2:206418/138‑1 (MQ=255)                                                                                                                 tCACGGATCAGCCAGCTGCTCGATCAGCAGGCG‑TGCA‑GCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCGCTGGCCAGCAGCgggg                                 <  6:180096/149‑1 (MQ=255)                                                                                                                                 tgctCGATCAGCAGGCG‑TGCA‑GCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGCCGCCCGCTGGCCAGCAGCGGGGCGACCCGCGACAACgg                 >  8:206469/1‑149 (MQ=255)                                                                                                                                   ctCGATCAGCAGGCG‑TGCA‑GCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCGCTGGCCAGCAGCGGGGCGACCCGCGACAACgg                 <  2:246362/147‑1 (MQ=255)                                                                                                                                   ctCGATCAGCAGGCG‑TGCA‑GCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCGCTGGCCAGCAGCGGGGCGACCCGCGACAACGGCa               <  6:124387/149‑1 (MQ=255)                                                                                                                                   ctCGATCAGCAGGCG‑TGCA‑GCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCGCTGGCCAGCAGCGGGGCGACCCGCGACAACGGCa               <  6:170184/149‑1 (MQ=255)                                                                                                                                   ctCGATCAGCAGGCG‑TGCA‑GCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCGCTGGCCAGCAGCGGGGCGACCCGCGACAACGGCa               <  6:366001/149‑1 (MQ=255)                                                                                                                                      gATCAGCAGGCG‑TGCA‑GCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCGCTg                                              >  1:96276/1‑115 (MQ=255)                                                                                                                                      gATCAGCAGGCG‑TGCA‑GCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCGCTg                                              <  2:96276/115‑1 (MQ=255)                                                                                                                                      gATCAGCAGGCG‑TGCA‑GCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCGCTGGCCAGCAGCGGGGCGACCCGCGACAACGGCAGAt            <  2:23180/149‑1 (MQ=255)                                                                                                                                       tgcagcagGC‑‑TGCACGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCGCTGGCCAGCAGCGGGGCGACCCGCGACAACGGCAGAt            >  3:254087/3‑148 (MQ=255)                                                                                                                                         cagcagGC‑‑TGCACGCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCCGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCGCTGGCCAGCAGCGGGGCGACCCGCGACAACGGCAGAtctc         <  4:254087/149‑1 (MQ=255)                                                                                                                                            cagGCG‑TGCA‑GCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCGCTGGCCAGCAGCGGGGCGACCCGCGACAACGGCAGATCTccc       >  2:20044/1‑148 (MQ=255)                                                                                                                                             agGCG‑TGCA‑GCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCGCTGGCCAGCAGCGGGGCGACCCGCGACAACGGCAGAt            >  1:3992/1‑142 (MQ=255)                                                                                                                                             agGCG‑TGCA‑GCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCGCTGGCCAGCAGCGGGGCGACCCGCGACAACGGCAGAt            <  2:3992/142‑1 (MQ=255)                                                                                                                                             agGCG‑TGCA‑GCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCGCTGGCCAGCAGCGGGGCGACCCGCGACAACGGCAGATCTCCCCa     <  8:30075/149‑1 (MQ=255)                                                                                                                                                 g‑TGCA‑GCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCGCTGGCCAGCAGCGGGGCGACCCGCGACAACGGCAGATCTCCCCAGcc  >  1:135339/1‑148 (MQ=255)                                                                                                                                                 |                                                                                                                                                      GCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCG‑TGCA‑GCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCGCTGGCCAGCAGCGGGGCGACCCGCGACAACGGCAGATCTCCCCAGCC  >  NC_002947/2566937‑2567228

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