Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,906,508	A→T	8.2%	L112Q (CTG→CAG)	PP_4315 ←	epimerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,906,508	0	A	T	8.2%	212.1 / 9.1	73	L112Q (CTG→CAG)	PP_4315	epimerase
Reads supporting (aligned to +/- strand):  ref base A (32/35);  new base T (3/3);  total (35/38)
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.89e-01

CAACGTTTCGAGCAGCATGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCACGGTGGGGGT  >  NC_002947/4906360‑4906656                                                                                                                                                     |                                                                                                                                                     cAACGTTTCGAGCAGCATGGCGTCATCGATAACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCCGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCAGGGCcag                                                                                                                                                     >  2:2310030/1‑150 (MQ=255)  aaCGTTTCGAGCAGCATGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCcagc                                                                                                                                                    <  1:21695/150‑1 (MQ=255)  aaCGTTTCGAGCAGCATGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGACCGCTGGCATGTCGACCACGGCGGGAGTTTTGTGGGTAAAGTCCATGGCcaga                                                                                                                                                    >  2:21695/1‑149 (MQ=255)      tttCGAGCAGCATGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGc                                                                                                                                                <  2:223326/150‑1 (MQ=255)      tttCGAGCAGCATGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGc                                                                                                                                                <  2:389824/150‑1 (MQ=255)          gAGCAGCATGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATc                                                                                                                                            >  2:184046/1‑150 (MQ=255)            gcagcaTGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCTGCTGGCCAtgga                                                                                                                                          <  1:979773/150‑4 (MQ=255)                   ggCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGTGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAgg                                                                                                                                                 <  2:1721795/136‑1 (MQ=255)                   ggCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGTGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAgg                                                                                                                                                 >  1:1721795/1‑136 (MQ=255)                   ggCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACg                                                                                                                                   <  1:200985/150‑1 (MQ=255)                    gCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACgg                                                                                                                                  >  2:1464306/1‑150 (MQ=255)                    gCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACgg                                                                                                                                  <  1:357914/150‑1 (MQ=255)                     cGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGc                                                                                                                                 <  1:1245457/150‑1 (MQ=255)                               caccacGACATAATCGTCGCTAAGGTACACCGCCTTTGTCTTGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCTATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCaa                                                                                                                       <  2:1354741/150‑1 (MQ=255)                                 ccacGACATAATCGTCGCTAAGGTAGACCGCATTAGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACt                                                                                                                     >  1:1332961/1‑150 (MQ=255)                                  cacGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTc                                                                                                                    >  1:1640473/1‑150 (MQ=255)                                  cacGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTc                                                                                                                    >  2:2025346/1‑150 (MQ=255)                                  cacGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCACTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCATGCCATCAGCGCCACGGCTGACGCGCAACTc                                                                                                                    <  2:486380/150‑1 (MQ=255)                                          aaTCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGc                                                                                                            >  2:367832/1‑150 (MQ=255)                                            tcgtcgCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCgg                                                                                                          <  1:367832/150‑1 (MQ=255)                                            tcgtcgCTAAGGTACAACGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCgg                                                                                                          >  1:580949/1‑150 (MQ=255)                                            tcgtccCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGAAAGATTCCGACCACGACCGCTACGGGTTGTTTGGCGTGGACGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCgg                                                                                                          <  1:1112485/150‑1 (MQ=255)                                                  cTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTg                                                                                                    >  1:244330/1‑150 (MQ=255)                                                   tAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTga                                                                                                   >  2:992178/1‑150 (MQ=255)                                                    aaGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGATGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGTCAGCAGGTCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTgaa                                                                                                  <  2:1166097/150‑1 (MQ=255)                                                       gTACACCGCCTTCGCGTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTgaagcg                                                                                               >  1:1678182/1‑150 (MQ=255)                                                         acacCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTgaagcgaa                                                                                             <  1:2275024/150‑1 (MQ=255)                                                         acacCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGCGTATTgaagcgaa                                                                                             >  2:1930457/1‑150 (MQ=255)                                                          cacCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTgaagcgaag                                                                                            <  2:244330/150‑1 (MQ=255)                                                          aacCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTgaagcgaag                                                                                            <  2:587825/149‑1 (MQ=255)                                                           acCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTgaagcgaagc                                                                                           >  1:1428215/1‑150 (MQ=255)                                                             cGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGa                                                                                                                              >  2:1677741/1‑113 (MQ=255)                                                             cGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGa                                                                                                                              <  1:1677741/113‑1 (MQ=255)                                                                 ggcGCCTGGCTCACGCGCCGTGCGCGCAAAAGCCCCTCTGGCATGTCGACCGCTACAGGTTGTTTGGCGGGGAAGTCCATCGCCAGCAGGCCATCAGCGCCACGGCTGACTCGCAACTCGCCACCGCGGGTATTGAAGCGAAGCACTTgc                                                                                     <  2:713240/148‑1 (MQ=255)                                                                   cGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTgcgc                                                                                   <  2:1640473/150‑1 (MQ=255)                                                                    gCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGcc                                                                                  >  1:103392/1‑150 (MQ=255)                                                                      cTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTg                                                                                >  2:1829591/1‑150 (MQ=255)                                                                       tGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGc                                                                               >  2:1555026/1‑150 (MQ=255)                                                                           tCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGGTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCAc                                                                           <  1:2180198/150‑1 (MQ=255)                                                                                cccAGCGCCTGCAAAAGGACCGTTGCGATGACACCCTCAACGCTGGGGGTGGGGTGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGt                                                                      <  2:1423808/150‑1 (MQ=255)                                                                                 ccAGCGCCTGCAAAAGGCCCGCTGTCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGtt                                                                     <  2:106166/150‑1 (MQ=255)                                                                                 ccAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCt                                                                              <  2:2066059/141‑1 (MQ=255)                                                                                 ccAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCt                                                                              >  1:2066059/1‑141 (MQ=255)                                                                                     cgcCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTc                                                                 <  1:1555026/150‑1 (MQ=255)                                                                                        cTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGaa                                                              <  2:1332961/150‑1 (MQ=255)                                                                                        cTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTAGCTCGaa                                                              <  1:1464306/150‑1 (MQ=255)                                                                                           cAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAg                                                           >  1:874116/1‑150 (MQ=255)                                                                                               aGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACg                                                       >  1:724173/1‑150 (MQ=255)                                                                                                   ccGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAg                                                   <  1:1160752/150‑1 (MQ=255)                                                                                                       tGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCCCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGa                                               <  1:2306852/150‑1 (MQ=255)                                                                                                         gCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTAtt                                                                                                     <  2:884350/94‑1 (MQ=255)                                                                                                         gCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTAtt                                                                                                     >  1:884350/1‑94 (MQ=255)                                                                                                           aTGTCGACCGCTACGGGTTGTTTGGCGGTGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTc                                                                 >  2:549476/1‑128 (MQ=255)                                                                                                           aTGTCGACCGCTACGGGTTGTTTGGCGGTGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTc                                                                 <  1:549476/128‑1 (MQ=255)                                                                                                             gTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCaa                                         >  1:1029454/1‑150 (MQ=255)                                                                                                                     cTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCAt                                 <  1:1829591/150‑1 (MQ=255)                                                                                                                               tttGGCGGTGAATTTCATGGTTAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGACACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGTTCACCCAGTTGCTCGAACAGCACGTAAGCGGAGGCCAAGATAGCATGGCCGCACAg                       >  2:923951/1‑150 (MQ=255)                                                                                                                                          aGTCCATGGCAAGCAGGCCATCAGCGCCACGGCTGCGGCGCAACTCACCACAACGGGAATTGCAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCAAACAGCACGTGAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCa            >  2:1853751/1‑150 (MQ=255)                                                                                                                                           gctgaTGGCCTGCTGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAg                                              >  1:1052986/5‑115 (MQ=255)                                                                                                                                           gctgaTGGCCTGCTGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAg                                              <  2:1052986/111‑1 (MQ=255)                                                                                                                                           gTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCAc           >  2:2299207/1‑150 (MQ=255)                                                                                                                                            tCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCACg          >  2:353989/1‑150 (MQ=255)                                                                                                                                            tCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCACg          >  1:853545/1‑150 (MQ=255)                                                                                                                                            tCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCACg          >  1:671040/1‑150 (MQ=255)                                                                                                                                            ctgaTGGCCTGCTGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACtt                                                                                       <  2:903953/70‑1 (MQ=255)                                                                                                                                            ctgaTGGCCTGCTGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACtt                                                                                       >  1:903953/4‑73 (MQ=255)                                                                                                                                             tgaTGGCCTGCTGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCACgg         >  2:1510639/3‑150 (MQ=255)                                                                                                                                              cATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAg                                              >  1:1672677/1‑112 (MQ=255)                                                                                                                                              cATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAg                                              <  2:1672677/112‑1 (MQ=255)                                                                                                                                              cATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCACGGt        >  1:1758390/1‑150 (MQ=255)                                                                                                                                              cATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCACGGt        <  1:184046/150‑1 (MQ=255)                                                                                                                                                  gCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCACGGTgggg    <  1:2310030/150‑1 (MQ=255)                                                                                                                                                  gCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCACGGTgggg    <  2:874116/150‑1 (MQ=255)                                                                                                                                                  gCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTACGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCACGGTgggg    <  2:846003/150‑1 (MQ=255)                                                                                                                                                  cccaccaGGCCATCAGCCCCACGGCTGACGCGCACCTCCCCACTGCCCGTATTGAACCGAAGCATTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTACTCCGAGGCCAAGGTAGCATGGCCCCACAGGTCGACTTCCACGGTgggg    <  1:923951/149‑1 (MQ=255)                                                                                                                                                   ccAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCACGGTg       >  2:848639/1‑146 (MQ=255)                                                                                                                                                   ccAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCACGGTg       <  1:848639/146‑1 (MQ=255)                                                                                                                                                    cagcagGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCACGGTGGGGGt  >  2:939319/1‑150 (MQ=255)                                                                                                                                                     agcagGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAg                                    <  2:516983/115‑1 (MQ=255)                                                                                                                                                     agcagGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAg                                    >  1:516983/1‑115 (MQ=255)                                                                                                                                                     |                                                                                                                                                     CAACGTTTCGAGCAGCATGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCACGGTGGGGGT  >  NC_002947/4906360‑4906656

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