Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	436,899	C→A	9.7%	W130L (TGG→TTG)	PP_0358 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	436,899	0	C	A	9.7%	221.2 / 12.5	72	W130L (TGG→TTG)	PP_0358	membrane protein
Reads supporting (aligned to +/- strand):  ref base C (34/31);  new base A (2/5);  total (36/36)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.65e-01

GGGCAGCGCAGCCACCGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCA‑CTTGCCACACC‑‑‑‑‑‑‑‑TGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGC‑‑CATTGCAGGCGG‑CGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGCCAGTTGCCGGCGATCAG  >  NC_002947/436757‑437044                                                                                                                                                          |                                                                                                                                                   gggCAGCGCAGCCACCGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCA‑CTTGCCACACC‑‑‑‑‑‑‑‑TGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGC‑‑CATTg                                                                                                                                                >  8:361187/1‑147 (MQ=255)   gcagcgcagcCACCGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCA‑CTTGCCACATC‑‑‑‑‑‑‑‑TGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGC‑‑CATTGCAgg                                                                                                                                            >  4:267715/1‑149 (MQ=255)     agcgcagcCACCGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCA‑CTTGCCACACC‑‑‑‑‑‑‑‑TGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGC‑‑CATTGCAggcg                                                                                                                                          <  5:338890/149‑1 (MQ=255)     agcgcagcCACCGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCA‑CTTGCCACACC‑‑‑‑‑‑‑‑TGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGC‑‑CATTGCAggcg                                                                                                                                          <  3:264432/149‑1 (MQ=255)      gcgcagcCACCGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCCGGCTGCCCA‑CTTGCCACACC‑‑‑‑‑‑‑‑TGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGC‑‑CATTGCAggcgg                                                                                                                                         <  5:343909/149‑1 (MQ=255)      gcgcagcCACCGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCA‑CTTGCCACACC‑‑‑‑‑‑‑‑TGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCTGC‑‑AATGGCTggcgg                                                                                                                                         >  3:139277/1‑149 (MQ=255)       cgcagcCACCGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCA‑CTTGCCACACC‑‑‑‑‑‑‑‑TGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCTGC‑‑AATGGCTggcggt                                                                                                                                        <  8:47251/149‑2 (MQ=255)       cgcagcCACCGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCA‑CTTGCCACACC‑‑‑‑‑‑‑‑TGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCTGC‑‑AATGGCTggcggt                                                                                                                                        <  2:290012/149‑2 (MQ=255)        gcagcCACCGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCA‑CTTGCCACACC‑‑‑‑‑‑‑‑TGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCTGC‑‑AATGGCTggcggtg                                                                                                                                       <  4:194202/149‑3 (MQ=255)         cagcCACCGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCA‑CTTGCCACATC‑‑‑‑‑‑‑‑TGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGC‑‑CATTGCAGGCGg‑cgc                                                                                                                                     <  3:267715/149‑1 (MQ=255)         cagcCACCGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCA‑CTTGCCACACC‑‑‑‑‑‑‑‑TGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCTGC‑‑AATGGCTggcggtgg                                                                                                                                      <  4:139277/149‑4 (MQ=255)         cagcCACCGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCA‑CTTGCCACACC‑‑‑‑‑‑‑‑TGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCTGC‑‑AATGGCTggcggtgg                                                                                                                                      <  1:64605/149‑4 (MQ=255)             cACCGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCA‑CTTGCCACACC‑‑‑‑‑‑‑‑TGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGC‑‑CATTGCAGGCGG‑CGCAGgc                                                                                                                                 <  6:267754/149‑1 (MQ=255)                  gcgcCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCA‑CTTGCCACACC‑‑‑‑‑‑‑‑TGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGC‑‑CATTGCACGCGG‑CGCAGTCTCTCa                                                                                                                            >  6:405151/1‑149 (MQ=255)                      cTGCTTACGGATCAGCCCCTAGAAGCCGAAGCTCAACGCCAGCACCAGAGCTACCCACGGCAGGCTGCCCA‑CTTGCCACACC‑‑‑‑‑‑‑‑TGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGC‑‑CATTGCAGGCGG‑CGCAGGCGCTCACCaa                                                                                                                        <  8:341716/149‑1 (MQ=255)                         cTTAGGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCA‑CTAGCCACACC‑‑‑‑‑‑‑‑TGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGC‑‑CATTGCAGGCGG‑CGCAGGCGCTCACCAAGa                                                                                                                      <  4:312008/148‑1 (MQ=255)                         cTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCA‑CTTGCCACACC‑‑‑‑‑‑‑‑TGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGC‑‑CATTGCAGGCGG‑CGCAGGCGCTCACCAAg                                                                                                                       <  7:361187/147‑1 (MQ=255)                         cTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCA‑CTTGCCACACC‑‑‑‑‑‑‑‑TGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGC‑‑CATTGCAGGCGG‑CGCAGGCGCTCACCAAGa                                                                                                                      <  4:359088/148‑1 (MQ=255)                            aCGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCA‑CTTGCCACACC‑‑‑‑‑‑‑‑TGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGC‑‑CATTGCAGGCGG‑CGCAGGCGCTCACCAAGAAt                                                                                                                    >  1:51713/1‑147 (MQ=255)                                      ccATAGAATCAGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCA‑CTTGCCACACC‑‑‑‑‑‑‑‑TGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGC‑‑CATTGCAGGCGG‑CGCAGGCGCTCACCAAGAATCAGCATGCCcag                                                                                                        <  5:405151/149‑1 (MQ=255)                                        aTAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCA‑CTTGCCACACC‑‑‑‑‑‑‑‑TGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGC‑‑CATTGCAGGCGG‑CGCAGGCGCTCACCAAGAATCAGCATGCCcag                                                                                                        >  6:166816/1‑147 (MQ=255)                                        aTAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCA‑CTTGCCACACC‑‑‑‑‑‑‑‑TGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGC‑‑CATTGCAGGCGG‑CGCAGGCGCTCACCAAGAATCAGCATGCCcag                                                                                                        >  8:223121/1‑147 (MQ=255)                                        aTAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCA‑CTTGCCACACC‑‑‑‑‑‑‑‑TGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGC‑‑CATTGCAGGCGG‑CGCAGGCGCTCACCAAGAATCAGCATGCCcag                                                                                                        >  6:304861/1‑147 (MQ=255)                                             aGCCGAAGCTAAACGCCATCACCAGTTATACCCACGGCTGGCTGCCCATTTTTCCTCACCCTCCCTTATCCCTCGCCACACCCACGGCAGCCATGCCCACCGCCAGC‑‑CATTGCAGGCGG‑CGCAGGCGCTCACCAAGAATCAGCATGcc                                                                                                           <  1:337503/148‑1 (MQ=255)                                               ccGAAGCGCAACGCCAGCACCAGTGATACCGAGGGCAGGCTTCCCT‑CTTGCCACACC‑‑‑‑‑‑‑‑TGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCATC‑‑CATTGCAGGCGG‑CGCAGGAGCTCTCCTAGAATCAGGCCGGCCAGCAGTACAtt                                                                                               <  4:126347/149‑1 (MQ=255)                                                      tCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCA‑CTTGCCACACC‑‑‑‑‑‑‑‑TGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGC‑‑CATTGCAGGCGG‑CGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGc                                                                                        >  4:310676/1‑149 (MQ=255)                                                      tCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCA‑CTTGCCACACC‑‑‑‑‑‑‑‑TGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGC‑‑CATTGCAGGCGG‑CGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGc                                                                                        >  1:298226/1‑149 (MQ=255)                                                        aaCGCCAGCACCAGTGATACCCACGGCAGGCTGCCCA‑CTTGCCACACC‑‑‑‑‑‑‑‑TGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGC‑‑CATTGCAGGCGG‑CGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCg                                                                                       <  5:304861/148‑1 (MQ=255)                                                        aaCGCCAGCACCAGTGATACCCACGGCAGGCTGCCCA‑CTTGCCACACC‑‑‑‑‑‑‑‑TGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGC‑‑CATTGCAGGCGG‑CGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACACTGATCAGCg                                                                                       <  7:223121/148‑1 (MQ=255)                                                         aCGCCAGCACCAGTGATACCCACGGCAGGCTGCCCA‑CTTGCCACACC‑‑‑‑‑‑‑‑TGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGC‑‑CATTGCAGGCGG‑CGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCggg                                                                                     <  8:315319/149‑1 (MQ=255)                                                              aGCACCAGTGATACCCACGGCAGGCTGCCCA‑CTTGCCACACC‑‑‑‑‑‑‑‑TGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGC‑‑CATTGCAGGCGG‑CGCAGGCGCTCACCAAGAATCAGCATGCCcagca                                                                                                      >  1:162058/1‑127 (MQ=255)                                                                cACCAGTGATACCCACGGCAGGCTGCCCA‑CTTGCCACACC‑‑‑‑‑‑‑‑TGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGC‑‑CATTGCAGGCGG‑CGCAGGCGCTCACCAAGAATCAg                                                                                                                 >  5:165198/1‑114 (MQ=255)                                                                cACCAGTGATACCCACGGCAGGCTGCCCA‑CTTGCCACACC‑‑‑‑‑‑‑‑TGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGC‑‑CATTGCAGGCGG‑CGCAGGCGCTCACCAAGAATCAg                                                                                                                 <  6:165198/114‑1 (MQ=255)                                                                 aCCAGTGATACCCACGGCAGGCTGCCCA‑CTTGCCACACC‑‑‑‑‑‑‑‑TGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGC‑‑CATTGCAGGCGG‑CGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGt                                                                              >  4:179445/1‑148 (MQ=255)                                                                 aCCAGTGATACCCACGGCAGGCTGCCCA‑CTTGCCACACC‑‑‑‑‑‑‑‑TGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGC‑‑CATTGCAGGCGG‑CGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTa                                                                             <  3:310676/149‑1 (MQ=255)                                                                   cAGTGATACCCACGGCAGGCTGCCCA‑CTTGCCACACC‑‑‑‑‑‑‑‑TGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGC‑‑CATTGCAGGCGG‑CGCCGGCGCTCACCAAGACTCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTaa                                                                            >  3:45990/1‑148 (MQ=255)                                                                   cAGTGATACCCACGGCAGGCTGCCCA‑CTTGCCACACC‑‑‑‑‑‑‑‑TGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGC‑‑CATTGCAGGCGG‑CGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTaa                                                                            >  1:167062/1‑148 (MQ=255)                                                                           cccACGGCAGGCTGCCCA‑CTTGCCACACC‑‑‑‑‑‑‑‑TGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGC‑‑CATTGCAGGCGG‑CGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAgg                                                                   >  5:196337/1‑149 (MQ=255)                                                                            ccACGGCAGGCTGCCCA‑CTTGCCACACC‑‑‑‑‑‑‑‑TGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGC‑‑CATTGCAGGCGG‑CGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATc                                                                                           <  4:344975/124‑1 (MQ=255)                                                                            ccACGGCAGGCTGCCCA‑CTTGCCACACC‑‑‑‑‑‑‑‑TGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGC‑‑CATTGCAGGCGG‑CGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATc                                                                                           >  3:344975/1‑124 (MQ=255)                                                                            ccACGGCAGGCTGCCCA‑CTTGCCACACC‑‑‑‑‑‑‑‑TGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGC‑‑CATTGCAGGCGG‑CGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGc                                                                  >  6:265705/1‑149 (MQ=255)                                                                            ccACGGCAGGCTGCCCA‑CTTGCCACACC‑‑‑‑‑‑‑‑TGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGC‑‑CATTGCAGGCGG‑CGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGc                                                                  >  2:365866/1‑149 (MQ=255)                                                                             cACGGCAGGCTGCCCA‑CTTGCCACACC‑‑‑‑‑‑‑‑TGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCGTAATTCCAGGCGG‑CGCAGGCGCGCCCCAAGATTATGCTTGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCAAgc                                                                   >  1:327958/1‑148 (MQ=255)                                                                              aCGGCAGGCTGCCCA‑CTTGCCACACC‑‑‑‑‑‑‑‑TGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGC‑‑CATTGCAGGCGg‑cgc                                                                                                                                     <  6:121720/80‑1 (MQ=255)                                                                              aCGGCAGGCTGCCCA‑CTTGCCACACC‑‑‑‑‑‑‑‑TGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGC‑‑CATTGCAGGCGg‑cgc                                                                                                                                     <  6:115914/80‑1 (MQ=255)                                                                              aCGGCAGGCTGCCCA‑CTTGCCACACC‑‑‑‑‑‑‑‑TGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGC‑‑CATTGCAGGCGg‑cgc                                                                                                                                     >  5:115914/1‑80 (MQ=255)                                                                              aCGGCAGGCTGCCCA‑CTTGCCACACC‑‑‑‑‑‑‑‑TGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGC‑‑CATTGCAGGCGg‑cgc                                                                                                                                     >  5:121720/1‑80 (MQ=255)                                                                                ggCAGGCTGCCCA‑CTTGCCACACC‑‑‑‑‑‑‑‑TGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGC‑‑CATTGCAGGCGG‑CGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGc                                                              <  2:167062/149‑1 (MQ=255)                                                                                  cAGGCTGCCCA‑CTTGCCACACC‑‑‑‑‑‑‑‑TGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGC‑‑CATTGCAGGCGG‑CGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCt                                                            <  8:119354/149‑1 (MQ=255)                                                                                  cAGGCTGCCCA‑CTTGCCACACC‑‑‑‑‑‑‑‑TGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGC‑‑CATTGCAGGCGG‑CGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCt                                                            <  6:196337/149‑1 (MQ=255)                                                                                     gCTGCCCA‑CTTGCCACACC‑‑‑‑‑‑‑‑TGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGC‑‑CATTGCAGGCGG‑CGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTAc                                                                                                  <  7:56708/108‑1 (MQ=255)                                                                                     gCTGCCCA‑CTTGCCACACC‑‑‑‑‑‑‑‑TGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGC‑‑CATTGCAGGCGG‑CGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTAc                                                                                                  >  8:56708/1‑108 (MQ=255)                                                                                                 gCCACACC‑‑‑‑‑‑‑‑TGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGC‑‑CATTGCAGGCGG‑CGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGCGCATAGGGCCg                                              >  8:366560/1‑149 (MQ=255)                                                                                                  ccACACC‑‑‑‑‑‑‑‑TGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGC‑‑CATTGCAGGCGG‑CGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCgt                                             >  4:398034/1‑149 (MQ=255)                                                                                                    acacC‑‑‑‑‑‑‑‑TGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGC‑‑CATTGCAGGCGG‑CGCAGGCGCTCACCAAGAATCAg                                                                                                                 <  2:359621/79‑1 (MQ=255)                                                                                                    acacC‑‑‑‑‑‑‑‑TGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGC‑‑CATTGCAGGCGG‑CGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCgttg                                           >  4:53634/1‑149 (MQ=255)                                                                                                        c‑‑‑‑‑‑‑‑tgctgCGCCACACCCACGGCAGCCATGCCCACCGCCAGC‑‑CATTGCAGGCGG‑CGCAGGCGCTCACCAAGACTCAGCATGCCCAGCAGTACAGTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACTGCCGTTGTTCa                                       >  1:344323/1‑149 (MQ=255)                                                                                                                  gctgcGCCACACCCACGGCAGCCATGCCCACCGCCAGC‑‑CATTGCAGGCGG‑CGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCAcc                                     >  5:293508/1‑149 (MQ=255)                                                                                                                   ctgcGCCACACCCACTGCAGCCATGCCCACCGCCAGC‑‑CATTGCAGGCGG‑CGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCg                                    <  7:51124/149‑1 (MQ=255)                                                                                                                    tgcGCCACACCCACGGCAGCCATGCCCACCGCCAGC‑‑CATTGCAGGCGG‑CGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGa                                   <  3:179445/149‑1 (MQ=255)                                                                                                                      cgcCACACCCACGGCAGCCATGCCCACCGCCAGC‑‑CATTGCAGGCGG‑CGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGAc                                  >  6:213411/1‑148 (MQ=255)                                                                                                                         cacaCCCACGGCAGCCATGCCCACCGCCAGC‑‑CATTGCAGGCGGCCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACca                                >  1:79349/1‑148 (MQ=255)                                                                                                                          acacCCACGGCAGCCATGCCCACCGCCAGC‑‑CATTGCAGGCGG‑CGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCacac                             <  5:265705/149‑1 (MQ=255)                                                                                                                            acCCACGGCAGCCATGCCCACCGCCAGC‑‑CATTGCAGGCGG‑CGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGcc                                                                       <  5:75015/105‑1 (MQ=255)                                                                                                                            acCCACGGCAGCCATGCCCACCGCCAGC‑‑CATTGCAGGCGG‑CGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGcc                                                                       >  6:75015/1‑105 (MQ=255)                                                                                                                              ccACGGCAGCCATGCCCACCGCCAGC‑‑CATTGCAGGCGG‑CGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAg                         >  3:253333/1‑149 (MQ=255)                                                                                                                              ccACGGCAGCCATGCCCACCGCCAGC‑‑CATTGCAGGCGG‑CGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAg                         <  6:293508/149‑1 (MQ=255)                                                                                                                                aCGGCAGCCATGCCCACCGCCAGC‑‑CATTGCAGGCGG‑CGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGt                       >  3:381186/1‑149 (MQ=255)                                                                                                                                         gtGCCCACCGCCAGC‑‑CATTGCAGGCGG‑CGCAGGAGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCAAGGCTTCCCTCTAGCATAAGTACGTTGTTCACCTACAACACTTAGGTCATCCAGtt              >  5:223567/2‑149 (MQ=255)                                                                                                                                         aTGCCTACCGCCAGC‑‑CATTGCAGGCGG‑CGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGGTCCCCGACCCCACGGAGGTCAGCGAGtt              >  8:79927/1‑149 (MQ=255)                                                                                                                                         aTGCCCACCGCCAGC‑‑CATTGCAGGCGG‑CGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGCCAGtt              <  4:381186/149‑1 (MQ=255)                                                                                                                                          tGCCCACCGCCAGC‑‑CATTGAAGGCGG‑CGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGCCAGTTg             >  4:95434/1‑149 (MQ=255)                                                                                                                                                     agc‑‑caTTGCAGGCGG‑CGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGCCAGTTGCCGGCGATCAg  >  3:102131/1‑149 (MQ=255)                                                                                                                                                          |                                                                                                                                                   GGGCAGCGCAGCCACCGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCA‑CTTGCCACACC‑‑‑‑‑‑‑‑TGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGC‑‑CATTGCAGGCGG‑CGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGCCAGTTGCCGGCGATCAG  >  NC_002947/436757‑437044

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