Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	626,559	T→G	7.5%	K213Q (AAG→CAG)	dedA ←	inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	626,559	0	T	G	7.5%	217.0 / 9.9	80	K213Q (AAG→CAG)	dedA	inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (35/39);  new base G (2/4);  total (37/43)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.44e-01

AGCGTCGCTGATGCCGTCCAGCATTGGCAGCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCAGGCCGCCGACC  >  NC_002947/626411‑626707                                                                                                                                                     |                                                                                                                                                     aGCGTCGCTGATGCCGTCCAGCATTGGCAGCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTg                                                                                                                                                     <  1:1505720/150‑1 (MQ=255)   cGTCGCTGATGCCGTCCAGCATTGGCAGCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGc                                                                                                                                                   <  1:344471/150‑1 (MQ=255)      cGCTGATGCCGTCCAGCATTGGCAGCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGt                                                                                                                                                >  2:102595/1‑150 (MQ=255)       gCTGATGCCGTCCAGCATTGGCAGCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTa                                                                                                                                               <  2:1549692/150‑1 (MQ=255)        cTGATGCCGTCCAGCATTGGCAGCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTAc                                                                                                                                              >  2:1678787/1‑150 (MQ=255)        cTGATGCCGTCCAGCATTGGCAGCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTAc                                                                                                                                              <  1:204720/150‑1 (MQ=255)        cTGATGCCGTCCAGCATTGGCAGCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCACCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTAc                                                                                                                                              >  2:1185384/1‑150 (MQ=255)              ccGTCCAGCATTGGCAGCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCa                                                                                                                                        <  2:1862240/150‑1 (MQ=255)               cGTCCAGCATTGGCAGCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAg                                                                                                                                       <  1:1678787/150‑1 (MQ=255)               cGTCCAGCATTGGCAGCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGGTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAg                                                                                                                                       <  1:1883347/150‑1 (MQ=255)                gTCCAGCATTGGCAGCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTGGGCCGTACGGc                                                                                                                                           >  1:743025/1‑145 (MQ=255)                gTCCAGCATTGGCAGCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTGGGCCGTACGGc                                                                                                                                           <  2:743025/145‑1 (MQ=255)                gTCCAGCATTGGCAGCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTGGGCCGTACGGCCAAgg                                                                                                                                      <  2:481886/150‑1 (MQ=255)                gTCCAGCATTGGCAGCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTGGGCCGTACGGCCAAgg                                                                                                                                      <  1:529710/150‑1 (MQ=255)                 tCCAGCATTGGCAGCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTCCCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGgc                                                                                                                                     <  2:859127/150‑1 (MQ=255)                   cAGCATTGGCAGCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGgcgc                                                                                                                                   >  1:131261/1‑150 (MQ=255)                   cAGCATTGGCAGCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGgcgc                                                                                                                                   >  2:547339/1‑150 (MQ=255)                       aTTGGCAGCCGCTAGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGTCGTACGGCCCAGGCGCCCAc                                                                                                                               <  2:368673/150‑1 (MQ=255)                          ggCAGCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCANgca                                                                                                                            <  2:1192536/150‑1 (MQ=255)                          ggCAGCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACgca                                                                                                                            <  1:952149/150‑1 (MQ=255)                                 gCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGcc                                                                                                                     <  1:925372/150‑1 (MQ=255)                                   tCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCga                                                                                                                   <  2:1868597/150‑1 (MQ=255)                                      cgTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGaa                                                                                                                <  2:393983/150‑1 (MQ=255)                                            gCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCg                                                                                                          <  2:1145232/150‑1 (MQ=255)                                               tGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCa                                                                                                       >  1:110678/1‑150 (MQ=255)                                               tGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGATTGCGTCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCa                                                                                                       <  1:408522/150‑1 (MQ=255)                                                gCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGCGCTTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCAc                                                                                                      <  1:1482696/150‑1 (MQ=255)                                                  cATTGCTGTGGCGGGATCGGGTAACGCGCCAGAGGACGATTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTAGCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGACGCCCACGCAGCAAGCCGAGAATCAACGGCACCa                                                                                                    >  2:1479518/1‑150 (MQ=255)                                                    ttGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAACAAGCTGTGAAGCATCGGCACCAGc                                                                                                  <  2:169098/150‑1 (MQ=255)                                                        tgtgGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACa                                                                                              <  2:773957/150‑1 (MQ=255)                                                         gtgGCGTGATCGGGTAACGCGCCAGAGTGCGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAg                                                                                             >  1:1793540/1‑150 (MQ=255)                                                             cGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGgatg                                                                                         >  2:1248905/1‑150 (MQ=255)                                                              gTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCCGCACCAGCGACAGgatga                                                                                        <  1:65442/150‑1 (MQ=255)                                                                 aTCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGgatgatga                                                                                     >  2:559848/1‑150 (MQ=255)                                                                  tCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGgatgatgat                                                                                    <  2:908211/150‑1 (MQ=255)                                                                   cGGGTAACGCGGCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATGATCGGCACCAGCGACAGGatgatgatg                                                                                   <  1:627882/150‑1 (MQ=255)                                                                    gggTAACGCTCCAGAGTACGCTTGCCCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGc                                                                                  <  1:1833618/150‑1 (MQ=255)                                                                    gggTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGc                                                                                  >  2:1427355/1‑150 (MQ=255)                                                                             gcCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCTCATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCccaccacca                                                                         >  1:1190792/1‑150 (MQ=255)                                                                             gcCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCccaccacca                                                                         >  1:1442745/1‑150 (MQ=255)                                                                                 gagTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAg                                                                     >  2:1837967/1‑150 (MQ=255)                                                                                 gagTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAg                                                                     <  2:1479851/150‑1 (MQ=255)                                                                                 gagTACGCTTGCACCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCCGCCAGCCGAGAATCATCGGCACCAGCTACAGGATGATGATGACCACCACCATCAg                                                                     <  1:1185384/150‑1 (MQ=255)                                                                                  agTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTAACTCTTCAGTGAGCCTTGGAAGGTGTTGGCAGTCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGc                                                                    >  1:1671454/1‑150 (MQ=255)                                                                                               cgccAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCCGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGGGAATGATGCCCACCACCATCAGCGACAGGTgctgct                                                       >  2:1095494/1‑150 (MQ=255)                                                                                                  cAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTtga                                                    <  1:1248789/150‑1 (MQ=255)                                                                                                     gcgcTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAACCGAGAATCATCGGCACCAGCGACAGGAAGATGATGCCCACCACCATCAGCGACAGGTGCTGCAtgatga                                                 >  1:466254/1‑150 (MQ=255)                                                                                                      cgcTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCTGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCCCCACCAGCGACAGGTGCTGCTTGATGaa                                                <  2:213777/150‑1 (MQ=255)                                                                                                       gcTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAAt                                               >  2:31682/1‑150 (MQ=255)                                                                                                       gcTGAATGACCACATGCCGTTCGTTCTTTACCTGTGCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAAt                                               <  1:26523/150‑1 (MQ=255)                                                                                                             tGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATTGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACg                                         >  1:1090245/1‑150 (MQ=255)                                                                                                             tGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACg                                         >  1:1775674/1‑150 (MQ=255)                                                                                                              gACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGt                                        >  2:413274/1‑150 (MQ=255)                                                                                                               aCCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGtt                                       >  1:981103/1‑150 (MQ=255)                                                                                                                 cacaTGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGc                                     >  2:621369/1‑150 (MQ=255)                                                                                                                   caTGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCg                                   <  2:981103/150‑1 (MQ=255)                                                                                                                    aTGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGACCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGAAGAATGGCACGTTGCCga                                  >  2:192851/1‑150 (MQ=255)                                                                                                                      gCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCgaag                                <  1:559848/150‑1 (MQ=255)                                                                                                                        cgttcgttCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCgaagaa                              <  1:621369/150‑1 (MQ=255)                                                                                                                          ttcgttcTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTtgat                                                   <  2:1622514/127‑1 (MQ=255)                                                                                                                          ttcgttcTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTtgat                                                   >  1:1622514/1‑127 (MQ=255)                                                                                                                             gttcTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGc                         <  1:1427355/150‑1 (MQ=255)                                                                                                                              ttcTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGcc                        >  1:1747441/1‑150 (MQ=255)                                                                                                                              ttcTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCAAAACCGACAGTATCTGCGTGAGGAATGGCACGTTGCCGAAGAAGTAGcc                        >  2:1754079/1‑150 (MQ=255)                                                                                                                              ttcTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGAAAGGATGATGATTCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGcc                        >  1:1771772/1‑150 (MQ=255)                                                                                                                               tcTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGccc                       >  2:1272077/1‑150 (MQ=255)                                                                                                                                    aCCTGTTCAGTGAGCCTTGGCCGTCAGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGt                  <  2:331395/150‑1 (MQ=255)                                                                                                                                    aCCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCATCGATAGGATGATGACGCCCACCCCCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGt                  >  1:1077217/1‑150 (MQ=255)                                                                                                                                      cTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGa                <  1:1015525/150‑1 (MQ=255)                                                                                                                                         ttCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATAAACGATGGCGACGGCGGGATGAAGGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCa             >  2:989884/1‑150 (MQ=255)                                                                                                                                           cAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCAgg           <  1:1837967/150‑1 (MQ=255)                                                                                                                                            aGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGGTGTGCGGGATGAGCGACAGTTGCTGCTTGATGGCTGGCCCCACGGTGACCAGGTATTCCAAGGTGACCAGgc          >  1:35420/1‑150 (MQ=255)                                                                                                                                             gTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGc                         >  2:209047/1‑134 (MQ=255)                                                                                                                                             gTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGc                         <  1:209047/134‑1 (MQ=255)                                                                                                                                              tGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCCGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGGGACCAGgccg        >  1:1271154/1‑150 (MQ=255)                                                                                                                                                 gCCTGGGCCGTACGGCCAAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGgat                                                                                          <  2:1498589/65‑1 (MQ=38)                                                                                                                                                 gCCTGGGCCGTACGGCCAAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGgat                                                                                          >  1:1498589/1‑65 (MQ=38)                                                                                                                                                  ccTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCAGGCCGCCGa    <  2:1873560/150‑1 (MQ=255)                                                                                                                                                   cTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCAGGCCGCCGAc   <  2:1442745/150‑1 (MQ=255)                                                                                                                                                   cTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCAGGCCGCCGAc   >  1:1766145/1‑150 (MQ=255)                                                                                                                                                    ttGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCAGGCCGCCGAcc  >  1:478020/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                                     AGCGTCGCTGATGCCGTCCAGCATTGGCAGCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCAGGCCGCCGACC  >  NC_002947/626411‑626707

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