Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,567,085:1	+C	7.6%	coding (843/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,085	1	.	C	7.6%	255.4 / 14.8	66	coding (843/879 nt)	PP_2250	LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base . (30/31);  new base C (2/3);  total (32/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 = 1.00e+00

GCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCG‑TGCA‑GCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCGCTGGCCAGCAGCGGGGCGACCCGCGACAACGGCAGATCTCCCCAGCCCAA  >  NC_002947/2566937‑2567231                                                                                                                                                       |                                                                                                                                                   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:68472/148‑1 (MQ=255)                                                             aaaGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCG‑TGCA‑GCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGtt                                                                                         <  4:172233/148‑1 (MQ=255)                                                             aaaGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCG‑TGCA‑GCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTg                                                                                        >  7:30075/1‑149 (MQ=255)                                                             aaaGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCG‑TGCA‑GCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTg                                                                                        >  3:342772/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:180096/1‑149 (MQ=255)                                                                             aGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCG‑TGCA‑GCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAAtt                                                                        >  5:124387/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                                                                                                 <  5:58463/119‑1 (MQ=255)                                                                                  aaaTTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCG‑TGCA‑GCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGc                                                                                                 >  6:58463/1‑119 (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                                                                                         <  5:149735/108‑1 (MQ=255)                                                                                                     gacggacgGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCG‑TGCA‑GCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGtt                                                                                         >  6:149735/1‑108 (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                                               <  2:206418/138‑1 (MQ=255)                                                                                                                 tCACGGATCAGCCAGCTGCTCGATCAGCAGGCG‑TGCA‑GCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCGCTGGc                                               >  1:206418/1‑138 (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:366001/149‑1 (MQ=255)                                                                                                                                   ctCGATCAGCAGGCG‑TGCA‑GCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCGCTGGCCAGCAGCGGGGCGACCCGCGACAACGGCa                  <  6:170184/149‑1 (MQ=255)                                                                                                                                   ctCGATCAGCAGGCG‑TGCA‑GCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCGCTGGCCAGCAGCGGGGCGACCCGCGACAACGGCa                  <  6:124387/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)                                                                                                                                                    gCA‑GCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCGCTGGCCAGCAGCGGGGCGACCCGCGACAACGGCAGATCTCCCCAGCCCaa  <  4:197044/149‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                   GCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCG‑TGCA‑GCCTGCTGCAGCGGCTGTTGTTTGCGCCACAGGATATGCACGGGCAGTTGCAATTCGTTGCGTGTATTACGGAATTCCAGGCGCACCAGACGCCCGCTGGCCAGCAGCGGGGCGACCCGCGACAACGGCAGATCTCCCCAGCCCAA  >  NC_002947/2566937‑2567231

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