Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,838,341	Δ1 bp	5.1%	coding (230/414 nt)	PP_1643 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,838,341	0	T	.	5.1%	249.8 / 6.7	78	coding (230/414 nt)	PP_1643	membrane protein
Reads supporting (aligned to +/- strand):  ref base T (37/37);  new base . (2/2);  total (39/39)
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.71e-01

AACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCAAGCCAGCCGAATACCGCCCGCGCCGGGTCGAACGCCGCCAGCACGCCCTTGATGAAATAGATATTCACCACGAAGCAGGTCCAGGCAT‑GCGCCCGGGCG‑C‑TGCCTTTCAGCATGCCTGGCAACAGCAACAGCAACGGCACCAGCTCGATCGCCAGGATCACCTCGACCCTGGCCCCATGCAGGTTGGCGAACCACAGGTTGTTCACCACCAGCAGGGCGATCAGGCCGAAGAAG  >  NC_002947/1838194‑1838487                                                                                                                                                    |                                                                                                                                                      aaCTTCCAGCGCACGTACAGTAGCCCGGCAATGCACAGCCCCAGGCTCACCAGCACTTCAAGCCAGCCGAATACCGCCCGCGCCGGGTCGAACGCCGCCAGCACGCCCTTGATGAAATAGATATTCACCACGAAGCAGGTCCAGGCAT‑g                                                                                                                                                     >  4:224753/1‑149 (MQ=255)        aGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCAAGCCAGCCGAATACCGCCCGCGCCGGGTCGAACGCCGCCAGCACGCCCTTGATGAAATAGATATTCACCACGAAGCAGGTCCAGGCAT‑GCGCCCgg                                                                                                                                              <  2:369110/149‑1 (MQ=255)             aCGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCAAGCCAGCCGAATACCGCCCGCGCCGGGTCGAACGCCGCCAGCACGCCCTTGATGAAATAGATATTCACCACGAAGCAGGTCCAGGCAT‑GCGCCCGGGCG‑C‑t                                                                                                                                       <  3:299661/149‑1 (MQ=255)                tACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCAAGCCAGCCGAATACCGCCCGCGCCCGGACGACCCCCGCCAGCACGCCCATGATGAAATAGATATTCACCACGAAGCAGGTCCAAGCAT‑GCGCCCCAGCG‑C‑AGc                                                                                                                                     >  5:280480/1‑148 (MQ=255)                     tAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCAAGCCAGCCGAATACCGCCCGCGCCGGGTCGAACGCCGCCAGCACGCCCTTGATGAAATAGATATTCACCACGAAGCAGGTCCAGGCAT‑GCGCCCGGGCG‑C‑TGCCTTTCa                                                                                                                               <  5:184062/149‑1 (MQ=255)                         ccGGCAATGAACAGCCCCAGGCTCACCAGCACTTCAAGCCAGCCGAATACCGCCCGCGCCGGGTCGAACGCCGCCAGCACGCCCTTGATGAAATAGATATTCACCACGAAGCAGGTCCAGGCAT‑GCGCCCGGGCG‑C‑TGCCTTTCAGCa                                                                                                                            >  8:204974/1‑148 (MQ=255)                          cGGCAATGAACAGCCCCAGGCTCACCAGCACTTCAAGCCAGCCGAATACCGCCCGCGCCGGGTCGAACGCCGCCAGCACGCCCTTGATGAAATAGATATTCACCACGAAGCAGGTCCAGGCAT‑GCGCCCGGGCG‑C‑TGCCTTTCAGCATg                                                                                                                          >  2:61927/1‑149 (MQ=255)                          cGGCAATGAACAGCCCCAGGCTC‑CCAGCACTTCAAGCCAGCCGAATACCGCCCGCGCCGGGTCGAACGCCGCCAGCACGCCCTTGATGAAATAGATATTCACCACGAAGCAGGTCCAGGCAT‑GCGCCCGGGCG‑C‑TGCCTTTCAGCATGc                                                                                                                         >  2:269397/1‑149 (MQ=255)                            gCAATGAACAGCCCCAGGCTCACCAGCACTTCAAGCCAGCCGAATACCGCCCGCGCCGGGTCGAACGCCGCCAGCACGCCCTTGATGAAATAGATATTCACCACGAAGCAGGTCCAGGCAT‑GCGCCCGGGCG‑C‑TGCCTTTCAGCATGcc                                                                                                                        <  5:30755/149‑1 (MQ=255)                              aaTGAAGAGCCCCAGGCTCACCAGCACTTCATGGCATCCGAATACCTCCCGCGCCGGGTCGAACGGCGCCAGGACGCGCTTGATGAAATAGATATTCACCACGATGCAGGTCCAGGCAT‑GCGCCCGGGCG‑C‑TGCCTTTCAGCATGCCTg                                                                                                                      <  6:280480/149‑1 (MQ=255)                                tGAACAGCCCCAGGCTCACCAGCACTTCAAGCCAGCCGAATACCGCCCGCGCCGGGTCGAACGCCGCCAGCACGCCCTTGATGAAATAGATATTCACCACGAAGCAGGTCCAGGCAT‑GCGCCCGGGCG‑C‑TGCCTTTCAGCATGCCTGgc                                                                                                                    <  6:75293/149‑1 (MQ=255)                                 gAACACCCCCAGCCTCACCAGCACTTCAAGCCAGCCGAATACCGCCCGCGCCGGGACGAACACCGCCAGCACGCCCTTGATGAAATAGAAATTCACCACGAAGCAGGACCAGGCAT‑GCACCCTCACG‑C‑TGCCTAACAGCAAGCCAGgc                                                                                                                    >  2:282264/1‑148 (MQ=255)                                   gcAGCCCCAGGCTCACCAGCACTTCAAGCCAGCCGAATACCGCCCGCGCCGGGTCGAACGCCGCCAGCACGCCCTTGATGAAATAGATATTCACGACGAAGCAGGTCCAGGCAT‑GCGCCCGGGCG‑C‑TGCCTTTCAGCATGCCTGgcaa                                                                                                                  >  3:8530/2‑148 (MQ=255)                                   aCAGCCCCAGGCTCACCAGCACTTCAAGCTAGCCGAATACCGCCCGCGCCGGGTCGAACGCCGCCAGCACGCCCTTGATGAAATAGATATTCACCACGAAGCAGGTCCAGGCAT‑GCGCCCGGGCG‑C‑TGCCTTTCAGCATGCCTGgcaac                                                                                                                 <  3:250446/149‑1 (MQ=255)                                   aCAGCCCCAGGCTCACCAGCACTTCAAGCCAGCCGAATACCGCCCGCGCCGGGTCGAACGCCGCCAGCACGCCCTTGATGAAATAGATATTCACCACGAAGCAGGTCCAGGCAT‑GCGCCCGGGCG‑C‑TGCCTTTCAGCATGCCTGgcaa                                                                                                                  >  4:188413/1‑148 (MQ=255)                                    cAGCCCCAGGCTCACCAGCACTTCAAGCCAGCCGAATACCGCCCGCGCCGGGTCGAACGCCGCCAGCACGCCCTTGATGAAATAGATATTCACCACGAAGCAGGTCCAGGCAT‑GCGCCCGGGCG‑C‑TGCCTTTCAGCATGCCTGgcaaca                                                                                                                <  2:4097/149‑1 (MQ=255)                                       ccccAGGCTCACCAGCACTTCAAGCCAGCCGAATACCGCCCGCGCCGGGTCGAACGCCGCCAGCACGCCCTTGATGAAATAGATATTCACCACGAAGCAGGTCCAGGCAT‑GCGCCCGGGCG‑C‑TGCCATTCAGCATGCCTGgcaacagc                                                                                                              >  6:197196/1‑148 (MQ=255)                                        cccAGGCTCACCAGCACTTCAAGCCAGCCGAATACCGCCCGCGCCGGGTCGAACGCCGCCAGCACGCCCTTGATGAAATAGATATTCACCACGAAGCAGGTCCAGGCAT‑GCGCCCGGGCG‑C‑TGCCTCTCAGCATGCCTGgcaacagcaa                                                                                                            <  7:57133/149‑1 (MQ=255)                                         ccAGGCTCACGAGCACTTCAAGCCAGCCGAATACCGCCCGCGCCGGGTCGAACGCCGCCAGCACGCCCTTGATGAAATAGATATTCACCACGAAGCAGGTCCAGGCAT‑GCGCCCGGGCG‑C‑TGCCTTTCAGCATGCCTGgcaacagcaac                                                                                                           <  3:224753/149‑1 (MQ=255)                                         ccAGGCTCACCAGCACTTCAAGCCAGCCGAATACCGCCCGCGCCGGGTCGAACGCCGCCAGCACGCCCTTGATGAAATAGATATTCACCACGAAGCAGGTCCAGGCAT‑GCGCCCGGGCG‑C‑TGCCTTTCAGCATGCCTGgcaacagcaa                                                                                                            >  5:154934/1‑148 (MQ=255)                                         ccAGGCTCACCAGCACTTCAAGCCAGCCGAATACCGCCCGCGCCGGGTCGAACGCCGCCAGCACGCCCTTGATGAAATAGATATTCACCACGAAGCAGGTCCAGGCA‑‑GCGCCCGGGCG‑CATGCCt                                                                                                                                   >  2:309762/1‑125 (MQ=255)                                         ccAGGCTCACCAGCACTTCAAGCCAGCCGAATACCGCCCGCGCCGGGTCGAACGCCGCCAGCACGCCCTTGATGAAATAGATATTCACCACGAAGCAGGTCCAGGCA‑‑GCGCCCGGGCG‑CATGCCt                                                                                                                                   <  1:309762/125‑1 (MQ=255)                                          cAGGCTCACCAGCACTTCAAGCCAGCCGAATACCGCCCGCGCCGGGTCGAACGCCGCCAGCACGCCCTTGATGAAATAGATATTCACCACGAAGCAGGTCCAGGCAT‑GCGCCCGGGCG‑C‑TGCCTTTCa                                                                                                                               <  8:310375/128‑1 (MQ=255)                                           aGGCTCACCAGCACTTCAAGCCAGCCGAATACCGCCCGCGCCGGGTCGAACGCCGCCAGCACGCCCTTGATGAAATAGATATTCACCACGAAGCAGGTCCAGGCAT‑GCGCCCGGgc                                                                                                                                            >  4:298127/1‑116 (MQ=255)                                           aGGCTCACCAGCACTTCAAGCCAGCCGAATACCGCCCGCGCCGGGTCGAACGCCGCCAGCACGCCCTTGATGAAATAGATATTCACCACGAAGCAGGTCCAGGCAT‑GCGCCCGGgc                                                                                                                                            <  3:298127/116‑1 (MQ=255)                                                 aCCAGCACTTCAAGCCAGCCGAATACCGCCCGCGCCGGGTCGAACGCCGCCAGCACGCCCTTGATGAAATAGATATTCACCACGAAGCAGGTCCAGGCAT‑GCGCCCGGGCG‑C‑TGCCTTTCAGCATGCCTGGCAACAGCAACAGCAACg                                                                                                    <  1:61927/148‑1 (MQ=255)                                                 aCCAGCACTTCAAGCCAGCCGAATACCGCCCGCGCCGGGTCGAACGCCGCCAGCACGCCCTTGATGAAATAGATATTCACCACGAAGCAGGTCCAGGCA‑‑GCGCCCGGGCgcc‑tgcct                                                                                                                                   >  1:44554/1‑117 (MQ=255)                                                 aCCAGCACTTCAAGCCAGCCGAATACCGCCCGCGCCGGGTCGAACGCCGCCAGCACGCCCTTGATGAAATAGATATTCACCACGAAGCAGGTCCAGGCA‑‑GCGCCCGGGCG‑CATGCCt                                                                                                                                   <  2:44554/117‑1 (MQ=255)                                                         ttCAAGCTAGCCGAATACCGCCCGCGCCGGGTCGAACGCCGCCAGCACGCCCTTGATGAAATAGATATTCACCACGAAGCAGGTCCAGGCAT‑GCGCCCGGGCG‑C‑TGCCTTTCAGCATGCCTGgcaacagcaa                                                                                                            >  1:26661/1‑132 (MQ=255)                                                         ttCAAGCTAGCCGAATACCGCCCGCGCCGGGTCGAACGCCGCCAGCACGCCCTTGATGAAATAGATATTCACCACGAAGCAGGTCCAGGCAT‑GCGCCCGGGCG‑C‑TGCCTTTCAGCATGCCTGgcaacagcaa                                                                                                            <  2:26661/132‑1 (MQ=255)                                                         ttCAAGCCAGCCGAATACCGCCCGCGCCGGGTCGAACGCCGCCAGCACGCCCTTGATGAAATAGATATTCACCACGAAGCAGGTCCAGGCAT‑GCGCCCGGGCG‑C‑TGCCTTTCAGCATGCCTGGCAACAGCAACAGCAACGGCACCAGCt                                                                                           >  3:60556/1‑149 (MQ=255)                                                            aaGCCAGCCGAATAGCGCCCGCGCCGGGTCGAACGCCGCCAGCACGCCCTTGATGAAATAGATATTCACCACGAAGCAGGTCCATGCAT‑GCGCCTGGGCG‑C‑TGCCTTTCAGCATGCCTGGCAACAGCAACAGCAACGGCACCAGCTCg                                                                                         <  1:282264/148‑1 (MQ=255)                                                            aaGCCAGCCGAATACCGCCCGCGCCGGGTCGAACGCCGCCAGCACGCCCTTGATGAAATAGATATTCACCACGAAGCAGGTCCAGGCAT‑GCGCCCGGGCG‑C‑TGCCTTTc                                                                                                                                >  6:29281/1‑109 (MQ=255)                                                            aaGCCAGCCGAATACCGCCCGCGCCGGGTCGAACGCCGCCAGCACGCCCTTGATGAAATAGATATTCACCACGAAGCAGGTCCAGGCAT‑GCGCCCGGGCG‑C‑TGCCTTTc                                                                                                                                <  5:29281/109‑1 (MQ=255)                                                             agccatccGAATACCGCCCGCGCCGGGTCGAACGCCGCCATCAAGCCCTTGATGAAATACAAATTCACCACGAAGCATGTCCAGGCATGGCGCCCGGGCG‑C‑TGCCTTTCAGCATGCCTGGCAACAGCAACAGCAACGGCACCAGCTCGa                                                                                        <  1:291048/149‑1 (MQ=255)                                                              gccagccGAATACCGCCCGCGCCGGGTCGAACGCCGCCAGCACGCCCTTGATGAAATAGATATTCACCACGAAGCAGGTCCAGGCAT‑GCGCCCGGGCG‑C‑TGCCTTTCAGCATGCCTGGCAACAGCAACAGCAACGGCACCAGCTCGATc                                                                                      <  5:197196/149‑1 (MQ=255)                                                                   ccGAATACCGCCCGCGCCGGGTCGAACGCCGCCAGCACGCCCTTGATGAAATAGATATTCACCACGAAGCAGGTCCAGGCAT‑GCGCCCGGGCG‑C‑TGCCTTTCAGCATGCCTGGCAACAGCAACAGCAACGGCACCAGCTCGATc                                                                                      >  5:226687/1‑144 (MQ=255)                                                                   ccGAATACCGCCCGCGCCGGGTCGAACGCCGCCAGCACGCCCTTGATGAAATAGATATTCACCACGAAGCAGGTCCAGGCAT‑GCGCCCGGGCG‑C‑TGCCTTTCAGCATGCCTGGCAACAGCAACAGCAACGGCACCAGCTCGATc                                                                                      <  6:226687/144‑1 (MQ=255)                                                                   ccGAATACCGCCCGCGCCGGGTCGAACGCCGCCAGCACGCCCTTGATGAAATAGATATTCACCACGAAGCAGGTCCAGGCAT‑GCGCCCGGGCG‑C‑TGCCTTTCAGCATGCCTGGCAACAGCAACAGCAACGGCACCAGCTCGATCGCCAg                                                                                 >  7:15905/1‑149 (MQ=255)                                                                    cGAATACCGCCCGCGCCGGGTCGAACGCCGCCAGCACGCCCTTGATGAAATAGATATTCACCACGAAGCAGGTCCAGGCAT‑GCGCCCGGGCG‑C‑TGCCTTTCAGCATGCCTGGCAACAGCAACAGCAACGGCACCAGCTCGATCGCCAg                                                                                 >  7:326153/1‑148 (MQ=255)                                                                          ccgcccgcGCCGGGTCGAACGCCGCCAGCACGCCCTTGATGAAATAGATATTCACCACGAAGCAGGTCCAGGCAT‑GCGCCCGGGCG‑C‑TGCCTTTCAGCATGCCTGGCAACAGCAACAGCAACGGCACCAGCTCGATCGCCAGGATCAcc                                                                          <  6:154934/149‑1 (MQ=255)                                                                           cgcccgcGCCGGGTCGAACGCCGCCAGCACGCCCTTGATGAAATAGATATTCACCCCGAAGCAGGTCCAGGCAT‑GCGCCCGGGCG‑C‑TGCCTTTCAGCATGCCTGGCAACAGCAACAGCAACGGCACCAGCTCGATCGCCAGGATCACCt                                                                         >  3:294938/1‑149 (MQ=255)                                                                           cgcccgcGCCGGGTCGAACGCCGCCAGCACGCCCTTGATGAAATAGATATTCACCACGAAGCAGGTCCAGGCAT‑GCGCCCGGGCG‑C‑TGCCTTTCAGCATGCCTGGCAACAGCAACAGCAACGGCACCAGCTCGATCGCCAGGATCACCt                                                                         >  8:217868/1‑149 (MQ=255)                                                                            gcccgcGCCGGGTCGAACGCCGCCAGCACGCCCTTGATGAAATAGATATTCACCACGAAGCAGGTCCAGGCAT‑GCGCCCGGGCG‑C‑TGCCTTTCAGCATGCCTGGCAACAGCAACAGCAACGGCACCAGCTCGATc                                                                                      >  3:168137/1‑135 (MQ=255)                                                                            gcccgcGCCGGGTCGAACGCCGCCAGCACGCCCTTGATGAAATAGATATTCACCACGAAGCAGGTCCAGGCAT‑GCGCCCGGGCG‑C‑TGCCTTTCAGCATGCCTGGCAACAGCAACAGCAACGGCACCAGCTCGATc                                                                                      <  4:168137/135‑1 (MQ=255)                                                                             cccgcGCCGGGTCGAACGCCGCCAGCACGCCCTTGATGAAATAGATATTCACCACGAAGCAGGTCCAGGCAT‑GCGCCCGGGCG‑C‑TGCCTTTCAGCATGCCTGgcaacagcaa                                                                                                            >  8:187034/1‑112 (MQ=255)                                                                             cccgcGCCGGGTCGAACGCCGCCAGCACGCCCTTGATGAAATAGATATTCACCACGAAGCAGGTCCAGGCAT‑GCGCCCGGGCG‑C‑TGCCTTTCAGCATGCCTGgcaacagcaa                                                                                                            <  7:187034/112‑1 (MQ=255)                                                                              ccgcGCCGGGTCGAACGCCGCCAGCACGCCCTTGATGAAATAGATATTCACCACGAAGCAGGTCCAGGCAT‑GCGCCCGGGCG‑C‑TGCCTTTCAGCATGCCTGGCAACAGCAACAGCAACGGCACCAGCTCGATCGCCAGGATCACCTCGa                                                                      >  1:152360/1‑149 (MQ=255)                                                                               cgcgCCGGGTCGAACGCCGCCAGCACGCCCTTGATGAAATAGATATTCACCACGAAGCAGGTCCAGGCAT‑GCGCCCGGGCG‑C‑TGCCTTTCAGCATGCCTGGCAACAGCAACAGCAACGGCACCAGCTCGATCGCCAGGATCACCTCGAc                                                                     <  4:60556/149‑1 (MQ=255)                                                                                 cgcCGGGTCGAACGCCGCCAGCACGCCCCTGATGAAATAGATATTCACCACGAAGCAGGTCCAGGCAT‑GCGCCCGGGCG‑C‑TGCCTTTCAGCATGACTGGCAACAGCAACAGCAACTGCAACAGCTCGATCGCCACTATAACCTCGa                                                                      >  6:360419/1‑146 (MQ=255)                                                                                 agcCTGGTCGAACGCCGCCAGCACGCCCATGATGAAATAGATATTCACCACGAAGCAGGTCCAGGCAT‑GCGCCCGGGCG‑C‑TGCCTTTCAGCATGCCTGGCAACAGCAACAGCAACGGCACCAGCTCGATCGGCAGGATCACCTCGa                                                                      <  5:360419/145‑1 (MQ=255)                                                                                           aaCGCCGCCAGCACGCCCTTGATGAAATAGATATTCACCACGAAGCAGGTCCAGGCAT‑GCGCCCGGGCG‑C‑TGCCTTTCAGCATGCCTGGCAACAGCAACAGCAACGGCACCAGCTCGATCGCCAGGATCACCTCGACCCTGGCCCCATg                                                         <  5:355324/149‑1 (MQ=255)                                                                                                  ccAGCACGCCCTTGATGAAATAGATATTCACCACGAAGCAGGTCCAGGCAT‑GCGCCCGGGCG‑C‑TGCCTTTCAGCATGCCTGGCAACAGCAACAGCAACGGCACCAGCTCGATCGCCAGGATCACCTCGACCCTGGCCCCATGCAGGTTg                                                  >  3:321925/1‑149 (MQ=255)                                                                                                   cAGCACGCCCTTGATGAAATAGATATTCACCACGAAGCAGGTCCAGGCAT‑GCGCCCGGGCG‑C‑TGCCTTTCAGCATGCCTGGCAACAGCAACAGCAACGGCACCAGCTCGATCGCCAGGATCACCTCGACCCTGGCCCCATGCAGGTTgg                                                 <  2:152360/149‑1 (MQ=255)                                                                                                     gCACGCCCTTGATGAAATAGATATTCACCACGAAGCAGGTCCAGGCAT‑GCGCCCGGGCG‑C‑TGCCTTTCAGCATGCCTGGCAACAGCAACAGCAACGGCACCAGCTCGATCGCCAGGATCACCTCGACCCTGGCCCCATGCAGGTTGGc                                                >  8:19187/1‑148 (MQ=255)                                                                                                       aCGCCCTTGATGAAATAGATATTCACCACGAAGCAGGTCCAGGCAT‑GCGCCCGGGCG‑C‑TGCCTTTCAGCATGCCTGGCAACAGCAACAGCAACGGCACCAGCTCGATCGCCAGGATCACCTCGACCCTGGCCCCATGCAGGTTGGCGaa                                             <  3:188413/149‑1 (MQ=255)                                                                                                         gCCCTTGATGAAATAGATATTCACCACGAAGCAGGTCCAGGCAT‑GCGCCCGGGCG‑C‑TGCCTTTCAGCATGCCTGGCAACAGCAACAGCAACGGCACCAGCTCGATCGCCAGGATCACCTCGACCCTGGCCCCATGCAGGTTGGCGAAcc                                           <  7:204974/149‑1 (MQ=255)                                                                                                                 tgaAATAGATATTCACCACGAAGCAGGTCCAGGCAT‑GCGCCCGGGCG‑C‑TGCCTTTCAGCATGCCTGGCAACAGCAACAGCAACGGCACCAGCTCGATCGCCAGGATCACCTCGACCCTGGCCCCATGCAGGTTGGCGAACCACAGgttg                                   <  8:15905/149‑1 (MQ=255)                                                                                                                  gaAATAGATATTCACCACGAAGCAGGTCCAGGCAT‑GCGCCCGGGCG‑C‑TGCCTTTCAGCATGCCTGGCAACAGCAACACCAACGGCACCAGCTCGATCGCCAGGCTCACCTCGCCCCTGGCCCCATGCAGGTTGGCGAACCACAGgtagt                                  >  5:10250/1‑149 (MQ=255)                                                                                                                   aaaTAGATATTCACCACGAAGCAGGTCCAGGCAT‑GCGCCCGGGCG‑C‑TGCCTTTCAGCATGCCTGGCAACAGCAACAGCAACGGCACCAGCTCGATCGCCAGGATCACCTCGACCCTGGCCCCATGCAGGTTGGCGAACCACAGgttgtt                                 <  7:217868/149‑1 (MQ=255)                                                                                                                    aaTAGATATTCACCACGAAGCAGGTCCAGGCAT‑GCGCCCGGGCG‑C‑TGCCTTTCAGCATGCCTGGCAACAGCAACAGCAACGGCACCAGCTCGATCGCCAGGATCACCTCGACCCTGGCCCCATGCAGGTTGGCGAACCACAGGTTGTTc                                <  7:19187/149‑1 (MQ=255)                                                                                                                               accacGAAGCAGGTCCAGGCAT‑GCGCCCGGGCG‑C‑TGCCTTTCAGCATGCCTGGCAACAGCAACAGCAACGGCACCAGCTCGATCGCCAGGATCACCTCGACCCTGGCCCCATGCAGGTTGGCGAACCACAGGTTGTTCACCACcagcag                     >  2:113247/1‑149 (MQ=255)                                                                                                                               accacGAAGCAGGTCCAGGCAT‑GCGCCCGGGCG‑C‑TGCCTTTCAGCATGCCTGGCAACAGCAACAGCAACGGCACCAGCTCGATCGCCAGGATCACCTCGACCCTGGCCCCATGCAGGTTGGCGAACCACAGGTTGTTCACCACcagcag                     >  8:215778/1‑149 (MQ=255)                                                                                                                               accacGAAGCAGGTCCAGGCAT‑GCGCCCGGGCG‑C‑TGCCTTTCAGCATGCCTGGCAACAGCAACAGCAACGGCACCAGCTCGATCGCCAGGATCACCTCGACCCTGGCCCCATGCAGGTTGGCGAACCACAGGTTGTTCACCACcagcag                     >  7:291121/1‑149 (MQ=255)                                                                                                                               accacGAAGCAGGTCCAGGCAT‑GCACCCGGGCG‑C‑TGCCTTTCAGCATGCCTGGCAACAGCAACAGCAACGGCACCAGCTCGATCGCCAGGATCACCTCGACCCTGGCCCCATGCAGGTTGGCGAACCACAGGTTGTTcac                              <  7:173348/140‑1 (MQ=255)                                                                                                                               accacGAAGCAGGTCCAGGCAT‑GCACCCGGGCG‑C‑TGCCTTTCAGCATGCCTGGCAACAGCAACAGCAACGGCACCAGCTCGATCGCCAGGATCACCTCGACCCTGGCCCCATGCAGGTTGGCGAACCACAGGTTGTTcac                              >  8:173348/1‑140 (MQ=255)                                                                                                                                 cacGAAGCAGGTCCAGGCAT‑GCGCCCGGGCG‑C‑TGCCTTTCAGCATGCCTGGCAACAGCAACAGCAACGGCACCAGCTCGATCGCCAGGATCACCTCGACCCTGGCCCCATGCAGGTTGGCGAACCACAGGTTGTTCACCACCAGCAggg                   <  4:8530/149‑1 (MQ=255)                                                                                                                                    gAAGCAGGTCCAGGCAT‑GCGCCCGGGCG‑C‑TGCCTTTCAGCATGCCTGGCAACAGCAACAGCAACGGCACCAGCTCGATCGCCAGGATCACCTCGACCCTGGCCCCATGCAGGTTGGCGAACCACAGGTTGTTCACCACCAGCAGGGCGa                >  3:163560/1‑149 (MQ=255)                                                                                                                                    gAAGCAGGTCCAGGCAT‑GCGCCCGGGCG‑C‑TGCCTTTCAGCATGCCTGGCAACAGCAACAGCAACGGCACCAGCTCGATCGCCAGGATCACCTCGACCCTGGCCCCATGCAGGTTGGCGAACCACAGGTTGTTCACCACCAGCAGGGCGa                >  1:356081/1‑149 (MQ=255)                                                                                                                                       gCAGGTCCAGGCAT‑GCGCCCGGGCG‑C‑TGCCTTTCAGCATGCCTGGCAACAGCAACAGCAACGGCACCAGCTCGATCGCCAGGATCACCTCGACCCTGGCCCCATGCAGGTTGGCGAACCACAGGTTGTTCACCACCAGCAGGGCGAt               >  4:21613/1‑147 (MQ=255)                                                                                                                                             ccAGGCAT‑GCGCCCGGGCG‑C‑TGCCTTTCAGCATGCCTGGTAACAGCAACAGCAACGGCACCAGCTCGATCGCCAGGATCACCTCGACCCTGGCCCCATGCAGGTTGGCGAACCACAGGTTGTTCACCACCAGCAGGGCGa                <  5:140232/140‑1 (MQ=255)                                                                                                                                             ccAGGCAT‑GCGCCCGGGCG‑C‑TGCCTTTCAGCATGCCTGGTAACAGCAACAGCAACGGCACCAGCTCGATCGCCAGGATCACCTCGACCCTGGCCCCATGCAGGTTGGCGAACCACAGGTTGTTCACCACCAGCAGGGCGa                >  6:140232/1‑140 (MQ=255)                                                                                                                                             ccAGGCAT‑GCGCCCGGGCG‑C‑TGCCTTTCAGCATGCCTGGCAACAGCAACAGCAACGGCACCAGCTCGATCGCCAGGATCACCTCGACCCTGGCCCCATGCAGGTTGGCGAACCACAGGTTGTTCACCACCAGCAGGGCGATCAGGCCg        >  6:279786/1‑148 (MQ=255)                                                                                                                                               aGGCAT‑GCGCCCGGGCG‑C‑TGCCTTTCAGCATGCCTGGCAACAGCAACAGCAACGGCACCAGCTCGATCGCCAGGATCACCTCGACCCTGGCCCCATGCAGGTTGGCGAACCACAGGTTGTTCACCACCAGCAGGGCGATCAGGCCg        >  7:105115/1‑146 (MQ=255)                                                                                                                                               aGGCAT‑GCGCCCGGGCG‑C‑TGCCTTTCAGCATGCCTGGCAACAGCAACAGCAACGGCACCAGCTCGATCGCCAGGATCACCTCGACCCTGGCCCCATGCAGGTTGGCGAACCACAGGTTGTTCACCACCAGCAGGGCGATCAGGCCg        >  5:374531/1‑146 (MQ=255)                                                                                                                                                 gCAT‑GCGCCCGGGCG‑C‑TGCCTTTCAGCATGCCTGGCAACAGCAACAGCAACGGCACCAGCTCGATCGCCAGGATCACCTCGACCCTGGCCCCATGCAGGTTGGCGAACCACAGGTTGTTCACCACCAGCAGGGCGATCAGGCCgaagaa   <  1:269397/149‑1 (MQ=255)                                                                                                                                                 gCAT‑GCGCCCGGGCG‑C‑TGCCTTTCAGCATGCCTGGCAACAGCAACAGCAACGGCACCAGCTCGATCGCCAGGATCACCTCGACCCTGGCCCCATGCAGGTTGGCGAACCACAGGTTGTTCACCACCAGCAGGGCGATCAGGCCgaa      <  8:105115/146‑1 (MQ=255)                                                                                                                                                  cAT‑GCGCCCGGGCG‑C‑TGCCTTTCAGCATGCCTGGCAACAGCAACAGCAACGGCACCAGCTCGATCGCCAGGATCACCTCGACCCTGGCCCCATGCAGGTTGGCGAACCACAGGTTGTTCACCACCAGCAGGGCGATCAGGCCGaagaag  <  8:326153/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                      AACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCAAGCCAGCCGAATACCGCCCGCGCCGGGTCGAACGCCGCCAGCACGCCCTTGATGAAATAGATATTCACCACGAAGCAGGTCCAGGCAT‑GCGCCCGGGCG‑C‑TGCCTTTCAGCATGCCTGGCAACAGCAACAGCAACGGCACCAGCTCGATCGCCAGGATCACCTCGACCCTGGCCCCATGCAGGTTGGCGAACCACAGGTTGTTCACCACCAGCAGGGCGATCAGGCCGAAGAAG  >  NC_002947/1838194‑1838487

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