Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,042,008	T→A	6.3%	D32E (GAT→GAA)	PP_3563 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,042,008	0	T	A	6.3%	190.3 / 3.1	63	D32E (GAT→GAA)	PP_3563	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (27/32);  new base A (2/2);  total (29/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.76e-01

GATGAGCAGCGCACACCACGACGCTGCCAGTCATAGGGAGATCCACATGCGCGTCCTGACATCCGTTGTCGCCTTGTCACTGGGCCTGGTCGTTTCGACCGGGGCGCTGGCAGCCACCGGCGAAGAGGCGCAGTTGATCGATTCGATCAATGTCTACCGCAGCAAGGCCCAGCCCTGCGGTGGCGAGGCGTCGTTGGAACTGCCGCCGCTCAACAGCGACACACGCCTTGCGCTGTCGCCAGAGGGCACCCGCGACCTGCAACAGGCCATGACCCGCGCCGCCTACCCCA  >  NC_002947/4041867‑4042156                                                                                                                                              |                                                                                                                                                     gATGAGCAGCGCACACCACGACGCTGCCAGTCATAGGGAGATCCACATGCGCGTCCTGACATCCGTTGTCGCCTTGTCACTGGGCCTGGTCGTTTCGACCGGGGCGCTGGCAGCCACCGGCGAAGAGGCGCAGTTGATCGATTCGATc                                                                                                                                                >  2:51513/1‑148 (MQ=255) gATGAGCAGCGCACACCACGACGCTGCCAGTCATAGGGAGATCCACATGCGCGTCCTGACATCCGTTGTCGCCTTGTCACTGGCCCTGGTCGTTTCGACCGGGGCGCTGGCAGCCACCGGCGAAGAGGCGCAGTTGATCGATTCGATc                                                                                                                                                >  8:39216/1‑148 (MQ=255)   tGAGCAGCGCACACCACGACGCTGCCAGTCATAGGGAGATCCACATGCGCGTCCTGACATCCGTTGTCGCCTTGTCACTGGGCCTGGTCGTTTCGACCGGGGCGCTGGCAGCCACCGGCGAAGAGGCGCAGTTGATCGATTCGATCAAt                                                                                                                                             >  8:5301/1‑149 (MQ=255)   tGAGCAGCGCACACCACGACGCTGCCAGTCATAGGGAGATCCACATGCGCGTCCTGACATCCGTTGTCGCCTTGTCACTGGGCCTGGTCGTTTCGACCGGGGCGCTGGCAGCCACCGGCGAAGAGGCGCAGTTGATCGAATCGATCaac                                                                                                                                             <  5:383732/149‑2 (MQ=255)     agcagcGCACACCACGACGCTGCCAGTCATAGGGAGATCCACATGCGCGTCCTGACATCCGTTGTCGCCTTGTCACTGGGCCTGGTCGTTTCGACCGGGGCGCTGGCAGCCACCGGCGAAGAGGCGCAGTTGATCGAATCGATCaactg                                                                                                                                           <  3:388426/149‑4 (MQ=255)      gcagcGCACACCACGACGCTGCCAGTCTTAGGGAGATCCACATGCGCGTCCTGACATCCGTTGTCGCCTTGTCACTGGGCCTGGTCGTTTCGACCGGGGCGCTGGCAGCCACCGGCGAAGAGGCGCAGTTGATCGATTCGATCAATGTc                                                                                                                                          >  5:180696/1‑149 (MQ=255)        agcGCACACCACGACGCTGCCAGTCATAGGGAGATCCACATGCGCGTCCTGACATCCGTTGTCGCCTTGTCACTGGGCCTGGTCGTTTCGACCGGGGCGCTGGCAGCCACCGGCGAAGAGGCGCAGTTGATCGATTCGATCAATGTCTa                                                                                                                                        <  2:74739/149‑1 (MQ=255)                 cacGACGCTGCCAGTCATAGGGAGATCCACATGCGCGTCCTGACATCCGTTGTCGCCTTGTCACTGGGCCTGGTCGTTTCGACCGGGGCGCTGGCAGCCACCGGCGAAGAGGCGCAGTTGATCGATTCGATCAATGTCTACCgcagc                                                                                                                                 >  7:288067/1‑147 (MQ=255)                   cgacgCTGCCAGTCATAGGGAGATCCACATGCGCGTCCTGACATCCGTTGTCGCCTTGTCACTGGGCCTGGTCGTTTCGACCGGGGCGCTGGCAGCCACCGGCGAAGAGGCGCAGTTGATCGATTCGATCAATGTCTACCGCAGCAAgg                                                                                                                             <  7:135825/149‑1 (MQ=255)                           ccAGTCATAGGGAGATCCACATGCGCGTCCTGACATCCGTTGTCGCCTTGTCACTGGGCCTGGTCGTTTCGACCGGGGCGCTGGCAGCCACCGGCGAAGAGGCGCAGTTGATCGATTCGATCAATGTCTACCGCAGCAAGGCCCAGccc                                                                                                                     <  4:370426/149‑1 (MQ=255)                            cAGTCATAGGGAGATCCACATGCGCGTCCTGACATCCGTTGTCGCCTTGTCACTGGGCCTGGTCGTTTCGACCGGGGCGCTGGCAGCCACCGGCGAAGAGGCGCAGTTGATCGATTCGATCAATGTCTACCGCAGCAAGGCCCAGCCCt                                                                                                                    >  7:131839/1‑149 (MQ=255)                              gTCATAGGGAGATCCACATGCGCGTCCTGACATCCGTTGTCGCCTTGTCACTGGGCCTGGTCGTTTCGACCGGGGCGCTGGCAGCCACCGGCGAAGAGGCGCAGTTGATCGATTCGATCAATGTCTACCGCAGCAAGGCCCAGCCCTGc                                                                                                                  <  8:288067/149‑1 (MQ=255)                                 aTAGGGAGATCCACATGCGCGTCCTGACATCCGTTGTCGCCTTGTCACTGGGCCTGGTCGTTTCGACCGGGGCGCTGGCAGCCACCGGCGAAGAGGCGCAGTTGATCGATTCGATCAATGTCTACCGCAGCAAGGCCCAGCCCTGCGGt                                                                                                               <  2:192731/149‑1 (MQ=255)                                 aTAGGGAGATCCACATGCGCGTCCTGACATCCGTTGTCGCCTTGTCACTGGGCCTGGTCGTTTCGACCGGGGCGCTGGCAGCCACCGGCGAAGAGGCGCAGTTGATCGATTCGATCAATGTCTACCGCAGCAAGGCCCAGCCCTGCGGt                                                                                                               >  6:205538/1‑149 (MQ=255)                                         aTCCACATGCGCGTCCTGACATCCGTTGTCGCCTTGTCACTGGGCCTGGTCGTTTCGACCGGGGCGCTGGCAGCCACCGGCGAAGAGGCGCAGTTGATCGATTCGATCAATGTCTACCGCAGCAAGGCCCAGCCCTGCGGTGGCGAGGc                                                                                                       <  2:23572/149‑1 (MQ=255)                                         aTCCACATGCGCGTCCTGACATCCGTTGTCGCCTTGTCACTGGGCCTGGTCGTTTCGACCGGGGCGCTGGCAGCCACCGGCGAAGAGGCGCAGTTGATCGATTCGATCAATGTCTACCGCAGCAAGGCCCAGCCCTGCGGTGGCGAGGc                                                                                                       <  1:59038/149‑1 (MQ=255)                                             acaTGCGCGTCCTGACATCCGTTGTCGCCTTGTCACTGGGCCTGGTCGTTTCGACCGGGGCGCTGGCAGCCACCGGCGAAGAGGCGCAGTTGATCGATTCGATCAATGTCTACCGCAGCAAGGCCCAGCCCTGCGGTGGCGAGGcgtcg                                                                                                   >  4:249508/1‑149 (MQ=255)                                                tGCGCGTCCTGACATCCGTTGTCGCCTTGTCACTGGGCCTGGTCGTTTCGACCGGGGCGCTGGCAGCCACCGGCGAAGAGGCGCAGTTGATCGATTCGATCAATGTCTACCGCAGCAAGGCCCAGCCCTGCGGTGGCGAGGCGTCGTTg                                                                                                <  3:201054/149‑1 (MQ=255)                                                  cgcgTCCTGACATCCGTTGTCGCCTTGTAACTGGGCCTGGTCGTTTCGACCGGGGCGCTGGCAGCCACCGGCGAAGAGGCGCAGTTGATCGATTCGATCAATGTCTACCGCAGCAAGGCCCAGCCCTGCGGTGGCGAGGCGTCGTTGGa                                                                                              <  3:109264/149‑1 (MQ=255)                                                   gcgTCCTGCCATCCGTTGTCGCCTTGTCACTGGGCCTGGTCGTTTCGACCGGGGCGCTGGCAGCCACCGGCGAAGAGGCGCAGTTGATCGATTCGATCAATGTCTACCGCAGCAAGGCCCAGCCCTGCGGTGGCGAGGCGTCGTTGGaa                                                                                             <  7:127778/149‑1 (MQ=255)                                                         tGACATCCGTTGTCGCCTTGTCACTGGGCCTGGTCGTTTCGACCGGGGCGCTGGCAGCCACCGGCGAAGAGGCGCAGTTGATCGATTCGATCAATGTCTACCGCAGCAAGGccta                                                                                                                         <  4:290602/115‑3 (MQ=255)                                                         tGACATCCGTTGTCGCCTTGTCACTGGGCCTGGTCGTTTCGACCGGGGCGCTGGCAGCCACCGGCGAAGAGGCGCAGTTGATCGATTCGATCAATGTCTACCGCAGCAAGGccga                                                                                                                         >  3:290602/1‑113 (MQ=255)                                                         tGACATCCGTTGTCGCCTTGTCACTGGGCCTGGTCGTTTCGACCGGGGCGCTGGCAGCCACCGGCGAAGAGGCGCAGTTGATCGATTCGATCAATGTCTACCGCAGCAAGGCCCAGCCCTGCGGTGGCGAGGCGTCGTTGGAACTgccg                                                                                       >  2:382884/1‑149 (MQ=255)                                                                  ttGTCGCCTTGTCACTGGGCCTGGTCGTTTCGACCGGGGCGCTGGCAGCCACCGGCGAAGAGGCGCAGTTGATCGATTCGATCAATGTCTACCGCATCAAGGCCCAGCCCTGCGGTGGCGAGGCGTCGTTGGAACTTCCGCCGCTCAAc                                                                              <  6:89904/149‑1 (MQ=255)                                                                  ttGTCGCCTTGTCACTGGGCCTGGTCGTTTCGACCGGGGCGCTGGCAGCCACCGGCGAAGAGGCGCAGTTGATCGATTCGATCAATGTCTACCGCAGCAAGGCCCAGCCCTGCGGTGGCGAGGCGTCGTTGGAACTGCCGCCGCTCaa                                                                               >  5:337315/1‑148 (MQ=255)                                                                     tCGCCTTGTCACTGGGCCTGGTCGTTTCGACCGGGGCGCTGGCAGCCACCTGCGAAGAGGCGCAGTTGATCGATTCGATCAATGTCTACCGCAGCAAGGCCCAGCCCTGCGGTGGCGAGGCGTCGTTGGAACTGCCGCCGCTCAACAGc                                                                           <  4:115865/149‑1 (MQ=255)                                                                     tCGCCTTGTCACTGGGCCTGGTCGTTTCGACCGGGGCGCTGGCAGCCACCGGCGAAGAGGCGCAGTTGATCGATTCGATCAATGTCTACCGCAGCAAGGCCCAGCCCTGCGGTGGCGAGGCGTCGTTGGAACTGCCGCCGCTCAACAGc                                                                           >  7:352762/1‑149 (MQ=255)                                                                       gCCTTGTCACTGGGCCTGGTCGTTTCGACCGGGGCGCTGGCAGCCACCGGCGAAGAGGCGCAGTTGATCGATTCGATCAATGTCTACCGCAGCAAGGCCCAGCCCTGCGGTGGCGAGGCGTCGTTGGAACTGCCGCCGCTCAACAGCGa                                                                         >  6:192652/1‑149 (MQ=255)                                                                        ccTTGTCACTGGGCCTGGTCGTTTCGACCGGGGCGCTGGCAGCCACCGGCGAAGAGGCGCAGTTGATCGATTCGATCGATGTCTACCGCAGCAAGGCCCAGCCCTGCGGTGGCGAGGCGTCGTTGGAACTGCCGCCGCTCAACAGCg                                                                          <  6:310178/147‑1 (MQ=255)                                                                        ccTTGTCACTGGGCCTGGTCGTTTCGACCGGGGCGCTGGCAGCCACCGGCGAAGAGGCGCAGTTGATCGATTCGATCGATGTCTACCGCAGCAAGGCCCAGCCCTGCGGTGGCGAGGCGTCGTTGGAACTGCCGCCGCTCAACAGCg                                                                          >  5:310178/1‑147 (MQ=255)                                                                               aCTGGGCCTGGTCGTTTCGACCGGGGCGCTGGCAGCCACCGGCGAAGAGGCGCAGTTGATCGATTCGATCAATGTCTACCGCAGCAAGGCCCAGCCCTGCGGTGGCGAGGCGTCGTTGGAACTGCCGCCGCTCAACAGCGACACACGcc                                                                 >  7:183447/1‑149 (MQ=255)                                                                               aCTGGGCCTGGTCGTTTCGACCGGGGCGCTGGCAGCCACCGGCGAAGAGGCGCAGTTGATCGATTCGATCAATGTCTACCGCAGCAAGGCCCAGCCCTGCGGTGGCGAGGCGTCGTTGGAACTGCCGCCGCTCAACAGCGACACACGcc                                                                 >  7:185290/1‑149 (MQ=255)                                                                                     ccTGGTCGTTTCGACCGGGGCGCTGGCAGCCACCGGCGAAGAGGCGCAGTTGATCGATTCGATCAATGTCTACCGCAGCAAGGCCCAGCCCTGCGGTGGCGAGGCGTCGTTGGAACTGCCGCCTCTCAACAGCGACACACGCCTTgcgc                                                           <  7:288258/149‑1 (MQ=255)                                                                                      cTGGTCGTTTCGACCGGGGCGCTGGCAGCCACCGGCGAAGAGGCGCAGTTGATCGATTCGATCAATGTCTACCGCAGCAAGGCCCAGCCCTGCGGTGGCGAGGCGTCGTTGGAACTGCCGCCGCTCAACAGCGACACACGCCTTGCGCt                                                          <  7:39216/149‑1 (MQ=255)                                                                                      cTGGTCGTTTCGACCGGGGCGCTGGCAGCCACCGGCGAAGAGGCGCAGTTGATCGATTCGATCAATGTCTACCGCAGCAAGGCCCAGCCCTGCGGTGGCGAGGCGTCGTTGGAACTGCCGCCGCTCAACAGCGACACACGCCTTGCGCt                                                          <  1:51513/149‑1 (MQ=255)                                                                                         gTCGTTTCGACCGGGGCGCTGGCAGCCACCGGCGAAGAGGCGCAGTTGATCGATTCGATCAATGTCTACCGCAGCAAGGCCCAGCCCTGCGGTGGCGAGGCGTCGTTGGAACTGCCGCCGCTCAACAGCGACACACGCCTTGCGCTGt                                                        <  2:267579/148‑1 (MQ=255)                                                                                         gTCGTTTCGACCGGGGCGCTGGCAGCCACCGGCGAAGAGGCGCAGTTGATCGATTCGATCAATGTCTACCGCAGCAAGGCCCAGCCCTGCGGTGGCGAGGCGTCGTTGGAACTGCCGCCGCTCAACAGCGACACACGCCTTGCGCTGTc                                                       <  1:72568/149‑1 (MQ=255)                                                                                             tttCGACCGGGGCGCTGGCAGCCACCGGCGAAGAGGCGCAGTTGATCGATTCGATCAATGTCTACCGCAGCAAGGCCCAGCCCTGCGGTGGCGAGGCGTCGTTGGAACTGCCGCCGCTCAACAGCGACACACGCCTTGCGCTGTCGc                                                     >  2:209639/1‑147 (MQ=255)                                                                                             tttCGACCGGGGCGCTGGCAGCCACCGGCGAAGAGGCGCAGTTGATCGATTCGATCAATGTCTACCGCAGCAAGGCCCAGCCCTGCGGTGGCGAGGCGTCGTTGGAACTGCCGCCGCTCAACAGCGACACACGCCTTGCGCTGTCGc                                                     >  2:52040/1‑147 (MQ=255)                                                                                                 gACCGGGGCGCTGGCAGCCACCGGCGAAGAGGCGCAGTTGATAGATTCGATCAATGTCTACCGCAGCAAGGCCCAGCCCTGCGGTGGCGAGGCGTCGTTGGAACTGCCGCCGCTCAACAGCGACACACGCCTTGCGCTGTCGCAAGAgg                                               >  4:37562/1‑149 (MQ=255)                                                                                                  aCCGGGGCGCTGGCAGCCACCGGCGAAGAGGCGCAGTTGATCGATTCGATCAATGTCTACCGCAGCAAGGCCCAGCCCTGCGGTGGCGAGGCGTCGTTGGAACTGCCGCCGCTCAACAGCGACACACGCCTTGCGCTGTCGCCAGAggg                                              >  2:145688/1‑149 (MQ=255)                                                                                                      gggCGCTGGCAGCCACCGGCGAAGAGGCGCAGTTGATCGATTCGATCAATGTCTACCGCAGCAAGGCCCAGCCCTGCGGTGGCGAGGCGTCGTTGGAACTGCCGCCGCTCAACAGCGACACACGCCTTGCGCTGTCGCCAGAGGGCAcc                                          <  7:319780/149‑1 (MQ=255)                                                                                                           cTGGCAGCCACCGGCGAAGAGGCGCAGTTGATCGATTCGATCAATGTCTACCGCAGCAAGGCCCAGCCCTGCGGTGGCGAGGCGTCGTTGGAACTGCCGCCGCTCAACAGCGACACACGCCTTGCGCTGTCGCCAGAGGGCACCCGCGa                                     <  5:192652/149‑1 (MQ=255)                                                                                                            tGGCAGCCACCGGCGAAGAGGCGCAGTTGATCGATTCGATCAATGTCTACCGCAGCAAGGCCCAGCCCTGCGGTGGCGAGGCGTCGTTGGAACTGCCGCCGCTCAACAGCGACACACGCCTTGCGCTGTCGCCAGAGGGCACCCGCGAc                                    <  8:203895/149‑1 (MQ=255)                                                                                                              gCAGCCACCGGCGAAGAGGCGCAGTTGATCGATTCGATCAATGTCTACCGCAGCAAGGCCCAGCCCTGCGGTGGCGAGGCGTCGTTGGAACTGCCGCCGCTCAACAGCGACACACGCCTTGCGCTGTCGCCAGAGGGCACCCGCGACCt                                  <  8:185290/149‑1 (MQ=255)                                                                                                                aGCCACCGGCGAAGAGGCGCAGTTGATCGATTCGATCAATGTCTACCGCAGCAAGGCCCAGCCCTGCGGTGGCGAGGCGTCGTTGGAACTGCCGCCGCTCAACAGCGACACACGCCTTGCGCTGTCGCCAGAGGGCACCCGCGACCTg                                 >  6:253702/1‑148 (MQ=255)                                                                                                                aGCCACCGGCGAAGAGGCGCAGTTGATCGATTCGATCAATGTCTACCGCAGCAAGGCCCAGCCCTGCGGTGGCGAGGCGTCGTTGGAACTGCCGCCGCTCAACAGCGACACACGCCTTGCGCTGTCGCCAGAGGGCACCCGCGACCTGc                                <  6:180696/149‑1 (MQ=255)                                                                                                                    aCCGGCGAAGAGGCGCAGTTGATCGATTCGATCAATGTCTACCGCAGCAAGGCCCAGCCCTGCGGTGGCGAGGCGTCGTTGGAACTGCCGCCGCTCAACAGCGACACACGCCTTGCGCTGTCGCCAGAGGGCACCCGCGACCTGCAACa                            <  1:382884/149‑1 (MQ=255)                                                                                                                     ccGGCGAAGAGGCGCAGTTGATCGATTCGATCAATGTCTACCGCAGCAAGGCCCAGCCCTGCGGTGGCGAGGCGTCGTTGGAACTGCCGCCGCTCAACAGCGACACACGCCTTGCGCTGTCGCCAGAGGGCACCCGCGACCTGCAACAg                           <  8:183447/149‑1 (MQ=255)                                                                                                                      cGGCGAAGAGGCGCAGTTGATCGATTCGATCAATGTCTACCGCAGCAAGGCCCAGCCCTGCGGTGGCGAGGCGTCGTTGGAACTGCCGCCGCTCAACAGCGACACACGCCTTGCGCTGTCGCCAGAGGGCACCCGCGACCTGCAACAgg                          <  6:337315/149‑1 (MQ=255)                                                                                                                      cGGCGAAGAGGCGCAGTTGATCGATTCGATCAATGTCTACCGCAGCAAGGCCCAGCCCTGCGGTGGCGAGGCGTCGTTGGAACTGCCGCCGCTCAACAGCGACACACGCCTTGCGCTGTCGCCAGAGGGCACCCGCGACCTGCAACAgg                          >  2:91925/1‑149 (MQ=255)                                                                                                                         cGAAGAGGCGCAGTTGATCGATTCGATCAATGTCTACCGCAGCAAGGCCCAGCCCTGCGGTGGCGAGGCGTCGTTGGAACTGCCGCCGCTCAACAGCGACACACGCCTTGCGCTGTCGCCAGAGGGCACCCGCGACCTGCAACAGGCCa                       <  3:341388/149‑1 (MQ=255)                                                                                                                         cGAAGAGGCGCAGTTGATCGATTCGATCAATGTCTACCGCAGCAAGGCCCAGCCCTGCGGTGGCGAGGCGTCGTTGGAACTGCCGCCGCTCAACAGCGACACACGCCTTGCGCTGTCGCCAGAGGGCACCCGCGACCTGCAACAGGCCa                       <  3:335903/149‑1 (MQ=255)                                                                                                                         cGAAGAGGCGCAGTTGATCGATTCGATCAATGTCTACCGCAGCAAGGCCCAGCCCTGCGGTGGCGAGGCGTCGTTGGAACTGCCGCCGCTCAACAGCGACACACGCCTTGCGCTGTCGCCAGAGGGCACCCGCGACCTGCAACAGGCCa                       <  7:5301/149‑1 (MQ=255)                                                                                                                            agagGCGCAGTTGATCGATTCGATCAATGTCTACCGCAGCAAGGCCCAGCCCTGCGGTGGCGAGGCGTCGTTGGAACTGCCGCCGCTCAACAGCGACACACGCCTTGCGCTGTCGCCAGAGGTCACCCGCGACCTGCAACAGGCCATGa                    <  1:209639/149‑1 (MQ=255)                                                                                                                            agagGCGCAGTTGATCGATTCGATCAATGTCTACCGCAGCAAGGCCCAGCCCTGCGGTGGCGAGGCGTCGTTGGAACTGCCGCCGCTCAACAGCGACACACGCCTTGCGCTGTCGCCAGAGGGCACCCGCGACCTGCAACAGGCCATGa                    <  5:253702/149‑1 (MQ=255)                                                                                                                                gcgcAGTTGATCGATTCGATCAATGTCTACCGCAGCAAGGCCCAGCCCTGCGGTGGCGAGGCGTCGTTGGAACTGCCGCCGCTCAACAGCGACACACGCCTTGCGCTGTCGCCAGAGGGCACCCGCGACCTGCAACAGGCCATGACCcg                >  2:84979/1‑149 (MQ=255)                                                                                                                                  gcAGTTGATCGATTCGATCAATGTCTACCGCAGCAAGGCCCAGCCCTGCGGTGGCGAGGCGTCGTTGGAACTGCCGCCGCTCAACAGCGACACACGCCTTGCGCTGTCGCCAGAGGGCACCCGCGACCTGCAACAGGCCATGACCcgcg              >  3:15045/1‑149 (MQ=255)                                                                                                                                   acattGATCGAATCGATCAATGTCTACCGCAGCAAGGCCCAGCCCTGCGGTGGCGAGGCGTCGTTGGAACTGCCGCCGCTCAACAGCGACACACGCCTTGCGCTGTCGCCAGAGGGCACCCGCGACCTGCAACAGGCCATGACCCgcgc             >  5:114072/4‑149 (MQ=255)                                                                                                                                   acattGATCGAATCGATCAATGTCTACCGCAGCAAGGCCCAGCCCTGCGGTGGCGAGGCGTCGTTGGAACTGCCGCCGCTCAACAGCGACACACGCCTTGCGCTGTCGCCAGAGGGCACCCGCGACCTGCAACAGGCCATGACCCgcgc             >  5:30318/4‑149 (MQ=255)                                                                                                                                      ttGATCGATTCGATCAATGTCTACCGCAGCAAGGCCCAGCCCTGCGGTGGCGAGGCGTCGTTGGAACTGCCGCCGCTCAACAGCGACACACGCCTTGCGCTGTCGCCAGAGGGCACCCGCGACCTGCAACAGGCCATGACCCGcgccgc          >  6:352891/1‑149 (MQ=255)                                                                                                                                        gATCGATTCGATCAATGTCTACCGCAGCAAGGCCCAGCCCTGCGGTGGCGAGGCGTCGTTGGAACTGCCGCCGCTCAACAGCGACACACGCCTTGCGCTGTCGCCAGAGGGCACCCGCGACCTGCAACAGGCCATGACCCGCgccgcc         >  2:339129/1‑148 (MQ=255)                                                                                                                                        gATCGATTCGATCAATGTCTACCGCAGCAAGGCCCAGCCCTGCGGTGGCGAGGCGTCGTTGGAACTGCCGCCGCTCAACAGCGACACACGCCTTGCGCTGTCGCCAGAGGGCACCCGCGACCTGCAACAGGCCATGACCCGCgccgcc         >  2:326664/1‑148 (MQ=255)                                                                                                                                              ttcgatCAATGTCTACCGCAGCAAGGCCCAGCCCTGCGGTGGCGAGGCGGCGTTGGAACTGCCGCCGCTCAACAGCGACACACGCCTTGCGCTGACGCCAGAGGGCACCCGCGACCTGCACCAGGCCATGACCCGCCCCGCCTACCCCa  >  5:71898/1‑149 (MQ=255)                                                                                                                                              |                                                                                                                                                     GATGAGCAGCGCACACCACGACGCTGCCAGTCATAGGGAGATCCACATGCGCGTCCTGACATCCGTTGTCGCCTTGTCACTGGGCCTGGTCGTTTCGACCGGGGCGCTGGCAGCCACCGGCGAAGAGGCGCAGTTGATCGATTCGATCAATGTCTACCGCAGCAAGGCCCAGCCCTGCGGTGGCGAGGCGTCGTTGGAACTGCCGCCGCTCAACAGCGACACACGCCTTGCGCTGTCGCCAGAGGGCACCCGCGACCTGCAACAGGCCATGACCCGCGCCGCCTACCCCA  >  NC_002947/4041867‑4042156

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