Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,082,901	G→C	8.3%	A111A (GCC→GCG)	PP_1861 ←	LysR family transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,082,901	0	G	C	8.3%	187.5 / 7.1	60	A111A (GCC→GCG)	PP_1861	LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (29/26);  new base C (2/3);  total (31/29)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.62e-01

CGACCAGGCGGGTGTCGGCCAAGGGGCCGGCGCGCAGCACCAGGTCGACCTCGCCCAGGTGGCTGCCCTGCAGGTCGACAAAGCTGTCGATCAACCGCAGTTGCACGTCCAGTCCCGGGTAGGCCACCAGGAAGT‑CGGCAATGGCCGGGGCCAGGTGGCGGCGGCCAAACGCGG‑CAGGGGCGTCGATGCGGATCAGGCCTTCCGGGGCATTGCTCAGTGACACCGCTTCAGCCCGGGCCAGGCGCAACTCCTCGATGATGCGCCTGGCCCGCTCGGCAAAGGCGTTGCCTGCAGGC  >  NC_002947/2082753‑2083048                                                                                                                                                      |                                                                                                                                                     cGACCAGGCGGGTGTCGGCCAAGGGGCCGGCGCGCAGCACCAGGTCGACCTCGCCCAGGTGGCTGCCCTGCAGGTCGACAAAGCTGTCGATCAACCGCAGTTGCACGTCCAGTCCCGGGTAGGCCACCAGGAAGT‑CGGCAATGGCCggg                                                                                                                                                      <  5:264493/149‑1 (MQ=255)        gCGGGTGTCGGCCAAGGGGCCGGCGCGCAGCACCAGGTCGACCTCGCCCAGGTGGCTGCCCTGCAGGTCGACAAAGCTGTCGATCAACCGCAGTTGCACGTCCAGTCCCGGGTAGGCCACCAGGAAGT‑CGGCAATGGCCGGGGCCAGGt                                                                                                                                               >  6:324668/1‑149 (MQ=255)             tgtCGGCCAAGGGGCCGGCGCGCAGCACCAGGTCGACCTCGCCCAGGTGGCTGCCCTGCAGGTCGACAAAGCTGTCGATCAACCGCAGTTGCACGTCCAGTCCCGGGTAGGCCACCAGGAAGT‑CGGCAATGGCCGGGGCCAGGTggcgg                                                                                                                                          >  1:60965/1‑149 (MQ=255)                cGGCCAAGGGGCCGGCGCGCAGCACCAGGTCGACCTCGCCCAGGTGGCTGCCCTGCAGGTCGACAAAGCTGTCGATCAACCGCAGTTGCACGTCCAGTCCCGGGTAGGCCACCAGGAAGT‑CGGCAATGGCCGGGGCCAGGTggcggcgg                                                                                                                                       <  6:114181/149‑1 (MQ=255)                 ggCCAAGGGGCCGGCGCGCAGCACCAGGTCGACCTCGCCCAGGTGGCTGCCCTGCAGGTCGACAAAGCTGTCGATCAACCGCAGTTGCACGTCCAGTCCCGGGTAGGCCACCAGGAAGT‑CGGCAATGGCCGGGGCCAGGTggcggcggc                                                                                                                                      <  5:80530/149‑1 (MQ=255)                    cAAGGGGCCGGCGCGCAGCACCAGGTCGACCTCGCCCAGGTGGCTGCCCTGCAGGTCGACAAAGCTGTCGATCAACCGCAGTTGCACGTCCAGTCCCGGGTAGGCCACCAGGAAGT‑CGGCAATGGCCGCCGCCACCTggccccggc                                                                                                                                      >  5:68678/1‑146 (MQ=255)                    cAAGGGGCCGGCGCGCAGCACCAGGTCGACCTCGCCCAGGTGGCTGCCCTGCAGGTCGACAAAGCTGTCGATCAACCGCAGTTGCACGTCCAGTCCCGGGTAGGCCACCAGGAAGT‑CGGCAATGGCCGCCGCCACCTggccccggc                                                                                                                                      <  6:68678/146‑1 (MQ=255)                           gcGGCGCGCAGCACCAGGTCGACCTCGCCCAGGTGGCTGCCCTGCAGGTCGACAAAGCTGTCGATCAACCGCAGTTGCACGTCCAGTCCCGGGTAGGCCACCAGGAAGT‑CGGCAATGGCCGGGGCCAGGTGGCGGCGGCCAAACGCGG‑c                                                                                                                           <  8:325364/148‑1 (MQ=255)                                     ccaCCTGGGCGAGGTCGACCTGGTGGCTGCCCTGCAGGTCGACAAAGCTGTCGATCAACCGCAGTTGCACGTCCAGTCCCGGGTAGGCCACCAGGAAGT‑CGGCAATGGCCGGGGCCAGGTGGCGGCGGCCAAACGCGG‑CAGGGGCGTc                                                                                                                  >  6:294451/2‑148 (MQ=255)                                        ccAGGTCGACCTCGCCCAGGTGGCTGCCCTGCAGGTCGACAAAGCTGTCGATCAACCGCAGTTGCACGTCCAGTCCCGGGTAGGCCACCAGGAAGT‑CGGCAATGGCCGGGGCCAGGTGGCGGCGGCCAAACGCGG‑CAGGGGCGTCGATg                                                                                                              >  8:28539/1‑149 (MQ=255)                                         cAGGTCGACCTCGCCCAGGTGGCTGCCCTGCAGGTCGACAAAGCTGTCGATCAACCGCAGTTGCACGTCCAGTCCCGGGTAGGCCACCAGGAAGT‑CGGCAATGGCCGGGGCCAGGTGGCGGCGGCCAAACGCGG‑CAGGGGCGt                                                                                                                   >  6:104874/1‑143 (MQ=255)                                         cAGGTCGACCTCGCCCAGGTGGCTGCCCTGCAGGTCGACAAAGCTGTCGATCAACCGCAGTTGCACGTCCAGTCCCGGGTAGGCCACCAGGAAGT‑CGGCAATGGCCGGGGCCAGGTGGCGGCGGCCAAACGCGG‑CAGGGGCGt                                                                                                                   <  5:104874/143‑1 (MQ=255)                                          aGGTCGACCTCGCCCAGGTGGCTGCCCTGCAGGTCGACAAAGCTGTCGATCAACCGCAGTTGCACGTCCAGTCCCGGGTAGGCCACCAGGAAGT‑CGGCAATGGCCGGGGCCAGGTGGCGGCGGCCAAACGCGG‑CAGGGGCGTCGATGCg                                                                                                            <  2:60965/149‑1 (MQ=255)                                                aCCTCGCCCAGGTGGCTGCCCTGCAGGTCGACAAAGCTGTCGATCAACCGCAGTTGCACGTCCAGTCCCGGGTAGGCCACCAGGAAGT‑CGGCAATGGCCGGGGCCAGGTGGCGGCGGCCAAAcgc                                                                                                                               <  3:131651/125‑1 (MQ=255)                                                aCCTCGCCCAGGTGGCTGCCCTGCAGGTCGACAAAGCTGTCGATCAACCGCAGTTGCACGTCCAGTCCCGGGTAGGCCACCAGGAAGT‑CGGCAATGGCCGGGGCCAGGTGGCGGCGGCCAAAcgc                                                                                                                               >  4:131651/1‑125 (MQ=255)                                                aCCTCGCCCAGGTGGCTGCCCTGCAGGTCGACAAAGCTGTCGATCAACCGCAGTTGCACGTCCAGTCCCGGGTAGGCCACCAGGAAGT‑CGGCAATGGCCGGGGCCAGGTGGCGGCGGCCAAACGCGG‑CAGGGGCGTCGATGCGGATCAg                                                                                                      <  5:271851/149‑1 (MQ=255)                                                aCCTCGCCCAGGTGGCTGCCCTGCAGGTCGACAAAGCTGTCGATCAACCGCAGTTGCACGTCCAGTCCCGGGTAGGCCACCAGGAAGT‑CGGCAATGGCCGCCGCCACCTGGCCCCGGCCatt                                                                                                                                  >  1:3637/1‑120 (MQ=255)                                                aCCTCGCCCAGGTGGCTGCCCTGCAGGTCGACAAAGCTGTCGATCAACCGCAGTTGCACGTCCAGTCCCGGGTAGGCCACCAGGAAGT‑CGGCAATGGCCGCCGCCACCTGGCCCCGGCCatt                                                                                                                                  <  2:3637/122‑3 (MQ=255)                                                 gctggCCCATGTGGCTGCCCTGCAGGTCGACAAAGCTGTCGATCAACCGCAGTTGCACGTCCAGTCCCGGGTAGGCCACCAGGAAGT‑CGGCAATGGCCGGGGCCAGGTGGCGGCGGCCAAACGCGG‑Ca                                                                                                                          <  4:64051/124‑1 (MQ=255)                                                 ccTCGCCCAGGTGGCTGCCCTGCAGGTCGACAAAGCTGTCGATCAACCGCAGTTGCACTTCCAGTCCCGGGTAGGCCAGCAGGAAGT‑CGGCAATGGCCGGGGCCAGGTGGCGGCGGCCAAACGCGG‑CAGGGGCGTCGATGCGGATCAgg                                                                                                     <  6:25698/149‑1 (MQ=255)                                                 ccTCGCCCAGGTGGCTGCCCTGCAGGTCGACAAAGCTGTCGATCAACCGCAGTTGCACGTCCAGTCCCGGGTAGGCCACCAGGAAGT‑CGGCAATGGCCGGGGCCAGGTGGCGGCGGCCAAACGCGG‑CAGGGGCGTCGATGCGGATCAgg                                                                                                     <  4:320634/149‑1 (MQ=255)                                                              gCTGCCCTGCAGGTCGACAAAGCTGTCGATCAACCGCAGTTGCACGTCCAGTCCCGGGTAGGCCACCAGGAAGT‑CGGCAATGGCCGGGGCCAGGTGGCGGCGGCCAAACGCGG‑CAGGGGCGTCGATGCGGATCAGGCCTTCCGGGGCAt                                                                                        >  5:295758/1‑149 (MQ=255)                                                                  cccTGCAGGTCGACAAAGCTGTCGATCAACCGCAGTTGCACGTCCAGTCCCGGGTAGGCCACCAGGAAGT‑CGGCAATGGCCGGGGCCAGGTGGCGGCGGCCAAACGCGG‑CAGGGGCGTCGATGCGGATCAGGCCTTCCGGGGCATTg                                                                                      <  3:84773/147‑1 (MQ=255)                                                                     tGCAGGTCGACAAAGCTGTCGATCAACCGCAGTTGCACGTCCAGTCCCGGGTAGGTCACCAGGAAGT‑CGGCAATGGCCGGGGCCAGGTGGCGGCGGCCAAACGCGG‑CAGGGGCGTCGATGCGGATCAGGCCTTCCGGGGCATTGCTCAg                                                                                 >  8:266527/1‑149 (MQ=255)                                                                       cAGGTCGACAAAGCTGTCGATCAACCGCAGTTGCACGTCCAGTCCCGGGTAGGTCACCAGGAAGT‑CGGCAATGGCCGGGGCCAGGTGGCGGCGGCCAAACGCGG‑CAGGGGCGTCGATGCGGATCAGGCCTTCCGGGGCATTGCTCAGt                                                                                <  7:266527/148‑1 (MQ=255)                                                                       cAGGTCGACAAAGCTGTCGATCAACCGCAGTTGCACGTCCAGTCCCGGGTAGGCCACCAGGAAGT‑CGGCAATGGCCGGGGCCAGGTGGCGGCGGCCAAACGCGGCCAGGGGCGTCGATGCGGATCAGGCCTTCCGGGGCATTGCTCAGt                                                                                <  1:245398/149‑1 (MQ=255)                                                                       cAGGTCGACAAAGCTGTCGATCAACCGCAGTTGCACGTCCAGTCCCGGGTAGGCCACCAGGAAGT‑CGGCAATGGCCGGGGCCAGGTGGCGGCGGCCAAACGCGG‑CAGGGGCGTCGATGCGGATCAGGCCTTCCGGGGCATTGCTCAGt                                                                                <  8:34367/148‑1 (MQ=255)                                                                       cAGGTCGACAAAGCTGTCGATCAACCGCAGTTGCACGTCCAGTCCCGGGTAGGCCACCAGGAAGT‑CGGCAATGGCCGGGGCCAGGTGGCGGCGGCCAAACGCGG‑CAGGGGCGTCGATGCGGATCAGGCCTTCCGGGGCATTGCTCAGt                                                                                <  5:324668/148‑1 (MQ=255)                                                                          gTCGACAAAGCTGTCGATCAACCGCAGTTGCACGTCCAGTCCCGGGTAGGCCACCAGGAAGT‑CGGCAATGGCCGGGGCCAGGTGGCGGCGGCCAAACGCGG‑CAGGGGCGTCGATGCGGATCAGGCCTTCCGGGGCATTGCTCAGTGaca                                                                            >  6:122315/1‑149 (MQ=255)                                                                          gTCGACAAAGCTGTCGATCAACCGCAGTTGCACGTCCAGTCCCGGGTAGGCCACCAGGAAGT‑CGGCAATGGCCGGGGCCAGGTGGCGGCGGCCAAACGCGG‑CAGGGGCGTCGATGCGGATCAGGCCTTCCGGGGCATTGCTCAGTGaca                                                                            <  7:28539/149‑1 (MQ=255)                                                                                     tGTCGATCAACCGCAGTTGCACGTCCAGTCCCGGGTAGGCCACCAGGAAGTACGGCAATGGCCGGGGCCAGGTGGCGGCGGCCAAACGCGG‑CAGGGGCGTCGATGCGGATCAGGCCTTCCGGGGCATTGCTCAGTGACACCGCTTCAGc                                                                  >  1:259567/1‑149 (MQ=255)                                                                                      gTCGATCAACCGCAGTTGCACGTCCAGTCCCGGGTAGGCCACCAGGAAGT‑CGGCAATGGCCGGGGCCAGGTGGCGGCGGCCAAACGCGG‑CAGGGGCGTCGATGCGGATCAGGCCTTCCGGGGCATTGCTCAGTGACACCGCTTCAGccc                                                                >  8:67527/1‑149 (MQ=255)                                                                                           tCAACCGCAGTTGCACGTCCAGTCCCGGGTAGGCCACCAGGAAGT‑CGGCAATGGCCGGGGCCAGGTGGCGGCGGCCAAACGCGG‑CAGGGGCGTCGATGCGGATCAGGCCTTCCGGGGCATTGCTCAGTGACACCGCTTCAGCCCGGGcc                                                           >  2:82814/1‑149 (MQ=255)                                                                                              aCCGCAGTTGCACGTCCAGTCCCGGGTAGGCCACCAGGAAGT‑CGGCAATGGCCGGGGCCAGGTGGCGGCGGCCAAACGCGG‑CAGGGGCGTCGCTGCGGATCAGGCCTTCCGGGGACTTGCTCAGTGACACCGCTTCAGCCCGGGCCAgg                                                        >  7:61543/1‑149 (MQ=255)                                                                                               ccGCAGTTGCACGTCCAGTCCCTGGTAGGCCACCAGGAAGT‑CGGCAATGGCCGGGGCCAGGTGGCGGCGGCCAAACGCGG‑CAGGGGCGTCGATGCGGATCAGGCCTTCCGGGGCATTGCTCAGTGACACCGc                                                                        <  7:237564/132‑1 (MQ=255)                                                                                               ccGCAGTTGCACGTCCAGTCCCGGGTAGGCCACCAGGAAGT‑CGGCAATGGCCGGGGCCAGGTGGCGGCGGCCAAACGCGG‑CAGGGGCGTCGATGCGGATCAGGCCTTCCGGGGCATTGCTCAGTGACACCGCTTCAGCCCGGGCCAGgc                                                       >  3:23212/1‑149 (MQ=255)                                                                                               ccGCAGTTGCACGTCCAGTCCCGGGTAGGCCACCAGGAAGT‑CGGCAATGGCCGGGGCCAGGTGGCGGCGGCCAAACGCGG‑CAGGGGCGTCGATGCGGATCAGGCCTTCCGGGGCATTGCTCAGTGACACCGCTTCAGCCCGGGCCAGgc                                                       >  4:122253/1‑149 (MQ=255)                                                                                               ccGCAGTTGCACGTCCAGTCCCGGGTAGGCCACCAGGAAGT‑CGGCAATGGCCGGGGCCAGGTGGCGGCGGCCAAACGCGG‑CAGGGGCGTCGAGGCGGATCAGGCCTGCCGGGGCCTTGCTCAGTGACACCGc                                                                        >  8:237564/1‑132 (MQ=255)                                                                                                  cAGTTGCACGTCCAGTCCCGGGTAGGCCACCAGGAAGT‑CGGCAATGGCCGGGGCCAGGTGGCGGCGGCCAAACGCGG‑CAGGGGCGTCGATGCGGATCAGGCCTTCCGGGGCAt                                                                                        <  8:323445/113‑1 (MQ=255)                                                                                                  cAGTTGCACGTCCAGTCCCGGGTAGGCCACCAGGAAGT‑CGGCAATGGCCGGGGCCAGGTGGCGGCGGCCAAACGCGG‑CAGGGGCGTCGATGCGGATCAGGCCTTCCGGGGCAt                                                                                        >  7:323445/1‑113 (MQ=255)                                                                                                  cAGTTGCACGTCCAGTCCCGGGTAGGCCACCAGGAAGT‑CGGCAATGGCCGGGGCCAGGTGGCGGCGGCCAAACGCGG‑CAGGGGCGTCGATGCGGATCAGGCCTTCCGGGGCATTGCTCAGTGACACCGCTTCAGCCCGGGCCAGgcgc                                                     >  5:80458/1‑148 (MQ=255)                                                                                                       gCACGTCCAGTCCCGGGTAGGCCACCAGGAAGT‑CGGCAATGGCCGGGGCCAGGTGGCGGCGGCCAAACGCGG‑CAGGGGCGTCGATGCGGATCAGGCCTTCCGGGGCATTGCTCAGTGACACCGCTTCAGCCCGGGCCAGGCGCatcatc                                               >  2:197809/1‑144 (MQ=255)                                                                                                           gTCCAGTCCCGGGTAGGCCACCAGGAAGT‑CGGCAATGGCCGGGGCCAGGTGGCGGCGGCCAAACGCGG‑CAGGGGCGTCGATGCGGATCAGGCCTTcc                                                                                               <  7:28379/97‑1 (MQ=255)                                                                                                           gTCCAGTCCCGGGTAGGCCACCAGGAAGT‑CGGCAATGGCCGGGGCCAGGTGGCGGCGGCCAAACGCGG‑CAGGGGCGTCGATGCGGATCAGGCCTTcc                                                                                               >  8:28379/1‑97 (MQ=255)                                                                                                                 tCCCGGGTAGGCCACCAGGAAGT‑CGGCAATGGCCGGGGCCAGGTGGCGGCGGCCAAACGCGG‑CAGGGGCGTCGATGCGGATCAGGCCTTCCGGGGCATTGCTCAGTGACACCGCTTCAGCCCGGGCCAGGCGCAACTCCTCGATGATgc                                     >  2:189775/1‑149 (MQ=255)                                                                                                                 tCCCGGGTAGGCCACCAGGAAGT‑CGGCAATGGCCGGGGCCAGGTGGCGGCGGCCAAACGCGG‑CAGGGGCGTCGATGCGGATCAGGCCTTCCGGGGCATTGCTCAGTGACACCGCTTCAGCCCGGGCCAGGCGCAACTCCTCGATGATgc                                     >  4:55883/1‑149 (MQ=255)                                                                                                                   ccGGGTAGGCCACCAGGAAGT‑CGGCAATGGCCGGGGCCAGGTGGCGGCGGCCAAACGCGG‑CAGGGGCGTCGATGCGGATCAGGCCTTCCGGGGCATTGCt                                                                                    >  5:64459/1‑100 (MQ=255)                                                                                                                   ccGGGTAGGCCACCAGGAAGT‑CGGCAATGGCCGGGGCCAGGTGGCGGCGGCCAAACGCGG‑CAGGGGCGTCGATGCGGATCAGGCCTTCCGGGGCATTGCt                                                                                    <  6:64459/100‑1 (MQ=255)                                                                                                                      ggTAGGCCACCAGGAAGT‑CGGCAATGGCCGGGGCCAGGTGGCGGCGGCCAAACGCGG‑CAGGGGCGTCGATGCGGATCAGGCCTTCCGGGGCATTGCTCAGTGACACCGCTTCAGCCCGGGCCAGGCGCAACTCCTCGATGATGCGCCTg                                >  7:130502/1‑149 (MQ=255)                                                                                                                      ggTAGGCCACCAGGAAGT‑CGGCAATGGCCGGGGCCAGGTGGCGGCGGCCAAACGCGG‑CAGGGGCGTCGATGCGGATCAGGCCTTCCGGGGCATTGCTCAGTGACACCGCTTCAGCCCGGGCCAGGCGCAACTCCTCGATGATGCGCCTg                                <  8:197675/149‑1 (MQ=255)                                                                                                                        tAGGCCACCAGGAAGT‑CGGCAATGGCCGGGGCCAGGTGGCGGCGGCCAAACGCGG‑CAGGGGCGTCGATGCGGATCAGGCCTTCCGGGGCATTGCTCAGTGACACCGCTTCAGCCCGGGCCAGGCGCAACTCCTCGATGATGCGCCTGGc                              <  5:122315/149‑1 (MQ=255)                                                                                                                            ccaccaGGAAGT‑CGGCAATGGCCGGGGCCAGGTGGCGGCGGCCAAACGCGG‑CAGGGGCGTCGATGCGGATCAGGCCTTCCGGGGCATTGCTCAGTGACACCGCTTCAGCCCGGGCCAGGCGCAACTCCTCGATGATGCGCCTGGCCCGc                          >  4:196921/1‑149 (MQ=255)                                                                                                                               ccaGGAAGT‑CGGCAATGGCCGGGGCCAGGTGGCGGCGGCCAAACGCGG‑CAGGGGCGTCGATGCGGATCAGGCCTTCCGGGGCATTGCTCAGTGACACCGCTTCAGCCCGGGCCAGGCGCAACTCCTCGATGATGCGCCTGGCCCGCTCg                       >  3:59980/1‑149 (MQ=255)                                                                                                                                     aGT‑CGGCAAGGGCCGGGGCCAGGGGGCGGCGGCAAAACGCGG‑CAGGGGCGTCGATGCGGATCAGGCCTTCCGGGGCATTGCTCAGTCACACCGCTTCAGCCCGGGCCAGGCGCAACTCCTCGATGATGCGCCTGGCCCGCTCGGCAAAg                 <  8:130502/149‑1 (MQ=255)                                                                                                                                      gT‑CGGCAATGGCCGGGGCCAGGTGGCGGCGGCCAAACGCGG‑CAGGGGCGTCGATGCGGATCAGGCCTTCCGGGGCATTGCTCAGTGACACCGCTTCAGCCCGGGCCAGGCGCAACTCCTCGATGATGCGCCTGGCCCGCTCGGCAAAgg                >  3:106444/1‑149 (MQ=255)                                                                                                                                          gcggttGGCCGCCGCCACCTGGCCCCGGCCAAACGCGG‑CAGGGGCGTCGATGCGGATCAGGCCTTCCGGGGCATTGCTCAGTGACACCGCTTCAGCCCGGGCCAGGCGCAACTCCTCGATGATGCGCCTGGCCCGCTCGGCAAAGGc               <  6:263621/142‑1 (MQ=255)                                                                                                                                           gCAATGGCCGGGGCCAGGTGGCGGCGGCCAAACGCGG‑CAGGGGCGTCGATGCGGATCAGGCCTTCCGGGGCATTGCTCAGTGACACCGCTTCAGCCCGGGCCAGGCGCAACTCCTCGATGATGCGCCTGGCCCGCTCGGCAAAGGCGtt            <  7:67527/149‑1 (MQ=255)                                                                                                                                              aTGGCCGGGGCCAGGTGGCGGCGGCCAAACGCGG‑CAGGGGCGTCGATGCGGATCAGGCCTTCCGGGGCATTGCTCAGTGACACCGCTTCAGCCCGGGCCAGGCGCAACTCCTCGATGATGCGCCTGGCCCGCTCGGCAAAGGCGTTGcc         <  4:106444/149‑1 (MQ=255)                                                                                                                                              aTGGCCGGGGCCAGGTGGCGGCGGCCAAACGCGG‑CAGGGGCGTCGATGCGGATCAGGCCTTCCGGGGCATTGCTCAGTGACACCGCTTCAGCCCGGGCCAGGCGCAACTCCTCGATGATGCGCCTGGCCCGCTCGGCAAAGGCGNTGcc         <  5:294451/149‑1 (MQ=255)                                                                                                                                                ggCCGGGGCCAGGTGGCGGCGGCCAAACGCGG‑CAGGGGCGTCGATGCGGATCAGGCCTTCCGGGGCATTGCTCAGTGACACCGCTTCAGCCCGGGCCAGGCGCAACTCCTCGATGATGCGCCTGGCCCGCTCGGCAAAGGCGTTGCCTg       >  5:235561/1‑149 (MQ=255)                                                                                                                                                ggCCGGGGCCAGGTGGAGGCGGCCAAACGCGG‑CAGGGGCGTCGATGCGGATCAGGCCTTCCGGGGCATTGCTCAGTGACACCGCTTCAGCCCGGGCCAGGCGCCACTCCTCGATGATGCGCCTGGCCCGCTCGGCAAAGGCGTTGCCTg       >  6:67591/1‑149 (MQ=255)                                                                                                                                                     gggCCAGGTGGCGGCGGCCAAACGCGG‑CAGGGGCGTCGATGCGGATCAGGCCTTCCGGGGCATTGCTCAGTGACACCGCTTCAGCCCGGGCCAGGCGCAACTCCTCGATGATGCGCCTGGCCCGCTCGGCAAAGGCGTTGCCTGCAGGc  <  1:82814/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                     CGACCAGGCGGGTGTCGGCCAAGGGGCCGGCGCGCAGCACCAGGTCGACCTCGCCCAGGTGGCTGCCCTGCAGGTCGACAAAGCTGTCGATCAACCGCAGTTGCACGTCCAGTCCCGGGTAGGCCACCAGGAAGT‑CGGCAATGGCCGGGGCCAGGTGGCGGCGGCCAAACGCGG‑CAGGGGCGTCGATGCGGATCAGGCCTTCCGGGGCATTGCTCAGTGACACCGCTTCAGCCCGGGCCAGGCGCAACTCCTCGATGATGCGCCTGGCCCGCTCGGCAAAGGCGTTGCCTGCAGGC  >  NC_002947/2082753‑2083048

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