Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,428,464	C→A	6.2%	G101V (GGC→GTC)	sohB ←	peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,428,464	0	C	A	6.2%	203.3 / 3.5	64	G101V (GGC→GTC)	sohB	peptidase
Reads supporting (aligned to +/- strand):  ref base C (24/36);  new base A (2/2);  total (26/38)
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.03e-01

GCCAGTTGCGACGCAGCCAAGCCGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGG‑CCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTGGGCCTTGTCGAGCAGGCCGGACTCCAGGCGCTCACGCAGGTC  >  NC_002947/4428317‑4428606                                                                                                                                                     |                                                                                                                                               gcCAGTTGCGACGCAGCCAAGCCGTAAATGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGG‑CCGAGGCTTTGATGTCGc                                                                                                                                                >  1:198010/1‑148 (MQ=255)   cAGTTGCGACGCAGCCAAGCCGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGG‑CCGAGGCTTTGATGTCGCCAt                                                                                                                                             >  5:70192/1‑149 (MQ=255)        gCGACGCAGCCAAGCCGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGG‑CCGAGGCTTTGATGTCGCCATCaaa                                                                                                                                         <  8:20134/148‑1 (MQ=255)            cGCAGCCAAGCCGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGG‑CCGAGGCTTTGATGTCGCCATCAAAGTCCa                                                                                                                                    <  3:292265/149‑1 (MQ=255)            cGCAGCCAAGCCGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTACGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTACGGTGG‑CCGAGGCTTTGATGTCGCCATCAAAGTCCa                                                                                                                                    <  5:157884/149‑1 (MQ=255)             gCAGCCAAGCCGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCCGAGGCTTTGATGTCGCCATCAAAGTCCa                                                                                                                                    <  1:32559/149‑1 (MQ=255)             gCAGCCAAGCCGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGG‑CCGAGGCTTTGATGTCGc                                                                                                                                                <  8:248508/136‑1 (MQ=255)             gCAGCCAAGCCGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGG‑CCGAGGCTTTGATGTCGc                                                                                                                                                >  7:248508/1‑136 (MQ=255)             gCAGCCAAGCCGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGG‑CCGAGGCTTTGATGTCGCCATCAAAGTCCAg                                                                                                                                   <  2:54728/149‑1 (MQ=255)               aGCCAAGCCGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGG‑CCGAGGCTTTGATGTCGCCATCAAAGTc                                                                                                                                      >  4:153984/1‑144 (MQ=255)               aGCCAAGCCGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGG‑CCGAGGCTTTGATGTCGCCATCAAAGTc                                                                                                                                      <  3:153984/144‑1 (MQ=255)               aGCCAAGCCGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGG‑CCGAGGCTTTGATGTCGCCATCAAAGTCCAGTa                                                                                                                                 <  5:245912/149‑1 (MQ=255)                   aaGCCGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGG‑CCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGaa                                                                                                                             >  1:95954/1‑149 (MQ=255)                    aGCCGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGG‑CCGAGGCTTTGATGGCGACATCAAAg                                                                                                                                        <  2:63301/137‑1 (MQ=255)                    aGCCGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGG‑CCGAGGCTTTGATGGCGACATCAAAg                                                                                                                                        >  1:63301/1‑137 (MQ=255)                       cGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGG‑CCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATAcc                                                                                                                         <  5:223203/149‑1 (MQ=255)                             tgtgCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGG‑CCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGctc                                                                                                                   >  2:149056/1‑149 (MQ=255)                               tgCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGG‑CCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGctc                                                                                                                   <  2:95954/147‑1 (MQ=255)                               tgCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGG‑CCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGctc                                                                                                                   <  2:243924/147‑1 (MQ=255)                                    cAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGG‑CCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCt                                                                                                            >  4:176995/1‑149 (MQ=255)                                     aGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGG‑CCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTc                                                                                                           >  3:212845/1‑149 (MQ=255)                                      ggCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGG‑CCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATAGCCGGCTCTTGTCCTCg                                                                                                          <  5:24702/149‑1 (MQ=255)                                       gCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGG‑CCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCgg                                                                                                         <  6:70192/149‑1 (MQ=255)                                       gCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGG‑CCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCgg                                                                                                         >  3:128552/1‑149 (MQ=255)                                          aCCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGG‑CCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCt                                                                                                      >  8:109391/1‑149 (MQ=255)                                            cgccgcTTTCCCGGCGCAGTACCACTTCGTCACGCGCGGTAGCGATCGTCAGCAGAGCGGTGATTTCGTTGAGCAGGCTTTAGGTGG‑CCGAGGATTTGATGTCGCCATCAAAGTCCAGTACGAATACCCTGCTCTTGTACTCGGCCTTg                                                                                                    >  1:348180/1‑149 (MQ=255)                                            cgccgcTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGG‑CCGAGGCTTTGATGTCGCCATCaaa                                                                                                                                         >  2:119821/1‑112 (MQ=255)                                             gccgcTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGG‑CCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGc                                                                                                   >  8:232871/1‑149 (MQ=255)                                               cgcTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGG‑CCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGccc                                                                                                 <  4:107695/149‑1 (MQ=255)                                                 cgttCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGG‑CCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCct                                                                                                >  7:236023/3‑148 (MQ=255)                                                    tCCATGCGCAGTACCAGTTCGTCACGCGCGGTAGCGAGCGTCAGAAGTGCGGTGATTTCGTTGCGCAGGGTTTCGGTGG‑CCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCcttc                                                                                              <  7:65557/147‑1 (MQ=255)                                                      cAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGG‑CCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCttcttc                                                                                           <  8:236023/148‑1 (MQ=255)                                                      cAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGG‑CCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTtcttct                                                                                          <  8:221585/149‑1 (MQ=255)                                                        ggCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGG‑CCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTctgc                                                                                        >  4:223237/1‑149 (MQ=255)                                                        ggCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGG‑CCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTctgc                                                                                        <  5:266704/149‑1 (MQ=255)                                                        ggCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGG‑CCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTctgc                                                                                        <  2:31590/149‑1 (MQ=255)                                                          cgcAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGG‑CCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTctgctg                                                                                      >  2:16208/1‑149 (MQ=255)                                                           gcAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGTAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGG‑CCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGt                                                                                     >  7:117596/1‑149 (MQ=255)                                                                  cACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGG‑CCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTcc                                                                               <  3:290189/148‑1 (MQ=255)                                                                  cACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGG‑CCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTcc                                                                               >  4:290189/1‑148 (MQ=255)                                                                   aCTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGG‑CCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACTCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGc                                                                             >  7:257149/1‑149 (MQ=255)                                                                            cgcgcgGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGG‑CCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCctg                                                                    >  1:352442/1‑149 (MQ=255)                                                                            cgcgcgGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGG‑CCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCctg                                                                    <  7:232871/149‑1 (MQ=255)                                                                             gcgcgGTAGCGAGCGTCAGCAGTTCGGTGATTTCGTTGCGCAGGCTTTCGGTGT‑CCGAGGCTTTGATGTCTCCATCAAAGTACAGTACGAATACCCGGCTCTTGTCCTAGTCCTTGCCCTTCTTCTGATGTTATTCAGCCTTGGCatga                                                                   >  4:39210/1‑148 (MQ=255)                                                                             gcgcgGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGG‑CCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCctgc                                                                   <  1:16208/149‑1 (MQ=255)                                                                                     gcgagcgTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGG‑CCGAGGCTTTGATGTCGCCATCAAAg                                                                                                                                        <  7:131536/72‑1 (MQ=255)                                                                                     gcgagcgTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGG‑CCGAGGCTTTGATGTCGCCATCAAAg                                                                                                                                        <  7:136494/72‑1 (MQ=255)                                                                                     gcgagcgTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGG‑CCGAGGCTTTGATGTCGCCATCAAAg                                                                                                                                        >  8:136494/1‑72 (MQ=255)                                                                                     gcgagcgTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGG‑CCGAGGCTTTGATGTCGCCATCAAAg                                                                                                                                        >  8:131536/1‑72 (MQ=255)                                                                                          cgTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGG‑CCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACTCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAggg                                                      <  8:257149/149‑1 (MQ=255)                                                                                              agcagTTCGGTGATTTCGTTGCGCAGGCTTTCGGTGG‑CCGAGGCTTTGATGTCGCCATCAAAGTCCATTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTATGCTGTTTTTCCGCCTTGTCCTTCTGCTTGCGCAGGGCCtt                                                  <  2:348180/149‑1 (MQ=255)                                                                                              agcagTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGG‑CCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCtt                                                  <  1:240560/149‑1 (MQ=255)                                                                                                cagTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGG‑CCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGa                                                <  3:223237/149‑1 (MQ=255)                                                                                                cagTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGG‑CCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGa                                                <  2:352442/149‑1 (MQ=255)                                                                                                  gTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGG‑CCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGc                                              <  3:168183/149‑1 (MQ=255)                                                                                                             tCGTTGCGCAGGCTTTCGGTGG‑CCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTGGGCCTTGt                                    <  1:149056/148‑1 (MQ=255)                                                                                                               gTTGCGCAGGCTTTCGGTGG‑CCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTGGGCCTTGTCGa                                 <  3:176995/149‑1 (MQ=255)                                                                                                                     cAGGCTTTCGGTGG‑CCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTGGGCCTTGTCGAGCAGGc                           <  4:212845/149‑1 (MQ=255)                                                                                                                       ggCTTTCGGTGG‑CCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTGGGCCTTGTCGAGCAGGCCg                         <  4:128552/149‑1 (MQ=255)                                                                                                                         cTTTCGGTGG‑CCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTGGGCCTTGTCGAGCAGGCCGGa                       <  2:198010/149‑1 (MQ=255)                                                                                                                                     cGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTGGGCCTTGTCGAGCAGGCCGGACTCCAGGCGCt            <  8:117596/149‑1 (MQ=255)                                                                                                                                      gAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTGGGCCTTGTCGAGCAGGCCGGACTCCAGGCGCt            >  2:201512/1‑148 (MQ=255)                                                                                                                                          cTTTGATGGCGACATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTGGGCCTTGTCGAGCAGGCCGGACTCCAGGCGCTCACg        >  2:357155/1‑148 (MQ=255)                                                                                                                                              gATGGCGACATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTGGGCCTTGTCGAGCAGGCCGGACTCCAGGCGCTCACGCAgg    <  1:357155/148‑1 (MQ=255)                                                                                                                                               aTGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTgg                                           >  1:270925/1‑108 (MQ=255)                                                                                                                                               aTGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTgg                                           <  2:270925/108‑1 (MQ=255)                                                                                                                                               aTGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTGGGCCTTGTCGAGCAGGCCGGACTCCAGGCGCTCACGCAGGTc  >  2:120323/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                               GCCAGTTGCGACGCAGCCAAGCCGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGG‑CCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTGGGCCTTGTCGAGCAGGCCGGACTCCAGGCGCTCACGCAGGTC  >  NC_002947/4428317‑4428606

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