Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	821,679	G→C	6.5%	M350I (ATG→ATC)	PP_0707 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	821,679	0	G	C	6.5%	188.6 / 3.1	61	M350I (ATG→ATC)	PP_0707	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (33/24);  new base C (2/2);  total (35/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-01

GTGTTCAGTGCCTTCATCACAGGCCTGGGCCGCGCCATGCTGATGCTGCAGCGTCCGTCCTGGCGGTTGCCGGCGATCCACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGTGGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCACCCTGCGCAAGCATGACCTGGA  >  NC_002947/821536‑821825                                                                                                                                                |                                                                                                                                                   gtgtTCAGTGCCTTCATCACAGGCCTGGGCCGCGCCATGCTGATGCTGCAGCGTCCGTCCTGGCGGTTGCCGGCGATCCACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCAtgatga                                                                                                                                                <  2:190015/148‑1 (MQ=255)  tgtTGAGTGCCTTCATCACAGGCCTGGGCCGCGCCTTGATGATGCTGCAGCGTCCGTCCTGGCGGTTGACGGCTATCCACGACTAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGAcc                                                                                                                                              <  2:60124/149‑1 (MQ=255)       aGTGCCTTCATCACAGGCCTGGGCCGCGCCATGCTGATGCTGCAGCGTCCGTCCTGGCGGTTGCCGGCGATCCACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATg                                                                                                                                           >  2:213325/1‑147 (MQ=255)           ccTTCATCACAGGCCTGGGCCGCGCCATGCTGATGCTGCAGCGTCCGTCCTGGCGGTTGCCGGCGATCCACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAgcg                                                                                                                                      >  6:100989/1‑148 (MQ=255)           ccTTCATCACAGGCCTGGGCCGCGCCATGCTGATGCTGCAGCGTCCGTCCTGGCGGTTGCCGGCGATCCACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAgcg                                                                                                                                      >  7:4234/1‑148 (MQ=255)            cTTCATCACAGGCCTGGGCCGCGCCATGCTGATGCTGCAGCGTCCGTCCTGGCGGTTGCCGGCGATCCACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCa                                                                                                                                    >  1:18021/1‑149 (MQ=255)                     aGGCCTGGGCCGCGCCATGCTGATGCTGCAGCGTCCGTCCTGGCGGTTGCCGGCGATCCACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCg                                                                                                                           <  7:258872/149‑1 (MQ=255)                         cTGGGCCGCGCCATGCTGATGCTGCAGCGTCCGTCCTGGCGGTTGCCGGCGATCCACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATg                                                                                                                                           >  7:178172/1‑129 (MQ=255)                         cTGGGCCGCGCCATGCTGATGCTGCAGCGTCCGTCCTGGCGGTTGCCGGCGATCCACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATg                                                                                                                                           <  8:178172/129‑1 (MQ=255)                         cTGGGCCGCGCCATGCTGATGCTGCAGCGTCCGTCCTGGCGGTTGCCGGCGATCCACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGAt                                                                                                                       >  2:297486/1‑149 (MQ=255)                                gcgcCATGCTGATGCTGCAGCGTCCGTCCTGGCGGTTGCCGGCGATCCACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGc                                                                                                                 >  1:297066/1‑148 (MQ=255)                                         tgatgctgCAGCGTCCGTCCTGGCGGTTGCCGGCGATCCACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGa                                                                                                                        <  2:186998/132‑1 (MQ=255)                                         tgatgctgCAGCGTCCGTCCTGGCGGTTGCCGGCGATCCACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGa                                                                                                                        >  1:186998/1‑132 (MQ=255)                                           atgctgCAGCGTCCGTCCTGGCGGTTGCCGGCGATCCACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCt                                                                                                     >  2:92263/1‑149 (MQ=255)                                                 gaGCGTCGGTGCTGGCGGTTGCCGGCGATCCACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACgg                                                                                                <  2:297066/147‑1 (MQ=255)                                                  aGCGTCCGTCCTGGCGGTTGCCGGCGATCCACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTg                                                                                              <  1:92263/149‑1 (MQ=255)                                                   gCGTCCGTCCTGGCGGTTGCCGGCGATCCACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTgg                                                                                             >  1:92627/1‑149 (MQ=255)                                                     gtccgtccTGGCGGTTGCCGGCGATCCACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCAt                                                                                                                                   >  7:86706/1‑109 (MQ=255)                                                     gtccgtccTGGCGGTTGCCGGCGATCCACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCAt                                                                                                                                   <  8:86706/109‑1 (MQ=255)                                                        cgtccTGGCGGTTGCCGGCGATCCACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACgg                                                                                                >  8:274793/1‑141 (MQ=255)                                                        cgtccTGGCGGTTGCCGGCGATCCACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACgg                                                                                                <  7:274793/141‑1 (MQ=255)                                                            cTGGCGGTTGCCGGCGATCCACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGgc                                                                                                        >  6:77098/1‑129 (MQ=255)                                                            cTGGCGGTTGCCGGCGATCCACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGgc                                                                                                        <  5:77098/129‑1 (MQ=255)                                                               gCGGTTGCCGGCGATCCACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGTTGACGGCCACCGTCAGCgcca                                                                                 <  7:27401/149‑4 (MQ=255)                                                                  gtggCCGGCGATCCACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGa                                                                              <  4:82228/146‑1 (MQ=255)                                                                           gATCCACGACGGAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACggg                                                                                  <  2:80592/136‑1 (MQ=255)                                                                           gATCCACGACGGAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACggg                                                                                  >  1:80592/1‑136 (MQ=255)                                                                           gATCCACGACGGAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACggg                                                                                  <  6:37942/136‑1 (MQ=255)                                                                           gATCCACGACGGAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACggg                                                                                  >  5:37942/1‑136 (MQ=255)                                                                                acgacgAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCtggtggtg                                                                >  2:197315/1‑149 (MQ=255)                                                                                        gTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGTGGCGCTCAc                                                        <  1:297486/149‑1 (MQ=255)                                                                                        gTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGTGGCGCTCAc                                                        >  3:159394/1‑149 (MQ=255)                                                                                            ccAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGTGGCGCTCACCttt                                                    <  1:197315/149‑1 (MQ=255)                                                                                             cAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGTGGCGCTCACCTTTg                                                   >  8:44707/1‑149 (MQ=255)                                                                                              aGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGTGGCGCTCACCTTTGc                                                  <  5:100989/149‑1 (MQ=255)                                                                                                    ctgggctggTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGTGGCGCTCACCTTTGCAGgcgc                                            >  5:21475/1‑149 (MQ=255)                                                                                                            gTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGTGGCGCTCACCTTTGCAGGCGCGCTGCtgt                                    >  8:60386/1‑149 (MQ=255)                                                                                                            gTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGTGGCGCTCACCTTTGCAGGCGCGCTGCtgt                                    <  8:4234/149‑1 (MQ=255)                                                                                                            gTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGTGGCGCTCACCTTTGCAGGCGCGCTGCtgt                                    >  7:245445/1‑149 (MQ=255)                                                                                                                ccGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGTGGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTa                                >  6:104861/1‑149 (MQ=255)                                                                                                                 cGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGTGGCGCTCACCTTTGCAGGCGCGctgctg                                     <  2:74307/143‑1 (MQ=255)                                                                                                                 cGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGTGGCGCTCACCTTTGCAGGCGCGctgctg                                     >  1:74307/1‑143 (MQ=255)                                                                                                                  gAAGGTGCTGGCATTGGCACTGATGGTGATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCtggtggtgg                                                               >  1:221015/1‑116 (MQ=255)                                                                                                                  gAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCtggtggtgg                                                               <  2:221015/116‑1 (MQ=255)                                                                                                                       tGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGTGGCGCt                                                           >  6:281119/1‑115 (MQ=255)                                                                                                                       tGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGTGGCGCt                                                           <  5:281119/115‑1 (MQ=255)                                                                                                                        gCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGTGGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCa                        >  3:202730/1‑149 (MQ=255)                                                                                                                         cTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGTGGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCAc                       <  4:202730/149‑1 (MQ=255)                                                                                                                             cATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGTGGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCACCCTg                   >  4:4631/1‑149 (MQ=255)                                                                                                                                  gCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGTGGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCACCCTgcgc                >  3:261750/1‑147 (MQ=255)                                                                                                                                  gCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACAGCGCTGGTGGTGGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCACCCTgcgc                >  4:131854/1‑147 (MQ=255)                                                                                                                                    aCTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGTGGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCACCCTGCGCAAg             >  2:193719/1‑148 (MQ=255)                                                                                                                                    aCTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGTGGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCACCCTGCGCAAGc            <  1:213325/149‑1 (MQ=255)                                                                                                                                     tccaTGGTCATCATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGTGGCGCTCACCTTTGCAGgc                                              <  2:49002/111‑1 (MQ=255)                                                                                                                                     tccaTGGTCATCATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGTGGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCACCCTGCGCAAGCa           >  2:268365/4‑149 (MQ=255)                                                                                                                                     tccaTGGTCATCATGACCATGGAGCGCATGAACAGCATGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGTGGCGCTCACCTTTGCAGgc                                              >  1:49002/4‑114 (MQ=255)                                                                                                                                        aTGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGTGGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCTCCCTGCGcatgcatg         <  2:18021/148‑1 (MQ=255)                                                                                                                                         tGGTCATCATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGTGGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCACCCTGCGCAAGCATGAc       <  1:268365/149‑1 (MQ=255)                                                                                                                                            tCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGACTGACCGCGCTGGTGGTGGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCACCCTGCGCAAGCATGACCTg    >  4:47031/1‑149 (MQ=255)                                                                                                                                             cATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGTGGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACcgccgc                         >  6:173397/1‑127 (MQ=255)                                                                                                                                             cATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGTGGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACcgccgc                         <  5:173397/127‑1 (MQ=255)                                                                                                                                             cATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGTGGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCACCCTGCGCAAGCATGACCTg    >  7:248424/1‑148 (MQ=255)                                                                                                                                              atgatgACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGTGGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCACCCTGCGCAAGCATGACCTgg   <  4:261750/148‑1 (MQ=255)                                                                                                                                              atgatgACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGTGGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCACCCTGCGCAAGCATGACCTGGa  <  7:60386/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                   GTGTTCAGTGCCTTCATCACAGGCCTGGGCCGCGCCATGCTGATGCTGCAGCGTCCGTCCTGGCGGTTGCCGGCGATCCACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGTGGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCACCCTGCGCAAGCATGACCTGGA  >  NC_002947/821536‑821825

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