Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,254,470	G→C	6.7%	A423G (GCC→GGC)	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,470	0	G	C	6.7%	180.7 / 3.9	60	A423G (GCC→GGC)	PP_2852	sulfatase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base G (27/29);  new base C (2/2);  total (29/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.85e-01

CAATGTAAGGTTGAAGCCCCAGGACGAGGCCAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAA‑GTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTGAGCAGCGTCAGCATCCATGGCTTGTCCTGCTTCTGCAACTGGCCCACATAGCCCAGGGCGCCTTCGAA  >  NC_002947/3254322‑3254618                                                                                                                                                      |                                                                                                                                                     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                                                                                                                                              <  6:303130/145‑1 (MQ=255)             gAAGCCCCAGGACGAGGCCAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAA‑GTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCt                                                                                                                                              >  5:303130/1‑145 (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                                                                                                                                       <  5:331295/142‑3 (MQ=255)                       gACGAGGCCAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAA‑GTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGCAGCCGTCGCCTACCTgga                                                                                                                                       <  5:336876/142‑3 (MQ=255)                       gACGAGGCCAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAA‑GTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGCAGCCGTCGCCTACCTgga                                                                                                                                       >  6:331295/1‑140 (MQ=255)                       gACGAGGCCAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAA‑GTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGCAGCCGTCGCCTACCTgga                                                                                                                                       >  6:336876/1‑140 (MQ=255)                          gAGGCCAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAA‑GTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTc                                                                                                                             <  4:286238/149‑1 (MQ=255)                          gAGGCCAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAA‑GTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTc                                                                                                                             <  1:344420/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                                                                                                    >  3:312004/1‑149 (MQ=255)                                                   cGTGAGACTCGTCGGAA‑GTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTcggccggc                                                                                                    <  7:77938/149‑1 (MQ=255)                                                          cTCGTCGGAA‑GTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAAt                                                                                             <  4:293009/149‑1 (MQ=255)                                                           tcgtcgGAA‑GTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGaa                                                                                              >  6:281178/1‑147 (MQ=255)                                                                        accaccaccAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTggtgggt                                                                                >  8:346012/1‑149 (MQ=255)                                                                        accaccaccAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTggtgggt                                                                                <  4:118325/149‑1 (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                                            >  1:317728/1‑149 (MQ=255)                                                                                                            tCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTGAGCAGCGTCAGCATCCATGGCTTGTc                                            >  7:242020/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                                                                  >  1:1648/1‑109 (MQ=255)                                                                                                                              gCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTGagca                                                                  <  2:1648/109‑1 (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)                                                                                                                                                     ggCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTGAGCAGCGTCAGCATCCATGGCTTGTCCTGCTTCTGCAACTGGCCCACATAGCCCAGGGCGCCTTCg    >  6:364514/1‑148 (MQ=255)                                                                                                                                                      gCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTGAGCAGCGTCAGCATCCATGGCTTGTCCTGCTTCTGCAACTGGCCCACATAGCCCAGGGCGCCTTCgaa  <  8:70327/149‑1 (MQ=255)                                                                                                                                                      gCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTGAGCAGCGTCAGCATCCATGGCTTGTCCTGCTTCTGCAACTGGCCCACATAGCCCAGGGCGCCTTCgaa  <  2:208357/149‑1 (MQ=255)                                                                                                                                                      gCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTGAGCAGCGTCAGCATCCATGGCTTGTCCTGCTTCTGCAACTGGCCCACATAGCCCAGGGCGCCTTCgaa  <  4:351223/149‑1 (MQ=255)                                                                                                                                                      gCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTGAGCAGCGTCAGCATCCATGGCTTGTCCTGCTTCTGCAACTGGCCCACATAGCCCAGGGCGCCTTCg    >  3:351223/1‑147 (MQ=255)                                                                                                                                                      |                                                                                                                                                     CAATGTAAGGTTGAAGCCCCAGGACGAGGCCAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAA‑GTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTGAGCAGCGTCAGCATCCATGGCTTGTCCTGCTTCTGCAACTGGCCCACATAGCCCAGGGCGCCTTCGAA  >  NC_002947/3254322‑3254618

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