Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	497,700	G→T	8.0%	L523L (CTG→CTT) ‡	PP_0409 →	two‑component system sensor histidine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	497,700	0	G	T	8.0%	172.4 / 2.1	59	L523L (CTG→CTT) ‡	PP_0409	two‑component system sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base G (30/22);  new base T (2/3);  total (33/26)
Fisher's exact test for biased strand distribution p-value = 6.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATCTGGAGTTCCCCTATATCAGTGAGGGTAGCGAGGCTTTGGTGGGTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTC‑CTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGCCAGGTAGTGTGGGACGGTGTAGTCTGGGACATCACACAGGGC  >  NC_002947/497554‑497844                                                                                                                                                   |                                                                                                                                                  gATCTGGAGTTCCCCTATATCAGTGAGGGTAGAGAGGCTTTGGTGGGTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACctggccctggcc                                                                                                                                                <  2:596606/150‑1 (MQ=255)   tCTGGAGTTCCCCTATATCAGTGAGGGTAGCGAGGCTTTGGTGGGTTATGCACCCAGCGAGATCCAGCACCCGCAGTTGGGCCTGCGCAACCTGGTACACACCGTGGAGCTGGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAg                                                                                                                                              >  1:1163941/1‑150 (MQ=255)   tCTGGAGTTCCCCTATATCAGTGAGGGTAGCGAGGCTTTGGTGGGTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAg                                                                                                                                              >  1:1412694/1‑150 (MQ=255)   tCTGGAGTTCCCCTATATCAGTGAGGGTAGCGAGGCTTTGGTGGGTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCCTGCCAg                                                                                                                                              >  2:568948/1‑150 (MQ=255)           tCCCCTATATCAGTGAGGGTAGCGAGGCTTTGGTGGGTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCCTGCCAG‑GACCAgtc                                                                                                                                     <  1:568948/150‑3 (MQ=255)                  tatCAGTGAGGGTAGCGAGGCTTTGGTGGGTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTc                                                                                                                               <  2:281720/150‑1 (MQ=255)                  tatCAGTGAGGGTAGCGAGGCTTTGGTGGGTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTc                                                                                                                               >  2:136912/1‑150 (MQ=255)                     cAGTGAGGGTAGCGAGGCTTTGGTGGGTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACCAGTC‑CTg                                                                                                                           <  1:914997/150‑1 (MQ=255)                      aGTGAGGGTAGCGAGGCTTTGGTGGGTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTC‑CTgg                                                                                                                          >  1:832720/1‑150 (MQ=255)                      aGTGAGGGTAGCGAGGCTTTGGTGGGTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTC‑CTgg                                                                                                                          <  2:1412694/150‑1 (MQ=255)                         gAGGGTAGCGAGGCTTTGGTGGGTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTC‑CTGGCAg                                                                                                                       <  1:280465/150‑1 (MQ=255)                         gAGGGTAGCGAGGCTTTGGTGGGTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACAACAAGACGGGCACTTAGTAGCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTC‑CTGGCAg                                                                                                                       >  1:168879/1‑150 (MQ=255)                          tgccctgtgtgatgtcccagactacaccgtcccacactacctggccgttgcccagatgccgggcgctggccttgatatcggcccagcgctgttcgccctgccGGGTG‑AGGATGCGCCCCTGGCCAGCGACCAGGACTGGTC‑CTGGCAggg                                                                                                                     >  1:1141817‑M1/103‑150 (MQ=255)                           gggTAGCGAGGCTTTGGTGGGTTATGCACCCAGCGAGATCTAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAAGGGGCAGGTCCGGGCCCTGGCCAGCGAGGAGGACTGGTA‑TTGGCAggg                                                                                                                     >  2:945087/1‑150 (MQ=255)                            ggTAGCGAGGCTTTGGTGGGTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCCGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGTCCCTGCCAG‑GACCAGTCCTGGTCGCTGGccaggg                                                                                                                    <  1:1182492/150‑6 (MQ=255)                             cctgtgtgatgtcccagactacaccgtcccacactacctggccgttgcccagatgccgggcgctggccttgatatcggcccagcgctgttcgccctgccGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTC‑CTGGCAGGGgc                                                                                                                   <  1:236776‑M1/51‑1 (MQ=255)                                 cGAGGCTTTGGTGGGTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTC‑CTGGCAGGGGCGCAt                                                                                                               >  1:1254714/1‑150 (MQ=255)                                 cGAGGCTTTGGTGGGTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTC‑CTGGCAGGGGCGCAt                                                                                                               <  1:262835/150‑1 (MQ=255)                                   tgatgtcccagactacaccgtcccacactacctggccgttgcccagatgccgggcgctggccttgatatcggcccagcgctgttcgccctgccGGGTGCAGGACCTGGTCCTGGCCAGCGACCAGGACTGGTC‑CTGGCAGGGGCGCATcc                                                                                                             >  1:844115‑M1/94‑150 (MQ=255)                                   aGGCTTTGGTGGGTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTC‑CTGGCAGGGGCGCATcc                                                                                                             >  2:602270/1‑150 (MQ=255)                                       tttGGTGGGTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTC‑CTGGCAGGGGCGCATCCTCAc                                                                                                         >  1:1306653/1‑150 (MQ=255)                                           gtggGTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTC‑CTGGCAGGGGCGCATCCTCACCCgg                                                                                                     >  1:1023683/1‑150 (MQ=255)                                               ctacaccgtcccacactacctggccgttgcccagatgccgggcgctggccttgatatcggcccagcgctgttcgccctgccGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTC‑CTGGCAGGGGCGCATCCTCACCCGGCAgg                                                                                                 >  2:236776‑M1/82‑150 (MQ=255)                                                ttATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTC‑CTGGCAGGGGCGCATCCTCACCCGGCAggg                                                                                                <  2:832720/150‑1 (MQ=255)                                                        ccAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTC‑Ct                                                                                                                            >  1:467509/1‑114 (MQ=255)                                                        ccAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTC‑Ct                                                                                                                            <  2:467509/114‑1 (MQ=255)                                                           gCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTC‑CTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAActccaga                                                                                     <  2:1163941/150‑7 (MQ=255)                                                              agaTCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTC‑CTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTggg                                                                                  >  2:1178698/1‑150 (MQ=255)                                                                 tCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACt                                                                                                                                   >  2:443540/1‑99 (MQ=255)                                                                 tCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACt                                                                                                                                   <  1:443540/99‑1 (MQ=255)                                                                 tCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCCTGCCAG‑GACCAGTCCTGGTCGCTGGccagg                                                                                                                     <  1:706438/112‑5 (MQ=25)                                                                 tCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCCTGCCAG‑GACCAGTCCTGGTCGCTGGccagg                                                                                                                     >  2:706438/1‑108 (MQ=25)                                                                  ccAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTC‑CTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGa                                                                              <  1:136912/150‑1 (MQ=255)                                                                   cAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTAGCCAGCGACCAGGACTGGTC‑CTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGat                                                                             >  1:25330/1‑150 (MQ=255)                                                                      ccgttgcccagatgccgggcgctggccttgatatcggcccagcgctgttcgccctgccGGGTG‑AGGATGCGCCCCTGGCCAGCGACCAGGACTGGTC‑CTGGCAg                                                                                                                       >  1:1081016‑M1/59‑104 (MQ=255)                                                                      ccgttgcccagatgccgggcgctggccttgatatcggcccagcgctgttcgccctgccGGGTG‑AGGATGCGCCCCTGGCCAGCGACCAGGACTGGTC‑CTGGCAg                                                                                                                       <  2:1081016‑M1/46‑1 (MQ=255)                                                                         ccGCAGATGGGCCTGCGCAACCTGGTGTACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAG‑GACCAGTCCTGGTcgctg                                                                                                                           >  1:669608/1‑94 (MQ=255)                                                                         ccGCAGATGGGCCTGCGCAACCTGGTGTACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAG‑GACCAGTCCTGGTcgctg                                                                                                                           <  2:669608/98‑5 (MQ=255)                                                                         ccGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTC‑CTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAg                                                                       <  2:1254714/150‑1 (MQ=255)                                                                         ccGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTC‑CTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAg                                                                       <  2:1392164/150‑1 (MQ=255)                                                                           gCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTC‑CTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACa                                                                                          >  1:1056554/1‑129 (MQ=255)                                                                           gCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTC‑CTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACa                                                                                          <  2:1056554/129‑1 (MQ=255)                                                                             aGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGCCCTGGCCAACGCCAGCGACCAGGACTGGTC‑CTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCa                                                                   <  1:602270/150‑1 (MQ=255)                                                                                         gcAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTC‑CTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATc                                                       >  2:522007/1‑150 (MQ=255)                                                                                              cTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTC‑CTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGc                                                  >  1:90599/1‑150 (MQ=255)                                                                                                     aCCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTC‑CTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGcc                                           >  2:897386/1‑150 (MQ=255)                                                                                                      ccccGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTC‑CTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGCCa                                          >  2:1376161/1‑150 (MQ=255)                                                                                                      ccccGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTC‑CTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGCCa                                          >  2:248584/1‑150 (MQ=255)                                                                                                                cGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTC‑CTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGCCAGGTAGTGTgg                                >  1:1117582/1‑150 (MQ=255)                                                                                                                  tGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTC‑CTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGat                                                                             <  1:1396067/103‑1 (MQ=255)                                                                                                                  tGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTC‑CTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGat                                                                             >  2:1396067/1‑103 (MQ=255)                                                                                                                           tCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTC‑CTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTg                                                     <  1:652481/118‑1 (MQ=255)                                                                                                                           tCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTC‑CTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTg                                                     >  2:652481/1‑118 (MQ=255)                                                                                                                           tCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTC‑CTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGCCAGGTAGTGTGGGACGGTGTAGt                     <  1:522007/150‑1 (MQ=255)                                                                                                                              caGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTC‑CTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGCCAGGTAGTGTGGGACGGTGTAGTCTg                  <  1:1178698/150‑1 (MQ=255)                                                                                                                                 ggTGCAGGACCTCGCCCTGCCCAGAGACCATGGCTGCTC‑CTGGCAGCGGCGCCTGATCACCCGCCAGTGCCAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGCGCCACGGCCAGTTAGTGTGGGGCGGAGACGTCCGGGa               >  2:1186560/1‑150 (MQ=255)                                                                                                                                 GGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTC‑CTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGCCAGGTAGTGTGGGACGGTGTAGTCTgg                 <  1:990585‑M1/148‑1 (MQ=255)                                                                                                                                 GGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTC‑CTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGCCAGGTAGTGTGGGACGGTGTAGTCTGGGa               >  1:1376161‑M1/1‑150 (MQ=255)                                                                                                                                 GGTGCAGGACCTGCCCCTGGCCAGCGACCAGGACTGGTC‑CTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGCCAGGTAGTGTGGGACGGTGTAGTCTgg                 >  2:990585‑M1/1‑148 (MQ=255)                                                                                                                                     cAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTC‑CTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGCCATGTAGTGTGGGACGGTGTAGTCTGGGACATc           <  2:168879/150‑1 (MQ=255)                                                                                                                                           ctggccctggccAGCGACCAGGACTGGTC‑CTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGCCAGGTAGTGTGGGACGGTGTAGTCTGGGACATCACACAg     >  2:200684/1‑150 (MQ=255)                                                                                                                                           ctggccctggccAGCGACCAGGACTGGTC‑CTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGCCAGGTAGTGTGGGACGGTGTAGTCTGGGACATCACACAg     <  1:945087/150‑1 (MQ=255)                                                                                                                                             ggccctggccAGCGACCAGGACTGGTC‑CTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGCCAGGTAGTGTGGGACGGTGTAGTCTGGGACATCACACAggg   <  1:897386/150‑1 (MQ=255)                                                                                                                                              gccctggccAGCGACCAGGACTGGTC‑CTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGCCa                                          <  2:1335969/110‑1 (MQ=255)                                                                                                                                              gccctggccAGCGACCAGGACTGGTC‑CTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGCCa                                          >  1:1335969/1‑110 (MQ=255)                                                                                                                                              gccctggccAGCGACCAGGACTGGTC‑CTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGCCAGGTAGTGTGGGACGGTGTAGTCTGGGACATCACACAggag  >  2:299813/1‑148 (MQ=255)                                                                                                                                                   |                                                                                                                                                  GATCTGGAGTTCCCCTATATCAGTGAGGGTAGCGAGGCTTTGGTGGGTTATGCACCCAGCGAGATCCAGCACCCGCAGATGGGCCTGCGCAACCTGGTGCACCCCGAGGACCGTGCTGATTATCACAGGGTGCAGGACCTGGCCCTGGCCAGCGACCAGGACTGGTC‑CTGGCAGGGGCGCATCCTCACCCGGCAGGGCGAACAGCGCTGGGCCGATATCAAGGCCAGCGCCCGGCATCTGGGCAACGGCCAGGTAGTGTGGGACGGTGTAGTCTGGGACATCACACAGGGC  >  NC_002947/497554‑497844

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑