Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,605,506	G→C	8.3%	E308Q (GAG→CAG) ‡	PP_4077 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,605,506	0	G	C	8.3%	179.7 / 7.1	60	E308Q (GAG→CAG) ‡	PP_4077	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (26/29);  new base C (2/3);  total (28/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.60e-01

CCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACTTTTCGCCCTTTCGAATCGTTC  >  NC_002947/4605367‑4605652                                                                                                                                             |                                                                                                                                                   cccACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGcc                                                                                                                                            >  6:297540/1‑148 (MQ=255)    acaGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTAGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GGCAGCCTCGCTGGCCAgg                                                                                                                                        <  6:264436/149‑1 (MQ=255)     caGGTGGGCATCGAGTTCCATTTCACGCGTCACCTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GGCAGCCTCGCTGGCCAGGt                                                                                                                                       <  4:234132/149‑1 (MQ=255)            gCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCt                                                                                                                                 >  6:259021/1‑148 (MQ=255)               tCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAg                                                                                                                             <  8:124229/149‑1 (MQ=255)                 gAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GGCAGCCTCGCTGGCCAGGt                                                                                                                                       <  4:224288/137‑1 (MQ=255)                 gAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GGCAGCCTCGCTGGCCAGGt                                                                                                                                       >  3:224288/1‑137 (MQ=255)                 gAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGcc                                                                                                                           <  5:259021/149‑1 (MQ=255)                       cATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCAGAGAAGATTTCAGGCTGTTCTGCAGCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGTCCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGaa                                                                                                                      <  1:55966/148‑1 (MQ=255)                         tttCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCg                                                                                                                   <  2:262042/149‑1 (MQ=255)                            cACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAg                                                                                                                 <  5:105633/148‑1 (MQ=255)                                          cAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTAc                                                                                                  <  5:297540/149‑1 (MQ=255)                                          cAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTAc                                                                                                  <  3:122129/149‑1 (MQ=255)                                           aGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGACGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGTCGCAGTCCGCTAAc                                                                                                  >  8:263412/1‑148 (MQ=255)                                            gACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCt                                                                                                     <  2:237030/144‑1 (MQ=255)                                            gACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCt                                                                                                     >  1:237030/1‑144 (MQ=255)                                            gACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACg                                                                                                 >  8:24035/1‑148 (MQ=255)                                             acccccTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGaa                                                                                               <  8:278372/146‑1 (MQ=255)                                              cacCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGaaa                                                                                              >  5:88636/1‑149 (MQ=255)                                                 ccTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCt                                                                                           >  5:306820/1‑149 (MQ=255)                                                 ccTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCt                                                                                           >  1:240711/1‑149 (MQ=255)                                                           ggCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTg                                                                                                                                >  8:289127/1‑102 (MQ=255)                                                           ggCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTg                                                                                                                                <  7:289127/102‑1 (MQ=255)                                                            gCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGAt                                                                                <  4:157111/149‑1 (MQ=255)                                                              acaGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGaa                                                                                                                      <  6:277060/109‑1 (MQ=255)                                                              acaGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGaa                                                                                                                      >  5:277060/1‑109 (MQ=255)                                                                agaagaTTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGa                                                                                                                                                   >  4:144118/1‑78 (MQ=255)                                                                agaagaTTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGa                                                                                                                                                   <  3:144118/78‑1 (MQ=255)                                                                agaagaTTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAg                                                                            >  6:99951/1‑149 (MQ=255)                                                                  aagaTTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAGGt                                                                           >  8:98419/1‑148 (MQ=255)                                                                  aagaTTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGtc                                                                          <  6:88636/149‑1 (MQ=255)                                                                     aTTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGtcatc                                                                       <  1:251464/149‑1 (MQ=255)                                                                      tttCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGtcatca                                                                      <  5:99951/149‑1 (MQ=255)                                                                         cAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAGGTCATCAGc                                                                    <  7:98419/148‑1 (MQ=255)                                                                            ccTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCAccc                                                                >  4:242076/1‑149 (MQ=255)                                                                            ccTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCAccc                                                                >  5:113088/1‑149 (MQ=255)                                                                                ttCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGa                                                              <  3:242076/147‑1 (MQ=255)                                                                                    gCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGaacaac                                                         >  1:264789/1‑148 (MQ=255)                                                                                      aCCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCAc                                                                  <  8:36987/137‑1 (MQ=255)                                                                                      aCCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCAc                                                                  >  7:36987/1‑137 (MQ=255)                                                                                        cccGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAg                                                                     >  2:252020/1‑132 (MQ=255)                                                                                        cccGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAg                                                                     <  1:252020/132‑1 (MQ=255)                                                                                                tCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTTTCGATGAAGTCATAAGCACCCGAACAACCACTGACGGaa                                              <  7:263412/147‑1 (MQ=255)                                                                                                tCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGaa                                              <  2:264789/147‑1 (MQ=255)                                                                                                    ccTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCg                                        >  1:117718/1‑149 (MQ=255)                                                                                                    ccTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCg                                        >  3:151470/1‑149 (MQ=255)                                                                                                    ccTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCg                                        >  2:165167/1‑149 (MQ=255)                                                                                                       gTTCAAGCACAGTGCCGAGCCCATCGACCTGCGCCAGCGAGGCTGCCCAGGTCGAGATGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGt                                      >  1:188828/1‑149 (MQ=255)                                                                                                                       gAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGATCAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACttt                     >  6:8708/1‑149 (MQ=255)                                                                                                                         gCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACTTTTc                   <  6:306820/149‑1 (MQ=255)                                                                                                                           ccATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTAAGTGTCGATGAAGTCATCCGACCCCGAACAACCACGGACGGAAGTACCGGGGCGCACCGTCGGACATTgcgc                 >  4:228774/1‑149 (MQ=255)                                                                                                                            cATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACTTTTCGcc                <  2:240711/149‑1 (MQ=255)                                                                                                                             aTCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGAc                                                  <  7:245857/114‑1 (MQ=255)                                                                                                                             aTCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGAc                                                  >  8:245857/1‑114 (MQ=255)                                                                                                                                gACCTG‑GCCAGCGAGGCTGCCCAGTTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACTTTTCGCCCttt            <  3:228774/149‑1 (MQ=255)                                                                                                                                gACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACTTTTCGCCCttt            <  7:24035/149‑1 (MQ=255)                                                                                                                                 aCCTG‑GGCAGCCTCGCTGGCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGCGCACCTGCGGACTTTTCGCCCTTTc           >  3:44188/1‑149 (MQ=38)                                                                                                                                 aCCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACTTTTCGCCCTTTc           <  6:113088/149‑1 (MQ=255)                                                                                                                                           gCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGATCAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACTTTTCGCCCTTTCGAATCGTTc  <  5:8708/149‑1 (MQ=255)                                                                                                                                           gCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACTTTTCGCCCTTTCGAATCGTTc  >  6:175994/1‑149 (MQ=255)                                                                                                                                             |                                                                                                                                                   CCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACTTTTCGCCCTTTCGAATCGTTC  >  NC_002947/4605367‑4605652

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