Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,633,806	A→T	5.5%	D252V (GAC→GTC)	lepB →	signal peptidase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,633,806	0	A	T	5.5%	222.8 / 2.3	72	D252V (GAC→GTC)	lepB	signal peptidase I
Reads supporting (aligned to +/- strand):  ref base A (37/31);  new base T (2/2);  total (39/33)
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.23e-01

TGAGCCGTTATCGCATGCCGCCGGACCAGCAGTGGACCGTGCCTGCGGGCCATTACTTCATGATGGGCGACAACCGCGACAACTCCAACGATAGCCGTTTCTGGGATGACCCGAACATTCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACAT‑TGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCT  >  NC_002947/1633659‑1633949                                                                                                                                                    |                                                                                                                                                 tGAGCCGTTATCGCATGCCGCCGGACCAGCAGTGGACCGTGCCTGCGGGCCATTACTTCATGATGGGCGACAACCGCGACAACTCCAACGATAGCCGTTTCTGGGATGACCCGAACATTCCCAAGGAACTGCACGGCATGGTTccggac                                                                                                                                                 <  2:319113/149‑1 (MQ=255)  gAGCCGTTATCGCATGCCGCCGGACCAGCAGTGGACCGTGCCTGCGGGCCATTACTTCATGATGGGCGACAACCGCGACAACTCCAACGATAGCCGTTTCTGGGATGACCCGAACATTCCCAAGGAACTGCACGGCATGGTTccggacc                                                                                                                                                >  7:224674/1‑149 (MQ=255)             gCATGCCGCCGGACCAGCAGTGGACCGTGCCTGCGGGCCATTACTTCATGATGGGCGACAACCGCGACAACTCCAACGATAGCCGTTTCTGGGATGACCCGAACATTCCTAAGGAACTGCACGGCATGGTTCCGGTCCGGAAccatgcc                                                                                                                                     <  1:108877/149‑7 (MQ=255)              cATGCCGCCGGACCAGCAGTGGACCGTGCCTGCGGGCCATTACTTCATGATGGGCGACAACCGCGACAACTCCAACGATAGCCGTTTCTGGGATGACCCGAACATTCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACAT‑TGTCg                                                                                                                                   <  1:357384/149‑1 (MQ=255)              cATGCCGCCGGACCAGCAGTGGACCGTGCCTGCGGGCCATTACTTCATGATGGGCGACAACCGCGACAACTCCAACGATAGCCGTTTCTGCGATGACCCGAACATTCCCAAGGAACTGCACGGCATGGTTCCGGTCCGGAAccat                                                                                                                                        <  5:223711/145‑4 (MQ=255)              cATGCCGCCGGACCAGCAGTGGACCGTGCCTGCGGGCCATTACTTCATGATGGGCGACAACCGCGACAACTCCAACGATAGCCGTTTCTGCGATGACCCGAACATTCCCAAGGAACTGCACGGCATGGTTCCGGTCCGGAAccat                                                                                                                                        >  6:223711/1‑142 (MQ=255)                    gccgTACCAGAATTGGACCGTGCCTGCGGGCCATTACTTCATGATGGGCGACAACCGCTACAACTCCAACGATAGCCGTTTCTGTCATGACCCGAACATTCCCAACGAACTGCACGGCATGGTTCCGGACCGGAACAT‑TGTCGGc                                                                                                                                 <  5:277736/145‑1 (MQ=255)                    gccgGACCAGCAGTGGACCGTGCCTGCGGGCAATTACTTCATGATGGGCGACAACCGCGACAACTCCAACGATAGCCGTTTCTGGGATGACCAGAACATTCCCAAGGAACTGAACGGCATGGTTCCGGAACGGAAGCT‑TGTCGGc                                                                                                                                 >  6:277736/1‑145 (MQ=255)                     ccgGACCAGCAGTGGACCGTGCCTGCGGGCCATTACTTCATGATGGGCGACAACCGCGACAACTCCAACGATAGCCGTTTCTGGGATGACCCTAACATTCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACAT‑TGTCGGCAAGGc                                                                                                                            <  8:243946/149‑1 (MQ=255)                        gACCAGCAGTGGACCGTGCCTGCGGGCCATTACTTCATGATGGGCGACAACCGCGACAACTCCAACGATAGCCGTTTCTGGGATGACCCGAACATTCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACAT‑TGTCGGCAAGGCCtt                                                                                                                         >  1:307413/1‑149 (MQ=255)                          ccAGCAGTGGACCGTGCCTGCGGGCCATTACTTCATGATGGGCGACAACCGCGACAACTCCAACGATAGCCGTTTCTGGGATGACCCGAACATTCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACAT‑TGTCGGCAAGGCCTTcg                                                                                                                       >  3:207794/1‑149 (MQ=255)                          ccAGCAGTGGACCGTGCCTGCGGGCCATTACTTCATGATGGGCGACAACCGCGACAACTCCAACGATAGCCGTTTCTGGGATGACCCGAACATTCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACAT‑TGTCGGCAAGGCCTTcg                                                                                                                       <  8:224674/149‑1 (MQ=255)                           cagcagTGGACCGTGCCTGCGGGCCATTACTTCATGATGGGCGACAACCGCGACAACTCCAACGATAGCCGTTTCTGGGATGACCCGAACATTCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACAT‑TGTCGGCAAGGCCTTcgc                                                                                                                      >  8:165055/1‑149 (MQ=255)                              cagTGGACCGTGCCTGCGGGCCATTACTTCATGATGGGCGACAACCGCGACAACTCCAACGATAGCCGTTTCTGGGATGACCCGAACATTCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACAT‑TGTCGGCAAGGCCTTCGCGgt                                                                                                                   >  4:331451/1‑149 (MQ=255)                                gTGGACCGTGCCTGCGGGCCATTACTTCATGATGGGCGACAACCGCGACAACTCCAACGATAGCCGTTTCTGGGATGACCCGAACATTCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACAT‑TGTCGGCAAGGCCTTCGCGgtgt                                                                                                                 <  8:241345/149‑1 (MQ=255)                                   gACCGTGCCTGCGGGCCATTACTTCATGATGGGCGACAACCGCGACAACTCCAACGATAGCCGTTTCTGGGATGACCCGAACATTCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACAT‑TGTCGGCAAGGCCTTCGCGGTGTGGa                                                                                                              >  2:201832/1‑149 (MQ=255)                                     ccGTGCCTGCGGGCCATTACTTCATGATGGGCGACAACCGCGACAACTCCAACGATAGCCGTTTCTGGGATGACCCGAACATTCCCAAGGAACTGCACGGCATGGTTCCGGACAGGAACAT‑TGTCGGCAAGGCCTTCGCGGTGTGGAt                                                                                                             >  7:307445/1‑148 (MQ=255)                                                 gCCATTACTTCATGATGGGCGACAACCGCGACAACTCCAAAGATAGCCGTTTCTGGGATGACCCTAACATTCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACAT‑TGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCagag                                                                                                <  4:207794/149‑1 (MQ=255)                                                        cTTCATGATGGGCGACAACCGCGACAACTCCAACGATAGCCGTTTCTGGGATGACCCGAACATTCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACAT‑TGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCaaa                                                                                          <  7:295598/148‑1 (MQ=255)                                                        cTTCATGATGGGCGACAACCGCGACAACTCCAACGATAGCCGTTTCTGGGATGACCCGAACATTCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACAT‑TGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAAc                                                                                         >  4:225057/1‑149 (MQ=255)                                                           cATGATGGGCGACAACCGCGACAACTCCAACGATAGCCGTTTCTGGGATGACCCGAACATTCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACAT‑TGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACt                                                                                        >  3:197359/1‑147 (MQ=255)                                                                  ggCGACAACCGCGACAACTCCAACGATAGCCGTTTCTGGGATGACCCGAACATTCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACAT‑TGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACtt                                                                               <  3:285054/149‑1 (MQ=255)                                                                   gCGACAACCGCGACAACTCCAACGATAGCCGTTTCTGGGATGACCCGAACATTCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACAT‑TGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTg                                                                              >  6:91217/1‑149 (MQ=255)                                                                    cGACAACCGCGACAACTCCAACGATAGCCGTTTCTGGGATGACCCGAACATTCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACAT‑TGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGc                                                                             >  1:358919/1‑149 (MQ=255)                                                                         aCCGCGACAACTCCAACGATAGCCGTTTCTGGGATGACCCGAACATTCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACAT‑TGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAAc                                                                        <  3:225057/149‑1 (MQ=255)                                                                              gACAACTCCAACGATAGCCGTTTCTGGGATGACCCGAACATTCCCAAGGAACTGCATGGCATGGTTCCGGACCGGAACAT‑TGTc                                                                                                                                    <  3:181430/84‑1 (MQ=255)                                                                              gACAACTCCAACGATAGCCGTTTCTGGGATGACCCGAACATTCCCAAGGAACTGCATGGCATGGTTCCGGACCGGAACAT‑TGTc                                                                                                                                    >  4:181430/1‑84 (MQ=255)                                                                                    tCCAACGATAGCCGTTTCTGGGATGACCCGAACATTCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACAT‑TGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCagag                                                                                                <  2:309652/114‑1 (MQ=255)                                                                                    tCCAACGATAGCCGTTTCTGGGATGACCCGAACATTCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACAT‑TGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGt                                                             >  7:224917/1‑149 (MQ=255)                                                                                    tCCAACGATAGCCGTTTCTGGGATGACCCGAACATTCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACAT‑TGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGt                                                             >  8:273863/1‑149 (MQ=255)                                                                                    tCCAACGATAGCCGTTTCTGGGATGACCCGAACATTCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACAT‑TGTCGGCAAGGACTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCACAACATGTCCCGGGt                                                             >  8:47745/1‑149 (MQ=255)                                                                                      cAACGATAGCCGTTTCTGGGATGACCCGAACATTCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACAT‑TGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTg                                                           >  4:178216/1‑149 (MQ=255)                                                                                      cAACGATAGCCGTTTCTGGGATGACCCGAACATTCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACAT‑TGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTg                                                           <  7:273863/149‑1 (MQ=255)                                                                                      cAACGATAGCCGTTTCTGGGATGACCCGAACATTCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACAT‑TGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTg                                                           >  2:262423/1‑149 (MQ=255)                                                                                      cAACGATAGCCGTTTCTGGGATGACCCGAACATTCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACAT‑TGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTg                                                           >  3:57655/1‑149 (MQ=255)                                                                                        aCGATAGCCGTTTCTGGGATGACCCGAACATTCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACAT‑TGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGc                                                         >  8:170869/1‑149 (MQ=255)                                                                                        aCGATAGCCGTTTCTGGGATGACCCGAACATTCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACAT‑TGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGc                                                         <  4:57655/149‑1 (MQ=255)                                                                                            tAGCCGTTTCTGGGATGACCCGAACATTCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACAT‑TGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGa                                                     <  7:170869/149‑1 (MQ=255)                                                                                            tAGCCGTTTCTGGGATGACCCGAACATTCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACAT‑TGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGa                                                     >  3:160351/1‑149 (MQ=255)                                                                                              gCCGTTTCTGGGATGACCCGAACATTCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACAT‑TGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATc                                                   <  2:307413/149‑1 (MQ=255)                                                                                               ccGTTTCTGGGATGACCCGAACATTCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACAT‑TGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATc                                                   >  4:265128/1‑148 (MQ=255)                                                                                                  tttCTGGGATGACCCGAACATTCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACAT‑TGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCAtt                                               <  1:201832/149‑1 (MQ=255)                                                                                                  tttCTGGGATGACCCGAACATTCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACAT‑TGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCAtt                                               <  3:48373/149‑1 (MQ=255)                                                                                                   ttCTTGGATGACCCGAACATTCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACAT‑TGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCGAACCTGTCGCGGGTTGGCCTGATCCAtt                                               <  4:160351/148‑1 (MQ=255)                                                                                                     cTGGGATGACCCGAACATTCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACAT‑TGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGa                                             <  5:91217/148‑1 (MQ=255)                                                                                                          aTGACCCGAACATTCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACAT‑TGACGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCa                                       >  5:193774/1‑149 (MQ=255)                                                                                                               ccGAACATTCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACAT‑TGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGgc                                  >  1:267079/1‑149 (MQ=255)                                                                                                                  aaCATTCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACAT‑TGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCgg                                                               >  4:86232/1‑117 (MQ=255)                                                                                                                  aaCATTCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACAT‑TGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCgg                                                               <  3:86232/117‑1 (MQ=255)                                                                                                                   aCATTCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACAT‑TGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGc                                                                                               <  4:141768/84‑1 (MQ=255)                                                                                                                   aCATTCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACAT‑TGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGc                                                                                               >  3:141768/1‑84 (MQ=255)                                                                                                                     aTTCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACAT‑TGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTc                            >  2:19334/1‑149 (MQ=255)                                                                                                                         ccAAGGAACTGCACGGCATGGTTCCGGACCGGAACAT‑TGTCGGCAAGGCCTTCGCGGTGTGGATGa                                                                                                           <  6:59641/66‑1 (MQ=255)                                                                                                                         ccAAGGAACTGCACGGCATGGTTCCGGACCGGAACAT‑TGTCGGCAAGGCCTTCGCGGTGTGGATGa                                                                                                           >  5:59641/1‑66 (MQ=255)                                                                                                                         ccAAGGAACTGCACGGCATGGTTCCGGACCGGAACAT‑TGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCTCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCt                               <  8:56006/142‑1 (MQ=255)                                                                                                                         ccAAGGAACTGCACGGCATGGTTCCGGACCGGAACAT‑TGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCt                               >  7:56006/1‑142 (MQ=255)                                                                                                                         ccAAGGAACTGCACGGCATGGTTCCGGACCGGAACAT‑TGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCcag                         <  2:358919/148‑1 (MQ=255)                                                                                                                         ccAAGGAACTGCACGGCATGGTTCCGGACCGGAACAT‑TGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCcag                         >  2:25984/1‑148 (MQ=255)                                                                                                                          catggAACTGGACGGCATGGTTCCGGACCGGAACAT‑TGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCcagca                       <  3:20991/146‑1 (MQ=255)                                                                                                                          cAAGGAACTGCACGGCATGGTTCCGGACCGGAACAT‑TGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCcagca                       <  1:262423/149‑1 (MQ=255)                                                                                                                          cAAGGAACTGCACGGCATGGTTCCGGACCGGAACAT‑TGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCcagca                       <  2:267079/149‑1 (MQ=255)                                                                                                                          cAAGGAACTGCACGGCATGGTTCCGGACCGGAACAT‑TGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCcag                         >  8:354935/1‑147 (MQ=255)                                                                                                                          cAAGGAACTGCACGGCATGGTTCCGGACCGGAACAT‑TGTCGGCAAGGCCTTCGCGATGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCcagca                       <  4:214250/149‑1 (MQ=255)                                                                                                                          cAAGGAACTGCACGGCATGGTTCCGGACCGGAACAT‑TGTCGGCAAGGCCTTCGCGATGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCcag                         >  3:214250/1‑147 (MQ=255)                                                                                                                          cAAGGAACTGCACGGCATGGTTCCGGACCGGAACAT‑TGTCGGCAAGCCATTCGCGGTGTGGATGAGATGGCCAGAGCCCAAACTAAGCAAATTGCCCAACATGTCGCGGGTTGGCCTGATCCATTGATACCCAACGGCCCTGACcagca                       >  6:76114/1‑149 (MQ=255)                                                                                                                           aaGGAACTGCACGGCATGGTTCCGGACCGGAACATGTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCcag                         >  1:242957/1‑147 (MQ=255)                                                                                                                                    cACGGCATGGTTCCGGACCGGAACAT‑TGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGcc              <  1:19334/148‑1 (MQ=255)                                                                                                                                       ggCATGGTTCCGGACCGGAACAT‑TGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCACt                               <  6:185698/128‑1 (MQ=255)                                                                                                                                       ggCATGGTTCCGGACCGGAACAT‑TGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCACt                               >  5:185698/1‑128 (MQ=255)                                                                                                                                       gacaatgTTCCGGTCCGGAACAT‑TGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAAt          >  3:259337/7‑149 (MQ=255)                                                                                                                                          aTGGTTCCGGACCGGAACAT‑TGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCa       >  4:14390/1‑149 (MQ=255)                                                                                                                                           tGGTTCCGGACCGGAACAT‑TGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCAt      <  4:197359/149‑1 (MQ=255)                                                                                                                                           tGGTTCCGGACCGGAACAT‑TGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCAt      >  2:73573/1‑149 (MQ=255)                                                                                                                                              ttCCGGACCGGAACAT‑TGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCAAGCAGACAGCGCCGAATGCATTTc   >  8:82469/1‑149 (MQ=255)                                                                                                                                                ccggaccggaACAT‑TGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCt  >  4:261637/1‑148 (MQ=255)                                                                                                                                                ccggaccggaACAT‑TGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCt  <  1:73573/148‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                 TGAGCCGTTATCGCATGCCGCCGGACCAGCAGTGGACCGTGCCTGCGGGCCATTACTTCATGATGGGCGACAACCGCGACAACTCCAACGATAGCCGTTTCTGGGATGACCCGAACATTCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACAT‑TGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCT  >  NC_002947/1633659‑1633949

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