Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,141,212	T→A	6.9%	L44Q (CTG→CAG)	PP_3645 →	L‑proline amide hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,141,212	0	T	A	6.9%	253.8 / 7.3	84	L44Q (CTG→CAG)	PP_3645	L‑proline amide hydrolase
Reads supporting (aligned to +/- strand):  ref base T (39/39);  new base A (4/2);  total (43/41)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-01

TTCAAGGAGCGCCCTCCATGATTCAGCCTGTCGCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGACCACCTGCAAATTGCACAAGCCACCGTGATCGGCT  >  NC_002947/4141065‑4141361                                                                                                                                                    |                                                                                                                                                        ttCAAGGAGCGCCCTCCATGATTCAGCCTGTCGCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTtggcctgg                                                                                                                                                       >  1:2241836/1‑150 (MQ=255)  tCAAGGAGCGCCCTCCATGATTCAGCCTGTCGCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTtggcctggc                                                                                                                                                      <  1:1619948/150‑1 (MQ=255)      ggAGCGCCCTCCATGATTCAGCCTGTCGCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAAc                                                                                                                                                  <  2:649363/150‑1 (MQ=255)        aGCGGCCTCCATGATTCAGCCTGTCGCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑a                                                                                                                                               <  2:1182880/150‑1 (MQ=255)         gcgcCCTCCATGATTCAGCCTGTCACTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACCGCGTGGACCAGAAATACAAAATAGGAGGTGATCAGGTCGTTGGCCAGGCCAACG‑aa                                                                                                                                              <  2:1207091/150‑2 (MQ=255)          cgcCCTCCATGATTCAGCCTGTCGCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACt                                                                                                                                             <  1:2143287/150‑1 (MQ=255)            cccTCCATGATTCAGCCTGTCGCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑Ac                                                                                                                                          >  2:2094885/1‑150 (MQ=255)            cccTCCATGATTCAGCCTGTCGCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑Ac                                                                                                                                          <  1:1317044/150‑1 (MQ=255)             ccTCCATGATTCAGCCTGTCGCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑Acc                                                                                                                                         <  1:1443497/150‑1 (MQ=255)               tCCATGATTCAGCCTGTCGCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACcgc                                                                                                                                       >  1:2015103/1‑150 (MQ=255)                    gATTCAGCCTGTCGCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCAt                                                                                                                                  >  1:361056/1‑150 (MQ=255)                       tCAGCCTGTCGCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGc                                                                                                                               <  2:1045418/150‑1 (MQ=255)                         aGCCTGTCGCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCt                                                                                                                             >  1:2195085/1‑150 (MQ=255)                           ccTGTCGCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAGTCTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTAc                                                                                                                           <  1:1996174/150‑1 (MQ=255)                           ccTGTCGCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTAc                                                                                                                           >  2:2084107/1‑150 (MQ=255)                           ccTGTCGCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTAc                                                                                                                           <  2:1821666/150‑1 (MQ=255)                             tGTCGCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGGGTTGGCCTGAACAAGGAAATGTGGGGCGGTCAGTGCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGa                                                                                                                         >  1:1345773/1‑150 (MQ=255)                                cGCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACAt                                                                                                                      <  1:1573228/150‑1 (MQ=255)                                    gatctcacaCCTGCTGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCTc                                                                                                                  <  2:46933/145‑1 (MQ=255)                                    gAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCTc                                                                                                                  <  1:418229/150‑1 (MQ=255)                                      aCGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCTCgg                                                                                                                <  1:330246/150‑1 (MQ=255)                                       cGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCTCGGc                                                                                                               <  1:1435545/150‑1 (MQ=255)                                        gCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑Ac                                                                                                                                          <  1:1711242/122‑1 (MQ=255)                                        gCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑Ac                                                                                                                                          >  2:1711242/1‑122 (MQ=255)                                          acacCTGCCGGCACCAGCTATCTGGATGTCGGCCCGGGCTGCTCGGTGCNGTTTCTCCACGGCGGGGGCCTGAACCAGGAAATGTGGGGCGGTCAGTTCGTTTGCCTGGCCAACG‑ACT‑ACCGCTGCATCGCCTACGACATGCTCgagatc                                                                                                            >  1:46933/1‑145 (MQ=255)                                              cTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCATGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCTCGGCCACGGTc                                                                                                        >  1:942959/1‑150 (MQ=255)                                              cTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGTGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCACGTCATCGCCTACGACATCCCCGGCCACGGTc                                                                                                        >  2:22783/1‑150 (MQ=255)                                              cTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑a                                                                                                                                               <  1:1991031/112‑1 (MQ=255)                                              cTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑a                                                                                                                                               >  2:1991031/1‑112 (MQ=255)                                              cTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCTCGGCCACGGTc                                                                                                        >  2:389144/1‑150 (MQ=255)                                               tGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCa                                                                                                       >  2:1607766/1‑150 (MQ=255)                                                   ggCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGGACAAGGAAATGTGGGGCAGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCTCGGGCACGGTCAAAGc                                                                                                   >  2:1326526/1‑150 (MQ=255)                                                    gCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGcc                                                                                                  >  2:633264/1‑150 (MQ=255)                                                      aCCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCcgc                                                                                                <  2:2195085/150‑1 (MQ=255)                                                       ccAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCaa                                                                                                                                                   <  1:660410/100‑1 (MQ=255)                                                       ccAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCaa                                                                                                                                                   >  2:660410/1‑100 (MQ=255)                                                       ccAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCcgcg                                                                                               >  1:1141657/1‑150 (MQ=255)                                                           cTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTAcc                                                                                           <  1:2084107/150‑1 (MQ=255)                                                                 ggATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGgccgc                                                                                     >  1:411587/1‑150 (MQ=255)                                                                  gATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGgccgcc                                                                                    <  1:546661/150‑1 (MQ=255)                                                                    tGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGa                                                                                  <  2:694844/150‑1 (MQ=255)                                                                          ccAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGGTGGCCTGGCCAACG‑ACT‑ACCGCGGCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGGCGACCGGCCGCCGACACAcc                                                                            >  1:1306051/1‑150 (MQ=255)                                                                           cAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCAGGCCAACG‑aa                                                                                                                                              <  1:237204/84‑2 (MQ=255)                                                                           cAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCAGGCCAACG‑aa                                                                                                                                              >  2:237204/1‑83 (MQ=255)                                                                            aGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCGTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCCGCTGCCGGCACACCGc                                                                          >  2:548204/1‑150 (MQ=255)                                                                             gggCCAACCCGTGGTGCTGATCCACGGCGTGTGCCTGCCCAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCt                                                                         <  1:22783/150‑1 (MQ=255)                                                                             gggCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCt                                                                         >  1:515334/1‑150 (MQ=255)                                                                               gCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGcc                                                                                                                              <  2:2096887/95‑1 (MQ=255)                                                                               gCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGcc                                                                                                                              >  1:2096887/1‑95 (MQ=255)                                                                                ccAACCCGTGGTGCTGATCGACGCGGTGGGCCTGAACAAGGAAATGGGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGa                                                                      >  1:655186/1‑150 (MQ=255)                                                                                ccAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACACGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGa                                                                      >  2:882548/1‑150 (MQ=255)                                                                                 cAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAAc                                                                                                                                                  >  2:1103851/1‑75 (MQ=255)                                                                                 cAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAAc                                                                                                                                                  <  1:1103851/75‑1 (MQ=255)                                                                                 cAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGaa                                                                     <  1:155498/150‑1 (MQ=255)                                                                                        gggtgCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATg                                                              <  1:2066692/149‑1 (MQ=255)                                                                                           tgCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCaaa                                                                                                     <  1:2145192/108‑1 (MQ=255)                                                                                           tgCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCaaa                                                                                                     >  2:2145192/1‑108 (MQ=255)                                                                                             cTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGAc                                                         <  1:1625705/150‑1 (MQ=255)                                                                                                   cACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGa                                                          <  1:2151645/143‑1 (MQ=255)                                                                                                   cACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGa                                                          >  2:2151645/1‑143 (MQ=255)                                                                                                     cGGCGTGGGCCGGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGc                                                 <  2:428600/150‑1 (MQ=255)                                                                                                           gggCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACt                                           <  1:2094885/150‑1 (MQ=255)                                                                                                             gCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCcgc                                                                                                <  2:1652475/95‑1 (MQ=255)                                                                                                             gCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCcgc                                                                                                >  1:1652475/1‑95 (MQ=255)                                                                                                             gCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCTCGGCCACGGTAAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGc                                         <  1:1705000/150‑1 (MQ=255)                                                                                                               cTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCtt                                       >  1:2092024/1‑150 (MQ=255)                                                                                                               cTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCtt                                       >  2:1075222/1‑150 (MQ=255)                                                                                                               cTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCtt                                       >  2:116910/1‑150 (MQ=255)                                                                                                                 gAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTgc                                     >  2:1714166/1‑149 (MQ=255)                                                                                                                 gAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGa                                     <  2:2015103/150‑1 (MQ=255)                                                                                                                  aaCAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTg                                                                        <  1:1394761/114‑1 (MQ=255)                                                                                                                  aaCAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTg                                                                        >  2:1394761/1‑114 (MQ=255)                                                                                                                    cAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGacca                                  >  2:1446705/1‑150 (MQ=255)                                                                                                                     aaGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGaccac                                 >  2:1567395/1‑150 (MQ=255)                                                                                                                        gAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGACCACCTg                              <  2:942959/150‑1 (MQ=255)                                                                                                                              tgGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGACCACCTGCAAAtt                        <  2:1475242/150‑1 (MQ=255)                                                                                                                                gggCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGACCACCTGCAAATTGc                      >  2:1308474/1‑150 (MQ=255)                                                                                                                                gggCCGTCAGTTCGGTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCGACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGACCACCTGCAAATTGc                      >  2:2311701/1‑150 (MQ=255)                                                                                                                                    ggTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGACCACCTGCAAATTGCACaa                  <  1:1567395/150‑1 (MQ=255)                                                                                                                                    ggTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGACCACCTGCAAATTGCACaa                  <  2:2211523/150‑1 (MQ=255)                                                                                                                                    gcggtagtCGTTGGCCAGGCCAACGAACTGACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGACCACCTGCAAATTGcac                    >  1:1656272/8‑150 (MQ=255)                                                                                                                                      tCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGACCACCTGCAAATTGCACAAGc                >  2:132644/1‑150 (MQ=255)                                                                                                                                           tCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGACCACCTGCAAATTGCATAAGCCACCg           <  1:334595/150‑1 (MQ=255)                                                                                                                                           tCGTTGGCCAGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGACCACCTGCAAATTGCACAAGCCACCg           >  2:859300/1‑150 (MQ=255)                                                                                                                                             gTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGACCACCTGCAAATTGCACAAGCCACCGTg         >  2:1610239/1‑150 (MQ=255)                                                                                                                                             gTTGGCCAGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGACCACCTGCAAATTGCACAAGCCACCGTg         >  2:59972/1‑150 (MQ=255)                                                                                                                                                ggcctggccAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGACCACCTGCAAATTGCACAAGCCACCGTGATc      <  1:389144/150‑1 (MQ=255)                                                                                                                                                ggcctggccAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGACCACCTGCAAATTGCACAAGCCAACGTGATc      <  2:2241836/150‑1 (MQ=255)                                                                                                                                                 gcctggccAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGACCACCTGCAAATTGCACAAGCCACCGTGATCg     >  1:1687490/1‑150 (MQ=255)                                                                                                                                                    tggccAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGACCACCTGCAAATTGCACAAGCCACCGTGATCGGCt  >  1:317660/1‑150 (MQ=255)                                                                                                                                                    |                                                                                                                                                        TTCAAGGAGCGCCCTCCATGATTCAGCCTGTCGCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGACCACCTGCAAATTGCACAAGCCACCGTGATCGGCT  >  NC_002947/4141065‑4141361

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