Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,240,646	A→G	6.3%	N70D (AAT→GAT)	PP_4618 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,240,646	0	A	G	6.3%	187.3 / 2.3	63	N70D (AAT→GAT)	PP_4618	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (33/26);  new base G (2/2);  total (35/28)
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.48e-01

GCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGA‑TGATGGACAAGGTACTGCAGCGCCGCCGCGGCCAACCACTGACCCTGGCCATCCTCGCTTTGGAACTGGCCAGGCGGCTGTCCATACCGCTGG  >  NC_002947/5240499‑5240793                                                                                                                                                    |                                                                                                                                                     gCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGaa                                                                                                                                                     <  3:242411/149‑1 (MQ=255)          tgGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCt                                                                                                                                               >  5:25511/1‑146 (MQ=255)          tgGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCt                                                                                                                                               <  6:25511/146‑1 (MQ=255)          tgGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTggg                                                                                                                                            >  3:243169/1‑149 (MQ=255)          tgGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTggg                                                                                                                                            <  6:82067/149‑1 (MQ=255)              ttGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCttt                                                                                                                                        <  6:77834/149‑1 (MQ=255)                 ctgcCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTcag                                                                                                                                     >  8:29068/1‑149 (MQ=255)                  tgcCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTcag                                                                                                                                     >  8:250710/1‑148 (MQ=255)                    ccGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTcagcag                                                                                                                                  <  7:250710/149‑1 (MQ=255)                      gAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGa                                                                                                                                >  1:319454/1‑149 (MQ=255)                           cGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAat                                                                                                                           <  7:245594/149‑1 (MQ=255)                              cccTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAatat                                                                                                                         >  6:25256/1‑148 (MQ=255)                               ccTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCa                                                                                                                       >  3:139397/1‑149 (MQ=255)                                 tGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCAcc                                                                                                                     >  4:226334/1‑149 (MQ=255)                                        gAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCAc                                                                                                                      >  4:39066/1‑141 (MQ=255)                                        gAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCAc                                                                                                                      <  3:39066/141‑1 (MQ=255)                                           cccGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGccc                                                                                                           <  2:319454/149‑1 (MQ=255)                                           cccGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGccc                                                                                                           >  2:155126/1‑149 (MQ=255)                                                 gCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGAAAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATgc                                                                                                     >  4:35152/1‑149 (MQ=255)                                                  ccAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATgcg                                                                                                    >  4:87850/1‑149 (MQ=255)                                                       cTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCga‑t                                                                                              >  7:47196/1‑149 (MQ=255)                                                         tgctgcGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTgc                                                                                                               >  6:249451/1‑131 (MQ=255)                                                         tgctgcGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTgc                                                                                                               <  5:249451/131‑1 (MQ=255)                                                           ctgcgctgcACGACCTGCCGCGCGAGACACGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCga‑tga                                                                                            <  4:243169/147‑1 (MQ=255)                                                           ctgcgctgcACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCga‑tgat                                                                                           <  1:155126/148‑1 (MQ=255)                                                           ctgcgctgcACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGa‑tgatg                                                                                          <  2:2623/149‑1 (MQ=255)                                                                ctgcACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGA‑TGATGGAc                                                                                       >  4:80351/1‑147 (MQ=255)                                                                    aCGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGA‑TGATGGACAAGGTa                                                                                 <  4:139397/149‑1 (MQ=255)                                                                    aCGACCAGCAGCCCGAGATAAGCGCCAACATGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGATTTCAGCAGGATGAATATCACCCCTTGCGCCACCAATCGTAGA‑TGATGGACAAGGt                                                                                  >  6:196742/1‑148 (MQ=255)                                                                        ccTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGA‑TGATGGACAAGGTACTGc                                                                             <  3:80351/149‑1 (MQ=255)                                                                          tGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGA‑TGATGGACAAGGTa                                                                                 >  3:269478/1‑143 (MQ=255)                                                                          tGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGA‑TGATGGACAAGGTa                                                                                 <  4:269478/143‑1 (MQ=255)                                                                           gCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGA‑TGATGGACAAGGTACTGCAgc                                                                          >  6:254116/1‑149 (MQ=255)                                                                            cAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGA‑TGATGGACAAGGTACTGCAgc                                                                          <  7:29068/148‑1 (MQ=255)                                                                            cAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGA‑TGATGGACAAGGTACTGCAgc                                                                          <  8:47196/148‑1 (MQ=255)                                                                                  gagaTAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGA‑TGATGGACAAGGTACTGCAGcgccgccg                                                                   >  4:92420/1‑149 (MQ=255)                                                                                     aTAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGA‑TGATGGACAAGGTACTGCAGCGCCGCCGCgg                                                                >  5:234735/1‑149 (MQ=255)                                                                                      tAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGA‑TGATGGACAAGGTACTGCAGCGCCGCCGCgg                                                                <  3:87850/148‑1 (MQ=255)                                                                                             cAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGA‑TGATGGACAAGGTACTGCAGCGCCGCCGCGGCCAACCAc                                                        >  7:36170/1‑149 (MQ=255)                                                                                                    ccGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGA‑TGATGGACAAGGTAc                                                                                <  4:154180/118‑1 (MQ=255)                                                                                                    ccGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGA‑TGATGGACAAGGTAc                                                                                >  3:154180/1‑118 (MQ=255)                                                                                                    ccGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGA‑TGATGGACAAGGTACTGCAGCGCCGCCGCGGCCAACCACTGAcc                                                   <  5:25256/147‑1 (MQ=255)                                                                                                      gATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGA‑TGATGGACAAGGTACTGCAGCGCCGCCGCGGCCAACCACTGACCCTg                                                >  4:30268/1‑148 (MQ=255)                                                                                                      gATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGA‑TGATGGACAAGGTACTGCAGCGCCGCCGCGGCCAACCACTGACCCTg                                                <  3:30268/148‑1 (MQ=255)                                                                                                       aTGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGA‑TGATGGACAAGGTACTGCAGCGCCGCCGCGGCCAACCACTGACCCTGGc                                              >  1:77068/1‑149 (MQ=255)                                                                                                           ttCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGA‑TGATGGACAAGGTACTGCAGCGCCGCCGCGGCCAACCACTGACCCTGGCCATc                                          >  2:110757/1‑149 (MQ=255)                                                                                                                cTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGA‑TGATGGACAAGGTACTGCAGCGCCGCCGCGGCCAACCACTGACCCTGGCCATCCTCGc                                     <  5:254116/149‑1 (MQ=255)                                                                                                                   tGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGA‑TGATGGACAAGGTACTGCAGCGCCGCCGCGGCCAACCACTGACCCTGGCCATCCTCGCtt                                   <  3:226334/148‑1 (MQ=255)                                                                                                                     cGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGA‑TGATGGACAAGGTACTGCAGCGCCGCCGCGGCCAACCACTGACCCTGGCCATCCTCGCTTTgg                                <  2:77068/149‑1 (MQ=255)                                                                                                                            gCCCAGCCACTGCTGCGCCAGCTGGATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGA‑TGATGGACAAGGTACTGCAGCGCCGCCGCGGCCAACCACTGACCCTGGCCATCCTCGCTTTGGAACTGGc                         >  5:114199/1‑149 (MQ=255)                                                                                                                            gCCCAGCCACTGCTGCGCCAGCTGGATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGA‑TGATGGACAAGGTACTGCAGCGCCGCCGCGGCCAACCACTGACCCTGGCCATCCTCGCTTTGGAACTGGc                         >  5:119233/1‑149 (MQ=255)                                                                                                                              ccAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGA‑TGATGGACAAGGTACTGCAGCGCCGCCGCGGCCAACCACTGACCCTGGCCATCCTCGCTTTGGAACTGGc                         >  3:19524/1‑147 (MQ=255)                                                                                                                              ccAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGA‑TGATGGACAAGGTACTGCAGCGCCGCCGCGGCCAACCACTGACCCTGGCCATCCTCGCTTTGGAACTGGc                         >  1:231717/1‑147 (MQ=255)                                                                                                                                 gCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGA‑TGATGGACAAGGTACTGCAGCGCCGCCGCGGCCAACCACTGACCCTGGCCATCCTCGCTTTGGAACTGGCCAggc                    >  6:298027/1‑149 (MQ=255)                                                                                                                                 gCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGA‑TGATGGACAAGGTACTGCAGCGCCGCCGCGGCCAACCACTGACCCTGGCCATCCTCGCTTTGGAACTGGCCAggc                    >  2:225697/1‑149 (MQ=255)                                                                                                                                 gCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGA‑TGATGGACAAGGTACTGCAGCGCCGCCGCGGCCAACCACTGACCCTGGCCATCCTCGCTTTGGAACTGGCCAggc                    <  5:298027/149‑1 (MQ=255)                                                                                                                                    aCTGCTGCGCCAGCTGGATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGA‑TGATGGACAAGGTACTGCAGCGCCGCCGCGGCCAACCACTGACCCTGGCCATCCTCGCTTTGGAACTGGCCAggcggc                 <  6:119233/149‑1 (MQ=255)                                                                                                                                    aCTGCTGCGCCAGCTGGATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGA‑TGATGGACAAGGTACTGCAGCGCCGCCGCGGCCAACCACTGACCCTGGCCATCCTCGCTTTGGAACTGGCCAggcggc                 <  6:114199/149‑1 (MQ=255)                                                                                                                                           cgcCAGCTTATTGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGA‑TGATGTACAAGGTACTGCAGCGCCGCCGCGGCCAACCACTGACCCTGGCCATCCTCGCTTTGGAACTGGCCAGGCGGCTGTc             <  5:196742/146‑1 (MQ=255)                                                                                                                                           cgcCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGA‑TGATGGACAAGGTACTGCAGCGCCGCCGCGGCCCACCACTGCCCCTGGCCATCCTCGCTTTGGAACTGTCCAGGCGGCTGTCCa           >  7:79943/1‑148 (MQ=255)                                                                                                                                           cgcCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGA‑TGATGGACAAGGTACTGCAGCGCCGCCGCGGCCAACCACTGACCCTGGCCATCCTCGCTTTGGAACTGGCCAGGCGGCTGTcc            <  2:231717/147‑1 (MQ=255)                                                                                                                                             ccAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGA‑TGATGGACAAGGTACTGCAGCGCCGCCGCGGCCAACCACTGACCCTGGCCATCCTCGCTTTGGAACTGGCCAGGCGGCTGTCCATAc        >  6:310437/1‑149 (MQ=255)                                                                                                                                              cAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGA‑TGATGGACAAGGTACTGCAGCGCCGCCGCGGCCAACCACTGACCCTGGCCATCCTCGCTTTGGAACTGGCCAGGCGGCTGTCCATAcc       >  3:23449/1‑149 (MQ=255)                                                                                                                                                   gAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGA‑TGATGGACAAGGTACTGCAGCGCCGCCGCGGCCAACCACTGACCCTGGCCATCCTCGCTTTGGAACTGGCCAGGCGGCTGTCCATACCGCTgg  >  3:11144/1‑149 (MQ=255)                                                                                                                                                    aaTGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGAGTGATGGACAAGGTACTGCAGCGCCGCCGCGGCCAACCACTGACCCTGGCCATCCTCGCTTTGGAACTGGCCAGGCGGCTGTCCATACCGCTg   >  1:88888/1‑148 (MQ=255)                                                                                                                                                    |                                                                                                                                                     GCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGA‑TGATGGACAAGGTACTGCAGCGCCGCCGCGGCCAACCACTGACCCTGGCCATCCTCGCTTTGGAACTGGCCAGGCGGCTGTCCATACCGCTGG  >  NC_002947/5240499‑5240793

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