Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,548,024	A→G	5.9%	L60L (TTG→CTG)	rhtA ←	threonine/homoserine exporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,548,024	0	A	G	5.9%	197.9 / 2.2	67	L60L (TTG→CTG)	rhtA	threonine/homoserine exporter
Reads supporting (aligned to +/- strand):  ref base A (32/31);  new base G (2/2);  total (34/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.26e-01

GTGAATTCCAGGGCGACGGCGATACCGATTGGTACCGTTTGCAGGGCCATATAGAAGAGGAAGTTCATGCCGCCCAGTGCCATGCCGTAGATGACCACATTGCGCAACGTGTTGGCGGTCATGCGTACCCGCCACGGGCGCAATATAAGCAGCATG‑ATGATGCTGGCGAAGATCAGGCGCAGGGTGGTAGTGCCTTGTGCGCCGACGATCGGGAACATGCTTTTGGCCAGGGAGGCGCCGGACTGAATCGATGCCATGGCGATAAGCAGTAGGCC  >  NC_002947/2547882‑2548156                                                                                                                                               |                                                                                                                                      gTGAATTCCAGGGCGACGGCGATACCGATTGGTACCGTTTGCAGGGCCATATAGAAGAGGAAGTTCATGCCGCCCAGTGCCATGCCGTAGATGACCACATTGCGCAACGTGTTGGCGGTCATGCGTACCCGCCACGGGCGCAATATAag                                                                                                                                 <  4:386820/149‑1 (MQ=255) gTGAATTCCAGGGCGACGGCGATACCGATTGGTACCGTTTGCAGGGCCATATAGAAGAGGAAGTTCATGCCGCCCAGTGCCATGCCGTAGATGACCACATTGCGCAACGTGTTGGCGGTCATGCGTACCCGCCACGGGCGCAATATAag                                                                                                                                 >  2:370973/1‑149 (MQ=255) gTGAATTCCAGGGCGACGGCGATACCGATTGGTACCGTTTGCAGGGCCATATAGAAGAGGAAGTTCATGCCGCCCAGTGCCATGCCGTAGATGACCACATTGCGCAACGTGTTGGCGGTCATGCGTACCCGCCACGGGCGCAATATAag                                                                                                                                 <  6:339218/149‑1 (MQ=255)        ccAGGGCGACGGCGATACCGATTGGTACCGTTTGCAGGGCCATATAGAAGAGGAAGTTCATGCCGCCCAGTGCCATGCCGTAGATGACCACATTGCGCAACGTGTTGGCGGTCATGCGTACCCGCCACGGGCGCAATATAAGCAGCat                                                                                                                           >  5:364230/1‑148 (MQ=255)         cAGGGCGACGGCGATACCGATTGGTACCGTTTGCAGGGCCATATAGAAGAGGAAGTTCATGCCGCCCAGTGCCATGCCGTAGATGACCACATTGCGCAACGTGTTGGCGGTCATGCGTACCCGCCACGGGCGCAATATAAGCAGCatg‑a                                                                                                                        >  5:63520/1‑149 (MQ=255)         cAGGGCGACGGCGATACCGATTGGTACCGTTTGCAGGGCCATATAGAAGAGGAAGTTCATGCCGCCCAGTGCCATGCCGTAGATGACCACATTGCGCAACGTGTTGGCGGTCATGCGTACCCGCCACGGGCGCAATATAAGCAGCatg‑a                                                                                                                        >  1:19302/1‑149 (MQ=255)         cAGGGCGACGGCGATACCGATTGGTACCGTTTGCAGGGCCATATAGAAGAGGAAGTTCATGCCGCCCAGTGCCATGCCGTAGATGACCACATTGCGCAACGTGTTGGCGGTCATGCGTACCCGCCACGGGCGCAATATAAGCAGCatg                                                                                                                          <  1:370973/148‑1 (MQ=255)          aGGGCGACGGCGATACCGATTGGTACCGTTTGCAGGGCCATATAGAAGAGGAAGTTCATGCCGCCCAGTGCCATGCCGTAGATGACCACATTGCGCAACGTGTTGGCGGTCATGCGTACCCGCCACGGGCGCAATATAAGCAGCatg‑at                                                                                                                       >  6:354053/1‑149 (MQ=255)               gACGGCGATACCGATTGGTACCGTTTGCAGGGCCATATAGAAGAGGAAGTTCATGCCGCCCAGTGCCATGCCGTAGATGACCACATTGCGCAACGTGTTGGCGGTCATGCGTACCCGCCACGGGCGCAATATAAGCAGCATG‑ATGATGc                                                                                                                  <  5:354053/149‑1 (MQ=255)                      atacgtaTTGGTACCGTTTTCATGGCCATATAGTAGAGGAAGTTCATGCCGCCCAGTGCCATGCCGTAGATGAACCCATTTCGGAACGTGTTGGCGGTCATGCGTACCCGCCACGGGCGCAATATAAGCAGCATG‑ATGATGCTGGCGaa                                                                                                           <  7:67431/143‑1 (MQ=255)                      aTACCGATTGGTACCGTTTGCAGGGCCATATAGAAGAGGAAGTTCATGCCGCCCAGTGCCATGCCGTAGATGACCACATTGCGCAACGTGTTGGCGGTCATGCGTACCCGCCACGGGCGCAATATAAGCAGCATG‑ATGATGCTGGCGaa                                                                                                           <  6:81076/149‑1 (MQ=255)                                gTACCGTTTGCAGGGCCATATAGAAGAGGAAGTTCATGCCGCCCAGTGCCATGCCGTAGATGACCACATTGCGCAACGTGTTGGCGGTCATGCGTACCCGCCACGGGCGCAATATAAGCAGCATG‑ATGATGCTGGCGAAGATCAGgcgc                                                                                                 >  2:215399/1‑149 (MQ=255)                                gTACCGTTTGCAGGGCCATATAGAAGAGGAAGTTCATGCCGCCCAGTGCCATGCCGTAGATGACCACATTGCGCAACGTGTTGGCGGTCATGCGTACCCGCCACGGGCGCAATATAAGCAGCATG‑ATGATGCTGGCGAAGATCAGgcg                                                                                                  >  1:203469/1‑148 (MQ=255)                                  aCCGTTTGCAGGGCCATATAGAAGAGGAAGTTCATGCCGCCCAGTGCCATGCCGTAGATGACCACATTGCGCAACGTGTTGGCGGTCATGCGTACCCGCCACGGGCGCAATATAAGCAGCATG‑ATGATGCTGGCGAAGATCAGGCGCAg                                                                                               >  7:122500/1‑149 (MQ=255)                                          cAGGGCCATATAGAAGAGGAAGTTCATGCCGCCCAGTGCCATGCCGTAGATGACCACATTGCGCAACGTGTTGGCGGTCATGCGTACCCGCCACGGGCGCAATATAAGCAGCATG‑ATGATGCTGGCGAAGATCAGGCGCAGGGTGGTa                                                                                        <  1:215399/148‑1 (MQ=255)                                             ggCCATATAGAAGAGGAAGTTCATGCCGCCCAGTGCCATGCCGTAGATGACCACATTGCGCAACGTGTTGGCGGTCATGCGTACCCGCCACGGGCTCAATATAAGCAGCATG‑ATGATGCTGGCGAAGATCAGGCGCAGGGTGGTAGTGc                                                                                    >  2:283298/1‑149 (MQ=255)                                             ggCCATATAGAAGAGGAAGTTCATGCCGCCCAGTGCCATGCCGTAGATGACCACATTGCGCAACGTGTTGGCGGTCATGCGTACCCGCCACGGGCGCAATATAAGCAGCATG‑ATGATGCTGGCGAAGATCAGGCGCAGGGTGGTAGTGc                                                                                    >  2:148312/1‑149 (MQ=255)                                               ccATATAGAAGAGGAAGTTCATGCCGCCCAGTGCCATGCCGTAGATGACCACATTGCGCAACGTGTTGGCGGTCATGCGTACCCGCCACGGGCGCAATATAAGCAGCATG‑ATGATGCTGGCGAAGATCAGGCGCAGGGTGGTAGTGCCt                                                                                  >  6:121994/1‑149 (MQ=255)                                                   ataGAAGAGGAAGTTCATGCCGCCCAGTGCCATGCCGTAGATGACCACATTGCGCAACGTGTTGGCGGTCATGCGTACCCGCCACGGGCTCAATATAAGCAGCATGCATGATGCTGGCGAAGATCAGGCGCAGGGTGGTAGTGCCTtgt                                                                               <  1:283298/149‑1 (MQ=255)                                                       aagaGGAAGTTCATGCCGCCCAGTGCCATGCCGTAGATGACCACATTGCGCAACGTGTTGGCGGTCATGCGTACCCGCCACGGGCGCAATATAAGCAGCATG‑ATGATGCTGGCGAAGATCAGGCGCAGGGTGGTAGTGCCTTGTGCGcc                                                                          >  1:166259/1‑149 (MQ=255)                                                         gagGAAGTTCATGCCGCCCAGTGCCATGCCGTAGATGACCACATTGCGCAACGTGTTGGCGGTCATGCGTATCCGCCACGGGCGCAATATAAGCAGCATG‑ATGATGCTGGCGAAGATCAGGCGCAGGGTGGTAGTGCCTTGTGCGCcga                                                                        >  2:170529/1‑149 (MQ=255)                                                         gagGAAGTTCATGCCGCCCAGTGCCATGCCGTAGATGACCACATTGCGCAACGTGTTGGCGGTCATGCGTACCCGCCACGGGCGCAATATAAGCAGCATG‑ATGATGCTGGCGAAGATCAGGCGCAGGGTGGTAGTGCCTTGTGCGcctc                                                                        >  1:262451/1‑147 (MQ=255)                                                         gagGAAGTTCATGCCGCCCAGTGCCATGCCGTAGATGACCACATTGCGCAACGTGTTGGCGGTCATGCGTACCCGCCACGGGCGCAATATAAGCAGCATG‑ATGATGCTGGCGAAGATCAGGCGCAGGGTGGTAGTGCCTTGTGCGCcga                                                                        >  3:224232/1‑149 (MQ=255)                                                         gagGAAGTTCATGCCGCCCAGTGCCATGCCGTAGATGACCACATTGCGCAACGTGTTGGCGGTCATGCGTACCCGCCACGGGCGCAATATAAGCAGCATG‑ATGATGCTGGCGAAGATCAGGCGCAGGGTGGTAGTGCCTTGTGCGCcga                                                                        >  8:4943/1‑149 (MQ=255)                                                            gAAGTTCATGCCGCCCAGTGCCATGCCGTAGATGACCACATTGCGCAACGTGTTGGCGGTCATGCGTACCCGCCACGGGCGCAGTATAAGCAGCATG‑ATGATGCTGGCGAAGATCAGGCGCAGGGTGGTAGTGCCTTGTGCGCcgacga                                                                     >  5:122445/1‑149 (MQ=255)                                                            gAAGTTCATGCCGCCCAGTGCCATGCCGTAGATGACCACATTGCGCAACGTGTTGGCGGTCATGCGTACCCGCCACGGGCGCAGTATAAGCAGCATG‑ATGATGCTGGCGAAGATCAGGCGCAGGGTGGTAGTGCCTTGTGCGCcgacga                                                                     >  5:116703/1‑149 (MQ=255)                                                            gAAGTTCATGCCGCCCAGTGCCATGCCGTAGATGACCACATTGCGCAACGTGTTGGCGGTCATGCGTACCCGCCACGGGCGCAATATAAGCAGCATG‑ATGATGCTGGCGAAGATCAGGCGCAGGGTGGTAGTGCCTTGTGCGCcgacga                                                                     <  2:19302/149‑1 (MQ=255)                                                            gAAGTTCATGCCGCCCAGTGCCATGCCGTAGATGACCACATTGCGCAACGTGTTGGCGGTCATGCGTACCCGCCACGGGCGCAATATAAGCAGCATG‑ATGATGCTGGCGAAGATCAGGCGCAGGGTGGTAGTGCCTTGTGCGCcgacga                                                                     >  1:85691/1‑149 (MQ=255)                                                            gAAGTTCATGCCGCCCAGTGCCATGCCGTAGATGACCACATTGCGCAACGTGTTGGCGGTCATGCGTACCCGCCACGGGCGCAATATAAGCAGCATG‑ATGATGCTGGCGAAGATCAGGCGCAGGGTGGTAGTGCCTTGTGCGCcgacga                                                                     <  6:63520/149‑1 (MQ=255)                                                            gAAGTTCATGCCGCCCAGTGCCATGCCGTAGATGACCACATTGCGCAACGTGTTGGCGGTCATGCGTACCCGCCACGGGCGCAATATAAGCAGCATG‑ATGATGCTGGCGAAGATCAGGCGCAGGGTGGTAGTGCCTTGTGCGCcgacga                                                                     >  1:69776/1‑149 (MQ=255)                                                             aaGTTCATGCCGCCCAGTGCCATGCCGTAGATGACCACATTGCGCAACGTGTTGGCGGTCATGCGTACCCGCCACGGGCGCAATATAAGCAGCATG‑ATGATGCTGGCGAAGATCAGGCGCAGGGTGGTAGTGCCTTGTGCGCCGACGAt                                                                    >  5:201390/1‑149 (MQ=255)                                                               gTTCATGCCGCCCAGTGCCATGCCGTAGATGACCACATTGCGCAACGTGTTGGCGGTCATGCGTACCCGCCACGGGCGCAATATAAGCAGCATG‑ATGATGCTGGCGAAGATCAGGCGCAGGGTGGTAGTGCCTTGTGCGCCGACGATCg                                                                  <  2:69776/149‑1 (MQ=255)                                                                  cATGCCGCCCAGTGCCATGCCGTAGATGACCACATTGCGCAACGTGTTGGCGGTCATGCGTACCCGCCACGGGCGCAATATAAGCAGCATG‑ATGATGCTGGCGAAGATCAGGCGCAGGGTGGTAGTGCCTTGTGCGCCGACGATCggg                                                                <  4:224232/148‑1 (MQ=255)                                                                  cATGCCGCCCAGTGCCATGCCGTAGATGACCACATTGCGCAACGTGTTGGCGGTCATGCGTACCCGCCACGGGCGCAATATAAGCAGCATG‑ATGATGCTGGCGAAGATCAGGCGCAGGGTGGTAGTGCCTTGTGCGCCGACGATCggg                                                                <  6:364230/148‑1 (MQ=255)                                                                      ccgccCAGTGCCATGCCGTAGATGACCACATTGCGCAACGTGTTGGCGGTCATGCGTACCCGCCACGGGCGCAATATAAGCAGCATG‑ATGATGCTGGCGAAGATCAGGCGCAGGGTGGTAGTGCCTTGTGCGCCGACGATCGGGAACa                                                            <  1:290971/148‑1 (MQ=255)                                                                      ccgccCAGTGCCATGCCGTAGATGACCACATTGCGCAACGTGTTGGCGGTCATGCGTACCCGCCACGGGCGCAATATAAGCAGCATG‑ATGATGCTGGCGAAGATCAGGCGCAGGGTGGTAGTGCCTTGTGCGCCGACGATCGGGAACAt                                                           >  5:29715/1‑149 (MQ=255)                                                                        gccCAGTGCCATGCCGTAGATGACCACATTGCGCAACGTGTTGGCGGTCATGCGTACCCGCCACGGGCGCAATATAAGCAGCATG‑ATGATGCTGGCGAAGATCAGGCGCAGGGTGGTAGTGCCTTGTGCGCCGACGATCGGGAACATGc                                                         <  7:4943/149‑1 (MQ=255)                                                                           cAGTGCCATGCCGTAGATGACCACATTGCGCAACGTGTTGGCGGTCATGCGTACCCGCCACGGGCGCAGTATAAGCAGCATG‑ATGATGCTGGCGAAGATCAGGCGCAGGGTGGTAGTGCCTTGTGCGCCGACGATCGGGAACATGCttt                                                      <  6:122445/149‑1 (MQ=255)                                                                           cAGTGCCATGCCGTAGATGACCACATTGCGCAACGTGTTGGCGGTCATGCGTACCCGCCACGGGCGCAGTATAAGCAGCATG‑ATGATGCTGGCGAAGATCAGGCGCAGGGTGGTAGTGCCTTGTGCGCCGACGATCGGGAACATGCttt                                                      <  6:116703/149‑1 (MQ=255)                                                                           cAGTGCCATGCCGTAGATGACCACATTGCGCAACGTGTTGGCGGTCATGCGTACCCGCCACGGGCGCAATATAAGCAGCATG‑ATGATGCTGGCGAAGATCAGGCGCAGGGTGGTAGTGCCTTGTGCGCCGACGCTCGGGAACATGCttt                                                      >  4:291224/1‑149 (MQ=255)                                                                                 cATGCCGTAGATGACCACATTGCGCAACGTGTTGGCGGTCATGCGTACCCGCCACGGGCGCAATATAAGCAGCATG‑ATGATGCTGGCGAAGATCAGGCGCAGGGTGGTAGTGCCTTGTGCGCcgac                                                                       >  5:180872/1‑126 (MQ=255)                                                                                 cATGCCGTAGATGACCACATTGCGCAACGTGTTGGCGGTCATGCGTACCCGCCACGGGCGCAATATAAGCAGCATG‑ATGATGCTGGCGAAGATCAGGCGCAGGGTGGTAGTGCCTTGTGCGCcgac                                                                       <  6:180872/126‑1 (MQ=255)                                                                                 cATGCCGTAGATGACCACATTGCGCAACGTGTTGGCGGTCATGCGTACCCGCCACGGGCGCAATATAAGCAGCATG‑ATGATGCTGGCGAAGATCAGGCGCAGGGTGGTAGTGCCTTGTGCGCCGACGATCGGGAACATGCTTTTGGc                                                  >  2:150405/1‑147 (MQ=255)                                                                                  aTGCCGTAGATGACCACATTGCGCAACGTGTTGGCGGTCATGCGTACCCGCCACGGGCGCAATATAAGCAGCATG‑ATGATGCTGGCGAAGATCAGGCGCAGGGTGGTAGTGCCTTGTGCGCCGACGATCGGGAACATGCTTTTGGCCAg                                               >  3:123203/1‑149 (MQ=255)                                                                                                cacaTTGCGCAACGTGTTGGCGGTCATGCGTACCCGCCACGGGCGCAATATAAGCAGCATG‑ATGATGCTGGCGAAGATCAGGCGCAGGGTGGTAGTGCCTTGTGCGCCGACGATCGGGAACATGCTTTTGGCCAGGGAGGCGCCGGACt                                 >  3:228755/1‑149 (MQ=255)                                                                                                  caTTGCGCAACGTGTTGGCGGTCATGCGTACCCGCCACGGGCGCAATATAAGCAGCATG‑ATGATGCTGGCGAAGATCAGGCGCAgg                                                                                              <  4:276247/86‑1 (MQ=255)                                                                                                  caTTGCGCAACGTGTTGGCGGTCATGCGTACCCGCCACGGGCGCAATATAAGCAGCATG‑ATGATGCTGGCGAAGATCAGGCGCAgg                                                                                              >  3:276247/1‑86 (MQ=255)                                                                                                  caTTGCGCAACGTGTTGGCGGTCATGCGTACCCGCCACGGGCGCAATATAAGCAGCATG‑ATGATGCTGGCGAAGATCAGGCGCAGGGTGGTAGTGCCTTGTGCGCCGACGATCGGGAACATGCTTTTGGc                                                  <  2:167103/130‑1 (MQ=255)                                                                                                  caTTGCGCAACGTGTTGGCGGTCATGCGTACCCGCCACGGGCGCAATATAAGCAGCATG‑ATGATGCTGGCGAAGATCAGGCGCAGGGTGGTAGTGCCTTGTGCGCCGACGATCGGGAACATGCTTTTGGc                                                  >  1:167103/1‑130 (MQ=255)                                                                                                  caTTGCGCAACGTGTTGGCGGTCATGCGTACCCGCCACGGGCGCAATATAAGCAGCATG‑ATGATGCTGGCGAAGATCAGGCGCAGGGTGGTAGTGCCTTGTGCGCCGACGATCGGGAACATGCTTTTGGCCAGGGAGGCGCCGGACTGa                               <  2:203469/149‑1 (MQ=255)                                                                                                  caTTGCGCAACGTGTTGGCGGTCATGCGTACCCGCCACGGGCGCAATATAAGCAGCAGG‑ATGATGCTGGCGAAGATCAGGCGCAGGGTGGTAGTGCCTTGTGCGCCGACGATCGGGAACATGCTTTTGGCCAGGGAGGCGCCGGACTg                                >  6:265019/1‑148 (MQ=255)                                                                                                    ttGCGCAACGTGTTGGCGGTCATGCGTACCCGCCACGGGCGCAATATAAGCAGCATG‑ATGATGCTGGCGAAGATCAGGCGCAGGGTGGTAGTGCCTTGTGCGCCGACGATCGGGAACATGCTTTTGGCCAGGGAGGCGCCGGACTGAAt                             <  6:29715/149‑1 (MQ=255)                                                                                                    ttGCGCAACGTGTTGGCGGTCATGCGTACCCGCCACGGGCGCAATATAAGCAGCATG‑ATGATGCTGGCGAAGATCAGGCGCAGGGTGGTAGTGCCTTGTGCGCCGACGATCGGGAACATGCTTTTGGCCAGGGAGGCGCCGGACTGAAt                             <  5:265019/149‑1 (MQ=255)                                                                                                       cgcAACGTGTTGGCGGTCATGCGTACCCGCCACGGGCGCAATATAAGCAGCATG‑ATGATGCTGGCGAAGATCAGGCGCAGGGTGGTAGTGCCTTGTGCGCCGACGATCGGGAACATGCTTTTGGCCAGGGAGGCGCCGGACTGAATCg                           <  4:123203/148‑1 (MQ=255)                                                                                                           aCGTGTTGGCGGTCATGCGTACCCGCCACGGGCGCAATATAAGCAGCATG‑ATGATGCTGGCGAAGATCAGGCGCAGGGTGGTAGTGCCTTGTGCGCCGACGATCGGGAACATGCTTTTGGCCAGGGAGGCGCCGGACTGAATCGATGcc                      <  2:85691/149‑1 (MQ=255)                                                                                                            cGTGTTGGCGGTCATGCGTACCCGCCACGGGCGCAATATAAGCAGCATG‑ATGATGCTGGCGAAGATCAGGCGCAGGGTGGTAGTGCCTTGTGCGCCGACGATCGGGAACATGCTTTTGGCCAGGGAGGCGCCGGACTGAATCGATGCCa                     <  1:150405/149‑1 (MQ=255)                                                                                                                ttGGCGGTCATGCGTACCCGCCACGGGCGCAATATAAGCAGCATG‑ATGATGCTGGCGAAGATCAGGCGCAGGGTGGTAGTGCCTTGTGCGCCGACGATCGGGAACATGCTTTTGGCCAGGGAGGCGCCGGACTGAATCGATGCCATGGc                 <  5:121994/149‑1 (MQ=255)                                                                                                                ttGGCGGTCATGCGTACCCGCCACGGGCGCAATATAAGCAGCATG‑ATGATGCTGGCGAAGATCAGGCGCAGGGTGGTAGTGCCTTGTGCGCCGACGATCGGGAACATGCTTTTGGCCAGGGAGGCGCCGGACTGAATCGATGCCATGGc                 <  8:122500/149‑1 (MQ=255)                                                                                                                      gTCATGCGTACCCGCCACGGGCGCAATATAAGCAGCATG‑ATGATGCTGGCGAAGATCAGGCGCAGGGTGGTAGTGCCTTGTGCGCCGACGATCGGGAACATGCTTTTGGCCAGGGAGGCGCCGGACTGAATCGATGCCATGGCGATAAg           <  6:201390/149‑1 (MQ=255)                                                                                                                      gTCATGCGTACCCGCCACGGGCGCAATATAAGCAGCATG‑ATGATGCTGGCGAAGATCAGGCGCAGGGTGGTAGTGCCTTGTGCGCCGACGATCGGGAACATGCTTTTGGCCAGGGAGGCGCCGGACTGAATCGATGCCATGGCGATAAg           <  1:148312/149‑1 (MQ=255)                                                                                                                           gCGTACCCGCCACGGGCGCAATATAAGCAGCATG‑ATGATGCTGGCGAAGATCAGGCGCAGGGTGGTAGTGCCTTGTGCGCCGACGATCGGGAACATGCTTTTGGCCAGGGAGGCGCCGGACTGAATCGATGCCATGGCGATAAGCAg        >  3:256511/1‑147 (MQ=255)                                                                                                                            cGTACCCGCCACGGGCGCAATATAAGCAGCATG‑ATGATGCTGGCGAAGATCAGGCGCAGGGTGGTAGTGCCTTGTGCGCCGACGATCGGGAACATGCTTTTGGCCAGGGAGGCGCCGGACTGAATCGATGCCATGGCGATAAGCAGTAg     <  4:256511/149‑1 (MQ=255)                                                                                                                               aCCCGCCACGGGCGCAATATAAGCAGCATG‑ATGATGCTGGCGAAGATCAGGCGCAGGGTGGTAGTGCCTTGTGCGCCGACGATCGGGAACATGCTTTTGGCCAGGGAGGCGCCGGACTGAATCGATGCCATGGCGATAAGCAGTAGGcc  >  1:362736/1‑149 (MQ=255)                                                                                                                                  cGCCACGGGCGCAATATAAGCAGCATG‑ATGATGCTGGCGAAGATCAGGCGCAGGGTGGTAGTGCCTTGTGCGCCGACGATCGGGAACATGCTTTTGGCCAGGGAGGCGCCGGACt                                 >  1:129576/1‑115 (MQ=255)                                                                                                                                  cGCCACGGGCGCAATATAAGCAGCATG‑ATGATGCTGGCGAAGATCAGGCGCAGGGTGGTAGTGCCTTGTGCGCCGACGATCGGGAACATGCTTTTGGCCAGGGAGGCGCCGGACt                                 <  2:129576/115‑1 (MQ=255)                                                                                                                                     cACGGGCGCAATATAAGCAGCATG‑ATGATGCTGGCGAAGATCAGGCGCAGGGTGGTAGTGCCTTGTGCGCCGACGATCGGGAACATGCTTTTGGCCAGGGAGGCGCCGGACTGAATCGATGCCATGGCGATAAGCAGTAGGc   <  4:141949/142‑1 (MQ=255)                                                                                                                                     cACGGGCGCAATATAAGCAGCATG‑ATGATGCTGGCGAAGATCAGGCGCAGGGTGGTAGTGCCTTGTGCGCCGACGATCGGGAACATGCTTTTGGCCAGGGAGGCGCCGGACTGAATCGATGCCATGGCGATAAGCAGTAGGc   >  3:141949/1‑142 (MQ=255)                                                                                                                                               |                                                                                                                                      GTGAATTCCAGGGCGACGGCGATACCGATTGGTACCGTTTGCAGGGCCATATAGAAGAGGAAGTTCATGCCGCCCAGTGCCATGCCGTAGATGACCACATTGCGCAACGTGTTGGCGGTCATGCGTACCCGCCACGGGCGCAATATAAGCAGCATG‑ATGATGCTGGCGAAGATCAGGCGCAGGGTGGTAGTGCCTTGTGCGCCGACGATCGGGAACATGCTTTTGGCCAGGGAGGCGCCGGACTGAATCGATGCCATGGCGATAAGCAGTAGGCC  >  NC_002947/2547882‑2548156

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