Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,266,698	A→T	11.0%	A432A (GCT→GCA)	quiP ←	acyl‑homoserine lactone acylase QuiP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,266,698	0	A	T	11.0%	180.4 / 13.2	64	A432A (GCT→GCA)	quiP	acyl‑homoserine lactone acylase QuiP
Reads supporting (aligned to +/- strand):  ref base A (26/31);  new base T (5/2);  total (31/33)
Fisher's exact test for biased strand distribution p-value = 2.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.23e-01

GCATCGGGTCGGCATAGCCGTCCCAGTCGTAACGCCCGTCCCAGCCCGGCGAGGGCAGCAGGCCCTGGCCTTCGCGGCGATTGGGGTAGCGGCCACTGACTTGCCAGCCGATGTGCTCAGGCTCGGCAAAAACGAAGTTGACGGCAGCAGCCGTC‑ACTTCACGTGTGCTGTCGAAGGCCCGCTCCACGTTCTTGGCCCGCGTCAGGTCGAACAGGGCATCCAGGCTGCGATCACCTTTGAACTGGGGCAGGTTCAGCGCCAGGCCCAGGCTGGCGTTTTCTGGCTGGGCCAGCAGG  >  NC_002947/1266550‑1266845                                                                                                                                                     |                                                                                                                                                     gcaTCGGGTCGGCATAGCCGTCCCAGTCGTAACGCCCGTCCCAGCCCGGCGAGGGCAGCAGGCCCTGGCCTTCGCGGCGATTGGGGTAGCGGCCACTGACTTGCCAGCCGATGTGCTCAGGCTGGGCAAAAACGAAGTTGACGGcagcag                                                                                                                                                     >  1:941473/1‑150 (MQ=255)        gtgggCATAGCCGTCCCAGTCGTAACGCCCGTCCCAGCCCGGCGAGGGCAGCAGGCCCTGGCCTTCGCGGCGATTGGGGTAGCGGCCACTGACTTGCCAGCCGATGTGCTCAGGCTCGGCAAAAACGAAGTTGACGGCAGCAGCCGTC‑Ac                                                                                                                                             <  2:1108124/147‑1 (MQ=255)        gTCGGCATAGCCGTCCCAGTCGTAACGCCCGTCCCAGCCCGGCGAGGGCAGCAGGCCCTGGCCTTCGCGGCGATTGGGGTAGCGGCCACTGACTTGCCAGCCGATGTGCTCAGGCTCGGCAAAAACGAAGTTGACGGCAGCAGCCGTC‑Ac                                                                                                                                             >  2:280542/1‑150 (MQ=255)             cATAGCCGTCCCAGTCGTAACGCCCGTCCCAGCCCGGCGAGGGCAGCAGGCCCTGGCCTTCGCGGCGATTGGGGTAGCGGCCACTGACTGGCCAGCCGATGTGCTCAGGCTCGGCAAAAACGAAGTTGACGGCAGCAGCCGTC‑ACTTCAc                                                                                                                                        <  1:342774/150‑1 (MQ=255)              aTAGCCGTCCCAGTCGTAACGCCCGTCCCAGCCCGGCGAGGGCAGCAGGCCCTGGCCTTCGCGGCGATTGGGGTAGCGGCCACTGACTTGCCAGCCGATGTGCTCAGGCTCGGCAAAAACGAAGTTGACGGCAGCAGCCGTC‑ACTTCACg                                                                                                                                       <  2:1872051/150‑1 (MQ=255)                  ccGTCCCAGTCGTAACGCCCGTCCCAGCCCGGCGAGGGCAGCAGGCCCTGGCCTTCGCGGCGATTGGGGTAGCGGCCACTGACTTGCCAGCCGATGTGCTCAGGCTCGGCAAAAACGAAGTTGACGGCAGCAGCCGTC‑ACTTCACGtgtg                                                                                                                                   <  1:280542/150‑1 (MQ=255)                         aGTCGTAACGCCCGTCCCAGCCCGGCGAGGGCAGCAGGCCCTGGCCTTCGCGGCGATTGGGGTAGCGGCCACTGACTTGCCAGCCGATGTGCTCAGGCTCGGCAAAAACGAAGTTGACGGCACCAGCCGTC‑ACTTCACGTGTGCTGTCGa                                                                                                                            <  1:51915/150‑1 (MQ=255)                          gTCGTAACGCCCGTCCCAGCCCGGCGAGGGCAGCAGGCCCTGGCCTTCGCGGCGATTGGGGTAGCGGCCACTGACTTGCCAGCCGATGTGCTCAGGCTCGGCAAAAACGAAGTTGACGGCAGCAGCCGTC‑ACTTCACGTGTGCTGTCGaa                                                                                                                           <  1:883958/150‑1 (MQ=255)                          gTCGTAACGCAAGTTCTAGCCCGGCGAGGGCAGCAGGCCGTGGCTTTCGCGGCGATTGGGGTAGCGGCTACTTACTTGCCAGCCGATGTGCTCAGTCTCGGCAAAAACGAAGTTGACGGCAGCAGCCGTC‑ACTTCACGTGTGCTGTCGaa                                                                                                                           <  1:1832057/150‑1 (MQ=255)                             gTAACGCCCGTCCCAGCCCGGCGAGGGCAGCAGGCCCTGGCCTTCGCGGCGATTGGGGTAGCGGCCACTGACTTGCCAGCCGATGTGCTCAGGCTCGGCAAAAACGAAGTTGACGGCAGCAGCCGTC‑ACTTCACGTGTGCTGTCGAAGGc                                                                                                                        >  1:164206/1‑150 (MQ=255)                                  gCCCGTCCCAGCCCGGCGAGGGCAGCAGGCCCTGGCCTTCGCGGCGATTGGGGTAGCGGCCACTGACTTGCCAGCCGATGTGCTCAGGCTCGGCAAAAACGAAGTTGACGGCAGCAGCCGTC‑ACTTCACGTGTGCTGTCGAAGGCCCGCt                                                                                                                   >  1:868178/1‑150 (MQ=255)                                  gCCCGTCCCAGCCCGGCGAGGGCAGCAGGCCCTGGCCTTCGCGGCGATTGGGGTAGCGGCCACTGACTTGCCAGCCGATGTGCTCAGGCTCGGCAAAAACGAAGTTGACGGCAGCAGCCGTC‑ACTTCACGTGTGCTGTCGAAGGCCCGCt                                                                                                                   >  1:1599284/1‑150 (MQ=255)                                  gCCCGTCCCAGCCCGGCGAGGGCAGCAGGCCCTGGCCTTCGCGGCGATTGGGGTAGCGGCCACTGACTTGCCAGCCGATGTGCTCAGGCTCGGCAAAAACGAAGTTGACGGCAGCAGCCGTC‑ACTTCACGTGTGCTGTCGAAGGCCCGCt                                                                                                                   >  2:1616054/1‑150 (MQ=255)                                          cAGCCCGGCGAGGGCAGCAGGCCCTGGCCTTCGCGGCGATTGGGGTAGCGGCCACTGACTTGCCAGCCGATGTGCTCAGGCTCGGCAAAAACGAAGTTGACGGCAGCAGCCGTC‑ACTTCACGTGTGCTGTCGAAGGCCCGCTCCACGTTc                                                                                                           >  2:1523337/1‑150 (MQ=255)                                          cAGCCCGGCGAGGGCAGCAGGCCCTGGCCTTCGCGGCGATTGGGGTAGCGGCCACTGACTTGCCAGCCGATGTGCTCAGGCTCGGCAAAAACGAAGTTGACGGCAGCAGCCGTC‑ACTTCACGTGTGCTGTCGAAGGCCCGCTCCACGTTc                                                                                                           <  2:135767/150‑1 (MQ=255)                                          cAGCCCGGCGAGGGCAGCAGGCCCTGGCCTTCGCGGCGATTGGGGTAGCGGCCACTGACTTGCCAGCCGATGTGCTCAGGCTCGGCAAAAACGAAGTTGACGGCAGCAGCCGTC‑ACTTCACGTGTGCTGTCGAAGGCCCGCTCCACGTTc                                                                                                           >  1:1677527/1‑150 (MQ=255)                                             cccGGCGAGGGCAGCAGGCCCTGGCCTTCGCGGCGATTGGGGTAGCGGCCACTGACTTGCCAGCCGATGTGCTCAGGCTCGGCAAAAACGAAGTTGACGGCAGCAGCCGTC‑ACTTCACGTGTGCTGTCGAAGGCCCGCTCCACGTTCTTg                                                                                                        <  2:788433/150‑1 (MQ=255)                                             cccGGCGAGGGCAGCAGGCCCTGGCCTTCGCGGCGATTGGGGTAGCGGCCACTGACTTGCCAGCCGATGTGCTCAGGCTCGGCAAAAACGAAGTTGACGGCAGCAGCCGTC‑ACTTCACGTGTGCTGTCGAAGGCCCGCTCCACGTTCTTg                                                                                                        <  2:1677527/150‑1 (MQ=255)                                              ccGGCGAGGGCAGCAGGCCCTGGCCTTCGCGGCGATTGGGGTAGCGGCCACTGACTTGCCAGCCGATGTGCTCAGGCTCGGCAAAAACGAAGTTGACGGCAGCAGCCGTC‑ACTTCACGTGTGCTGTCGAAGGCCCGCTCCACGTTCTTgg                                                                                                       >  1:1219693/1‑150 (MQ=255)                                                        cagcagGCCCTGGCCTTCGCGGCGATTGGGGTAGCGGCCACTGACTTGCCAGCCGATGTGCTCAGGCTCGGCAAAAACGAAGTTGACGGCAGCAGCCGTC‑ACTTCACGTGTGCTGTCGAAGGCCCGCTCCACGTTCTTGGCCCGCGTCAg                                                                                             <  1:1523337/150‑1 (MQ=255)                                                           cagGCCCTGGCCTTCGCGGCGATTGGGGTAGCGGCCACTGACTTGCCAGCCGATGTGCTCAGGCTCGGCAAAAACGAAGTTGACGGCAGCAGCCGTC‑ACTTCACGTGTGCTGTCGAAGGCCCGCTCCACGTTCTTGGCCCGCGTCAGGTc                                                                                          <  2:164206/150‑1 (MQ=255)                                                                ccTGGCCTTCGCGGCGATTGGGGTAGCGGCCACTGACTTGCCAGCCGATGTGCTCAGGCTCGGCAAAAACGAAGTTGACGGCAGCAGCCGTC‑ACTTCACGTGTGCTGTCGAAGGCCCGCTCCACGTTCTTGGCCCGCGTCAGGTCGAACa                                                                                     >  2:1497070/1‑150 (MQ=255)                                                                ccTGGCCTTCGCGGCGATTGGGGTAGCGGCCACTGACTTGCCAGCCGATGTGCTCAGGCTCGGCAAAAACGAAGTTGACGGCAGCAGCCGTC‑ACTTCACGTGTGCTGTCGAAGGCCCGCTCCACGTTCTTGGCCCGCGTCAGGTCGAACa                                                                                     >  2:1612215/1‑150 (MQ=255)                                                                   ggCCTTCGCGGCGATTGGGGTAGCGGCCAGTGACTTGCCAGCCGATGTGCGCAGGCTCGGCAAAAACGAAGTTGACGGCAGCAGCCGTC‑ACTTCACGTGTGCTGTCGAAGGCCCGCTCCACGTTCTTGGCCCGCGTCAGGTCGAACAggg                                                                                  <  1:1648212/150‑1 (MQ=255)                                                                       ttCGCGGCGATTGGGGTAGCGGCCACTGACTTGCCAGCCGATGTGCTCAGGCTCGGCAAAAACGAAGTTGACGGCAGCAGCCGTC‑ACTTCACGTGTGCTGTCGAAGGCCCGCTCCACGTTCTTGGCCCGCGTCAGGTCGAACAGGGCATc                                                                              <  2:1219693/150‑1 (MQ=255)                                                                           cggcgATTGGGGTAGCGGCCACTGACTTGCCAGCCGATGTGCTCAGGCTCGGCAAAAACGAAGTTGACGGCAGCAGCCGTC‑ACTTCACGTGTGCTGTCGAAGGCCCGCTCCACGTTCTTGGCCCGCGTCAGGTCGAACAGGGCATCCAgg                                                                          >  1:458981/1‑150 (MQ=255)                                                                             gcgATTGGGGTAGCGGCCACTGACTTGCCAGCCGATGTGCTCAGGCTCGGCAAAAACGAAGTTGACGGCAGCAGCCGTC‑ACTTCACGTGTGCTGTCGAAGGCCCGCTCCACGTTCTTGGCCCGCGTCAGGTCGAACAGGGCATCCAGGCt                                                                        <  2:458981/150‑1 (MQ=255)                                                                             gcgATTGGGGTAGCGGCCACTGACTTGCCAGCCGATGTGCTCAGGCTCAGCAAAAACGAAGTTGACCGCAGCAGGCGCC‑GCTTCACGTGTGCTGTCGAAGGCCCGCTCCACGTTCTTGATCCGCGCCAGGTCGAAATGGGCAACCATGcg                                                                        >  2:407249/1‑149 (MQ=255)                                                                               gATTGGGGTAGCGGCCACTGACTTGCCAGCCGATGTGCTCAGGCTCGGCAAAAACGAAGTTGACGGCAGCAGCCGTC‑ACTTCACGTGTGCTGTCGAAGGCCCGCTCCACGTTCTTGGCCCGCGTCAGGTCGAACAgg                                                                                   >  2:1489119/1‑137 (MQ=255)                                                                               gATTGGGGTAGCGGCCACTGACTTGCCAGCCGATGTGCTCAGGCTCGGCAAAAACGAAGTTGACGGCAGCAGCCGTC‑ACTTCACGTGTGCTGTCGAAGGCCCGCTCCACGTTCTTGGCCCGCGTCAGGTCGAACAgg                                                                                   <  1:1489119/137‑1 (MQ=255)                                                                                aTTGGGGTAGCGGCCACTGACTTGCCAGCCGATGTGCTCAGGCTCGGCAAAAACGAAGTTGACGGCAGCAGCCGTC‑ACTTCACGTGTGCTGTCGAAGGCCCGCTCCACGTTCTTg                                                                                                        >  1:325517/1‑115 (MQ=255)                                                                                aTTGGGGTAGCGGCCACTGACTTGCCAGCCGATGTGCTCAGGCTCGGCAAAAACGAAGTTGACGGCAGCAGCCGTC‑ACTTCACGTGTGCTGTCGAAGGCCCGCTCCACGTTCTTg                                                                                                        <  2:325517/115‑1 (MQ=255)                                                                                aTTGGGGTAGCGGCCACTGACTTGCCAGCCGATGTGCTCAGGCTCGGCAAAAACGAAGTTGACGGCAGCAGCCGTC‑ACTTCACGTGTGCTGTCGAAGGCCCGCTCCACGTTCTTg                                                                                                        <  2:326490/115‑1 (MQ=255)                                                                                aTTGGGGTAGCGGCCACTGACTTGCCAGCCGATGTGCTCAGGCTCGGCAAAAACGAAGTTGACGGCAGCAGCCGTC‑ACTTCACGTGTGCTGTCGAAGGCCCGCTCCACGTTCTTg                                                                                                        >  1:326490/1‑115 (MQ=255)                                                                                         gCGGCCACTGACTTGCCAGCCGATGTGCTCAGGCTCGGCAAAAACGAAGTTGACGGCAGCAGCCGTC‑ACTTCACGTGTGCTGTCGAAGGCCCGCTCCACGTTCTTGGCCCGCGTCAGGTCGAACAGGGCATCCAGGCTGCGATCACCttt                                                            <  1:1497070/150‑1 (MQ=255)                                                                                           ggCCACTGACTTGCCAGCCGATGTGCTCAGGCTCGGCAAAAACGAAGTTGACGGCAGCAGCCGTC‑ACTTCACGTGTGCTGTCGAAGGCCCGCTCCACGTTCTTGGCCCGCGTCAGGTCGAACAGGGCATCCAGGCTGCGATCACCTTTGa                                                          >  2:1480618/1‑150 (MQ=255)                                                                                              cACTGCCTCGCCCGCCGATGTGCTCAGGCTCGGCAAAGACGAAGTTGACGGCAGCAGCCGTC‑ACTTCACGTGTGCTGTCGACGGGCCGCTCCCCGTTCTTGCCCCACGTGAGGTCGAACAGGGCGTCCAGGCTGGGATCAGCTTTGGACt                                                       >  2:1504952/1‑150 (MQ=255)                                                                                               aCTGACTTGCCAGCCGATGTGCTCAGGCTCGGCAAAAACGAAGTTGACGGCAGCAGCCGTC‑ACTTCACGTGTGCTGTCGAAGGCCCGCTCCACGTTCTTGGCCCGCGTCAGGTCGAACAGGGCATCCAGGCTGCGATCACCTTTGAACTg                                                      <  2:1599284/150‑1 (MQ=255)                                                                                                cTGACTTGCCAGCCGATGTGCTCAGGCTCGGCAAAAACGAAGTTGACATTTGCAGCCGTC‑ACTTCACGTGTGCTGTCGAAGGCCCGCTCCACGTTCTTGGCCCGCGTCAGGTCGAACAGGGCATCCAGGCTGCGATCACCTTTGAACTgg                                                     <  1:117178/150‑1 (MQ=255)                                                                                                 tGACTTGCCAGCCGATGTGCTCAGGCTCGGCAAAAACGAAGTTGACGGCAGCAGCCGTC‑ACTTCACGTGTGCTGTCGAAGGCCCGCTCCACGTTCTTGGCCCGCGTCAGGTCGAACAGGGCATCCAGGCTGCGATCACCTTTGAACTggg                                                    >  1:1197334/1‑150 (MQ=255)                                                                                                 tGACTTGCCAGCCGATGTGCTCAGGCTCGGCAAAAACGAAGTTGACGGCAGCAGCCGTC‑ACTTCACGTGTGCTGTCGAAGGCCCGCTCCACGTTCTTGGCCCGCGTCAGGTCGAACAGGGCATCCAGGCTGCGATCACCTTTGAACTggg                                                    >  1:248154/1‑150 (MQ=255)                                                                                                 ggACTTGCCAGCCGATGTGCTCAGGCTCGGCAAAAACGAAGTTGACGGCAGCAGCCGTC‑ACTTCACGTGTGCTGTCGAAGGCCCGCTCCACGTTCTTGGCCCGCGTCAGGTCGAACAGGGCATCCAGGCTGCGATCACCTTTGAACTggg                                                    <  2:677938/149‑1 (MQ=255)                                                                                                     ttGCCAGCCGATGTGCTCAGGCTCGGCAAAAACGAAGTTGAGGGTAGCAGCCGTC‑ACTTCACGTGTGCTGTCGAAGGCCCGCTCCACGTTCTTGGCCCGCGTCAGGTCGAACAGGGCATCCAGGCTGCGATCACCTTTGAACTGGGGCAg                                                <  1:407249/150‑1 (MQ=255)                                                                                                     ttGCCAGCCGATGTGCTCAGGCTCGGCAAAAACGAAGTTGACGGCAGCAGCCGTC‑ACTTCACGTGTGCTGTCGAAGGCCCGCTCCACGTTCTTGGCCCGCGTCAGGTCGAACAGGGCATCCAGGCTGCGATCACCTTTGAACTGGGGCAg                                                >  2:932126/1‑150 (MQ=255)                                                                                                        ccAGCCGATGTGCTCAGGCTCGGCAAAAACGAAGTTGACGGCAGCAGCCGTC‑ACTTCACGTGTGCTGTCGAAGGCCCGCTCCACGTTCTTGGCCCGCGTCAGGTCGAACAGGGCATCCAGGCTGCGATCACCTTTGAACTGGGGCAGGtt                                             >  1:489224/1‑150 (MQ=255)                                                                                                         cAGCCGATGTGCTCAGGCTCGGCAAAAACGAAGTTGACGGCAGCAGCCGTC‑ACTTCACGTGTGCTGTCGAAGGCCCGCTCCACGTTCTTGGCCCGCGTCAGGTCGAACAGGGCATCCAGGCTGCGATCACCTTTGAACTGGGGCAGGTTc                                            >  2:1431795/1‑150 (MQ=255)                                                                                                          aGCCGATGTGCTCAGGCTCGGCAAAAACGAAGTTGACGGCAGCAGCCGTC‑ACTTCACGTGTGCTGTCGAAGGCCCGCTCCACGTTCTTGGCCCGCGTCAGGTCGAACAGGGCATCCAGGCTGCGATCACCTTTGAACTGGGGCAGGTTCa                                           >  2:1478690/1‑150 (MQ=255)                                                                                                          aGCCGATGTGCTCAGGCTCGGCAAAAACGAAGTTGACGGCAGCAGCCGTC‑ACTTCACGTGTGCTGTCGAAGGCCCGCTCCACGTTCTTGGCCCGCGTCAGGTCGAACAGGGCATCCAGGCTGCGATCACCTTTGAACTGGGGCAGGTTCa                                           >  1:453814/1‑150 (MQ=255)                                                                                                            ccGATGTGCTCAGGCTCGGCAAAAACGAAGTTGACGGCAGCAGCCGTC‑ACTTCACGTGTGCTGTCGAAGGCCCGCTCCACGTTCTTGGCCCGCGTCAGGTCGAACAGGGCATCCAGGCTGCGATCACCTTTGAACTGGGGCAGGTTCAgc                                         <  2:489224/150‑1 (MQ=255)                                                                                                                    cTCAGGCTCGGCAAAAACGAAGTTGACGGCAGCAGCCGTC‑ACTTCACGTGTGCTGTCGAAGGCCCGCTCCACGTTCTTGGCCCGCGTCAGGTCGAACAGGGCATCCAGGCTGCGATCACCTTTGAACTGGGGCAGGTTCAGCGccaggcc                                 >  2:315103/1‑150 (MQ=255)                                                                                                                      cAGGCTCGGCAAAAACGAAGTTGACGGCAGCAGCCGTC‑ACTTCACGTGTGCTGTCGAAGGCCCGCTCCACGTTCTTGGCCCGCGTCAGGTCGAACAGGGCATCCAGGCTGCGATCACCTTTGAACTGGGGCAGGTTCAGCGccaggccca                               <  1:360157/150‑1 (MQ=255)                                                                                                                        ggCTCGGCAAAAACGAAGTTGACGGCGGCAGCCGTC‑ACTTCACGTGTGCTGTCGAAGGCCCGCTCCACGTTCTTGGCCCGCGTCAGGTCGAACAGGGCATCCAGGCTGCGATCACCTTTGAACTGGGGCAGGTTCAGCGccaggcccagg                             <  2:69717/150‑1 (MQ=255)                                                                                                                         gCTCGGCAAAAACGAAGTTGACGGCAGCAGCCGTC‑ACTTCACGTGTGCTGTCGAAGGCCCGCTCCACGTTCTTGGCCCGCGTCAGGTCGAACAGGGCATCCAGGCTGCGATCACCTTTGAACTGGGGCAGGTTCAGCGccaggcccag                              <  1:1901866/148‑1 (MQ=255)                                                                                                                         gCTCGGCAAAAACGAAGTTGACGGCAGCAGCCGTC‑ACTTCACGTGTGCTGTCGAAGGCCCGCTCCACGTTCTTGGCCCGCGTCAGGTCGAACAGGGCATCCAGGCTGCGATCACCTTTGAACTGGGGCAGGTTCAGCGccaggcccag                              >  2:1901866/1‑148 (MQ=255)                                                                                                                          cTCGGCAAAAACGAAGTTCACGGCAGCAGCCGTC‑ACTTCACGTGTGCTGTCGAAGGCCCTCTCCACGTTCTTGGCCCGCGTCAGGTCGAACAGGGCATCCAGGCTGCGATCACCTTTGAACTGGGGCAGGTTCAGCGCCAGGCCCAGGCt                           <  1:934613/150‑1 (MQ=255)                                                                                                                           tCGGCAAAAACGAAGTTGACGGCAGCAGCCGTC‑ACTTCACGTGTGCTGTCGAAGGCCCGCTCCACGTTCTTGGCCCGCGTCAGGTCGAACAGGGCATCCAGGCTGCGATCACCTTTGAACTGGGGCAGGTTCAGCGCCAGGCCCAGGCTg                          <  1:1612215/150‑1 (MQ=255)                                                                                                                           tCGGCAAAAACGAAGTTGACGGCAGCAGCCGTC‑ACTTCACGTGTGCTGTCGAAGGCCCGCTCCACGTTCTTGGCCCGCGTCAGGTCGAACAGGGCATCCAGGCTGCGATCACCTTTGAACTGGGGCAGGTTCAGCGCCAGGCCCAGGCTg                          <  1:1616054/150‑1 (MQ=255)                                                                                                                                     gtgaagtGACGGCTGCTGCCGTC‑ACTTCACGTGTGCTGTCGAAGGCCCGCTCCACGTTCTTGGCCCGCGTCAGGTCGAACAGGGCATCCAGGCTGCGATCACCTTTGAACTGGGGCAGGTTCAGCGCCAGGCCCAGGCTGGCGTTTTctg                >  1:1028186/7‑150 (MQ=255)                                                                                                                                         gTTGACGGCAGCAGCCGTC‑ACTTCACGTGTGCTGTCGAAGGCCCGCTCCACGTTCTTGGCCCGCGTCAGGTCGAACAGGGCATCCAGGCTGCGATCACCTTTGAACTGGGGCAGGTTCAGCGCCAGGCCCAGGCTGGCGTTTTctggctg            <  2:102168/150‑1 (MQ=255)                                                                                                                                          gtGACGGCTGCTGCCGTCAACTTCACGTGTGCTGTCGAAGGCCCGCTCCACGTTCTTGGCCcgc                                                                                                  >  2:210412/2‑64 (MQ=25)                                                                                                                                          gtGACGGCTGCTGCCGTCAACTTCACGTGTGCTGTCGAAGGCCCGCTCCACGTTCTTGGCCcgc                                                                                                  <  1:210412/63‑1 (MQ=25)                                                                                                                                          gtGACGGCTGCTGCCGTCAACTTCACGTGTGCTGTCGAAGGCCCGCTCCACGTTCTTGGCCCGCGTCAGGTCGAACAGGGCATCCAGGCTGCGATCACCTTTGAACTGGGGCAGGTTCAGCGCCAGGCCCAGGCTGGCGTTTTctggctg            >  2:522546/2‑150 (MQ=255)                                                                                                                                           tGACGGCTGCTGCCGTCAACTTCGCGTGTGCTGTCGAAGGCCCGCTCCACGTTCTTGGCCCGCGTCAGGTCGAACAGGGCATCCAGGCTGCGATCACCTTTGAACTGGGGCAGGTTCAGCGCCAGGCCCAGGCTGGCGTTTTctggctgg           >  2:376116/1‑150 (MQ=255)                                                                                                                                             aCGGCTGCTGCCGTC‑ACTTCACGTGTGCTGTCGAAGGCCCGCTCCACGTTCTTGGCCCGCGTCAGGTCGAACAGGGCATc                                                                              >  2:1373866/1‑80 (MQ=255)                                                                                                                                             aCGGCTGCTGCCGTC‑ACTTCACGTGTGCTGTCGAAGGCCCGCTCCACGTTCTTGGCCCGCGTCAGGTCGAACAGGGCATc                                                                              <  1:1373866/80‑1 (MQ=255)                                                                                                                                                gcagcaGCCGTC‑ACTTCACGTGTGCTGTCGAAGGCCCGCTCCACGTTCTTGGCCCGCGTCAGGTCGAACAGGGCATCCAGGc                                                                         >  1:1798130/1‑82 (MQ=255)                                                                                                                                                gcagcaGCCGTC‑ACTTCACGTGTGCTGTCGAAGGCCCGCTCCACGTTCTTGGCCCGCGTCAGGTCGAACAGGGCATCCAGGc                                                                         <  2:1798130/82‑1 (MQ=255)                                                                                                                                                  agcagcCGTC‑ACTTCACGTGTGCTGTCGAAGGCCCGCTCCACGTTCTTGGCCCGCGTCAGGTCGAACAGGGCATCCAGGCTGCGATCACCTTTGAACTGGGGCAGGTTCAGCGCCAGGCCCAGGCTGGCGTTTTCTGGCTGGGCcagcag   >  1:308958/1‑150 (MQ=255)                                                                                                                                                  agcagcCGTC‑ACTTCACGTGTGCTGTCGAAGGCCCGCTCCACGTTCTTGGCCCGCGTCAGGTCGAACAGGGCATCCAGGCTGCGATCACCTTTGAACTGGGGCAGGTTCAGCGCCAGGCCCAGGCTGGCGTTTTCTGGCTGGGCcagcag   >  1:42067/1‑150 (MQ=255)                                                                                                                                                   gcagcCGTC‑ACTTCACGTGTGCTGTCGAAGGCCCGCTCCACGTTCTTGGCCCGCGTCAGGTCGAACAGGGCATCCAGGCTGCGATCACCTTTGAACTGGGGCAGGTTCAGCGCCAGGCCCAGGCTGGCGTTTTCTGGCTGGGCCAGCAgg  >  1:1444225/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                                     GCATCGGGTCGGCATAGCCGTCCCAGTCGTAACGCCCGTCCCAGCCCGGCGAGGGCAGCAGGCCCTGGCCTTCGCGGCGATTGGGGTAGCGGCCACTGACTTGCCAGCCGATGTGCTCAGGCTCGGCAAAAACGAAGTTGACGGCAGCAGCCGTC‑ACTTCACGTGTGCTGTCGAAGGCCCGCTCCACGTTCTTGGCCCGCGTCAGGTCGAACAGGGCATCCAGGCTGCGATCACCTTTGAACTGGGGCAGGTTCAGCGCCAGGCCCAGGCTGGCGTTTTCTGGCTGGGCCAGCAGG  >  NC_002947/1266550‑1266845

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