Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,155,468	A→C	8.0%	Q209P (CAG→CCG) ‡	PP_1013 →	sensor histidine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,155,468	0	A	C	8.0%	230.3 / 7.4	74	Q209P (CAG→CCG) ‡	PP_1013	sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base A (35/32);  new base C (3/3);  total (38/36)
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.98e-01

CCATTATTCGCTGACCTTGGAGCCGCTGAAACCCCCGGTACTGGTCACGCAGATCCGCTTGGGCGAGGGCGAATGGTTGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTGCTGGTGCACTGGCAAAGCTGGCCACTGAAGCGCCTGGCGCGGGCCGCCCGAGAGATGTCGCTGGGTGCCGATGTGGCGCCAGTGGCTGAAGCAGGC  >  NC_002947/1155325‑1155616                                                                                                                                                |                                                                                                                                                     ccATTATTCGCTGACCTTGGAGCCGCTGAAACCCCCGGTACTGGTCACGCAGATCCGCTTGGGCGAGGGCGAATGGTTGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCg                                                                                                                                                 <  3:361390/149‑1 (MQ=255)    ttattCGCTGACCTTGGAGCCGCTGAAACCCCCGGTACTGGTCACGCAGATCCGCTTGGGCGAGGGCGAATGGTTGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGATCTGCTGGCGCGGCAGGCGTTCGGCt                                                                                                                                              >  7:297711/1‑149 (MQ=255)      attCGCTGACCTTGGAGCCGCTGAAACCCCCGGTACTGGTCACGCAGATCCGCTTGGGCGAGGGCGAATGGTTGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGATCTGCTGGCGCGGCAGGCGTTCGGCtt                                                                                                                                             <  8:297711/148‑1 (MQ=38)       ttCGCTGACCTTGGAGCCGCTGAAACCCCCGGTACTGGTCACGCAGATCCGCTTGGGCGAGGGCGAATGGTTGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGATCTGCTGGCGCGGCAGGCGTTCGGCTTc                                                                                                                                            <  1:36521/148‑1 (MQ=38)        tCGCTGACGATGGATCCGCTGATACCCCCGGAACTGGTCACTCTGATCCGCTTTGGCGAGGTCGAATGGTTGTACCTCGCCTCTCTGTTGCCCGAACCTTTCACGGGGCTGGTAGCCGAACGCCTGCCGCGCCAGCTGATCGGCttga                                                                                                                                           <  6:373737/148‑3 (MQ=255)          gCTGACCTTGGAGCCGCTGAAACCCCCGGTACTGGTCACGCAGATCCGCTTGGGCGAGGGCGAATGGTTGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCg                                                                                                                                        >  6:40697/1‑149 (MQ=255)          gCTGACCTTGGAGCCGCTGAAACCCCCGGTACTGGTCACGCAGATCCGCTTGGGCGAGGGCGAATGGTTGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCg                                                                                                                                        >  4:372858/1‑149 (MQ=255)              aCCTTGGAGCCGCTGAAACCCCCGGTACTGGTCACGCAGATCCGCTTGGGCGAGGGCGAATGGTTGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCt                                                                                                                                    >  6:314708/1‑149 (MQ=255)                     aGCCGCTGAAACCCCCGGTACTGGTCACGCAGATCCGCTTGGGCGAGGGCGAATGGTGGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACg                                                                                                                             <  1:261361/149‑1 (MQ=255)                       ccGCTGAAACCCCCGGTACTGGTCACGCAGATCCGCTTGGGCGAGGGCGAATGGTTGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGc                                                                                                                           >  3:45979/1‑149 (MQ=255)                         gCTGAAACCCCCGGTACTGGTCACGCAGATCCGCTTGGGCGAGGGCGAATGGTTGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCt                                                                                                                         >  8:56662/1‑149 (MQ=255)                         gCTGAAACCCCCGGTACTGGTCACGCAGATCCGCTTGGGCGAGGGCGAATGGTTGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCt                                                                                                                         >  8:230817/1‑149 (MQ=255)                           tGAAACCCCCGGTACTGGTCACGCAGATCCGCTTGGGCGAGGGCGAATGGTTGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGAt                                                                                                                                                   >  1:355335/1‑121 (MQ=255)                           tGAAACCCCCGGTACTGGTCACGCAGATCCGCTTGGGCGAGGGCGAATGGTTGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGAt                                                                                                                                                   <  2:355335/121‑1 (MQ=255)                                cccccGGTACTGGTCACGCAGATCCGCTTGGGCGAGGGCGAATGGTTGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCttgctgtt                                                                                                                  >  3:224234/1‑149 (MQ=255)                                 ccccGGTACTGGTCACGCAGATCCGCTTGGGCGAGGGCGAATGGTTGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCttgctgttg                                                                                                                 <  8:83752/149‑1 (MQ=255)                                 ccccGGTACTGGGCACGCAGATCCGCTTGGGCGAGGGCGAATGGTTGTACATCGCCTCGCGGTTGCCCGAACCGTACACGGGGCGGGAAGCCGAACGCCTGCCGCGCCAGCAGAGCGGCTTCATCGTGCTGACCACGGCCttgctgttg                                                                                                                 <  3:42831/149‑1 (MQ=255)                                  cccGGTACTGGTCACGCAGATCCGCTTGGGCGAGGGCGAATGGTTGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCttgctgttgc                                                                                                                >  8:231800/1‑149 (MQ=255)                                   ccGGTACTGGTCACGCAGATCCGCTTGGGCGAGGGCGAATGGTTGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCttgctgttgct                                                                                                               <  5:303397/149‑1 (MQ=255)                                    cGGTACTGGTCACGCAGATCCGCTTGGGCGAGGGCGAATGGTTGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCttgctgttgctg                                                                                                              <  7:230817/149‑1 (MQ=255)                                          tGGTCACGCAGATCCGCTTGGGCGAGGGCGAATGGTTGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATc                                                                                                        >  8:267671/1‑149 (MQ=255)                                            gTCACGCAGATCCGCTTGGGCGAGGGCGAATGGTTGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCgg                                                                                                      >  7:36529/1‑149 (MQ=255)                                            gTCACGCAGATCCGCTTGGGCGAGGGCGAATGGTTGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCgg                                                                                                      >  3:170024/1‑149 (MQ=255)                                            gTCACGCAGATCCGCTTGGGCGAGGGCGAATGGTTGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCg                                                                                                       <  5:40697/148‑1 (MQ=255)                                             tCACGCAGATCCGCTTGGGCGAGGGCGAATGGTTGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGc                                                                                                     <  3:216009/149‑1 (MQ=255)                                                 gCAGATCCGCTTGGGCGAGGGCGAATGGTTGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGcc                                                                                                                          <  8:321360/124‑1 (MQ=255)                                                 gCAGATCCGCTTGGGCGAGGGCGAATGGTTGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGcc                                                                                                                          >  7:321360/1‑124 (MQ=255)                                                 gCAGATCCGCTTGGGCGAGGGCGAATGGTTGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCttgctgttg                                                                                                                 <  6:12965/133‑1 (MQ=255)                                                 gCAGATCCGCTTGGGCGAGGGCGAATGGTTGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCttgctgttg                                                                                                                 >  5:12965/1‑133 (MQ=255)                                                         gCTTGGGCGAGGGCGAATGGTTGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTTCAGTTCATCGGCTTGCTGGTGCAc                                                                                         <  7:56662/149‑1 (MQ=255)                                                         gCTTGGGCGAGGGCGAATGGTTGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTGCTGGTg                                                                                            <  7:191915/146‑1 (MQ=255)                                                         ccTTGGGCGAGGGCGAATCGATGTACATCGACTCGCTGTTGCCAGAACCGTACACGCGGCTCGAAGCCTAAAGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGATGACAACGGCCTTGATGTTGATGTTCATCGCCTTGCTGGTg                                                                                            >  8:191915/2‑146 (MQ=255)                                                              ggcgagggcgaATGGTTGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTGCTGGTGCACTGGCa                                                                                    <  4:45979/149‑1 (MQ=255)                                                              ggcgagggcgaATGGTTGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTGCTGGTGCACTGGCa                                                                                    <  7:267671/149‑1 (MQ=255)                                                                cgagggcgaATGCTTGTACATCGCCTCCCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCATGCCGCGCCAGCAGATCGGCTTCATCGAGCTGACCACGGCCTTGCTGTTGCTGTTCAACGGCTTGCTGGTGCACTGGCaaa                                                                                  >  3:108030/1‑149 (MQ=255)                                                                    ggcgaATGGTTGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTGCTGGTGCACTGGCAAAGCTg                                                                              <  4:170024/149‑1 (MQ=255)                                                                     gcgaATGGTTGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCttg                                                                                                                       >  7:343320/1‑107 (MQ=255)                                                                     gcgaATGGTTGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCttg                                                                                                                       <  8:343320/107‑1 (MQ=255)                                                                                tACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTg                                                                                                  <  8:245149/117‑1 (MQ=255)                                                                                tACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTg                                                                                                  >  7:245149/1‑117 (MQ=255)                                                                                tACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTGCTGGTGCACTGGCAAAGCTGGCCACTGAAgc                                                                   <  7:231800/148‑1 (MQ=255)                                                                                  cATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATc                                                                                                                                         <  3:21601/76‑1 (MQ=255)                                                                                  cATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATc                                                                                                                                         >  4:21601/1‑76 (MQ=255)                                                                                   aTCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCATGGCCTTGCTGTTGCTGTTCATCGGCTTGCTg                                                                                               <  5:94387/117‑1 (MQ=255)                                                                                   aTCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCATGGCCTTGCTGTTGCTGTTCATCGGCTTGCTg                                                                                               >  6:94387/1‑117 (MQ=255)                                                                                   aTCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTGCTGGTGCACTGGCAAAGCTGGCCACTGAAGCGCCt                                                               >  2:256132/1‑149 (MQ=255)                                                                                   aTCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTGCTGGTGCACTGGCAAAGCTGGCCACTGAAGCGCCt                                                               >  4:317159/1‑149 (MQ=255)                                                                                     cGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGGTGTTGCTGTTCATCGGCTTGCTGGTGCACTGGCAAAGCTGGCCACTGAAGCGCCTgg                                                             >  6:342223/1‑149 (MQ=255)                                                                                     cGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTGCTGGTGCACTGGCAAAGCTGGc                                                                            >  7:12712/1‑134 (MQ=255)                                                                                     cGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTGCTGGTGCACTGGCAAAGCTGGc                                                                            <  8:12712/134‑1 (MQ=255)                                                                                                 gCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTGCTGGTGCACTGGCAAAGCTGGCCACTGAAGCGCCTg                                                              <  4:217976/136‑1 (MQ=255)                                                                                                 gCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTGCTGGTGCACTGGCAAAGCTGGCCACTGAAGCGCCTg                                                              >  3:217976/1‑136 (MQ=255)                                                                                                  cccGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCgg                                                                                                      <  5:373540/95‑1 (MQ=255)                                                                                                  cccGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCgg                                                                                                      >  6:373540/1‑95 (MQ=255)                                                                                                    cGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTGCTGGTGCACTGGCAAAGCTGGCCACTGAAGCGCCTGGCGCGGGCCGCCCgag                                              <  3:317159/149‑1 (MQ=255)                                                                                                     gAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTGCTGGTGCACTGGCAAAGCTGGCCACTGAAGCGCCTGGCGCGGGCCGCCCgaga                                             <  5:314708/149‑1 (MQ=255)                                                                                                      aaCCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTGCTGGTGCACTGGCAAAGCTGGc                                                                            >  4:166295/1‑117 (MQ=255)                                                                                                              acGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTGCTGGTGCACTGGCAAAGCTGGCCACTGAAGCGCCTGGCGCGGGCCGCCCGAGAGATGTCGCt                                    >  6:329156/1‑149 (MQ=255)                                                                                                              acGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTGCTGGTGCACTGGCAAAGCTGGCCACTGAAGCGCCTGGCGCGGGCCGCCCGAGAGATGTCGCt                                    >  4:352904/1‑149 (MQ=255)                                                                                                                ggggCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTGCTGGTGCACTGGCAAAGCTGGCCACTGAAGCGCCTGGCGCGGGCCGCCCGAGAGATGTCGCTgg                                  >  4:300806/1‑149 (MQ=255)                                                                                                                   cgatgAAGCCGATCTGCTGGCGCGGCAGGCGTTCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTGCTGGTGCACTGGCAAAGCTGGCCACTGAAGCGCCTGGCGCGGGCCGCCCGAGAGATGTCGCTGGGTg                               >  3:367442/5‑149 (MQ=38)                                                                                                                    cTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTGCTGGTGCACTGGCAAAGCTAGCCACTGAAGCGCCTGGCGCGGGCCGCCCGAGAGATGTCGCTGGGTGc                              <  8:36529/149‑1 (MQ=255)                                                                                                                     tGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTGCTGGTGCACTGGCAAAGCTGGCCACTGAAGCGCCTGGCGCGGGCCGCCCGAGAGATGTCGCTGGGTGcc                             >  5:345677/1‑149 (MQ=255)                                                                                                                     atgAAGCCGATCTGCTGGCGCGGCAGGCGTTCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTGCTGGTGCACTGGCAAAGCTGGCCACTGAAGCGCCTGGCGCGGGCCGCCCGAGAGATGTCGCt                                    <  6:44469/140‑1 (MQ=38)                                                                                                                     atgAAGCCGATCTGCTGGCGCGGCAGGCGTTCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTGCTGGTGCACTGGCAAAGCTGGCCACTGAAGCGCCTGGCGCGGGCCGCCCGAGAGATGTCGCt                                    >  5:44469/3‑142 (MQ=38)                                                                                                                       gAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTGCTGGTGCACTGGCAAAGCTGGCCACTGAAGCGCCTGGCGCGGGCCGCCCGAGAGATGTCGCTGGGTGCCGa                           >  4:301402/1‑149 (MQ=255)                                                                                                                       gAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTGCTGGTGCACTGGCAAAGCTGGCCACTGAAGCGCCTGGCGCGGGCCGCCCGAGAGATGTCGCTGGGTGCCGa                           <  4:224234/149‑1 (MQ=255)                                                                                                                       gAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTGCTGGTGCACTGGCAAAGCTGGCCACTGAAGCGCCTGGCGCGGGCCGCCCGAGAGATGTCGCTGGGTGCCGa                           >  1:193765/1‑149 (MQ=255)                                                                                                                             gAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTGCTGGTGCACTGGCAAAGCTGGCCACTGAAGCGCCTGGCGCGGGCCGCCCGAGAGATGTCGCTGGGTGCCGATGTGgc                     >  4:363892/1‑149 (MQ=255)                                                                                                                                   cTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTGCTGGTGCACTGGCAAAGCTGGCCACTGAAGCGCCTGGCGCGGGCCGCCCGAGAGATGTCGCTGGGTGCCGATGTGGCGCCAGt               <  1:256132/149‑1 (MQ=255)                                                                                                                                          gcCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTGCTGGTGCACTGGCAAAGCTGGCCACTGAAGCGCCTGGCGCGGGCCGCCCGAGAGATGTCGCTGGGTGCCGATGTGGCGCCAGTGGCTg          >  6:68114/1‑147 (MQ=255)                                                                                                                                          gcCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTGCTGGTGCACTGGCAAAGCTGGCCACTGAAGCGCCTGGCGCGGGCCGCCCGAGAGATGTCGCTGGGTGCCGATGTGGCGCCAGTGGCTGaa        <  7:46702/149‑1 (MQ=255)                                                                                                                                          gcCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTGCTGGTGCACTGGCAAAGCTGGCCACTGAAGCGCCTGGCGCGGGCCGCCCGAGAGATGTCGCTGGGTGCCGATGTGGCGCCAGTGGCTGaa        <  5:329156/149‑1 (MQ=255)                                                                                                                                          gcCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTGCTGGTGCACTGGCAAAGCTGGCCACTGAAGCGCCTGGCGCGGGCCGCCCGAGAGATGTCGCTGGGTGCCGATGTGGCGCCAGTGGCTGaa        <  2:193765/149‑1 (MQ=255)                                                                                                                                                agATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTGCTGGTGCACTGGCAAAGCTGGCCACTGAAGCGCCTGGCGCGGGCCGCCCGAGAGATGTCGCTGGGTGCCGATGTGGCGCCAGTGGCTGAAGCAggc  >  3:153971/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                                     CCATTATTCGCTGACCTTGGAGCCGCTGAAACCCCCGGTACTGGTCACGCAGATCCGCTTGGGCGAGGGCGAATGGTTGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTGCTGGTGCACTGGCAAAGCTGGCCACTGAAGCGCCTGGCGCGGGCCGCCCGAGAGATGTCGCTGGGTGCCGATGTGGCGCCAGTGGCTGAAGCAGGC  >  NC_002947/1155325‑1155616

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