Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	972,269	C→G	6.0%	R187R (CGC→CGG)	secD →	protein translocase subunit SecD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	972,269	0	C	G	6.0%	215.4 / 2.6	66	R187R (CGC→CGG)	secD	protein translocase subunit SecD
Reads supporting (aligned to +/- strand):  ref base C (34/28);  new base G (2/2);  total (36/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.96e-01

GGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGACCACCACCGATCGCAACGAGCTTGCGGTGCTGCGTCTGGCGCTGACCCAGGCGAAAGTCGCCGAGATCCGCGAATACTCGATCAAGCAGAACCTCAC  >  NC_002947/972125‑972415                                                                                                                                                 |                                                                                                                                                   ggacatggacaAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGaa                                                                                                                                                 >  8:24702/1‑148 (MQ=255)     atggacaAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGc                                                                                                                                            <  2:204659/149‑1 (MQ=255)        gccttgtccATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCAc                                                                                                                                          <  7:238225/141‑1 (MQ=255)             ggCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGccc                                                                                                                                    <  7:24702/149‑1 (MQ=255)               ccATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCt                                                                                                                                   >  5:193483/1‑148 (MQ=255)               ccATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCATGCCTCAGCAGGATGGCGGCATCATGCTGGGCTACGCTGACGATGCCACCAGCGAACAGGCACGAGccat                                                                                                                                   >  7:361459/1‑146 (MQ=255)                            cATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAAtt                                                                                                                     <  1:151406/149‑1 (MQ=255)                                   gTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGa                                                                                                              >  5:244735/1‑149 (MQ=255)                                    tCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGAt                                                                                                             >  3:268633/1‑149 (MQ=255)                                    tCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGAt                                                                                                             <  6:295771/149‑1 (MQ=255)                                     cTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCCCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGAtt                                                                                                            >  1:154452/1‑149 (MQ=255)                                        cGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTc                                                                                                          >  5:192594/1‑148 (MQ=255)                                        cGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTc                                                                                                          >  5:110045/1‑148 (MQ=255)                                         gAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGGGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGa                                                                                                        >  5:73841/1‑149 (MQ=255)                                         gAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGGGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGa                                                                                                        >  1:24513/1‑149 (MQ=255)                                             gCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCt                                                                                                     >  4:103120/1‑148 (MQ=255)                                              cGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGa                                                                                                   <  6:192594/149‑1 (MQ=255)                                              cGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGa                                                                                                   <  6:193483/149‑1 (MQ=255)                                                     aaGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGaccacccc                                                                                            >  4:81558/1‑147 (MQ=255)                                                     aaGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGaccaccac                                                                                            >  5:180732/1‑149 (MQ=255)                                                         ccTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGACCACCACCGAt                                                                                        <  6:110045/149‑1 (MQ=255)                                                          cTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGt                                                                                                                                        <  1:293106/100‑1 (MQ=255)                                                          cTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGt                                                                                                                                        >  2:293106/1‑100 (MQ=255)                                                          cTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGACCACCACCGATc                                                                                       >  4:95465/1‑149 (MQ=255)                                                            tgctgcGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAAt                                                                                                                >  3:86218/1‑122 (MQ=255)                                                            tgctgcGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAAt                                                                                                                <  4:86218/122‑1 (MQ=255)                                                                 cgcAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGACCACCACCGATCGCAAc                                                                                  >  3:168861/1‑147 (MQ=255)                                                                    aaaGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGCCGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGACCACCACCGATc                                                                                       >  3:110509/1‑139 (MQ=255)                                                                    aaaGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGACCACCACCGATc                                                                                       <  4:110509/139‑1 (MQ=255)                                                                     aaGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGACCACCACCGATCGCAACGAGCtt                                                                            >  6:166733/1‑149 (MQ=255)                                                                       gagCTCGTCCGCTACCTCAGCCTGCCTCAGAAGGATGGCGTCATCATGCTGGTATTCGCTGACGATGCCACCCGCTAACAGGCACGTGCCCTGATCCTCAAGAATTTCAATGATTTCGACCTGACCACCACCGATCGCAACGAGCTTGc                                                                          <  5:353130/149‑1 (MQ=255)                                                                        agCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGACCACCACCGATCGCAACGAGCTTGCg                                                                         <  8:116054/149‑1 (MQ=255)                                                                          cgcgTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGACCACCACCGATCGCAACGAGCTTGCGGt                                                                       >  1:168489/1‑149 (MQ=255)                                                                             gTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGACCACCACCGATCGCAACGAGCTTGCGGtgct                                                                    <  5:221956/149‑1 (MQ=255)                                                                                 gCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGACCACCACCGATCGCAACGAGCTTGCGGTGCTGCGt                                                                >  7:25902/1‑149 (MQ=255)                                                                                 gCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGACCACCACCGATCGCAACGAGCTTGCGGTGCTGCGt                                                                >  6:274923/1‑149 (MQ=255)                                                                                 gCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGACCACCACCGATCGCAACGAGCTTGCGGTGCTGCGt                                                                <  8:8669/149‑1 (MQ=255)                                                                                   tACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGACCACCACCGATCGCAACGAGCTTg                                                                           <  1:337938/136‑1 (MQ=255)                                                                                   tACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGACCACCACCGATCGCAACGAGCTTg                                                                           >  2:337938/1‑136 (MQ=255)                                                                                     ccGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGACCACCACCGATCGCAACGAGCTTGCGGTGCTGCGTCt                                                              <  3:103120/147‑1 (MQ=255)                                                                                     ccGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGACCACCACCGATCGCAACGAGCTTGCGGTGCTGCGTCTg                                                             <  8:25902/148‑1 (MQ=255)                                                                                     ccGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGACCACCACCGATCGCAACGAGCTTGCGGAGCTGCGTCTg                                                             <  5:166733/148‑1 (MQ=255)                                                                                       gCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGACCACCACCGATCGCAACGAGCTTGCGGTGCTGCGTCTGgcg                                                          <  4:168861/149‑1 (MQ=255)                                                                                           cgaagcccagcatGATGCCGCCATCCTGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGACCACCACCGATCGCAACGAGCTTGCGGTGCTGCGTCTGGCGCTGa                                                      >  5:259379/8‑149 (MQ=255)                                                                                               cctCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGGGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGACCACCACCGATCGCAACGAGCTTGGGGTGCTGCGTCTGGCGCTGAcc                                                    <  2:24513/147‑1 (MQ=255)                                                                                               cctCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGGGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGACCACCACCGATCGCAACGAGCTTGGGGTGCTGCGTCTGGCGCTGAcc                                                    <  6:73841/147‑1 (MQ=255)                                                                                                   agcagGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGACCACCACCGATCGCAACGAGCTTGCGGTGCTGCGTCTGGCGCTGACCCAgg                                                <  5:274923/147‑1 (MQ=255)                                                                                                   agcagGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGACACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGACCACCACCGATCGCAACGAGCTTGCGGTGCTGCGTCTGGCGCTGACCCAGGc                                               >  1:204989/1‑148 (MQ=255)                                                                                                    gcagGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCAGATCCGCAAGAATTTCAATGATTTCGACCTGACCACCACCGATCGCAACGAGCTTGCGGTGCTGCGTCTGGCGCTGACCCAGGCGa                                             <  3:10102/149‑1 (MQ=255)                                                                                                        gATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGACCACCACCGATCGCAACGAGCTTGCGGTGATGCGTCTGGCGCTGCCCCCGGCGAAAGt                                         >  4:230663/1‑149 (MQ=255)                                                                                                              ggcATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGACCACCACCGATCGCAACGAGCTTGCGGTGCTGCGTCTGGCGCTGACCCAGGCGAAAGTCGcc                                     >  3:267169/1‑147 (MQ=255)                                                                                                                catcatGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGACCACCACCGATCGCAACGAGCTTGCGGTGCTGCGTCTGGCGCTGACCCAGGCGAAAGTCGcc                                     >  2:35923/1‑145 (MQ=255)                                                                                                                     tGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGACCACCACCGATCGCAACGAGCTTGCGGTGCTGCGTCTGGCGCTGACCCAGGCGAAAGTCGCCGAGATCcgc                            <  6:244735/149‑1 (MQ=255)                                                                                                                     tGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGACCACCACCGATCGCAACGAGCTTGCGGTGCTGCGTCTGGCGCTGACCCAGGCGAAAGTCGCCGAGATCcgc                            <  3:230663/149‑1 (MQ=255)                                                                                                                        tGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGACCACCACCGATCGCAACGAGCTTGCGGTGCTGCGTCTGGCGCTGACCCAGGCGAAAGTCGCCGAGATCCGCGaa                         >  4:167427/1‑149 (MQ=255)                                                                                                                           gCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGACCACCACCGATCGCAACGAGCTTGCGGTGCTGCGTCTGGCGCTGACCCAGGCGAAAGTCGCCGAGATCCGCGAATAc                      >  6:164189/1‑149 (MQ=255)                                                                                                                             ctCGCTGACGCTGCCACCCGCGACCCGGCACGTGCCCTGATCCCCAAGAATTTCAATGAATTCGACCTGACCACCACCGATCGCAACGAGCTTGCGGTGCTGCGTCTGTCGCTGACCCAGCCGAAAGTCGCCGAGATCCGCGAATACTc                    <  8:361459/148‑1 (MQ=255)                                                                                                                              tCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGACCACCACCGATCGCAACGAGCTTGCGGTGCTGCGTCTGGCGCTGACCCAGGCGAAAGTCGCCGAGATCCGCGAATACTc                    >  5:83842/1‑148 (MQ=255)                                                                                                                                   gACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGACCACCACCGATCGCAACGAGCTTGCGGTGCTGCGTCTGGCGCTGACCCAGGCGAAAGTCGCCGAGATCCGCGAATACTCGATc                >  6:174653/1‑147 (MQ=255)                                                                                                                                      gATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGACCACCACCGATCGCAACGAGCTTGCGGTGCTGCGTCTGGCGCTGACCCAGGCGAAAGTCGCCGAGATCCGCGAATACTCGa                  >  7:202436/1‑142 (MQ=255)                                                                                                                                      gATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGACCACCACCGATCGCAACGAGCTTGCGGTGCTGCGTCTGGCGCTGACCCAGGCGAAAGTCGCCGAGATCCGCGAATACTCGa                  <  8:202436/142‑1 (MQ=255)                                                                                                                                           cACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGACCACCACCGATCGCAACGAGCTTGCGGTGCTGCGTCTGGCGCTGACCCAGGCGAAAGTCGCCGAGATCCGCGAATACTCGATCaa              >  3:330851/1‑141 (MQ=255)                                                                                                                                           cACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGACCACCACCGATCGCAACGAGCTTGCGGTGCTGCGTCTGGCGCTGACCCAGGCGAAAGTCGCCGAGATCCGCGAATACTCGATCaa              <  4:330851/141‑1 (MQ=255)                                                                                                                                            aCCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGACCACCACCGATCGCAACGAGCTTGCGGTGCTGCGTCTGGCGCTGACCCAGGCGAAAGTCGCCGAGATCCGCGAATACTCGATCAAGCAGAAc       <  6:82150/147‑1 (MQ=255)                                                                                                                                            aCCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGACCACCACCGATCGCAACGAGCTTGCGGTGCTGCGTCTGGCGCTGACCCAGGCGAAAGTCGCCGAGATCCGCGAATACTCGATCAAGCAGAAc       >  5:82150/1‑147 (MQ=255)                                                                                                                                             cccGCGAACAGGCACGTGCCCTGTTCCGCAAGAATTTCAATGATTTCGACCTGACCACCACCGATCGCAACGAGCTTGCGGTGCTGCGTCTGGCGCTGACCAAGGCGAAAGTCGCCGAGATCCGCGAATACTCGATCAAGCAGAACCt     >  2:338851/1‑148 (MQ=255)                                                                                                                                               cgcgAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGACCACCACCGATCGCAACGAGCTTGCGGTGCTGCGTCTGGCGCTGACCCAGGCGAAAGTCGCCGAGATCCGCGAATACTCGATCAAGCAGAACCTcac  <  2:168489/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                   GGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGACCACCACCGATCGCAACGAGCTTGCGGTGCTGCGTCTGGCGCTGACCCAGGCGAAAGTCGCCGAGATCCGCGAATACTCGATCAAGCAGAACCTCAC  >  NC_002947/972125‑972415

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