Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,428,461	T→G	5.8%	D102A (GAC→GCC)	sohB ←	peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,428,461	0	T	G	5.8%	215.3 / 4.0	69	D102A (GAC→GCC)	sohB	peptidase
Reads supporting (aligned to +/- strand):  ref base T (27/38);  new base G (2/2);  total (29/40)
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.65e-01

ACGCGCCAGTTGCGACGCAGCCAAGCCGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGG‑CCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTGGGCCTTGTCGAGCAGGCCGGACTCCAGGCGCTCACGCAGGTC  >  NC_002947/4428313‑4428606                                                                                                                                                      |                                                                                                                                                  aCGCGCCAGTTGCGACGCAGCCAAGCCGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGG‑CCGAGGCTTTGATGt                                                                                                                                                   >  2:297452/1‑149 (MQ=255)  cgcgCCAGTTGCGACGCAGCCAAGCCGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGG‑CCGAGGCTTTGATGTc                                                                                                                                                  <  7:172482/149‑1 (MQ=255)   gcgcCAGTTGCGACGCAGCCAAGCCGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGG‑CCGAGGCTTTGATGTCg                                                                                                                                                 >  1:54728/1‑149 (MQ=255)     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                                                                                                                                        >  1:63301/1‑137 (MQ=255)                        aGCCGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGG‑CCGAGGCTTTGATGGCGACATCAAAg                                                                                                                                        <  2:63301/137‑1 (MQ=255)                           cGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGG‑CCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATAcc                                                                                                                         <  5:223203/149‑1 (MQ=255)                                 tgtgCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGG‑CCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGctc                                                                                                                   >  2:149056/1‑149 (MQ=255)                                   tgCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGG‑CCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGctc                                                                                                                   <  2:243924/147‑1 (MQ=255)                                   tgCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGG‑CCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGctc                                                                                                                   <  2:95954/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                                                                                                         >  3:128552/1‑149 (MQ=255)                                           gCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGG‑CCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCgg                                                                                                         <  6:70192/149‑1 (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                                                                                        <  5:266704/149‑1 (MQ=255)                                                            ggCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGG‑CCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTctgc                                                                                        >  4:223237/1‑149 (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)                                                                                                                                                      |                                                                                                                                                  ACGCGCCAGTTGCGACGCAGCCAAGCCGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGG‑CCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTGGGCCTTGTCGAGCAGGCCGGACTCCAGGCGCTCACGCAGGTC  >  NC_002947/4428313‑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: ‑