Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,254,467	T→C	6.5%	Y424C (TAC→TGC)	PP_2852 ←	sulfatase domain‑containing protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,254,467	0	T	C	6.5%	178.8 / 2.4	61	Y424C (TAC→TGC)	PP_2852	sulfatase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (28/29);  new base C (2/2);  total (30/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.96e-01

CGCCAATGTAAGGTTGAAGCCCCAGGACGAGGCCAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAA‑GTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTGAGCAGCGTCAGCATCCATGGCTTGTCCTGCTTCTGCAACTGGCCCACATAGCCCAGGGCGCCT  >  NC_002947/3254319‑3254613                                                                                                                                                      |                                                                                                                                                   cgcCAATGTAAGGTTGAAGCCCCAGGACGAGGCCAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAA‑GTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGt                                                                                                                                                    <  3:166326/149‑1 (MQ=255)    cAATGTAAGGTTGAAGCCCCAGGACGAGGCCAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAA‑GTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAgg                                                                                                                                                 >  5:224275/1‑149 (MQ=255)     aaTGTAAGGTTGAAGCCCCAGGACGAGGCCAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAA‑GTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGc                                                                                                                                                >  4:48078/1‑149 (MQ=255)     aaTGTAAGGTTGAAGCCCCAGGACGAGGCCAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAA‑GTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGc                                                                                                                                                >  2:344420/1‑149 (MQ=255)       tGTAAGGTTGAAGCCCCAGGACGAGGCCAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAA‑GTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGa                                                                                                                                              >  1:183648/1‑149 (MQ=255)                gAAGCCCCAGGACGAGGCCAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAA‑GTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCt                                                                                                                                         >  5:303130/1‑145 (MQ=255)                gAAGCCCCAGGACGAGGCCAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAA‑GTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCt                                                                                                                                         <  6:303130/145‑1 (MQ=255)                gAAGCCCCAGGACGAGGCCAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAA‑GTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCt                                                                                                                                     >  2:27784/1‑149 (MQ=255)                 aaGCCCCAGTACGAGGCCAGGGGCACGCCATCGATGCCGTGAGACTCGTCGGAA‑GTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTTCGTCATCCAGGTAGGCGACGGCTGCCt                                                                                                                                     <  6:311590/148‑1 (MQ=255)                   gCCCCAGGACGAGGCCAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAA‑GTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGtt                                                                                                                                  >  3:293009/1‑149 (MQ=255)                        aGGACGAGGCCAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAA‑GTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGgt                                                                                                                             <  6:224275/149‑1 (MQ=255)                         ggACGAGGCCAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAA‑GTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGgtg                                                                                                                            <  6:26388/149‑1 (MQ=255)                          gACGAGGCCAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAA‑GTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGCAGCCGTCGCCTACCTgga                                                                                                                                  >  6:331295/1‑140 (MQ=255)                          gACGAGGCCAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAA‑GTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGCAGCCGTCGCCTACCTgga                                                                                                                                  <  5:336876/142‑3 (MQ=255)                          gACGAGGCCAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAA‑GTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGCAGCCGTCGCCTACCTgga                                                                                                                                  >  6:336876/1‑140 (MQ=255)                          gACGAGGCCAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAA‑GTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGCAGCCGTCGCCTACCTgga                                                                                                                                  <  5:331295/142‑3 (MQ=255)                             gAGGCCAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAA‑GTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTc                                                                                                                        <  1:344420/149‑1 (MQ=255)                             gAGGCCAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAA‑GTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTc                                                                                                                        <  4:286238/149‑1 (MQ=255)                             gAGGCCAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAA‑GTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCACGTAGGCGACGGCTGCCTGTTTCGGTGTGTc                                                                                                                        <  1:29819/149‑1 (MQ=255)                                gCCAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAA‑GTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCg                                                                                                                       >  8:102851/1‑147 (MQ=255)                                  cAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAA‑GTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAgc                                                                                                                   >  6:154612/1‑149 (MQ=255)                                  cAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAA‑GTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAgc                                                                                                                   >  3:118325/1‑149 (MQ=255)                                         aCGCCATCGATGCCGTGAGACTCGTCGGAA‑GTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAg                                                                                                            <  3:76025/149‑1 (MQ=255)                                          cGCCATCGATGCCGTGAGACTCGTCGGAA‑GTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAgg                                                                                                           <  7:54585/149‑1 (MQ=255)                                            ccATCGATGCCGTGAGACTCGTCGGAAGGTCACCACCACCTGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAgg                                                                                                           >  1:22570/1‑148 (MQ=255)                                                 gATGCCGTGAGACTCGTCGGAA‑GTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTcgg                                                                                                    >  1:361201/1‑149 (MQ=255)                                                      cGTGAGACTCGTCGGAA‑GTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTcggccggc                                                                                               <  7:77938/149‑1 (MQ=255)                                                      cGTGAGACTCGTCGGAA‑GTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTcggccggc                                                                                               >  3:312004/1‑149 (MQ=255)                                                             cTCGTCGGAA‑GTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAAt                                                                                        <  4:293009/149‑1 (MQ=255)                                                              tcgtcgGAA‑GTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGaa                                                                                         >  6:281178/1‑147 (MQ=255)                                                                           accaccaccAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTggtgggt                                                                           <  4:118325/149‑1 (MQ=255)                                                                           accaccaccAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTggtgggt                                                                           >  8:346012/1‑149 (MQ=255)                                                                           accaccacaAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGACCGGCGCCGAATAGGGCTggtgggt                                                                           >  2:80890/1‑149 (MQ=255)                                                                             caccaccaGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGc                                                                         <  3:48078/149‑1 (MQ=255)                                                                                caccaGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCa                                                                      >  7:317035/1‑149 (MQ=255)                                                                                 accaGTGTGTCGTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCAc                                                                     <  2:361201/149‑1 (MQ=255)                                                                                   caGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGGCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACgg                                                                   <  1:27784/149‑1 (MQ=255)                                                                                     gtgtgtCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTg                                                                 <  1:313573/149‑1 (MQ=255)                                                                                               gAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGACCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTGAGCAGCGt                                                         <  1:80890/147‑1 (MQ=255)                                                                                                 ggACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTGAGCAGCGTCAg                                                      >  3:314429/1‑148 (MQ=255)                                                                                                            cgcTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTGAGCAGCGTCAGCATCCATGGCtt                                          <  7:346012/149‑1 (MQ=255)                                                                                                               tCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTGAGCAGCGTCAGCATCCATGGCTTGTc                                       >  7:242020/1‑149 (MQ=255)                                                                                                               tCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTGAGCAGCGTCAGCATCCATGGCTTGTc                                       >  1:317728/1‑149 (MQ=255)                                                                                                                 gagaTTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTGAGCAGCGTCAGCATCCAt                                               >  2:328264/1‑139 (MQ=255)                                                                                                                 gagaTTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTGAGCAGCGTCAGCATCCAt                                               <  1:328264/139‑1 (MQ=255)                                                                                                                 gagaTTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTGAGCAGCGTCAGCATCCATGGCTTGTcc                                      <  2:317728/148‑1 (MQ=255)                                                                                                                 gagaTTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTGAGCAGCGTCAGCATCCATGGCTTGTCCt                                     >  1:208357/1‑149 (MQ=255)                                                                                                                 gagaTTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGGCGAAGCGCTCGAGGTACTCGGCCGTCGCCGAATAGGGCTGGTGGGTGCCCACGGTGAGCAGCGTCAGCATCCATGGCTTGTc                                       <  2:22570/147‑1 (MQ=255)                                                                                                                       gTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTGAGCAGCGTCAGCATCCATGGCTTGTCCTGCTTc                                <  4:314429/148‑1 (MQ=255)                                                                                                                        tCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTGAGCAGCGTCAGCATCCATGGCTTGTCCTGCTTCt                               <  2:183648/148‑1 (MQ=255)                                                                                                                            gAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTGAGCAGCGTCAGCATCCATGGCTTGTCCTGCTTCTGCAAc                          >  3:279005/1‑149 (MQ=255)                                                                                                                               aTGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTGAGCAGCGTCAGCATCCATGGCTTGTCCTGCTTCTGCAACTgg                       <  5:154612/149‑1 (MQ=255)                                                                                                                                tGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTGAGCAGCGTCAGCATCCATGGCTTGTCCTGCTTCTGCAACTGGc                      <  5:281178/149‑1 (MQ=255)                                                                                                                                tGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTGAGCAGCGTCAGCATCCATGGCTTGTCCTGCTTCTGCAACTGGc                      <  8:317035/149‑1 (MQ=255)                                                                                                                                 gCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTGagca                                                             <  2:1648/109‑1 (MQ=255)                                                                                                                                 gCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTGagca                                                             >  1:1648/1‑109 (MQ=255)                                                                                                                                   cccGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTGAGCAGCGTCAGCATCCATGGCTTGTCCTGCTTCTGCAACTGGCCca                   >  8:157411/1‑149 (MQ=255)                                                                                                                                    ccGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTGAGCAGCGTCAGCATCCATGGCTTGTCCTGCTTCTGCAACTGGCCca                   >  6:375096/1‑148 (MQ=255)                                                                                                                                    ccGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTGAGCAGCGTCAGCATCCATGGCTTGTCCTGCTTCTGCAACTGGCCca                   <  5:375096/148‑1 (MQ=255)                                                                                                                                                  cAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTGAGCAGCGTCAGCATCCATGGCTTGTCCTGCTTCTGCAACTGGCCCACATAGCCCAGGgcgc    >  7:70327/1‑149 (MQ=255)                                                                                                                                                    ggTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTGAGCAGCGTCAGCATCCATGGCTTGTCCTGCTTCTGCAACTGGCCCACATAGCCCAGGGCGCCt  <  7:157411/149‑1 (MQ=255)                                                                                                                                                     gTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTGAGCAGCGTCAGCATCCATGGCTTGTCCTGCTTCTGCAACTGGCCCACATAGACCAGGGCGCCt  >  7:202810/1‑148 (MQ=255)                                                                                                                                                      |                                                                                                                                                   CGCCAATGTAAGGTTGAAGCCCCAGGACGAGGCCAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAA‑GTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTGAGCAGCGTCAGCATCCATGGCTTGTCCTGCTTCTGCAACTGGCCCACATAGCCCAGGGCGCCT  >  NC_002947/3254319‑3254613

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