Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,659,684	T→A	5.9%	L243Q (CTG→CAG)	metK →	methionine adenosyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,659,684	0	T	A	5.9%	204.7 / 3.1	68	L243Q (CTG→CAG)	metK	methionine adenosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (37/27);  new base A (2/2);  total (39/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.79e-01

GTTTCGCAGAAAGACCTGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCA‑TGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCGCTACGTGGCCAAGAACATCGTCGCCGCCGGT  >  NC_002947/5659539‑5659829                                                                                                                                                  |                                                                                                                                                   gTTTCGCAGAAAGACCTGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGACGGCCGAACTACTGCACAAAGGCACCCAGTACCACAACAACCCGACCGGCAACTTCATCATCGGTGGACCGGTGGGTGACTGCGGCCTgcc                                                                                                                                                 >  3:280662/1‑147 (MQ=255)    tCGCAGAAAGACCTGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTgg                                                                                                                                              >  8:249897/1‑149 (MQ=255)     cGCAGAAAGACCTGCAAGAAGCCGTCATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTgg                                                                                                                                              <  3:392295/148‑1 (MQ=255)      gCAGAAAGACCTGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGcc                                                                                                                                            >  6:319877/1‑149 (MQ=255)        aGAAAGACCTGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCg                                                                                                                                           >  4:179692/1‑148 (MQ=255)        aGAAAGACCTGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCg                                                                                                                                           >  5:97440/1‑148 (MQ=255)        aGAAAGACCTGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCAGTCAGGCCg                                                                                                                                           >  4:79851/1‑148 (MQ=255)          aaaGACCTGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCAGTCAGGCCGCag                                                                                                                                        <  3:154548/149‑2 (MQ=255)             gACCTGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCAGTCAGGCCGCagtca                                                                                                                                     <  3:79851/149‑5 (MQ=255)                cTGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGatcat                                                                                                                                  >  2:137127/1‑149 (MQ=255)                   cAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCAtcgt                                                                                                                               <  1:137127/149‑1 (MQ=255)                    aagaagCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCAtcgtc                                                                                                                              <  4:356352/149‑1 (MQ=255)                               aTGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACgg                                                                                                                   <  3:179692/149‑1 (MQ=255)                                 ggAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTg                                                                                                                 >  4:303399/1‑149 (MQ=255)                                  gAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTgg                                                                                                                >  2:282517/1‑149 (MQ=255)                                  gAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTgg                                                                                                                >  4:127697/1‑149 (MQ=255)                                  gAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTg                                                                                                                 <  6:97440/148‑1 (MQ=255)                                   aGCTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGAACGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGc                                                                                                                                             >  8:170441/1‑119 (MQ=255)                                       gATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCA‑Tgg                                                                                                          >  5:111954/1‑149 (MQ=255)                                            tCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCA‑TGGCCCGt                                                                                                     >  6:385396/1‑149 (MQ=255)                                            tCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCA‑TGGCCCGt                                                                                                     >  1:55920/1‑149 (MQ=255)                                             cAAGCACACGCTGCCGGCCGAACTACTGAACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGCGGCA‑TGGCCCGTc                                                                                                    >  5:33561/1‑149 (MQ=255)                                                 cacaCGCTGCCAGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCA‑TGGCCCGTCACgg                                                                                                <  5:48462/149‑1 (MQ=255)                                                          ccggcctAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCA‑TGGCCCGTCACGGGGGTGgcgc                                                                                       <  3:257337/149‑1 (MQ=255)                                                          ccggccgAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCA‑TGGCCCGTCACGGGGGTGgc                                                                                         <  6:111954/147‑1 (MQ=255)                                                           cggccgAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCA‑TGGCCCGTCACGGGGGTGGcgcg                                                                                      <  3:303399/149‑1 (MQ=255)                                                            ggccgAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGatc                                                                                                                                    <  5:216901/103‑1 (MQ=255)                                                            ggccgAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGatc                                                                                                                                    >  6:216901/1‑103 (MQ=255)                                                            ggccgAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTACATCATCGGTGGCCCGGTGGGTGACTGCGCCCTGACTGCCCGCAAGATCATCGTCGACTCCTACGGAGGCA‑TGGCCCGTCACGGGGGTGGCGCGt                                                                                     >  5:8448/1‑149 (MQ=255)                                                             gccgAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCA‑TGGCCCGTCACGGGGGTGGCGCGtt                                                                                    <  2:55920/149‑1 (MQ=255)                                                                 aaCTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCA‑TGGCCCGTCACGGGGGTGGCGCGTtctc                                                                                 <  8:237263/148‑1 (MQ=255)                                                                 aaCTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCA‑TGGCCCGTCACGGGGGTGGCGCGTtctc                                                                                 >  7:237263/1‑148 (MQ=255)                                                                  actactGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCA‑TGGCCCGTCACGGGGGTGGCGCGTTCTCCg                                                                               >  4:149459/1‑149 (MQ=255)                                                                   ctactGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCA‑TGGCCCGTCACGGGGGTGGCGCGTTCTCCgg                                                                              <  7:249897/149‑1 (MQ=255)                                                                      ctGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCA‑TGGCCCGTCACGGGGGTGGCGCGTTCTCCGGc                                                                             <  3:127697/147‑1 (MQ=255)                                                                         cacaAAGGCACCCAGTACAACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACTGTGGCA‑TGGCCCGTAACGGGGGTGGCGCGTTCTCCGGCAAGGa                                                                        >  5:141229/1‑149 (MQ=255)                                                                            aaaGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCA‑TGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGAc                                                                       <  2:226670/147‑1 (MQ=255)                                                                            aaaGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCA‑TGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGAc                                                                       >  1:226670/1‑147 (MQ=255)                                                                              aGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCA‑TGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGt                                                                   >  7:181184/1‑149 (MQ=255)                                                                                  aCCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCA‑TGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTcc                                                                 >  6:318147/1‑147 (MQ=255)                                                                                   cccAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCA‑TGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAg                                                              >  2:116824/1‑149 (MQ=255)                                                                                    ccAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCAtcgt                                                                                                                               >  7:343218/1‑84 (MQ=255)                                                                                    ccAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCAtcgt                                                                                                                               <  8:343218/84‑1 (MQ=255)                                                                                    ccAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGCGGCA‑TGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAgg                                                             >  2:9905/1‑149 (MQ=255)                                                                                       gTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCA‑TGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTc                                                           >  2:100978/1‑148 (MQ=255)                                                                                          ccACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCA‑TGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGAcc                                                       >  3:85908/1‑149 (MQ=255)                                                                                             caTCACCCCGACCGGGAACTGCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACACCTACGGTGGCA‑TGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGGCGACCGtt                                                    <  4:280662/149‑1 (MQ=255)                                                                                              aTCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCA‑TGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTc                                                   <  8:181184/149‑1 (MQ=255)                                                                                               tCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCACTGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTc                                                   <  1:116824/149‑1 (MQ=255)                                                                                                    ccGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCA‑TGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGc                                             <  1:282517/149‑1 (MQ=255)                                                                                                         cGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCA‑TGGccc                                                                                                       <  8:362370/86‑1 (MQ=255)                                                                                                         cGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCA‑TGGccc                                                                                                       >  7:362370/1‑86 (MQ=255)                                                                                                         cGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCA‑TGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACg                                        >  1:220140/1‑149 (MQ=255)                                                                                                             aaCTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCA‑TGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACgccgg                                    >  6:115649/1‑149 (MQ=255)                                                                                                              aCTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCA‑TGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACgccggc                                   <  3:149459/149‑1 (MQ=255)                                                                                                                ttCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCA‑TGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACgccggccg                                 >  6:355300/1‑149 (MQ=255)                                                                                                                ttCATCATCGGTGGCCCGGTGGGTCACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCA‑TGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACgccggccg                                 >  1:310601/1‑149 (MQ=255)                                                                                                                   atcatcGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCTCAAGATCATCGTCGACTCCTACGGCGGCA‑TGGCCCGTCACGGGTGTGTCGCGTTATCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCGCt                               <  6:33561/148‑1 (MQ=255)                                                                                                                   atcatcGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCA‑TGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCGCTa                              >  5:106893/1‑149 (MQ=255)                                                                                                                          gTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCA‑TGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCGCTACGTGGc                        >  8:38334/1‑148 (MQ=255)                                                                                                                          gTGGCCCGGTGGGTCACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCA‑TGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCGCTACGTGGcc                       <  2:310601/149‑1 (MQ=255)                                                                                                                           tGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCA‑TGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCGCTACGTGGcc                       >  1:213294/1‑148 (MQ=255)                                                                                                                                cGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCA‑TGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCGCTACGTGGCCAAGAAc                 <  5:318147/149‑1 (MQ=255)                                                                                                                                          ttGCGGCCAGTCAGGCCGCAAGATCATCGTCGACTCCTACGGTGGCA‑TGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCGCTACGTGGCCAAGAACATCGTcgccg       >  7:132364/2‑149 (MQ=255)                                                                                                                                             aggCCTGACTGGCCGCAAGATCATCGTTGACTCCTACGGTGGCA‑TGGCCCGTGACGGGGGTGGCGCGTTCTCCGGAAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGGCGGCCGCTACGTGGCCAATAACATCGTCgccgcc     <  6:8448/147‑1 (MQ=255)                                                                                                                                               gCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCA‑TGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCAGGCCGCTACGTGGCCAAGAACATCGTCGCCGCCGGt  >  5:232037/1‑149 (MQ=255)                                                                                                                                                ccTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCA‑TGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGc                                             >  6:377118/1‑105 (MQ=255)                                                                                                                                                ccTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCA‑TGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGc                                             <  5:377118/105‑1 (MQ=255)                                                                                                                                                  tGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCA‑TGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCGCTACGTGGCCAAGAACAt               >  6:255636/1‑133 (MQ=255)                                                                                                                                                  tGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCA‑TGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCGCTACGTGGCCAAGAACAt               <  5:255636/133‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                   GTTTCGCAGAAAGACCTGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCA‑TGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCGCTACGTGGCCAAGAACATCGTCGCCGCCGGT  >  NC_002947/5659539‑5659829

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