Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,875,844	G→A	6.9%	G109G (GGC→GGT)	ilvA‑II ←	threonine deaminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,875,844	0	G	A	6.9%	213.2 / 3.6	72	G109G (GGC→GGT)	ilvA‑II	threonine deaminase
Reads supporting (aligned to +/- strand):  ref base G (34/33);  new base A (3/2);  total (37/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.78e-01

CGGTGCCCTGCCCGGCGATGGTGTGCGGATCGTCGTAAGGGTGAACGTAGACGAAGCCCTTTTCATCGACCAGTTTCAGCGAATAGGCCAGTGCCTCAGGGAACGAGTCGCCATGCAGCACCACCTTGCCACCGCGCGAGCGCACGCCTTCGACCTTGATCTCCGGAGTGGTCTTGGGCATCACGATGGTGGCCTTGATGCCCATTTCACGCGCGGCCAGGGCCAGGCCTTGGGCGTGGTTGCCGGCAGAGGCGGTCACCACCCCGCGGGCCAGTTCTTCCGGGCTCAGCTGG  >  NC_002947/5875699‑5875991                                                                                                                                                  |                                                                                                                                                    cGGTGCCCTGCCCGGCGATGGTGTGCGGATCGTCGTAAGGGTGAACGTAGACGAAGCCCTTTTCATCGACCAGTTTCAGCGAATAGGCCAGTGCCTCAGGGAACGAGTCGCCATGCAGCACCACCTTGCCACCGCGCGAGCGCACGCCt                                                                                                                                                  >  2:47045/1‑149 (MQ=255)  ggTGCCCTGCCCGGCGATGGTGTGCGGATCGTCGTAAGGGTGAACGTAGACGAAGCCCTTTTCATCGACCAGTTTCAGCGAATAGGCCAGTGCCTCAGGGAACGAGTCGCCATGCAGCACCACCTTGCCACCGCGCGAGCGCACGCCtt                                                                                                                                                 <  4:110462/149‑1 (MQ=255)      ccctgcccGGCGATGGTGTGCGGATCGTCGTAAGGGTGAACGGAGAAGAAGCCCTTTTCATCGACCAGTTTCAGCGAATAGGCCAGTGCCTCAGGGAACGAGTCGCCATGCAGCACCACCTTGCCACCGCGCGAGCGCACGCCTTCGaa                                                                                                                                             <  6:69312/149‑2 (MQ=255)       cctgcccGGCGATGGTGTGCGGATCGTCGTAAGGGTGAACGTAGACGAAGCCCTTTTCATCGACCAGTTTCAGCGAATAGGCCAGTGCCTCAGGGAACGAGTCGCCATGCAGCACCACCTTGCCACCGCGCGAGCGCACGCCTTCGAcc                                                                                                                                            <  5:63298/149‑1 (MQ=255)        ctgcccGGCGATGGTGTGCGGATCGTCGTAAGGGTGAACGTAGACGAAGCCCTTTTCATCGACCAGTTTCAGCGAATAGGCCAGTGCCTCAGGGAACGAGTCGCCATGCAGCACCACCTTGCCACCGCGCGAGCGCACGCCTTCGACCt                                                                                                                                           >  5:340013/1‑149 (MQ=255)          gcccGGCGATGGTGTGCGGATCGTCGTAAGGGTGAACGTTGACGAAGCCGTTTTCATCGACCAGTTTCAGCGAATAGGCCAGTGCCTCAGGGAACGAGTCGGCATGCAGCACCACCTTGCCCCCGCGGGAGCGCACGCCGTCGTCCtt                                                                                                                                          >  6:6321/1‑148 (MQ=255)            ccGGCGATGGTGTGCGGATCGTCGTAAGGGTGAACGTAGACGAAGCCCTTTTCATCGACCAGTTTCAGCGCATAGGCCAGTGCCTCAGGGAACGAGTCGCCATGCAGCACCACCTTGCCACCGCGCGAGCGCACGCCTTCGACCTTGAt                                                                                                                                       >  4:115910/1‑149 (MQ=255)                 gATGGTGTGCGGATCGTCGTAAGGGTGAACGTAGACGAAGCCCTTTTCATCGACCAGTTTCAGCGAATAGGCCAGTGCCTCAGGGAACGAGTCGCCATGCAGCACCACCTTGCCACCGCGCGAGCGCACGCCTTCGACCTTGATCTCCg                                                                                                                                  <  6:191726/149‑1 (MQ=255)                    ggTGTGCGGATCGTCGTAAGGGTGAACGTAGACGAAGCCCTTTTCATCGACCAGTTTCAGCGAATAGGCCAGTGCCTCAGGGAACGAGTCGCCATGCAGCACCACCTTGCCACCGCGCGAGCGCACGCCTTCGACCTTGATCTCCGGa                                                                                                                                <  3:316559/148‑1 (MQ=255)                     gtgtGCGGATCGTCGTAAGGGTGAACGTAGACGAAGCCCTTTTCATCGACCAGTTTCAGCGAATAGGCCAGTGCCTCAGGGAACGAGTCGCCATGCAGCACCACCTTGCCACCGCGCGAGCGCACGCCTTCGACCTTGATCTCCGGAGt                                                                                                                              <  1:65299/149‑1 (MQ=255)                       gtgCGGATCGTCGTAAGGGTGAACGTAGACGAAGCCCTTTTCATCGACCAGTTTCAGCGAATAGGCCAGTGCCTCAGGGAACGAGTCGCCATGCAGCACCACCTTGCCACCGCGCGAGCGCACGCCTTCGACCTTGATCTCCGGAGTgg                                                                                                                            >  1:334066/1‑149 (MQ=255)                           ggATCGTCGTAAGGGTGAACGTAGACGAAGCCCTTTTCATCGACCAGTTTCAGCGAATAGGCCAGTGCCTCAGGGAACGAGTCGCCATGCAGCACCACCTTGCCACCGCGCGAGCGCACGCCTTCGACCTTGATCTCCGGAGTGGTCtt                                                                                                                        <  7:24374/149‑1 (MQ=255)                             aTCGTCGTAAGGGTGAACGTAGACGAAGCCCTTTTCATCGACCAGTTTCAGCGAATAGGCCAGTGCCTCAGGGAACGAGTCGCCATGCAGCACCACCTTGCCACCGCGCGAGCGCACGCCTTCGACCTTGATCTCCGGAGTGGTCTTgg                                                                                                                      >  5:36737/1‑149 (MQ=255)                             aTCGTCGTAAGGGTGAACGTAGACGAAGCCCTTTTCATCGACCAGTTGCAGCGAATAGGCCAGTGCCTCAGGGAACGAGGCGCCATGCAGCACCACCTTGCCACCGCGCGAGCGCACGCCTTCGACCTTGATCTCCGGAGTGGTCTTgg                                                                                                                      <  2:177292/149‑1 (MQ=255)                                   gtAAGGGTGAACGTAGACGAAGCCCTTTTCATCGACCAGTTTCAGCGAATAGGCCAGTGCCTCAGGGAACGAGTCGCCATGCAGCACCACCTTGCCACCGCGCGAGCGCACGCCTTCGACCTTGATCTCCGGAGTGGTCTTGGGCATCa                                                                                                                <  1:47045/149‑1 (MQ=255)                                    tAAGGGTGAACGTAGACGAAGCCCTTTTCATCGCCCAGTTTCAGCGAATAGGCCAGTGCCTCAGGGAACGAGTCGCCATGCAGCACCACCTTGCCACCGCGCGAGCGCACGCCTTCGACCTTGATCTCCGGAGTGGTCTTGGGCATCAc                                                                                                               >  5:112783/1‑149 (MQ=255)                                    tAAGGGTGAACGTAGACGAAGCCCTTTTCATCGACCAGTTTCAGCGAATAGGCCAGTGCCTCAGGGAACGAGTCGCCATGCAGCACCACCTTGCCACCGCGCGAGCGCACACCTTCGACCTTGATCTCCGGAGTGGTCTTGGGCATCAc                                                                                                               >  3:364009/1‑149 (MQ=255)                                    tAAGGGTGAACGTAGACGAAGCCCTTTTCATCGACCAGTTTCAGCGAATAGGCCAGTGCCTCAGGGAACGAGTCGCCATGCAGCACCACCTTGCCACCGCGCGAGCGCACACCTTCGACCTTGATCTCCGGAGTGGTCTTGGGCATCAc                                                                                                               >  3:369740/1‑149 (MQ=255)                                         gTGAACGTAGACGAAGCCCTTTTCATCGACCAGTTTCAGCGAATAGGCCAGTGCCTCAGGGAACGAGTCGCCATGCAGCACCACCTTGCCACCGCGCGAGCGCACGCCTTCGACCTTGATCTCCGGAGTGGTCTTGGGCATCACGAtgg                                                                                                          >  8:277092/1‑149 (MQ=255)                                              cGTAGACGAAGCCCTTTTCATCGACCAGTTTCAGCGAATAGGCCAGTGCCTCAGGGAACGAGTCGCCATGCAGCACCACCTTGCCACCGCGCGAGCGCACGCCTTCGACCTTGATCTCCGGAGTGGTCTTGGGCATCACGATGGTGGcc                                                                                                     <  2:363454/149‑1 (MQ=255)                                              cGTAGACGAAGCCCTTTTCATCGACCAGTTTCAGCGAATAGGCCAGTGCCTCAGGGAACGAGTCGCCATGCAGCACCACCTTGCCACCGCGCGAGCGCACGCCTTCGACCTTGATCTCCGGAGTGGTCTTGGGCATCACGATGGTGGcc                                                                                                     >  8:67200/1‑149 (MQ=255)                                              cGTAGACGAAGCCCTTTTCATCGACCAGTTTCAGCGAATAGGCCAGTGCCTCAGGGAACGAGTCGCCATGCAGCACCACCTTGCCACCGCGCGAGCGCACGCCTTCGACCTTGATCTCCGGAGTGGTCTTGGGCATCACGATGGTGGcc                                                                                                     <  6:80194/149‑1 (MQ=255)                                              cGTAGACGAAGCCCTTTTCATCGACCAGTTTCAGCGAATAGGCCAGTGCCTCAGGGAACGAGTCGCCATGCAGCACCACCTTGCCACCGCGCGAGCGCACGCCTTCGACCATGATCTCCGGAGTGGTCTTGGGCATCACGATGGTGGcc                                                                                                     >  7:154764/1‑149 (MQ=255)                                                     gAAGCCCTTTTCATCGACCAGTTTCAGCGAATAGGCCAGTGCCTCAGGGAACGAGTCGCCATGCAGCACCACCTTGCCACCGCGCGAGCGCACGCCTTCGACCTTGATCTCCGGAGTGGTCTTGGGCATCACGATGGTGGCCTTGATGc                                                                                              >  3:166986/1‑149 (MQ=255)                                                     gAAGCCCTTTTCATCGACCAGTTTCAGCGAATAGGCCAGTGCCTCAGGGAAAGAGTCGCCATGCAGCACCACCTTGCCACCGCGCGAGCGCACGCCTTCGACCTTGATCTCCGGAGTGGTCTTGGGAATCACGATTGTGGCCTTGATGc                                                                                              >  6:273229/1‑149 (MQ=255)                                                           cTTTTCATCGACCAGTTTCAGCGAATAGGCCAGTGCCTCAGGGAACGAGTCGCCATGCAGCACCACCTTGCCACCGCGCGAGCGCACGCCTTCGACCTTGATCTCCGGAGTGGTCTTGGGCATCACGATGGTGGCCTTGATGCCCAttt                                                                                        <  7:232481/149‑1 (MQ=255)                                                             tttCATCGACCAGTTTCAGCGAATAGGCAAGTGCCTCAGGGAACGAGTCGCCATGCAGCACCACATTTCCACCGCGCCAGCGCACGCCATCGACCTTGATCTCCGGAGTGGACTTGGGCATCACCATGGTGGCCTAGATCCCCATTTCa                                                                                      >  5:9722/1‑149 (MQ=255)                                                                  tCGACCAGTTTCAGCGAATAGGCCAGTGCCTCAGGGAACGAGTCGCCATGCAGCACCACCTTGCCACCGCGCGAGCGCACGCCTTCGACCTTGATCTCCGGAGTGGTCTTGGGCATCACGATGGTGGCCTTGATGCCCATTTCAcgcgc                                                                                 >  1:94620/1‑149 (MQ=255)                                                                       cAGTTTCAGCGAATAGGCCAGTGCCTCAGGGAACGAGTCGCCATGCAGCACCACCTTGCCACCGCGCGAGCGCACACCTTCGACCTTGATCTCCGGAGTGGTCTTGGGCATCACGATGGTGGCCTTGATGCCCATTTCACGCGCggcca                                                                            <  4:369740/149‑1 (MQ=255)                                                                       cAGTTTCAGCGAATAGGCCAGTGCCTCAGGGAACGAGTCGCCATGCAGCACCACCTTGCCACCGCGCGAGCGCACACCTTCGACCTTGATCTCCGGAGTGGTCTTGGGCATCACGATGGTGGCCTTGATGCCCATTTCACGCGCggcca                                                                            <  4:364009/149‑1 (MQ=255)                                                                         gTTTCAGCGAATAGGCCAGTGCCTCAGGGAACGAGTCGCCATGCAGCACCACCTTGCCACCGCGCGAGCGCACGCCTTCGACCTTGATCTCCGGAGTGGTCTTGGGCATCACGATGGTGGCCTTGATGCCCATTTCACGCGCggccagg                                                                          <  6:112783/149‑1 (MQ=255)                                                                         gTTTCAGCGAATAGGCCAGTGCCTCAGGGAACGAGTCGCCATGCAGCACCACCTTGCCACCGCGCGAGCGCACGCCTTCGACCTTGATCTCCGGAGTGGTCTTGGGCATCACGATGGTGGCCTTGATGCCCATTTCACGCGCggccagg                                                                          <  6:36737/149‑1 (MQ=255)                                                                              aGCGAATTGGCCAGTCCCTGAGTGAAGGAGTAGCCATGGAGCACCCCCTTGCCCCGGCGCGAGGGCACGCCATCGCCCTTGATCTCCGGACTGGTCTTTGGCATCTCGATGGTGGCCTTGATGCCCATTTCACGCGCggccagggcca                                                                      <  6:9722/148‑1 (MQ=255)                                                                              aGCGAATAGGCCAGTGCCTCAGGGAACGAGTCGCCATGCAGCACCACCTTGCCACCGCGCGAGCGCACGCCTTCGACCTTGATCTCCGGAGTGGTCTTGGGCATCACGATGGTGGCCTTGATGCCCATTTCACGCGCggccagggccag                                                                     >  5:383730/1‑149 (MQ=255)                                                                                cGAATAGGCCAGTGCCTCAGGGAACGAGTCGCCATGCAGCACCACCTTGCCACCGCGCGAGCGCACGCCTTCGACCTTGATCTCCGGAGTGGTCTTGGGCATCACGATGGTGGCCTTGATGCCCATTTCACGCGCGGCCAGGGCCAGGc                                                                   <  6:119430/149‑1 (MQ=255)                                                                                  aaTAGGCCAGTGCCTCAGGGAACGAGTCGCCATGCAGCACCACCTTGCCACCGCGCGAGCGCACGCCTTCGACCTTGATCTCCGGAGTGGTCTTGGGCATCACGATGGTGGCCTTGATGCCCATTTCACGCGCGGCCAGGGCCAGGCCt                                                                 <  4:166986/149‑1 (MQ=255)                                                                                    tAGGCCAGTGCCTCAGGGAACGAGTCGCCATGCAGCACCACCTTGCCACCGCGCGAGCGCACGCCTTCGACCTTGATCTCCGGAGTGGTCTTGGGCATCACGATGGTGGCCTTGATGCCCATTTCACGCGCGGCCAGGGCCAGGCCTTg                                                               <  7:277092/149‑1 (MQ=255)                                                                                     aGGCCAGTGCCTCACGGAACGAGTCGCCATGCAGCACCACCTTGCCACCGCGCGAGCGCACGCCATCGACCTTGATCTCCGGAGTCGTCTTGGGCATCACGATGGGGGCCTTGATGCCCATTTCACGCGCGGCCAGGGCCAGGCCTTgg                                                              <  5:6321/149‑1 (MQ=255)                                                                                       gCCAGTGCCTCAGGGAACGAGTCGCCATGCAGCACCACCTTGCCACCGCGCGAGCGCACGCCTTCGACCTTGATCTCCGGAGTGGTc                                                                                                                          >  4:201873/1‑87 (MQ=255)                                                                                       gCCAGTGCCTCAGGGAACGAGTCGCCATGCAGCACCACCTTGCCACCGCGCGAGCGCACGCCTTCGACCTTGATCTCCGGAGTGGTc                                                                                                                          <  3:201873/87‑1 (MQ=255)                                                                                        ccAGTGCCTCAGGGAACGAGTCGCCATGCAGCACCACCTTGCCACCGCGCGAGCGCACGCCTTCGACCTTGATCTCCGGAGTGGTCTTGGGCATCACGATGGTGGCCTTGATGCCCATTTCACGCGCGGCCAGGGCCAGGCCTTGGGCg                                                           <  4:259859/149‑1 (MQ=255)                                                                                         cAGTGCCTCAGGGAACGAGTCGCCATGCAGCACCACCTTGCCACCGCGCGAGCGCACGCCTTCGACCTTGATCTCCGGAGTGGTCTTGGGCATCACGATGGTGGCCTTGATGCCCATTTCACGCGCGGCCAGGGCCAGGCCTTGGGCGt                                                          >  3:283308/1‑149 (MQ=255)                                                                                          aGTGCCTCAGGGAACGAGTCGCCATGCAGCACCACCTTGCCACCGCGCGAGCGCACGCCTTCGACCTTGATCTCCGGAGTGGTCTTGGGCATCACGATGGTGGCCTTGATGCCCATTTCACGCGCGGCCAGGGCCAGGCCTTGGGCGTg                                                         <  3:115910/149‑1 (MQ=255)                                                                                            tGCCTCAGGGAACGAGTCGCCATGCAGCACCACCTTGCCACCGCGCGAGCGCACGCCTTCGACCTTGATCTCCGGAGTGGTCTTGGGCATCACGATGGTGGCCTTGATGCCCATTTCACGCGCGGCCAGGGCCAGGCCTTGGGCGTGGt                                                       >  8:336476/1‑149 (MQ=255)                                                                                            ctcatCAGGGAACGAGTCGCCATGCAGCACCACCTTGCCACCGCGCGAGCGCACGCCTTCGACCTTGATCTCCGGAGTGGTCTTGGGCATCACGATGGTGGCCTTGATGCCCATTTCACGCGCGGCCAGGGCCAGGCCTTGGGCGTGGt                                                       >  6:18250/5‑149 (MQ=255)                                                                                               cTCAGGGAACGAGTCGCCATGCAGCACCACCTTGCCACCGCGCGAGCGCACGCCTTCGACCTTGATCTCCGGAGTGGTCTTGGGCATCACGATGGTGGCCTTGATGCCCATTTCACGCGCGGCCAGGGCCAGGCCTTGGGCGTGGTTGc                                                    <  7:67200/149‑1 (MQ=255)                                                                                               cTCAGGGAACGAGTCGCCATGCAGCACCACCTTGCCACCGCGCGAGCGCACACCTTCGACCTTGATCTCCGGAGTGGTCTTGGGCATCACGATGGTGGCCTTGATGCCCATTTCACGCGCGGCCAGGGCCAGGCCTTGGGCGTGGTTGc                                                    >  1:98080/1‑149 (MQ=255)                                                                                               cTCAGGGAACGAGTCGCCATGCAGAACCACCTTGCCACCGCGCGAGCGCACGCCTTCGACCTTGATCTCCGGAGAGGTATTGGGCATCACGATGGTGGCCTTGATGCCCATTTAACGAGCGGCCAGGGCAAGGCCTTGGGAGTGGTTGc                                                    >  8:199897/1‑149 (MQ=255)                                                                                                  aGGGAACGAGTCGCCATGCAGCACCACCTTGCCACCGCGCGAGCGCACGCCTTCGACCTTGATCTCCGGAGTGGTCTTGGGCATCACGATGGTGGCCTTGATGCCCATTTCACGCGCGGCCAGGGCCAGGCCTTGGGCGTGGTTGCCgg                                                 <  4:283308/149‑1 (MQ=255)                                                                                                  aGGGAACGAGTCGCCATGCAGCACCACCTTGCCACCGCGCGAGCGCACGCCTTCGACCTTGATCTCCGGAGTGGTCTTGGGCATCACGATGGTGGCCTTGATGCCCATTTCACGCGCGGCCAGGGCCAGGCCTTGGGCGTGGTTGCCgg                                                 >  6:339401/1‑149 (MQ=255)                                                                                                     gAACGAGTCGCCATGCAGCACCACCTTGCCACCGCGCGAGCGCACGCCTTCGACCTTGATCTCCGGAGTGGTCTTGGGCATCACGATGGTGGCCTTGATGCCCATTTCACGCGCGGCCAGGGCCAGGCCTTGGGCGTGGTTGCCgg                                                 >  6:30542/1‑146 (MQ=255)                                                                                                        cGAGTCGCCATGCAGCACCACCTTGCCACCGCGCGAGCGCACGCCTTCGACCTTGATCTCCGGAGTGGTCTTGGGCATCACGATGGTGGCCTTGATGCCCATTTCACGCGCGGCCAGGGCCAGGCCTTGGGCGTGGTTGCCGGCAGAgg                                           <  7:143588/149‑1 (MQ=255)                                                                                                         gAGTCGCCATGCAGCACCACCTTGCCACCGCGCGAGCGCACGCCTTCGACCTTGATCTCCGGAGTGGTCTTGGGCATCACGATGGTGGCCTTGATGCCCATTTCACGCGCGGCCAGGGCCAGGCCTTGGGCGTGGTTGCCGGCAGAgg                                           >  8:143588/1‑148 (MQ=255)                                                                                                           gTCGCCATGCAGCACCACCTTGCCACCGCGCGAGCGCACGCCTTCGACCTTGATCTCCGGAGTGGTCTTGGGCATCACGATGGTGGCCTTGATGCCCATTTCACGCGCGGCCAGGGCCAGGCCTTGGGCGTGGTTGCCGGCAGAGGCgg                                        >  1:190101/1‑149 (MQ=255)                                                                                                           gTCGCCATGCAGCACCACCTTGCCACCGCGCGAGCGCACGCCTTCGACCTTGATCTCCGGAGTGGTCTTGGGCATCACGATGGTGGCCTTGATGCCCATTTCACGCGCGGCCAGGGCCAGGCCTTGGGCGGGGTTGCCGGCAGAGGCgg                                        <  1:258240/149‑1 (MQ=255)                                                                                                            tCGCCATGCAGCACCACCTTGCCACCGCGCGAGCGCACGCCTTCGACCTTGATCTCCGGAGTGGTCTTGGGCATCACGATGGTGGCCTTGATGCCCATTTCACGCGCGGCCAGGGCCAGGCCTTGGGCGTGGTTGCCGGCAGAGGCGGt                                       <  3:243867/149‑1 (MQ=255)                                                                                                              gCCATGCAGCACCACCTTGCCACCGCGCGAGCGCACGCCTTCGACCTTGATCTCCGGAGTGGTCTTGGGCATCACGATGGTGGCCTTGATGCCCATTTCACGCGCGGGCAGGGCCAGGCCTTGGGCGTGGTTGCCGGCAGAGGCGGTca                                     >  2:285149/1‑149 (MQ=255)                                                                                                              gCCATGCAGCACCACCTTGCCACCGCGCGAGCGCACGCCTTCGACCTTGATCTCCGGAGTGGTCTTGGGCATCACGATGGTGGCCTTGATGCCCATTTCACGCGCGGCCAGGGCCAGGCCTTGGGCGTGGTTGCCGGCAGAGGCGGTca                                     >  6:129688/1‑149 (MQ=255)                                                                                                              gCCATGCAGCACCACCTTGCCACCGCGCGAGCGCACGCCTTCGACCTTGATCTCCGGAGTGGTCTTGGGCATCACGATGGTGGCCTTGATGCCCATTTCACGCGCGGCCAGGGCCAGGCCTTGGGCGTGGTTGCCGGCAGAGGCGGTca                                     <  5:18250/149‑1 (MQ=255)                                                                                                                  tGCAGCACCACCTTGTCACCGCGCGAGCGCACGCCTTCGACCTTGATCTCCGGAGTGGTCTTGGGCATCACGATGGTGGCCTTGATGCCCATTTCACGCGCGGCCAGGGCCAGGCCTTGGGCGTGGTTGCCGGCag                                              >  7:244662/1‑136 (MQ=255)                                                                                                                  tGCAGCACCACCTTGTCACCGCGCGAGCGCACGCCTTCGACCTTGATCTCCGGAGTGGTCTTGGGCATCACGATGGTGGCCTTGATGCCCATTTCACGCGCGGCCAGGGCCAGGCCTTGGGCGTGGTTGCCGGCag                                              <  8:244662/136‑1 (MQ=255)                                                                                                                    cagcaCCACCTTGCCACCGCGCGAGCGCACGCCTTCGACCTTGATCTCCGGAGTGGTCTTGGGCATCACGATGGTGGCCTTGATGCCCATTTCACGCGCGGCCAGGGCCAGGCCTTGGGCGTGGTTGCCGGCAGAGGCGGTCACCAccc                               >  4:347111/1‑149 (MQ=255)                                                                                                                         ccaccTTGCCACCGCGCGAGCGCACGCCTTCGACCTTGATCTCCGGAGTGGTCTTGGGCATCACGATGGTGGCCTTGATGCCCATTTCACGCGCGGCCAGGGCCAGGCCTTGGGCGTGGTTGCCGGCAGAGGCGGTCACCACCCCGCgg                          >  2:330637/1‑149 (MQ=255)                                                                                                                             cTTGCCACCGCGCGAGCGCACGCCTTCGACCTTGATCTCCGGAGTGGTCTTGGGCATCACGATGGTGGCCTTGATGCCCATTTCACGCGCGGCCAGGGCCAGGCCTTGGGCGTGGTTGCCGGCAGAGGCGGTCACCACCCCGCGGGc                        >  1:243849/1‑147 (MQ=255)                                                                                                                                  cACCGCGCGAGCGCACGCCTTCGACCTTGATCTCCGGAGTGGTCTTGGGCATCACGATGGTGGCCTTGATGCCCATTTCACGCGCGGCCAGGGCCAGGCCTTGGGCGTGGTTGCCGGCAGAGGCGGTCACCACCCCGCGGGCCAGttct                 >  2:78229/1‑149 (MQ=255)                                                                                                                                  cACCGCGCGAGCGCACGCCTTCGACCTTGATCTCCGGAGTGGTCTTGGGCATCACGATGGTGGCCTTGATGCCCATTTCACGCGCGGCCAGGGCCAGGCCTTGGGCGTGGTTGCCGGCAGAGGCGGTCACCACCCCGCGGGCCAGttct                 >  8:82665/1‑149 (MQ=255)                                                                                                                                   aCCGCGCGAGCGCACGCCTTCGACCTTGATCTCCGGAGTGGTCTTGGGCATCACGATGGTGGCCTTGATGCCCATTTCACGCGCGGCCAGGGCCAGGCCTTGGGCGTGGTTGCCGGCAGAGGCGGTCACCACCCCGCGGGCCAGttctt                <  6:383730/149‑1 (MQ=255)                                                                                                                                      gcgcgAGCGCACGCCTTCGACCTTGATCTCCGGAGTGGTCTTGGGCATCACGATGGTGGCCTTGATGCCCATTTCACGCGCGGCCAGGGCCAGGCCTTGGGCGTGGTTGCCGGCAGAGGCGGTCACCACCCCGCGGGCCAGTTCTTCCg             <  6:340013/149‑1 (MQ=255)                                                                                                                                          gAGCGCACGCCTTCGACCTTGATCTCCGGAGTGGTCTAGGGCATCACGATGGTGGCCTTGATGCCCATTTCACGCGCGGCCAGGGCCAGGCCTTGGGCGTGGTTGCCGGCAGAGGCGGTCACCACCCCGCGGGCCAGTTCTTCCGGGCt         >  2:56277/1‑149 (MQ=255)                                                                                                                                              gcACGCCTTCGACCTTGATCTCCGGAGTGGTCTTGGGCATCACGATGGTGGCCTTGATGCCCATTTCACGCGCGGCCAGGGCCAGGCCTTGGGCGTGGTTGCCGGCAGAGGCGGTCACCACCCCGCGGGCCAGTTCTTCCGGGCTCAGc     >  6:182687/1‑149 (MQ=255)                                                                                                                                                 cGCCTTCTAGCTTGATCGCCGGAGTGGTCTTGGTCATCACGATGGTGGCCTTGATGCCCATTTCCCGCGCGGCCAGGGCCAGGCCTTGGGCGTGGTTGCCGGCAGAGGCGGTCACCACCCGGCGGGCCAGTTCTTCCGGGCTCAGCTgg  <  8:154764/149‑1 (MQ=255)                                                                                                                                                 cGCCTTCGACCTTGATCTCCGGAGTGGTCTTGGGCATCACGATGGTGGCCTTGATGCCCATTTCACGCGCGGCCAGGGCCAGGCCTTGGGCGTGGTTGCCGGCAGAGGCGGTCACCACCCCGCGGGCCAGTTCTTCCGGGCTCAGCTgg  <  2:334066/149‑1 (MQ=255)                                                                                                                                                  gCCTTCGACCTTGATCTCCGGAGTGGTCTTGGGCATCACGATGGTGGCCTTGATGCCCATTTCACGCGCGGCCAGGGCCAGGCCTTGGGCGTGGTTGCCGGCAGAGGCGGTc                                      >  4:103358/1‑112 (MQ=255)                                                                                                                                                  gCCTTCGACCTTGATCTCCGGAGTGGTCTTGGGCATCACGATGGTGGCCTTGATGCCCATTTCACGCGCGGCCAGGGCCAGGCCTTGGGCGTGGTTGCCGGCAGAGGCGGTc                                      <  3:103358/112‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                    CGGTGCCCTGCCCGGCGATGGTGTGCGGATCGTCGTAAGGGTGAACGTAGACGAAGCCCTTTTCATCGACCAGTTTCAGCGAATAGGCCAGTGCCTCAGGGAACGAGTCGCCATGCAGCACCACCTTGCCACCGCGCGAGCGCACGCCTTCGACCTTGATCTCCGGAGTGGTCTTGGGCATCACGATGGTGGCCTTGATGCCCATTTCACGCGCGGCCAGGGCCAGGCCTTGGGCGTGGTTGCCGGCAGAGGCGGTCACCACCCCGCGGGCCAGTTCTTCCGGGCTCAGCTGG  >  NC_002947/5875699‑5875991

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