Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,277,853	T→G	7.0%	D104A (GAC→GCC)	ampC ←	beta‑lactamase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,277,853	0	T	G	7.0%	212.6 / 5.5	80	D104A (GAC→GCC)	ampC	beta‑lactamase
Reads supporting (aligned to +/- strand):  ref base T (36/38);  new base G (4/2);  total (40/40)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

ACTACCTGTTGCGGTGTTTGTACGGTGTCGGGGAACTGCAGTGGCAGGCAGTCGGCGCTATAGGCGCCCAGTTCCAGCACCGTCGCGTCGCCGATGGGGACGCCGGCCAATGCCGGGTGATAGCGCTTGGCCGGGGCCTCGAGGTCAAGCTTGCCCTCGGCG‑CTGGCAAGCGCGG‑CCAGGGTGGCCGTGTAGGTCTTGCTCAGCGAACCGATCTCGAACAGCGTATCGCGGCTGACCGGTACGTTGTCGGCCTTGCTGGCCACACCGTAGCTGAAGTAATGCGCGTGGTCGCGT  >  NC_002947/3277709‑3278002                                                                                                                                                 |                                                                                                                                                        aCTACCTGTTGCGGTGTTTGTACGGTGTCGGGGAACTGCAGTGGCAGGCAGTCGGCGCTATAGGCGCCCAGTTCCAGCACCGTCGCGTCGCCGATGGGGACGCCGGCCAATGCCGGGTGATAGCGCTTGGCCGGGGCCTCGAGGTCAAGc                                                                                                                                                    <  1:616921/150‑1 (MQ=255)   tACCTGTTGCGGTGTTTGTACGGTGTCGGGGAACTGCAGTGGCAGGCAGTCGGCGCTATAGGCGCCCAGTTCCAGCACCGTCGCGTCGCCGATGGGGACGCCGGCCAATGCCGGGTGATAGCGCTTGGCCGGGGCCTCGAGGTCAAGCtt                                                                                                                                                  >  2:1942777/1‑150 (MQ=255)   tACCTGTTGCGGTGTTTGTACGGTGTCGGGGAACTGCAGTGGCAGGCAGTCGGCGCTATAGGCGCCCAGTTCCAGCACCGTCGCGTCGCCGATGGGGACGCCGGCCAATGCCGGGTGATAGCGCTTGGCCGGGGCCTCGAGGTCAAGCtt                                                                                                                                                  >  2:2027087/1‑150 (MQ=255)   tACCTGTTGCGGTGTTTGTACGGTGTAGGGGAACTGCAGTGGCAGGCAGTCGGCGCTATAGGCGCCCAGTTCCAGCACCGTCGCGTCGCCGATGGGGACGCCGGCCAATGCCGGGTGATAGCGCTTGGCCGGGGCCTCGAGGTCAAGCtt                                                                                                                                                  >  1:427819/1‑150 (MQ=255)     ccTGTTGCGGTGTTTGTACGGTGTCGGGGAACTGCAGTGGCAGGCAGTCGGCGCTATAGGCGCCCAGTTCCAGCACCGTCGCGTCGCCGATGGGGACGCCGGCCAATGCCGGGTGATAGCGCTTGGCCGGGGCCTCGAGGTCAAGCTTGc                                                                                                                                                >  1:1857672/1‑150 (MQ=255)      cTGTTGCGGTGTTTGTACGGTGTCGGGGAACTGCAGTGGCAGGCAGTCGGCGCTATAGGCGCCCAGTTCCAGCACCGTCACTTGGACGATGGGGACGCCGGCCAATGCCGGGTGATCGCGCTTGGCCGGGGCCTCGAGGTCAAGCTTGcc                                                                                                                                               >  1:986876/1‑150 (MQ=255)             ggTGTTTGTACGGTGTGGGGGAACTGCAGTGGCAGGCAGTCGGCGCTATAGGCGCCCAGTTCCAGCACCGTCGCGTCGCCGATGGGGACGCCGGCCAATGCCGGGTGATAGCGCTTGGCCGGGGCCTCGAGGTCAAGCTTGCCCTCGgcg                                                                                                                                        >  1:1152304/1‑150 (MQ=255)              gtgtTTGTACGGTGTCGGGGAACTGCAGTGGCAGGCAGTCGGCGCTATAGGCGCCCAGTTCCAGCACCGTCGCGTCGCCGATGGGGACGCCGGCCAATGCCGGGTGATAGCGCTTGGCCGGGGCCTCGAGGTCAAGCTTGCCCTCGgcg‑c                                                                                                                                      >  2:775058/1‑150 (MQ=255)                gtTTGTACGGTGTCGGGGAACTGCAGTGGCAGGGAGTCGGGGGTATAGTCGCCCAGTGCCAGCACCGTCGCGTGGCCGCTGGGGACGCCGGCCAATGCCTGGTGATAGCGCTTGGCCGGGGGCTCGAGGTCAAGCTTGCCCTCGGCG‑CTg                                                                                                                                    <  2:175013/150‑1 (MQ=255)                gtTTGTACGGTGTCGGGGAACTGCAGTGGCAGGCAGTCGGCGCTATAGGCGCCCAGTTCCAGCACCGTCGCGTCGCCGATGGGGACGCCGGCCAATGCCGGGTGATAGCGCTTGGCCGGGGCCTCGAGGTCAAGCTTGCCCTCGGCG‑CTg                                                                                                                                    <  1:89655/150‑1 (MQ=255)                    gTACGGTGTCGGGGAACTGCAGTGGCAGGCAGTCGGCGCTATAGGCGCCCAGTTCCAGCACCGTCGCGTCGCCGATGGGGACGCCGGCCAATGCCGGGTGATAGCGCTTGGCCGGGGCCTCGAGGTCAAGCTTGCCCTCGGCG‑CTGGCaa                                                                                                                                >  2:2285948/1‑150 (MQ=255)                            tCGAGGAACCGCAGTGGCAGGCAGTCGCCGCCATAGGCGCCCAGTTCCAGCACCGTCGCGTCGCCGATGGGGACGCCGGCCAATGCCGGGTGATAGCGCTTGGCCGGGGCCTCGAGCTCAAGCTTGCCCTCGCCG‑CTTGCAAGCgccg‑cc                                                                                                                       >  2:1100156/1‑150 (MQ=255)                                 gAACTGCAGTGGCAGGCAGTCGGCGCTATAGGCGCCCAGTTCCAGCACCGTCGCGTCGCCGATGGGGACGCCGGCCAATGCCGGGTGATAGCGCTTGGCCGGGGCCTCGAGGGCAAGCTTg                                                                                                                                                 <  1:1460588/121‑1 (MQ=255)                                 gAACTGCAGTGGCAGGCAGTCGGCGCTATAGGCGCCCAGTTCCAGCACCGTCGCGTCGCCGATGGGGACGCCGGCCAATGCCGGGTGATAGCGCTTGGCCGGGGCCTCGAGGGCAAGCTTg                                                                                                                                                 >  2:1460588/1‑121 (MQ=255)                                     tGCAGTGGCAGGCAGTCGGCGCTATAGGCGCCCAGTTCCAGCACCGTTGCGTCGCCGATGGGGACGCCGGCCAATGCCGGGTGATAGCGCTTGGCCGGGGCCTCGAGGTCAAGCTTGCCCTCGGCG‑CTGGCAAGCGCGG‑CCAGGGTGGcc                                                                                                              <  1:1789485/150‑1 (MQ=255)                                     tGCAGTGGCAGGCAGTCGGCGCTATAGGCGCCCAGTTCCAGCACCGTCGCGTCGCCGATGGGGACGCCGGCCAATGCCGGGTGATAGCGCTTGGCCGGGGCCTCGAGGTCAAGCTTGCCCTCGGCG‑CTGGCAAGCGCGG‑CCAGGGTGGcc                                                                                                              >  1:369249/1‑150 (MQ=255)                                        aGTGGCAGGCAGTCGGCGCTATAGGCGCCCAGTTCCAGCACCGTCGCGTCGCCGATGGGGACGCCGGCCAATGCCGGGTGATAGCGCTTGGCCGGGGCCTCGAGGTCAAGCTTGCCCTCGGCG‑CTGGCAAGCGCGG‑CCAGGGTGGCCgtg                                                                                                           >  2:1394181/1‑150 (MQ=255)                                         gTGGCAGGCAGTCGGCGCTATAGGCGCCCAGTTCCAGCACCGTCGCGTCGCCGATGGGGACGCCGGCCAATGCCGGGTGATAGCGCTTGGCCGGGGCCTCGAGGTCAAGCTTGCCCTCGGCG‑CTGGCAAGCGCGG‑CCAGGGTGGCCgtgt                                                                                                          <  1:1394181/150‑1 (MQ=255)                                         gTGGCAGGCAGTCGCCGCTATAGGCGCCCAGTTCCAGCACCGCCGCGTGGGACTGGAGTTCAGCCGTGTGTTCCGGGTGATAGCGCTTGGCCGGGGCCTCGAGGTCAAGCTTGCCCTCGGCG‑CTGGCAAGCGCGG‑CCAGGGTGGCCgtgt                                                                                                          <  1:2027087/150‑1 (MQ=255)                                          tGGCAGGCAGTCGGCGCTATAGGCGCGCAGTTCCAGCACCGTCGCGTCGCCGATGGGGACGCCGGCCAATGCCGGGTGATAGCGCTTGGCCGGGGCCTCGAGGTCAAGCTTGCCCTCGGCG‑CTGGCAAGCGCGG‑CCAGGGTGGCCGTGTa                                                                                                         <  2:567397/150‑1 (MQ=255)                                           ggcaggcaGTCGGCGCTATAGGCGCCCAGTTCCAGCACCGTCGCGTCGCCGATGGGGACGCCGGCCAATGCCGGGTGATAGCGCTTGGCCGGGGCCTCGAGGTCAAGCTTGCCCTCGGCG‑CTGGCAAGCGCGG‑CCAGGGTGGCCGTGTAg                                                                                                        <  1:141381/150‑1 (MQ=255)                                            gcaggcagTCGGCGCTATAGGCGCCCAGTTCCAGCACCGTCGCGTCGCCGATGGGGACGCCGGCCAATGCCGGGTGATAGCGCTTGGCCGGGGCCTCGAGGTCAAGCTTGCCCTCGGCG‑CTGGCAAGCGCGG‑CCAGGGTGGCCGTGTAgg                                                                                                       <  2:2060560/150‑1 (MQ=255)                                            gcaggcagTCGGCGCTATAGGCGCCCAGTTCCAGCACCGTCGCGTCGCCGATGGGGACGCCGGCCAATGCCGGGTGATAGCGCTTGGCCGGGGCCTCGAGGTCAAGCTTGCCCTCGGCG‑CTGGCAAGCGCGG‑CCAGGGTGGCCGTGTAgg                                                                                                       >  1:273636/1‑150 (MQ=255)                                            gcaggcagTCGGCGCTATAGGCGCCCAGTTCCAGCACCGTCGCGTCGCCGATGGGGACGCCGGCCAATGCCGGGTGATAGCGCTTGGCCGGGGCCCCGAGGGCAAGCTTGACCTCga                                                                                                                                          <  1:2322767/117‑2 (MQ=255)                                            gcaggcagTCGGCGCTATAAGCGCCCAGTTCCAGCACCGTCGCGTCGCCGATGGGGACGCCGGCCAATGCCGGGTGATAGCGCTTGGCCGGGGCCCCGAGGGCAAGCTTGACCTCga                                                                                                                                          >  2:2322767/1‑116 (MQ=255)                                                 cagTCGGCGCTATAGGCGCCCAGTTCCAGCACCGTCGCGTCGCCGATGGGGACGCCGGCCAATGCCGGGTGATAGCGCTTGGCCGGGGCCTCGAGGTCAAGCTTGCCCTCGGCG‑CTGGCAAGCGCGG‑CCAGGGTGGCCGTGTGGGTCTTg                                                                                                  <  2:995060/150‑1 (MQ=255)                                                 cagTCGGCGCTATAGGCGCCCAGTTCCAGCACCGTCGCGTCGCCGATGGGGACGCCGGCCAATGCCGGGTGATAGCGCTTGGCCGGGGCCTCGAGGTCAAGCTTGCCCTCGGCG‑CTGGCAAGCGCGG‑CCAGGGTGGCCGTGTAGGTCTTg                                                                                                  <  1:468821/150‑1 (MQ=255)                                                 cagTCGGCGCTATAGGCGCCCAGTTCCAGCACCGTCGCGTCGCCGATGGGGACGCCGGCCAATGCCGGGTGATAGCGCTTGGCCGGGGCCTCGAGGTCAAGCTTGCCCTCGGCG‑CTGGCAAGCGCGG‑CCAGGGTGGCCGTGTAGGTCTTg                                                                                                  <  1:1108461/150‑1 (MQ=255)                                                   gTCGGCGCTATAGGCGCCCAGTTCCAGCACCGTCGCGTCGCCGATGGGGACGCCGGCCAATGCCGGGTGATAGCGCTTGGCCGGGGCCTCGAGGTCAAGCTTGACCTCGGCG‑CTGGCAAGCGCGG‑CCAGGGGGGGCGTGTAGGTCTTGCt                                                                                                >  1:2059113/1‑150 (MQ=255)                                                    tCGGCGCTATAGGCGCCCAGTTCCAGCACCGTCGCGTCGCCGATGGGGACGCCGGCCAATGCCGGGTGATAGCGCTTGGCCGGGGCCTCGAGGTCAAGCTTGCCCTCGGCG‑CTGGCAAGCGCGG‑CCAGGGTGGCCGTGTAGGTCTTGCTc                                                                                               <  2:748389/150‑1 (MQ=255)                                                     cGGCGCTATAGGCGCCCAGTTCCAGCACCGTCGCGTCGCCGATGGGGACGCCGGCCAATGCCGGGTGATAGCGCTTGGCCGGGGCCTCGAGGTCAAGCTTGCCCTCGGCG‑CCGGCAAGCGCGG‑CCAGGGTGGCCGTGTAGGTCTTGCTCa                                                                                              <  1:1942777/150‑1 (MQ=255)                                                      ggCGCTATAGGCGCCCAGTTCCAGCACCGTCGCGTCGCCGATGGGGACGCCGGCCAATGCCGGGTGATAGCGCTTGGCCGGGGCCTCGAGGTCAAGCTTGCCCTCGGCG‑CTGGCAAGCGCGG‑CCAGGGTGGCCGTGTAGGTCTTGCTCAg                                                                                             <  2:1494472/150‑1 (MQ=255)                                                      ggCGCTATAGGCGCCCAGTTCCAGCACCGTCGCGTCGCCGATGGGGACGCCGGCCAATGCCGGGTGATAGCGCTTGGCCGGGGCCTCGAGGTCAAGCTTGCCCTCGCCG‑CTGGCAAGCGCGG‑CCAGGGTGGCCGTGTAGGTCTTGCTCAg                                                                                             >  2:1429770/1‑150 (MQ=255)                                                         gcTATAGGCGCCCAGTTCCAGCACCGTCGCGTCGCCGATGGGGACGCCGGCCAATGCCGTGTGATAGCGCTTGGCCGGGGCCTCGAGGTCAAGCTTGCCCTCGGCG‑CTGGCAAGCGCGG‑CCAGGGTGGCCGTGTAGGTCt                                                                                                    <  1:404101/140‑1 (MQ=255)                                                         gcTATAGGCGCCCAGTTCCAGCACCGTCGCGTCGCCGATGGGGACGCCGGCCAATGCCGGGTGATAGCGCTTGGCCGGGGCCTCGAGGTCAAGCTTGCCCTCGGCG‑CTGGCAAGCGCGG‑CCAGGGTGGCCGTGTAGGTCt                                                                                                    >  2:404101/1‑140 (MQ=255)                                                         gcTATAGGCGCCCAGTTCCAGCACCGTCGCGTCGCCGATGGGGACGCAGGCCAATGGCCGGTGAGAGCGATTGCTCGGGGCCGCGAGGGCAAGCTTGCCCTCGGCG‑CTGGCAAGCGCGG‑CCAGGGTGGCCGTGTAGGTCTTGCTCAGCGa                                                                                          >  1:2037308/1‑150 (MQ=255)                                                           tataGGCGCCCAGTTCCAGCACCGTCGCGTCGCCGATGGGGACGCCGGCCAATGCCGGGTGATAGCGCTTGGCCGGGGCCTCGAGGTCAAGCTTGCCCTCGGCG‑CTGGCAAGCGCGG‑CCAGGGTGGCCGTGTAGGTCTTGCTCAGCGAAc                                                                                        <  2:2305971/150‑1 (MQ=255)                                                            ataGGCGCCCAGTTCCAGCACCGTCGCGTCGCCGATGGGGACGCCGGCCAATGCCGGGTGATAGCGCTTGGCCGGGGCCTCGAGGTCAAGCTTGCCCTCGGCG‑CTGGCAGGCCCGA‑CAAGGGGAGCCGTGTAGGTTTTGCTCGGCGAAcc                                                                                       >  2:1464644/1‑150 (MQ=255)                                                              aGGCGCCCAGTTCCAGCACCGTCGCGTCGCCGATGGGGACGCCGGCCAATGCCGGGTGATAGCGCTTGGCCGGGGCCTCGAGGTCAAGCTTGCCATCGGCG‑CTAGCAAGCGCGGTCCAAGG‑GGGCGAATAGGGCGTTCTCAGCGAACCGa                                                                                     >  2:1423700/1‑150 (MQ=255)                                                                  gcCCAGTTCCAGCACCGTCGCGTCGCCGATGGGGACGCCGGCCAATGCCGGGTGATAGCGCTTGGCCGGGGCCTCGAGGTCAAGCTTGCCCTCGGCG‑CTGGTAAGCGCGG‑CCAGGGTGGCCGTGTAGGTCTTGCTCAGCGAACCGAtctc                                                                                 <  1:775058/150‑1 (MQ=255)                                                                  gcCCAGTTCCAGCACCGTCGCGTCGCCGATGGGGACGCCGGCCAATGCCGGGTGATAGCGCTTGGCCGGGGCCTCGAGGTCAAGCTTGCCCTCGGCG‑CTGGCaa                                                                                                                                <  2:1799078/104‑1 (MQ=255)                                                                  gcCCAGTTCCAGCACCGTCGCGTCGCCGATGGGGACGCCGGCCAATGCCGGGTGATAGCGCTTGGCCGGGGCCTCGAGGTCAAGCTTGCCCTCGGCG‑CTGGCaa                                                                                                                                >  1:1799078/1‑104 (MQ=255)                                                                    ccAGTTCCAGCACCGTCGCGTCGCCGATGGGGACGCCGGCCAATGCCGGGTGATAGCGCTTGGCCGGGGCCTCGAGGTCAAGCTTGCCCTCGGCG‑CTGGCAAGCGCGG‑CCAGGGTGGCCGTGTAGGTCTTGCTCAGCGAACCGATCTCGa                                                                               >  1:1238469/1‑150 (MQ=255)                                                                      aGTTCCAGCACCGTCGCGTCGCCGATGGGGACGCCGGCCAATGCCGGGTGATAGCGCTTGGCCGGGGCCTCGAGGTCAAGCTTGCCCTCGGCG‑CTg                                                                                                                                    <  2:1973101/96‑1 (MQ=255)                                                                      aGTTCCAGCACCGTCGCGTCGCCGATGGGGACGCCGGCCAATGCCGGGTGATAGCGCTTGGCCGGGGCCTCGAGGTCAAGCTTGCCCTCGGCG‑CTg                                                                                                                                    >  1:1973101/1‑96 (MQ=255)                                                                       gTTCCAGCACCGTCGCGTCGCCGATGGGGACGCCGGCCAATGCCGGGTGATAGCGCTTGGCCGGGGCCTCGAGGTCAAGCTTGCCCTCGGCG‑CTCGCAAGCGCGG‑CCAGGGTGGCCGTGTAGGTCTTGCTCAGCGAACCGATCTCGAACa                                                                            <  2:1152304/150‑1 (MQ=255)                                                                           cAGCACCGTCGCGTCGCCGATGGGGACGCCGGCCAATGCCGGGTGATAGCGCTTGGCCGGGGCCTCGAGGTCAAGCTTGCCCTCGGCG‑CTGGCAAGCGCGG‑CCAGGGTGGCCGTGTAGGTCTTGCTCAGCGAACCGATCTCGAACAGCGt                                                                        >  1:1050504/1‑150 (MQ=255)                                                                             gCACCGTCGCGTCGCCGATGGGGACGCCGGCCAATGCCGGGTGATAGCGCTTGGCCGGGGCCTCGAGGTCAAGCTTGCCCTCGGCG‑CTGGCAAGCGCGG‑CCAGGGTGGCCGTGTAGGTCTTGCTCAGCGAACCGATCTCGAACAGCGTAt                                                                      <  1:2285948/150‑1 (MQ=255)                                                                               aCCGTCGCGTCGCCGATGGGGACGCCGGCCAATGCCGGGTGATAGCGCTTGGCCGGGGCCTCGAGGTCAAGCTTGCCCTCGGCG‑CTGGCAAGCGCGG‑CCAGGGTGGCCGTGTAGGTCTTGCTCAGCGAACCGATCTCGAACAGCGTATcg                                                                    >  2:1727947/1‑150 (MQ=255)                                                                                ccGTCGCGTCGCCGATGGGGACGCCGGCCAATGCCGGGTGATAGCGCTTGGCCGGGGCCTCGAGGTCAAGCTTGCCCTCGGCG‑CTGGCAAGCGCGG‑CCAGGGTGGCCGTGTAGGTCTTGCTCAGCGAACCGATCTCGAACAGCGTATcgc                                                                   <  1:1423700/150‑1 (MQ=255)                                                                                 cgtcgcgtcgCCTATGGGGACGCCGGCCAATGCCGGGTGATATCGCTTGGCCGGGGCCTCGAGGTCAAGCTTGCCCTCGGCG‑CTGGCAAGCGCGG‑CCAGGGTGGCCGTGTAGGTCTTGCTCAGCGAACCGATCTCGAACAGCGTATcgcg                                                                  <  1:68067/150‑1 (MQ=255)                                                                                  gtcgcgtcgcCGATGGGGACGCCGGCCAATGCCGGGTGATAGCGCTTGGCCGGGGCCTCGAGGTCAAGCTTGCCCTCGGCG‑CTGGCAAGCGCGG‑CCAGGGTGGCCGTGTAGGTCTTGCTCAGCGAACCGATCTCGAACAGCGTATCGCgg                                                                 <  2:369249/150‑1 (MQ=255)                                                                                  gtcgcgtcgcCGATGGGGACGCCGGCCAATGCCGGGTGATAGCGCTTGGCCGGGGCCTCGAGGTCAAGCTTGCCCTCGGCG‑CTGGCAAGCGCGG‑CCAGGGTGGCCGTGTAGGTCTTGCTCAGCGAACCGATCTCGAACAGCGTATCGCgg                                                                 <  1:1727947/150‑1 (MQ=255)                                                                                   tcgcgtcgcCGATGGGGACGCCGGCCAATGCCGGGTGATAGCGCTTGGCCGGGGCCTCGAGGTCAAGCTTGCCCTCGGCG‑CTGGCAAGCGCGG‑CCAGGGTGGCCGTGTAGGTCTTGCTCAGc                                                                                            <  1:40426/122‑1 (MQ=255)                                                                                   tcgcgtcgcCGATGGGGACGCCGGCCAATGCCGGGTGATAGCGCTTGGCCGGGGCCTCGAGGTCAAGCTTGCCCTCGGCG‑CTGGCAAGCGCGG‑CCAGGGTGGCCGTGTAGGTCTTGCTCAGc                                                                                            >  2:40426/1‑122 (MQ=255)                                                                                       gtcgcCGATGGGGACGCCGGCCAATGCCGGGTGATAGCGCTTGGCCGGGGCCTCGAGGTCAAGCTTGCCCTCGGCG‑CTGGCAAGCGCGG‑CCAGGGTGGCCGTGTAGGTCTTGCTCAGCGAACCGATCTCGAACAGCGTATCGCGGCTGAc                                                            >  1:467964/1‑150 (MQ=255)                                                                                             gATGGGGACGCCGGCCAATGCCGGGTGATAGCGCTTGGCCGGGGCCTCGAGGTCAAGCTTGCCCTCGGCG‑CTGGCAAGCGCGG‑CCAGGGTGGCCGTGTAGGTCTTGCTCAGCGAACCGATCTCGAACAGCGTATCGCGGCTGACCGGTAc                                                      <  2:427819/150‑1 (MQ=255)                                                                                                ggggATGTCGGCCAATGCAGGGTGATAGCGCTTGGCCGGGGCGTTGAGGTCAAGCTTGCCCTCGGCG‑CTGGCAAGCGCGG‑CCAGGGTGGCCGTGTAGGTTTTTCTCAGCCGACGGATCTCGAACTGCGTATGGCGGCGGACCGGTaaatt                                                   <  2:84596/150‑5 (MQ=255)                                                                                                  ggACGCCGGCCAATGCCGGGTGATAGCGCTTGGCCGGGGCCTCGAGGTCAAGCTTGCCCTCGGCG‑CTGGCAAGCGCGG‑CCAGGGTGGCCGTGTAGGTCTTGCTCAGCGAACCGATCTCGAACAGCGTATCGCGGCTGACCGGTACGTTg                                                  <  1:2131846/149‑1 (MQ=255)                                                                                                  ggACGCCGGCCAATGCCGGGTGATAGCGCTTGGCCGGGGCCTCGAGGTCAAGCTTGCCCTCGGCG‑CTGGCAAGCGCGG‑CCAGGGTGGCCGTGTAGGTCTTGCTCAGCGAACCGATCTCGAACAGCGTATCGCGGCTGACCGGTACGTTg                                                  >  2:2131846/1‑149 (MQ=255)                                                                                                   gACGCCGGCCAATGCCGGGTGATAGCGCTTGGCCGGGGCCTCGAGGTCAAGCTTGCCCTCGGCG‑CTGGCAAGCGCGG‑CCAGGGTGGCCGTGTAgg                                                                                                       >  2:1997579/1‑95 (MQ=255)                                                                                                   gACGCCGGCCAATGCCGGGTGATAGCGCTTGGCCGGGGCCTCGAGGTCAAGCTTGCCCTCGGCG‑CTGGCAAGCGCGG‑CCAGGGTGGCCGTGTAgg                                                                                                       <  1:1997579/95‑1 (MQ=255)                                                                                                        cGGCCAATGCCGGGTGATAGCGCTTGGCCGGGGCCTCGAGGTCAAGCTTGCCCTCGGCG‑CTGGCAAGCGCGG‑CCAGGGTGGCCGTGTAGGTCTTGCTCAGCGAACCGATCTCGaa                                                                              <  1:1112085/115‑1 (MQ=255)                                                                                                        cGGCCAATGCCGGGTGATAGCGCTTGGCCGGGGCCTCGAGGTCAAGCTTGCCCTCGGCG‑CTGGCAAGCGCGG‑CCAGGGTGGCCGTGTAGGTCTTGCTCAGCGAACCGATCTCGaa                                                                              >  2:1112085/1‑115 (MQ=255)                                                                                                         ggCCAATGCCGGGTGATAGCGCTTGGCCGGGGCCTCGAGGTCAAGCTTGCCCTCGGCG‑CTGGCAAGCGCGG‑CCAGGGTGGCCGTGTAGGTCTTGCTCAGCGAACCGAt                                                                                    <  2:187764/108‑1 (MQ=255)                                                                                                         ggCCAATGCCGGGTGATAGCGCTTGGACGGGGCCTCGAGGTCAAGCTTGCCCTCGGCG‑CTGGCAAGCGCGG‑CCAGGGTGGCCGTGTAGGTCTTGCTCAGCGAACCGAt                                                                                    >  1:187764/1‑108 (MQ=255)                                                                                                           ccAATGCCGGGTGATAGCGCTTGGCCGGGGCCTCGAGGTCAAGCTTGCCCTCGGCG‑CTGGCAAGCGCGG‑CCAGGGTGGCCGTGTAGGTCTTGCTCAGCGAACCGATCTCGAACAGCGTATCGCGGCTGACCGGTACGTTGTCGGCCTTGc                                        >  1:1191016/1‑150 (MQ=255)                                                                                                                gCCGGGTGATAGCGCTTGGCCGGGGCCTCGAGGTCAAGCTTGCCCTCGGCG‑CTGGCAAGCGCGG‑CCAGGGTGGCCGTGTAGGTCTTGCTCAGCGAACCGATCTCGAACAGCGTATCGCGGCTGACCGGTACGTTGTCGGCCTTGCTGGcc                                   >  1:1554229/1‑150 (MQ=255)                                                                                                                  cGGGTGATAGCGCTTGGCCGGGGCCTCGAGGTCAAGCTTGCCCTCGTCG‑CTGGCAAGCGCGG‑CCAGGGTGGCCGTGTAGGTCTTGCTCAGCGAACCGAGCTCGAACAGCGTATCGCGGCTGACCGGTACGTTGTCGGCCTTGCTGGCcac                                 <  1:841834/150‑1 (MQ=255)                                                                                                                  cGGGTGATAGCGCTTGGCCGGGGCCTCGAGGTCAAGCTTGCCCTCGGCG‑CTGGCAAGCGCGG‑CCAGGGTGGCCGTGTAGGTCTTGCTCAGCGAACCGATCTCGAACAGCGTATCGCGGCTGACCGGTACGTTGTCGGCCTTGCTGGCcac                                 >  1:1968464/1‑150 (MQ=255)                                                                                                                        aTAGCGCTTGGTCGGGGCCTAGAGGTCACGCTTGCCCTCGGCG‑CTGGCAAGCGCGG‑CCAGGGTGGCCGTGTAGGTCTTGCTCAGCGAACCGATCTCGAACAGCGTATCGCGGCTGACCGGTACGTTGTCGGCCTTGCTGGCCACACCGTa                           <  1:1100156/150‑1 (MQ=255)                                                                                                                                        gcgccGAGGGCAAGCTTGACCTCGGCG‑CTGGCAAGCGCGG‑CCAGGGTGGCCGTGTAGGTCTTGCTCAGCGAACCGATCTCGAACAGCGTATCGCGGCTGACCGGTACGTTGTCGGCCTTGCTGGCCACACCGTAGCTGAAGTAATGcgcg           >  2:2212391/5‑150 (MQ=255)                                                                                                                                           tCGAGGTCAAGCTTGCCCTCGGCG‑CTGGCAAGCGCGG‑CCAGGGTGGCCGTGTAGGTCTTGCTCAGCGAACCGATCTCGAACAGCGTATCGCGGc                                                                >  2:297125/1‑94 (MQ=255)                                                                                                                                           tCGAGGTCAAGCTTGCCCTCGGCG‑CTGGCAAGCGCGG‑CCAGGGTGGCCGTGTAGGTCTTGCTCAGCGAACCGATCTCGAACAGCGTATCGCGGc                                                                <  1:297125/94‑1 (MQ=255)                                                                                                                                           tCGAGGTCAAGCTTGCCCTCGGCG‑CTGGCAAGCGCGG‑CCAGGGTGGCCGTGTAGGTCTTGCTCAGCGAACCGATCTCGAACAGCGTATCGCGGCTGACCGGTACGTTGTCGGc                                             <  1:1402288/113‑1 (MQ=255)                                                                                                                                           tCGAGGTCAAGCTTGCCCTCGGCG‑CTGGCAAGCGCGG‑CCAGGGTGGCCGTGTAGGTCTTGCTCAGCGAACCGATCTCGAACAGCGTATCGCGGCTGACCGGTACGTTGTCGGc                                             >  2:1402288/1‑113 (MQ=255)                                                                                                                                           tCGAGGTCAAGCTTGCCCTCGGCG‑CTGGCAAGCGCGG‑CCAGGGTGGCCGTGTAGGTCTTGCTCAGCGAACCGATCTCGAACAGCGTATCGCGGCTGACCGGTACGTTGTCGGCCTTGCTGGCCACACCGTAGCTGaa                     >  2:225175/1‑137 (MQ=255)                                                                                                                                           tCGAGGTCAAGCTTGCCCTCGGCG‑CTGGCAAGCGCGG‑CCAGGGTGGCCGTGTAGGTCTTGCTCAGCGAACCGATCTCGAACAGCGTATCGCGGCTGACCGGTACGTTGTCGGCCTTGCTGGCCACACCGTAGCTGaa                     <  1:225175/137‑1 (MQ=255)                                                                                                                                                gTCAAGCTTGCCCTCGGCG‑CTGGCAAGCGCGG‑CCAGGGTGGCCGTGTAGGTCTTGCTCAGCGAACCGATCTCGAACAGCGTATCGCGGCTGACCGGTACGTTGTCGGCCTTGCTGGCCACACCGTAGCTGAAGTAATGCGCGTGGTcgcg   >  2:57065/1‑150 (MQ=255)                                                                                                                                                gTCAAGCTTGCCAT‑GTCGACGGGCAAGCGCGG‑CCAGGGTTGCCAGGTAGGTCTTGCTCAGCGAACCGATCTCGAACAGCGTATCGCGGCTGACCCGTACGTTGTCGGCCTTGCTGGCCACACCGTAGCTGAAGTAATGCGCGTGGTcgcg   >  1:1248234/1‑150 (MQ=255)                                                                                                                                                 tCAAGCTTGCCCTCGGCG‑CTGGCAAGCGCGG‑CCAGGGTGGCCGTGTAGGTCTTGCTCAGCGAACCGATCTCGAACAGCGTATCGCGGCTGACCGGTACGTTGTCGGCCTTGCTGGCCACACCGTAGCTGAAGTAATGCGCGTGGTCGCGt  >  1:2004543/1‑150 (MQ=255)                                                                                                                                                 |                                                                                                                                                        ACTACCTGTTGCGGTGTTTGTACGGTGTCGGGGAACTGCAGTGGCAGGCAGTCGGCGCTATAGGCGCCCAGTTCCAGCACCGTCGCGTCGCCGATGGGGACGCCGGCCAATGCCGGGTGATAGCGCTTGGCCGGGGCCTCGAGGTCAAGCTTGCCCTCGGCG‑CTGGCAAGCGCGG‑CCAGGGTGGCCGTGTAGGTCTTGCTCAGCGAACCGATCTCGAACAGCGTATCGCGGCTGACCGGTACGTTGTCGGCCTTGCTGGCCACACCGTAGCTGAAGTAATGCGCGTGGTCGCGT  >  NC_002947/3277709‑3278002

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑