Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,637,741	G→T	6.2%	E107D (GAG→GAT)	PP_3204 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,637,741	0	G	T	6.2%	257.0 / 5.6	81	E107D (GAG→GAT)	PP_3204	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (37/39);  new base T (3/2);  total (40/41)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.63e-01

CCCTGGACAACGGCGATGATCCGACAGCGCAGCCGCTGCAATTGAGCCCTGGCCCGCTGGGCAGCGTGATTCGCGTGGTGGACATACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAA‑CCCTGGACTACGGCATCGTTACCGAGGGCGAGGTCTGGCTGGTGCTGGACGATGGCGAAGTCCACCTCAAACGAGGC  >  NC_002947/3637596‑3637885                                                                                                                                                  |                                                                                                                                                  cccTGGACAACGGCGATGATCCGACAGCGCAGCCGCTGCAATTGAGCCCTGGCCCGCTGGGCAGCGTGATTCGCGTGGTGGACATACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATc                                                                                                                                                <  7:142901/149‑1 (MQ=255)  ccTGGACAACGGCGATGATCCGACAGCGCAGCCGCTGCAATTGAGCCCTGGCCCGCTGGGCAGCGTGATTCGCGTGGTGGACATACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCg                                                                                                                                               <  6:294840/149‑1 (MQ=255)      gACAACGGCGATGATCCGACAGCGCAGCCGCTGCAATTGAGCCCTGGCCCGCTGGGCAGCGTGATTCGCGTGGTGGACATACCACCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCggc                                                                                                                                             <  5:52661/147‑1 (MQ=255)        cAACGGCGATGATCCGACAGCGCAGCCGCTGCAATTGAGCCCTGGCCCGCTGGGCAGCGTGATTCGCGTGGTGGACATACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCggcgagg                                                                                                                                         <  8:226315/149‑1 (MQ=255)           cGGCGATGATCCGACAGCGCAGCCGCTGCAATTGAGCCCTGGCCCGCTGGGCAGCGTGATTCGCGTGGTGGACATACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCggcgaggc                                                                                                                                        >  5:315853/1‑147 (MQ=255)              cGATGATCCGACAGCGCAGCCGCTGCAATTGAGCCCTGGCCCGCTGGGCAGCGTGATTCGCGTGGTGGACATACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTg                                                                                                                                   >  1:120540/1‑149 (MQ=255)               gatgatCCGACAGCGCAGCCGCTGCAATTGAGCCCTGGCCCGCTGGGCAGCGTGATTCGCGTGGTGGACATACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGc                                                                                                                                  >  2:249197/1‑149 (MQ=255)                 tgatCCGAACGCGCATCCGCTTCAATTGAGCCCTGGCCCGCTGGGCAGCGTGATTCGCGTTGTGGACATACCGCCCGACAGCGTGCAGAACCTGTTCAGTTCCGAAGATGCTGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCg                                                                                                                                <  2:348389/149‑1 (MQ=255)                  gatCCGACAGCGCAGCCGCTGCAATTGAGCGCTGGCCCGCTGGGCAGCGTGATTCGCGTGGTGGACATACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCgg                                                                                                                               >  1:31734/1‑149 (MQ=255)                     ccGACAGCGCAGCCGCTGCAATTGAGCCCTGGCCCGCTGGGCAGCGTGATTCGCGTGGTGGACATACCGCCCGACAGCGTGCAGAACCAGGGCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCAc                                                                                                                            <  2:111784/149‑1 (MQ=255)                       gACAGCGCAGCCGCTGCAATTGAGCCCTGGCCCGCTGGGCAGCGTGATTCGCGTGGTGGACATACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACgg                                                                                                                          >  3:82435/1‑149 (MQ=255)                            cgcAGCCGCTGCAATTGAGCCCTGGCCCGCTGGGCAGCGTGATTCGCGTGGTGGACATACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCa                                                                                                                       >  2:28396/1‑147 (MQ=255)                             gcAGCCGCTGCAATTGAGCCCTGGCCCGCTGGGCAGCGTGATTCGCGTGGTGGACATACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAg                                                                                                                    >  2:165792/1‑149 (MQ=255)                                 ccGCTGCAATTGAGCCCTGGCCCGCTGGGCAGCGTGATTCGCGTGGTGGACATACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGa                                                                                                                >  5:264964/1‑149 (MQ=255)                                   gCTGCAATTGAGCCCTGGCCCGCTGGGCAGCGTGATTCGCGTGGTGGACATACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGAtt                                                                                                              <  6:315853/149‑1 (MQ=255)                                     tGCAATTGAGCCCTGGCCCGCTGGGCAGCGTGATTCGCGTGGTGGACATACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGATCTCggcgaa                                                                                                                                          >  3:332460/1‑118 (MQ=255)                                     tGCAATTGAGCCCTGGCCCGCTGGGCAGCGTGATTCGCGTGGTGGACATACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGATCTCggcgaa                                                                                                                                          <  4:332460/119‑2 (MQ=255)                                      gCAATTGAGCCCTGGCCCGCTGGGCAGCGTGATTCGCGTGGTGGACATACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCg                                                                                                            >  3:251366/1‑148 (MQ=255)                                       cAATTGAGCCCTGGCCCGCTGGGCAGCGTGATTCGCGTGGTGGACATACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCg                                                                                                            >  6:238895/1‑147 (MQ=255)                                       cAATTGAGCCCTGGCCCGCTGGGCAGCGTGATTCGCGTGGTGGACATACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGaa                                                                                                          <  1:249197/149‑1 (MQ=255)                                           tGAGCCCTGGCCCGCTGGGCAGCGTGATTCGCGTGGTGGACATACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGca                                                                                                       >  1:303450/1‑148 (MQ=255)                                           tGAGCCCTGGCCCGCTGGGCAGCGTGATTCGCGTGGTGGACATACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGca                                                                                                       >  3:183700/1‑148 (MQ=255)                                                    gCCCGCTGGGCAGCGTGATTCGCGTGGTGGACATACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGa                                                                                               <  7:333527/147‑1 (MQ=255)                                                       cGCTGGGCAGCGTGATTCGCGTGGTGGACATACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATgcac                                                                                          >  1:103866/1‑149 (MQ=255)                                                         cTGGGCAGCGTGATTCGCGTGGTGGACATACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCGGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGc                                                                                                                   >  2:123803/1‑122 (MQ=255)                                                              cAGCGTGATTCGCGTGGTGGACATACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGcaccgcaccg                                                                                   <  5:238895/149‑1 (MQ=255)                                                                    gATTCGCGTGGTGGACATACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAA‑cc                                                                             <  8:211798/148‑1 (MQ=255)                                                                        cgcgTGGTGGACATACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAA‑CCCTgg                                                                         <  1:15224/148‑1 (MQ=255)                                                                           gtggtgGACATACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAA‑CCCTGGACTa                                                                     <  7:103137/149‑1 (MQ=255)                                                                             ggtggACATACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAA‑CCCTGGACTACg                                                                   >  3:336133/1‑149 (MQ=255)                                                                              gtggACATACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAA‑CCCTGGACTACg                                                                   <  4:82435/148‑1 (MQ=255)                                                                               tggACATACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGCAA‑CCCTGGACTACgg                                                                  >  7:336252/1‑148 (MQ=255)                                                                                 gACATACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCggcga                                                                                                                                           >  7:204264/1‑74 (MQ=255)                                                                                 gACATACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCggcga                                                                                                                                           <  8:204264/74‑1 (MQ=255)                                                                                     tACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCTCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTc                                                                                                             <  2:75211/100‑1 (MQ=255)                                                                                     tACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAA‑CCCTg                                                                          <  7:78406/134‑1 (MQ=255)                                                                                     tACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAA‑CCCTg                                                                          >  8:78406/1‑134 (MQ=255)                                                                                     tACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAA‑CCCTGGACTCCGGCATCGt                                                            >  3:221809/1‑148 (MQ=255)                                                                                     tACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCACCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTc                                                                                                             >  1:75211/1‑100 (MQ=255)                                                                                      aCCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAA‑CCCTGGACTACGGCATCGTTa                                                          >  6:372217/1‑149 (MQ=255)                                                                                      aCCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAA‑CCCTGGACTACGGCATCGTTa                                                          >  8:139554/1‑149 (MQ=255)                                                                                        cGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAA‑CCCTGGACTACGGCATCGTTAcc                                                        >  5:223166/1‑149 (MQ=255)                                                                                          cccGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAA‑CCCTGGACTACGGCATCGTTACCGa                                                      <  5:372217/149‑1 (MQ=255)                                                                                            cGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAA‑CCCTGGACTACGGCATCGTTACCGAgg                                                    <  7:148936/149‑1 (MQ=255)                                                                                            cGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAA‑CCCTGGACTACGGCATCGTTACCGAgg                                                    >  6:100515/1‑149 (MQ=255)                                                                                              aCAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAA‑CCCTGGACTACGGCATCGTTACCGAGGGc                                                  <  1:28396/149‑1 (MQ=255)                                                                                                  cGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAA‑CCCTGGACTACGGCATCGTTACCGAGGGCGAgg                                              <  4:183700/149‑1 (MQ=255)                                                                                                      cAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAA‑CCCTGGACTACGGCATCGTTACCGAGGGCGAGGTctg                                          >  5:287552/1‑149 (MQ=255)                                                                                                        gAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAA‑CCCTGGACTACGGCATCGTTACCGAGGGCGAGGTctggc                                        <  2:31734/149‑1 (MQ=255)                                                                                                         aaCCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAA‑CCCTGGACTACGGCATCGTTACCGAg                                                     <  3:354581/135‑1 (MQ=255)                                                                                                         aaCCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAA‑CCCTGGACTACGGCATCGTTACCGAg                                                     >  4:354581/1‑135 (MQ=255)                                                                                                         aaCCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAA‑CCCTGGACTACGGCATCGTTACCGAGGGCGAGGTctggct                                       >  6:187776/1‑149 (MQ=255)                                                                                                          aCCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAA‑CCCTGGACTACGGCATCGTTACCGAGGGCGAGGTctggctg                                      >  2:301494/1‑149 (MQ=255)                                                                                                           ccAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAA‑CCCTGGACTACGGCATCGTTACCGAGGGCGAGGTctggctgg                                     >  6:286884/1‑149 (MQ=255)                                                                                                            cAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCAACGAAA‑CCCTGGACTACGGCATCGTTACCGAGGGCGAGGTCTGGCTGGt                                    <  6:264964/149‑1 (MQ=255)                                                                                                              ggTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAA‑CCCTGGACTACGGCATCGTTACCGAGGGCGAGGTCTGGCAGGTGc                                  <  2:303450/149‑1 (MQ=255)                                                                                                                tCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAA‑CCCTGGACTACGGCATCGTTACCGAGGGCGAGGTCTGGCTGGTGCTg                                >  3:57219/1‑149 (MQ=255)                                                                                                                tCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAA‑CCCTGGACTACGGCATCGTTACCGAGGGCGAGGTCTGGCTGGTGCTg                                <  5:100515/149‑1 (MQ=255)                                                                                                                   gTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAA‑CCCTGGACTACGGCATCGTTACCGAGGGCGAGGTCTGGCTGGTGCTGGAc                             <  4:336133/149‑1 (MQ=255)                                                                                                                    tGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAA‑CCCTGGACTACGGCATCGTTACCGAGGGCGAGGTCTGGCTGGTGCTGGAc                             >  6:258752/1‑148 (MQ=255)                                                                                                                       cGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAA‑CCCTGGACTACGGCATCGTTACCGAGGGCGAGGTCTGGCTGGTGCTGGACGATg                         >  7:231364/1‑149 (MQ=255)                                                                                                                       cGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAA‑CCCTGGACTACGGCATCGTTACCGAGGGCGAGGGCTGGCTGGTGCTGGACGATg                         >  6:120236/1‑149 (MQ=255)                                                                                                                       cGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAA‑CCCTGGACTACGGCAGCGTTACCGAGGGCGAGGTCTGGCTGGTGCTGGACGATg                         >  7:225479/1‑149 (MQ=255)                                                                                                                        gAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAA‑CCCTGGACTACGGCATCGTTACCGAGGGCGAGGTCTGGCTGGTGCTGGACGATgg                        <  5:187776/149‑1 (MQ=255)                                                                                                                             tGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAACCCCTGGACTACGGCATCGTTACCGAGGGCGAGGTCTGGCTGGTGCTGGACGATGGCGaa                    <  1:301494/149‑1 (MQ=255)                                                                                                                               cAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAA‑CCCTGGACTACGGCATCGTTACCGAGGGCGAGGTCTGGCTGGTGCTGGACGATGGCGAAGt                  <  1:165792/148‑1 (MQ=255)                                                                                                                               cAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAA‑CCCTGGACTACGGCATCGTTACCGAGGGCGAGGTCTGGCTGGTGCTGGACGATGGCGAAGTc                 <  6:147984/149‑1 (MQ=255)                                                                                                                               cAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAA‑CCCTGGACTACGGCATCGTTACCGAGGGCGAGGTCTGGCTGGTGCTGGACGATGGCGAAGTc                 <  6:287552/149‑1 (MQ=255)                                                                                                                                aGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAA‑CCCTGGACTACGGCATCGTTACCGAGGGCGAGGTCTGGCTGGTGCTGGACGATGGCGAAGTcc                <  2:103866/149‑1 (MQ=255)                                                                                                                                  ccGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAA‑CCCTGGACTACGGCATCGTTACCGAGGGCGAGGTCTGGCTGGTGCTGGACGATGGCGAAGTCCAc              <  3:12580/149‑1 (MQ=255)                                                                                                                                  ccGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAA‑CCCTGGACTACGGCATCGTTACCGAGGGCGAGGTCTGGCTGGTGCTGGACGATGGCGAAGTCCAc              <  4:46898/149‑1 (MQ=255)                                                                                                                                    actcgcctCGCCGATCTCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAA‑CCCTGGACTACGGCATCGTTACCGAGGGCGAGGTCTGGCTGGTGCTGGACGATGGCGAAGTc                 >  7:285927/8‑144 (MQ=255)                                                                                                                                    actcgcctCGCCGATCTCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAA‑CCCTGGACTACGGCATCGTTACCGAGGGCGAGGTCTGGCTGGTGCTGGACGATGGCGAAGTc                 <  8:285927/137‑1 (MQ=255)                                                                                                                                     cTGTGTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCAGATTCGAAGCACAAGCTGATGCACCGCACCGAAA‑CCCTGGACTACGGCATCGTTACCGAGGGCGAGGTCTGGCTGGTGCTGGACGATGGCGAAGTCCACCTc           <  6:223166/149‑1 (MQ=255)                                                                                                                                      tcgcctCGCCGATCTCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAA‑CCCTGGACTACGGCATCGTTACCGAGGGCGAGGTCTGGCTGGTGCTGGACGATGGCGAAGTCCACCTc           >  4:218535/6‑148 (MQ=255)                                                                                                                                       gTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAA‑CCCTGGACTACGGCATCGTTACCGAGGGCGAGGTCTGGCTGGTGCTGGACGATGGCGAAGTCCACCTCa          >  2:43558/1‑148 (MQ=255)                                                                                                                                          ttCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAA‑CCCTGGACTACGGCATCGTTACCGAGGGCGAGGTCTGGCTGGTGCTGGACGATGGCGAAGTCCACCTCAAACg      >  2:365785/1‑149 (MQ=255)                                                                                                                                          ttCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAA‑CCCTGCACTAAGGCATCGTTACCGAGGGCGAGGTCTGGCTGGTGCTGGCGGATGGCGAAGTCAACATCAAACg      >  7:143179/1‑149 (MQ=255)                                                                                                                                             gCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAA‑CCCTGGACTACGGCATCGTTACCGCGGGCGAGGTCTGGCTGGTGCTGGACGATGGCGAAGTCCACCTCAAACGAg    <  1:193669/148‑1 (MQ=255)                                                                                                                                             gCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAA‑CCCTGGACTACGGCATCGTTACCGAGGGCGAGGTCTGGCTGGTGCTGGACGATGGCGAAGTCCACCTCAAACGAgg   <  1:365785/149‑1 (MQ=255)                                                                                                                                              ccGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAA‑CCCTGGACTACGGCATCGTTACCGAGGGCGAGGTCTGGCTGGTTCTGGACGATGGCGAAGTCCACCTCAAACGAGGc  <  4:251366/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                  CCCTGGACAACGGCGATGATCCGACAGCGCAGCCGCTGCAATTGAGCCCTGGCCCGCTGGGCAGCGTGATTCGCGTGGTGGACATACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAA‑CCCTGGACTACGGCATCGTTACCGAGGGCGAGGTCTGGCTGGTGCTGGACGATGGCGAAGTCCACCTCAAACGAGGC  >  NC_002947/3637596‑3637885

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