Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,605,418	A→C	5.9%	Q263P (CAG→CCG)	PP_4077 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,605,418	0	A	C	5.9%	202.3 / 4.8	69	Q263P (CAG→CCG)	PP_4077	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (33/31);  new base C (2/2);  total (35/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.67e-01

AGCCCTGAGGAGGCCTTGCTGCCCTACCCGCAGAACGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCG  >  NC_002947/4605273‑4605567                                                                                                                                                  |                                                                                                                                                      aGCCCTGAGGAGGCCTTGCTGCCCTACCCGCAGAACGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTg                                                                                                                                                   >  1:755362/1‑150 (MQ=255)  gCCCTGAGGAGGCCTTGCTGCCCTACCCGCAGAACGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTgg                                                                                                                                                  <  2:1715811/150‑1 (MQ=255)    ccTGAGGAGGCCTTGCTGCCCTACCCGCAGAACGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGc                                                                                                                                                >  1:499147/1‑150 (MQ=255)         ggaggCCTTGCTGCCCTACCCGCAGAACGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGa                                                                                                                                           <  2:1359463/150‑1 (MQ=255)          gaggCCTTGCTGCCCTACCCGCAGAACGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAg                                                                                                                                          <  1:662946/150‑1 (MQ=255)                      cccTACCCGCAGAACGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCAcg                                                                                                                              >  1:640029/1‑150 (MQ=255)                       ccTACCCGCAGAACGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGt                                                                                                                                         >  2:531393/1‑138 (MQ=255)                       ccTACCCGCAGAACGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGt                                                                                                                                         <  1:531393/138‑1 (MQ=255)                         tACCCGCAGAACGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGt                                                                                                                           >  1:969258/1‑150 (MQ=255)                            ccGCAGAACGTGTTAGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCaa                                                                                                                        >  2:573584/1‑150 (MQ=255)                              gCAGAACGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGCGTTCCATTTCACGCGTCAACt                                                                                                                      >  2:1499624/1‑150 (MQ=255)                              gCAGAACGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTc                                                                                                                                 >  1:1824320/1‑139 (MQ=255)                              gCAGAACGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTc                                                                                                                                 <  2:1824320/139‑1 (MQ=255)                              gCAGAACGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACt                                                                                                                      >  1:1553592/1‑150 (MQ=255)                                aGAAGGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTGGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGc                                                                                                                    <  2:1119280/150‑1 (MQ=255)                                    cGTGTTTGATTGCTACCGTATTCTGCAGGAATACTTCGTGTTCACGAAGCCCTTCCACTTCTTCATCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGTGCATCGAGGTCCATTTCACGCGTCAACTGCCAGa                                                                                                                <  2:457165/150‑1 (MQ=255)                                                 tACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCAAGCGTAAAATGCAAGACACCCtccccggt                                                                                                   >  2:1566119/1‑143 (MQ=255)                                                   ccGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTgg                                                                                                 <  1:1566119/150‑1 (MQ=255)                                                           tGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCAcg                                                                                                                              <  2:582266/113‑1 (MQ=255)                                                           tGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCAcg                                                                                                                              >  1:582266/1‑113 (MQ=255)                                                            gCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGc                                                                                                                    <  1:181055/122‑1 (MQ=255)                                                            gCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGc                                                                                                                    >  2:181055/1‑122 (MQ=255)                                                            gCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACagaaga                                                                                        <  2:755362/150‑1 (MQ=255)                                                             cAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGAt                                                                                       >  1:720591/1‑150 (MQ=255)                                                              aGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGAtt                                                                                      >  1:1284676/1‑150 (MQ=255)                                                                gAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAAc                                                                                                                       <  1:267300/115‑1 (MQ=255)                                                                gAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAAc                                                                                                                       >  2:267300/1‑115 (MQ=255)                                                                gAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTc                                                                                    <  1:1709715/150‑1 (MQ=255)                                                                gAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTc                                                                                    >  2:745483/1‑150 (MQ=255)                                                                   tACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGc                                                                                 <  1:1764295/150‑1 (MQ=255)                                                                   tACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGc                                                                                 <  2:479297/150‑1 (MQ=255)                                                                     cTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCt                                                                               >  1:552009/1‑150 (MQ=255)                                                                            ttCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGaca                                                                                                              <  1:330403/112‑1 (MQ=255)                                                                            ttCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGaca                                                                                                              >  2:330403/1‑112 (MQ=255)                                                                               ccGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAg                                                                                  >  2:1733228/1‑137 (MQ=255)                                                                               ccGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAg                                                                                  <  1:1733228/137‑1 (MQ=255)                                                                                  aaGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCg                                                                  <  2:552009/150‑1 (MQ=255)                                                                                     cgcTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCTGACCCCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCt                                                               <  1:573584/150‑1 (MQ=255)                                                                                        ttCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTc                                                            >  1:1531317/1‑150 (MQ=255)                                                                                        ttCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTc                                                            >  1:1445268/1‑150 (MQ=255)                                                                                        ttCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTCCCCACAGGTGGGAAAAGAGTCTTAGTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTc                                                            <  2:1911389/150‑1 (MQ=255)                                                                                           cACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAAc                                                         <  2:1445268/150‑1 (MQ=255)                                                                                           cACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAAc                                                         <  2:1531317/150‑1 (MQ=255)                                                                                             cttcttCAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAg                                                                                  <  2:300868/123‑1 (MQ=255)                                                                                             cttcttCAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAg                                                                                  >  1:300868/1‑123 (MQ=255)                                                                                              ttcttcAGCATGACCGGCCTGGGAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTg                                                      >  2:1575959/1‑150 (MQ=255)                                                                                                 ttcAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCATGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTc                                                   >  1:444259/1‑150 (MQ=255)                                                                                                 ttcAGAATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTc                                                   >  1:662705/1‑150 (MQ=255)                                                                                                  tcAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCa                                                  <  2:720591/150‑1 (MQ=255)                                                                                                   cAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCaa                                                 <  2:444259/150‑1 (MQ=255)                                                                                                      cATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCa                                                                                   >  1:1856511/1‑113 (MQ=255)                                                                                                      cATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCa                                                                                   <  2:1856511/113‑1 (MQ=255)                                                                                                         tACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAg                                           <  2:969258/150‑1 (MQ=255)                                                                                                           ccGGACTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTg                                         >  1:1014854/1‑150 (MQ=255)                                                                                                            cGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCACACGGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGc                                        <  1:1070680/150‑1 (MQ=255)                                                                                                                 tGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGTGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGc                                   >  2:401106/1‑150 (MQ=255)                                                                                                                   gAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGAACCCCGGCTGTTAACATGTTAAAGCACAGTGCCGAGccc                                 >  1:949968/1‑150 (MQ=255)                                                                                                                         cTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGAc                           >  1:1129234/1‑150 (MQ=255)                                                                                                                         cTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGAc                           <  1:1806861/150‑1 (MQ=255)                                                                                                                          tgtgGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGAcc                          >  2:1244580/1‑150 (MQ=255)                                                                                                                            tgGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTg                        <  2:1553592/150‑1 (MQ=255)                                                                                                                              gCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTggg                      >  2:44714/1‑149 (MQ=255)                                                                                                                                cGGCACAAGCCTGCCCACAGGTGTGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCa                    <  1:401106/150‑1 (MQ=255)                                                                                                                                    acaAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGa                <  2:1885271/150‑1 (MQ=255)                                                                                                                                         gatGCCCACCTGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAg               <  1:664179/144‑1 (MQ=255)                                                                                                                                         gatGCCCACCTGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAg               >  2:664179/3‑146 (MQ=255)                                                                                                                                          atGCCCACCTGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGt                                                     <  1:1546846/106‑1 (MQ=255)                                                                                                                                          atGCCCACCTGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGt                                                     >  2:1546846/2‑107 (MQ=255)                                                                                                                                              ccACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAg      >  2:314965/1‑150 (MQ=255)                                                                                                                                               cacaGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGt                                                             <  2:400011/94‑1 (MQ=255)                                                                                                                                               cacaGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGt                                                             >  1:400011/1‑94 (MQ=255)                                                                                                                                                  aGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCg  >  1:1823478/1‑150 (MQ=255)                                                                                                                                                  aGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCg  >  2:33680/1‑150 (MQ=255)                                                                                                                                                  |                                                                                                                                                      AGCCCTGAGGAGGCCTTGCTGCCCTACCCGCAGAACGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCG  >  NC_002947/4605273‑4605567

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