Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,539,279	C→G	6.6%	K569N (AAG→AAC)	PP_4872 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,539,279	0	C	G	6.6%	193.1 / 4.4	61	K569N (AAG→AAC)	PP_4872	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (29/28);  new base G (2/2);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.63e-01

ACGGGTAGAACGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCAGGGGGCCACGCTCGGTGTGCTCCGGGGCAATCTTCAGGTAGCCGCCAACGTGGTGGGTGACCAGTTCCTTGACGTACTCCGG  >  NC_002947/5539131‑5539422                                                                                                                                                     |                                                                                                                                                aCGGGTAGAACGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAAc                                                                                                                                                 <  7:66449/149‑1 (MQ=255)    ggTATAACGCCTGTTCGTGGTCGGCGCGGAAGAGGTTGAGCTTCAGACTCTGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCTCGGCAATGAAGTACGGGATCTGGTACTGGTCGTTGACCGCTTCTTTCGAGAACTtc                                                                                                                                              <  5:116554/149‑1 (MQ=255)     gTAGAACGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTtct                                                                                                                                             >  7:160990/1‑149 (MQ=255)      tAGAACGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAAGTtctc                                                                                                                                            >  3:54869/1‑149 (MQ=255)      tAGAACGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTtctc                                                                                                                                            >  8:53341/1‑149 (MQ=255)      tAGAACGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTtctc                                                                                                                                            <  2:287060/149‑1 (MQ=255)          aCGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTTCCCGCTTCTTTCGAGAACTTAACGaa                                                                                                                                         >  7:331175/1‑148 (MQ=255)          aCGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAAGTTCTCGaaa                                                                                                                                        <  8:48614/149‑2 (MQ=255)          aCGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAAc                                                                                                                                        <  8:160990/149‑1 (MQ=255)            gCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAAGTTCTCGaaaga                                                                                                                                      <  4:54869/149‑4 (MQ=255)              cTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATgc                                                                                                                                    >  7:143419/1‑149 (MQ=255)                   ccTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCtt                                                                                                                                >  6:180181/1‑148 (MQ=255)                    cTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTg                                                                                                                               >  2:248506/1‑148 (MQ=255)                       gTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCTAGAACTTCTCGAACATGCGCTTGAAGc                                                                                                                           >  4:308452/1‑149 (MQ=255)                       gTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAg                                                                                                                            <  8:143419/148‑1 (MQ=255)                       gTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGc                                                                                                                           >  6:192355/1‑149 (MQ=255)                        tCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCg                                                                                                                          >  6:190837/1‑149 (MQ=255)                            cgcgGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCg                                                                                                                      <  1:248506/149‑1 (MQ=255)                               ggAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAg                                                                                                                   <  6:18932/149‑1 (MQ=255)                                gAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAgg                                                                                                                  >  5:181114/1‑149 (MQ=255)                                    acGTGGCCCTTCAGCCACAGGTCCAGGTTCCTCTTGTCTTGGTCGGTGGTAGCGGGGTGCGCGGCAATGAAGTAGGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGATGCGGTCGTAGGTGc                                                                                                               <  8:331175/147‑1 (MQ=255)                                       ttGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGa                                                                                                            >  8:36068/1‑148 (MQ=255)                                          cccTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATAc                                                                                                         <  6:181114/148‑1 (MQ=255)                                           gcTTCAGCCACAGGGCCAGGTTCATGATGTCTTCGTCGGTTGTACCGGTGTGCGCGGCAATGAAGTAGGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTGTCGATCATGCGCTTGAAGCg                                                                                                                          <  4:260852/129‑1 (MQ=255)                                           ccTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCg                                                                                                                          >  3:260852/1‑130 (MQ=255)                                                 gCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCtt                                                                                                  >  5:64359/1‑148 (MQ=255)                                                  ccACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTtca                                                                                                <  6:64359/149‑1 (MQ=255)                                                    acaGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTtcatc                                                                                              <  7:36068/149‑1 (MQ=255)                                                    acaGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTtcatc                                                                                              <  5:180181/149‑1 (MQ=255)                                                     caGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGc                                                                                                    <  6:49392/142‑1 (MQ=255)                                                     caGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGc                                                                                                    >  5:49392/1‑142 (MQ=255)                                                     caGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTtcatca                                                                                             >  2:89892/1‑149 (MQ=255)                                                          ccAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTg                                                                                        <  3:148366/149‑1 (MQ=255)                                                                     atGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCAGGGGGcc                                                                             <  5:190837/149‑1 (MQ=255)                                                                          ttCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCAGGGGGCCACGCt                                                                        >  4:191560/1‑149 (MQ=255)                                                                             gtcgGTGGTACCGGGGTGCGTGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCAGGGGGCCACGCTCgg                                                                     >  8:81736/1‑149 (MQ=255)                                                                              tcgGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCAGGGGGCCACGCTCGgt                                                                    <  3:191560/149‑1 (MQ=255)                                                                                   ggtACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCAGGGGGCCACGCTCGGTGTGCt                                                               >  8:51882/1‑149 (MQ=255)                                                                                                 cgGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCAGGGGGCCACGCTCGGTGTGCTCCGGGGCAATCTTc                                                 <  5:192355/149‑1 (MQ=255)                                                                                                  ggCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCatcatc                                                                                           <  1:266526/106‑1 (MQ=255)                                                                                                  ggCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCatcatc                                                                                           >  2:266526/1‑106 (MQ=255)                                                                                                      aTGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCAGGGGGCCACGCTCGGTGTGCTCCGGGGCAATCTTCAGGTa                                            <  7:51882/149‑1 (MQ=255)                                                                                                                     aGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCAGGGGGCCACGCTCGGTGTGCTCCGGGGCAATCTTCAGGTAGCCGCCAACgtggt                              <  7:53341/148‑1 (MQ=255)                                                                                                                          cTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCAGGGGGCCACGCTCGGTGTGCTCCGGGGCAATCTTCAGGTAGCCGCCAACGTGGTGGGTGa                        <  1:89892/149‑1 (MQ=255)                                                                                                                           tGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCAGGGGGCCACGCTCGGTGTGCTCCGGGGCAATCTTCAGGTAGCCGCCAACGTGGTGGGTGAc                       >  2:332927/1‑149 (MQ=255)                                                                                                                            gCTCCTTGCCCGCTTCTTTCTAGAACTTCTCGAACATGCGATTGAAGCTGTCGTAGTTGCCGATACAAGGCTTCATGATCTTGTCCAGGGGGCAACGCTCGGTGTGATCCGGGTCAATCTACAGTTAGACGACAACGAGGTGGGTGAc                       >  6:376795/1‑148 (MQ=255)                                                                                                                             cTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCAGGGGGCCACGCTCGGTGTGCTCCGGGGCAATCTTCAGGTa                                            >  1:332203/1‑126 (MQ=255)                                                                                                                             cTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCAGGGGGCCACGCTCGGTGTGCTCCGGGGCAATCTTCAGGTa                                            <  2:332203/126‑1 (MQ=255)                                                                                                                                cTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCAGGGGGCCACGCTCGGTGTGCTCCGGGGCAATCTTCAGGTAGCCGCCAACGTGGTGGGTGACCAGtt                  >  4:303864/1‑149 (MQ=255)                                                                                                                                cTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCAGGGGGCCACGCTCGGTGTGCTCCGGGGCAATCTTCAGGTAGCCGCCAACGTGGTGGGTGACCAGtt                  >  7:334585/1‑149 (MQ=255)                                                                                                                                 ttGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCAGGGGGCCACGCTCGGTGTGCTCCGGGGCAATCTTCAGGTAGCCGCCAACGtggtgg                            <  2:74487/138‑1 (MQ=255)                                                                                                                                 ttGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCAGGGGGCCACGCTCGGTGTGCTCCGGGGCAATCTTCAGGTAGCCGCCAACGtggtgg                            >  1:74487/1‑138 (MQ=255)                                                                                                                                 ttGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCAGGGGGCCACGCTCGGTGTGCTCCGGGGCAATCTTCAGGTAGCCGCCAACGTGGTGGGTGACCAGTTc                 >  5:124298/1‑149 (MQ=255)                                                                                                                                   gCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCAggggg                                                                               >  2:6170/1‑85 (MQ=255)                                                                                                                                    cccGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCAGGGGGCCACGCTCGGTGTGCTCCGGGGCAATCTTCAGGTAGCCGCCAACGTGGTGGGTGACCAGTTCCtt              <  8:334585/149‑1 (MQ=255)                                                                                                                                         ttcttTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCAGGGGGCCACGCTCGGTGTGCTCCGGGGCAATCTTCAGGTAGCCGCCAACGTggtgggt                          >  5:40553/1‑132 (MQ=255)                                                                                                                                         ttcttTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCAGGGGGCCACGCTCGGTGTGCTCCGGGGCAATCTTCAGGTAGCCGCCAACGTggtgggt                          <  6:40553/132‑1 (MQ=255)                                                                                                                                         catgttCGAGAAGTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCAGGGGGCCACGCTCGGTGTGCTCCGGGGCAATCTTCAGGTAGCCGCCAACGTGGTGGGTGACCAGTTCCTTGACg          >  7:301624/5‑148 (MQ=255)                                                                                                                                           cttTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCAGGGGGCCACGCTCGGTGTGCTCCGGGGCAATCTTCAGGTAGCCGCCAACGTGGTGGGTGACCAGTTCCTTGACGTa        <  1:332927/148‑1 (MQ=255)                                                                                                                                                gagaACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCAGGGGGCCACGCTCGGTGTGCTCCGGGGCAATCTTCAGGTAGCCGCCAACGTGGTGGGTGACCAGTTCCTTGACGTACTCCgg  <  3:308452/149‑1 (MQ=255)                                                                                                                                                 agaACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCAGGGGGCCACGCTCGGTGTGCTCCGGGGCAATCTTCAGGTAGCCGCCAACgtggt                              <  3:236318/120‑1 (MQ=255)                                                                                                                                                 agaACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCAGGGGGCCACGCTCGGTGTGCTCCGGGGCAATCTTCAGGTAGCCGCCAACgtggt                              >  4:236318/1‑120 (MQ=255)                                                                                                                                                     |                                                                                                                                                ACGGGTAGAACGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGCTTCATCATCTTGTCCAGGGGGCCACGCTCGGTGTGCTCCGGGGCAATCTTCAGGTAGCCGCCAACGTGGTGGGTGACCAGTTCCTTGACGTACTCCGG  >  NC_002947/5539131‑5539422

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