Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,704,637	T→A	11.8%	M1* (TTG→TAG) †‡	PP_1499 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,704,637	0	T	A	11.8%	170.3 / 13.3	68	M1* (TTG→TAG) †‡	PP_1499	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (28/32);  new base A (5/3);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.68e-01

GAGGACAACGGCGACGGCGTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAAC‑CATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGTGAAGTGGGAC  >  NC_002947/1704491‑1704779                                                                                                                                                   |                                                                                                                                                gaggaCAACGGCGACGGCGTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGc                                                                                                                                               <  2:104872/149‑1 (MQ=255) gaggaCAACGGCGACGGCGTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTAGCGGGCTTgg                                                                                                                                               >  7:304600/1‑148 (MQ=255)    gaCAACGGCGACGGCGTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCcg                                                                                                                                            >  7:178027/1‑149 (MQ=255)       aaCGTCGACGGCGTGCTGATCGACAGACAGGCGGCTCGGGTGCTGGCCTACCTGGACGACCCGGCCAAATCAGACCTGTTCGTCGAAGTGGCGCGTTTCCAGCGCTCGGGAAAATGAACGTGTTGGCGTCTTCTCGGGCTTTCAcgcg                                                                                                                                          >  2:231952/1‑148 (MQ=255)       aaCGGCGACGGCGTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCcgcg                                                                                                                                          >  6:88776/1‑148 (MQ=255)        acggcgacggcgTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGaa                                                                                                                                        <  4:86016/149‑1 (MQ=255)        acggcgacggcgTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCAAGCCCGCGaa                                                                                                                                        <  4:243819/149‑1 (MQ=255)        acggcgacggcgTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCCTCTTCGCGGGCAAGCCCGCGaa                                                                                                                                        <  7:50521/149‑1 (MQ=255)          ggcgacggcgTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGaa                                                                                                                                        >  2:58178/1‑147 (MQ=255)          ggcgacggcgTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAaga                                                                                                                                      <  4:128427/149‑1 (MQ=255)             gacggcgTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGc                                                                                                                                   >  7:256412/1‑149 (MQ=255)                   gTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTAc                                                                                                                             <  5:170296/149‑1 (MQ=255)                        gATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGc                                                                                                                                               <  1:118983/126‑1 (MQ=255)                        gATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGc                                                                                                                                               >  2:118983/1‑126 (MQ=255)                             aCAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCgtg                                                                                                                   >  1:24286/1‑149 (MQ=255)                             aCAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCgtg                                                                                                                   <  8:178027/149‑1 (MQ=255)                             aCAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTAGCGGGCTTGACCGAGAAGAGGCCGGTACACGTATCgtg                                                                                                                   >  6:277579/1‑149 (MQ=255)                                gCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTGGCGGGCGTGGCCGCGCAGGGGCCGGGTCAGGGTTGgtggt                                                                                                                 >  3:192052/1‑146 (MQ=255)                                ccccGGCGGCGGGGGTGCTGGCCTACCTGGACCACCCGGCAAACTCCGACCTGTTCGGCGATGGGGCGCGGTTCCAGCGCTCGGGAAAAGGCCCGTGTTGGGGTCTTCGGGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTa                                                                                                                 <  4:192052/147‑1 (MQ=255)                                  ccGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTa                                                                                                              <  2:292084/149‑1 (MQ=255)                                      cggcgCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGaaa                                                                                                          <  2:27865/149‑1 (MQ=255)                                       ggcgCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAg                                                                                                         <  7:94804/149‑1 (MQ=255)                                         cgcgGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGa                                                                                                       <  5:175845/149‑1 (MQ=255)                                         cgcgGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGa                                                                                                       <  2:197393/149‑1 (MQ=255)                                                gCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCg                                                                                                >  3:314312/1‑149 (MQ=255)                                                   ggCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTg                                                                                             >  3:149304/1‑149 (MQ=255)                                                    gCCTACCTGGACGACCCGTCCAACTCCGACGTGTTCGTCGAAGTGGCGCGGTTTCAGCGCTCGGGCAAATTACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCTAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTg                                                                                             <  8:154798/148‑1 (MQ=255)                                                    gCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTGGGGAAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGa                                                                                                       <  5:245270/138‑1 (MQ=255)                                                    gCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGGTCCCGGTgg                                                                                            >  8:104593/1‑149 (MQ=255)                                                      cTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGAc                                                                                          >  4:215984/1‑149 (MQ=255)                                                      cTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGAc                                                                                          <  4:7772/149‑1 (MQ=255)                                                       tACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCATCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGAc                                                                                          >  7:90499/1‑148 (MQ=255)                                                         ccTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGAGGAGCGGCCCGTTCCCGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGGCgc                                                                                        >  3:346118/1‑147 (MQ=255)                                                           tGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGACAg                                                                                                                           >  5:220481/1‑111 (MQ=255)                                                              acgacCGGGCCAACTCGGAGAGGTTCGTCTAAGTGGCGCGGTTCCAGCGCTCGGTCAAATTACGGTGTTGGGTTCTTCGCGGGGTTGCCCGCTATGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCa                                                                                  <  2:159190/149‑1 (MQ=255)                                                              acgacCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCa                                                                                  <  3:141849/149‑1 (MQ=255)                                                                gacCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGc                                                                                <  3:215984/149‑1 (MQ=255)                                                                      gCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAAC‑Ca                                                                         <  8:256412/149‑1 (MQ=255)                                                                         aaCTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAAC‑CATcc                                                                      <  1:58178/149‑1 (MQ=255)                                                                              cGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCaaa                                                                             >  7:345433/1‑138 (MQ=255)                                                                              cGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGACCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCaaa                                                                             <  8:345433/138‑1 (MQ=255)                                                                                aCCTGTTGGTCGATGGGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAAC‑CATCCGCGCCCt                                                               <  8:304600/149‑1 (MQ=255)                                                                                aCCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAAC‑CATCCGCGCCCt                                                               <  5:88776/149‑1 (MQ=255)                                                                                        gtcgAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAAC‑CATACGCGCCCTGTCCGAcc                                                       >  1:8121/1‑149 (MQ=255)                                                                                                   gcgGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAAC‑CATCCGCGCCCTGTCCGACCGCATCGTcgcg                                            >  6:364985/1‑149 (MQ=255)                                                                                                   gcgGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCATCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAAC‑CATCCGCGCCCTGTCCGACCGCATCGTcgcg                                            <  8:90499/149‑1 (MQ=255)                                                                                                     ggTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACGCATCCGCGCCCTGTCCGACCGCATCGTCGCgg                                           <  1:231952/149‑1 (MQ=255)                                                                                                           aGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGATGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAAC‑CATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGAc                                    <  4:149304/149‑1 (MQ=255)                                                                                                           aGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAAC‑CATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGAc                                    >  1:307030/1‑149 (MQ=255)                                                                                                            gcgcTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAAC‑CATCCGCGCCCTGTCCGACCGCATCGTCGCgg                                           <  4:41359/141‑1 (MQ=255)                                                                                                            gcgcTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAAC‑CATCCGCGCCCTGTCCGACCGCATCGTCGCgg                                           >  3:41359/1‑141 (MQ=255)                                                                                                               cTCGGGCAAATGACCGTGTTGGCGTCTTCCCCGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAAC‑CATCCCCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGcc                                 >  8:31367/1‑148 (MQ=255)                                                                                                                                  cggCCTCTTCGCGGGCAAGCCCGCGAAGACGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAAC‑CATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCgg             >  2:309182/2‑149 (MQ=255)                                                                                                                                    gcctCTTCGCGGGCAAGCCCGCGAAGACGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAAC‑CATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGt            >  4:150229/4‑148 (MQ=255)                                                                                                                                    gCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAAC‑CATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGt            <  4:314312/148‑1 (MQ=255)                                                                                                                                    gCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAAC‑CATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGt            <  2:24286/148‑1 (MQ=255)                                                                                                                                      gTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAAC‑CATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGTGaa         <  5:364985/149‑1 (MQ=255)                                                                                                                                       tCTTCGCGGGCAAGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAAC‑CATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGTGAAg        >  4:151509/1‑149 (MQ=255)                                                                                                                                       tCTTCGCGGGCAAGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAAC‑CATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGTGAAg        >  3:215392/1‑149 (MQ=255)                                                                                                                                        cTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAAC‑CATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGTGAAGt       >  6:326453/1‑149 (MQ=255)                                                                                                                                        cTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGGGGACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAAC‑CCTCCGCGCCCTGTCCGACCGCATCGGCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGTGAAGt       <  7:255985/149‑1 (MQ=255)                                                                                                                                         ttCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAAC‑CATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGTGAAGTg      >  1:144352/1‑149 (MQ=255)                                                                                                                                         ttCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAAC‑CATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGTGAAGTg      >  2:234079/1‑149 (MQ=255)                                                                                                                                         ttCGCGGGCAAGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAAC‑CATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCtt                    >  3:290615/1‑135 (MQ=255)                                                                                                                                         ttCGCGGGCAAGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAAC‑CATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCtt                    <  4:290615/135‑1 (MQ=255)                                                                                                                                            gcgGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAAC‑CATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGTGAAGTGGga   >  8:11628/1‑149 (MQ=255)                                                                                                                                             cgGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAAC‑CATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGTGAAGTGGgac  <  4:346118/149‑1 (MQ=255)                                                                                                                                              gggCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAAC‑CATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGTGAAGTGGgac  >  2:108397/1‑148 (MQ=255)                                                                                                                                                   |                                                                                                                                                GAGGACAACGGCGACGGCGTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAAC‑CATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGTGAAGTGGGAC  >  NC_002947/1704491‑1704779

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