Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	284,588	T→A	13.2%	intergenic (+34/+29)	betT‑II → / ← atsK	choline/carnitine/betaine transporter/alpha‑ketoglutarate‑dependent sulfonate dioxygenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	284,588	0	T	A	13.2%	182.5 / 19.3	68	intergenic (+34/+29)	betT‑II/atsK	choline/carnitine/betaine transporter/alpha‑ketoglutarate‑dependent sulfonate dioxygenase
Reads supporting (aligned to +/- strand):  ref base T (30/29);  new base A (4/5);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.43e-01

CCTGAGCGAAGGCAGCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCA‑GACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCACGTCATTCTCCTG  >  NC_002947/284440‑284736                                                                                                                                                     |                                                                                                                                                      ccTGAGCAAAGGCAGCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTg                                                                                                                                                      >  1:1077764/1‑150 (MQ=255)          aGGCAGCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGT‑CTGGCCTCTTCGCGGGCAAGCCCGCGAaga                                                                                                                                            <  2:465495/150‑1 (MQ=255)          aGGCAGCCAGAACGATGCCCTGGTGGGATACAGCAAGGCGCAGATCATCACCGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGAAAACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTCCCCGCGAag                                                                                                                                             <  2:181442/150‑1 (MQ=255)           ggCAGCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAaga                                                                                                                                            <  1:532154/150‑1 (MQ=255)              aGCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGc                                                                                                                                         <  2:129721/150‑1 (MQ=255)                   aaaatATGCCCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCA‑GAc                                                                                                                                   <  2:447073/146‑1 (MQ=255)                    aaCTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCA‑GACa                                                                                                                                  <  1:5674/150‑1 (MQ=255)                           aCCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCA‑GACAGGCAAcc                                                                                                                           >  2:552634/1‑150 (MQ=255)                            ccTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCA‑GACAGGCAACCg                                                                                                                          <  1:552634/150‑1 (MQ=255)                             cTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCGCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCAAGCCCGCGAAGAGGCCAGGACAggagagct                                                                                                                          <  1:1162727/150‑7 (MQ=38)                                      tACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCA‑GACAGGCAACCGCTAGAAAGGt                                                                                                                >  1:338814/1‑150 (MQ=255)                                      tACAGCAAGGAGCAGATCATCAACGACAGCCTAGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGTTCAAGCTCTCGTGTCCTGGCCTCTTCGCGGGGTTGCCCGCGAAGAGGCCA‑GACAGGCAACCGCTAGAAAGGt                                                                                                                <  2:216070/150‑1 (MQ=255)                                          gCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCGCTTCGCGGGCTTGCCCGCGAAGAGGCCA‑GACAGGCAACCGCTAGAAAGGTGCAt                                                                                                            <  2:175465/150‑1 (MQ=255)                                           cAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCA‑GACAGGCAACCGCTAGAAAGGTGCATc                                                                                                           >  1:640448/1‑150 (MQ=255)                                           cAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCA‑GACAGGCAACCGCTAGAAAGGTGCATc                                                                                                           >  2:1252644/1‑150 (MQ=255)                                           cAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCA‑GACAGGCAACCGCTAGAAAGGTGCATc                                                                                                           <  1:1252644/150‑1 (MQ=255)                                                   aGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCA‑GACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGa                                                                                                   >  2:20279/1‑150 (MQ=255)                                                         tcaACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCATGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCA‑GACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCg                                                                                             >  1:896664/1‑150 (MQ=255)                                                                     tCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGaa                                                                                                                                              <  2:737940/90‑1 (MQ=255)                                                                     tCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGaa                                                                                                                                              >  1:737940/1‑90 (MQ=255)                                                                      cGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCA‑GACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCAc                                                                                <  2:704573/150‑1 (MQ=255)                                                                         ccAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCA‑GACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGt                                                                                       <  2:350840/140‑1 (MQ=255)                                                                         ccAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCA‑GACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGt                                                                                       >  1:350840/1‑140 (MQ=255)                                                                          cAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCA‑GACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACccgc                                                                            >  2:945157/1‑150 (MQ=255)                                                                            gTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCA‑GACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACccgccg                                                                          >  2:1003898/1‑150 (MQ=255)                                                                             tACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCA‑GACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCgg                                                                         <  1:226955/150‑1 (MQ=255)                                                                             tACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCA‑GACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCgg                                                                         >  1:1220996/1‑150 (MQ=255)                                                                                    ggCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCA‑GACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTtcggtc                                                                  <  1:58879/150‑1 (MQ=255)                                                                                       acaTGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCA‑GACAGGCAACCGCTAGAAAGGTGCATCCCCaa                                                                                                      >  2:151069/1‑111 (MQ=255)                                                                                       acaTGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCA‑GACAGGCAACCGCTAGAAAGGTGCATCCCCaa                                                                                                      <  1:151069/111‑1 (MQ=255)                                                                                           gCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCA‑GACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAatc                                                           <  1:20279/150‑1 (MQ=255)                                                                                           gCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCA‑GACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAatc                                                           >  1:876320/1‑150 (MQ=255)                                                                                              gTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCA‑GACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAatcatc                                                        <  1:829867/150‑1 (MQ=255)                                                                                              gTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCA‑GACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAatcatc                                                        <  2:640448/150‑1 (MQ=255)                                                                                               tACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCA‑GACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATcc                                                       >  1:914178/1‑150 (MQ=255)                                                                                                 ccTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCA‑GACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAatcatc                                                        >  2:491817/1‑147 (MQ=255)                                                                                                 ccTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCA‑GACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAatcatc                                                        <  1:491817/147‑1 (MQ=255)                                                                                                   tGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCA‑GACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAAGCATCCACCg                                                   >  1:953458/1‑150 (MQ=255)                                                                                                       ccTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCA‑GACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCaaa                                               >  2:567363/1‑150 (MQ=255)                                                                                                              cGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCA‑GACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTgg                                        >  1:1261275/1‑150 (MQ=255)                                                                                                              cGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCA‑GACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTgg                                        >  2:663373/1‑150 (MQ=255)                                                                                                                cTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCA‑GACAGGTAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGt                                      >  2:1222831/1‑150 (MQ=255)                                                                                                                 tGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCA‑GACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTc                                     >  1:408988/1‑150 (MQ=255)                                                                                                                      aaGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCA‑GACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGTTAATCATCCACCGCAAAATGCTGGGTCAcgcg                                <  1:621790/150‑1 (MQ=255)                                                                                                                       aGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCA‑GACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCgg                               <  2:338814/150‑1 (MQ=255)                                                                                                                         ctctCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCA‑GACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGtt                             <  2:876320/150‑1 (MQ=255)                                                                                                                         ctctCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCA‑GACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCCGTGCGTCTCCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGtt                             <  2:1077764/150‑1 (MQ=255)                                                                                                                         ctctCCTGTCCTGGCCTCTGCGCGGGCTTGCCCGCGAAGAGGCCA‑GACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCCTCCACCGCAAAATGCTGGGTCACGCGGtt                             <  2:115123/150‑1 (MQ=255)                                                                                                                             cctgtcctgGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCA‑GACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTcc                         <  2:953458/150‑1 (MQ=255)                                                                                                                              ctgtcctgGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCA‑GACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTccc                        >  1:1231299/1‑150 (MQ=255)                                                                                                                              ctgtcctgGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCA‑GACAGGCAACCGCTAGAAAGGTGAATCCCCAAGGCTGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCg                                                   >  1:777833/1‑123 (MQ=255)                                                                                                                              ctgtcctgGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCA‑GACAGGCAACCGCTAGAAAGGTGAATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCg                                                   <  2:777833/123‑1 (MQ=255)                                                                                                                               tgtcctgGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCA‑GACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGACTGTCGCCCGGTGCATCAACCGCCGGGTCGCCCGCTAAGCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCa                       >  1:287823/1‑150 (MQ=255)                                                                                                                               ctgtctgGCCTCTTCGCGGGCAAGCCCGCGAAGAGGCCA‑GACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCgg                               <  2:189650/138‑1 (MQ=255)                                                                                                                               ctgtctgGCCTCTTCGCGGGCAAGCCCGCGAAGAGGCCA‑GACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCgg                               >  1:189650/5‑142 (MQ=255)                                                                                                                                  tctgGCCTCTTCGCGGGCAAGCCCGCGAAGAGGCCA‑GGCAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGaa                    >  1:901762/2‑150 (MQ=255)                                                                                                                                   ctgGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCA‑GACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAAc                   <  1:663373/150‑1 (MQ=255)                                                                                                                                       ccTCTTCGCGGGCTTGCCCGCGAAGAGGCCA‑GACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCa               <  2:408988/150‑1 (MQ=255)                                                                                                                                       ccTCTTCGCGGGCTTGCCCGCGAAGAGGCCA‑GACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCa               >  1:378904/1‑150 (MQ=255)                                                                                                                                        ctctTCGCGGGCTTGCCCGCGAAGAGGCCA‑GACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCAc              >  2:253118/1‑150 (MQ=255)                                                                                                                                        ctctTCGCGGGCTTGCCCGCGAAGAGGCCA‑GACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCAc              >  2:1121855/1‑150 (MQ=255)                                                                                                                                         tctTCGCGGGCTTGCCCGCGAAGAGGCCA‑GACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCACg             >  1:645252/1‑150 (MQ=255)                                                                                                                                         tctTCGCGGGCAAGCCCGCGAAGAGGCCA‑GACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCAc                                                     >  2:494541/1‑110 (MQ=255)                                                                                                                                         tctTCGCGGGCAAGCCCGCGAAGAGGCCA‑GACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCAc                                                     <  1:494541/110‑1 (MQ=255)                                                                                                                                          ctTCGCGGGCTTGCCCGCGAAGAGGCCA‑GACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCACGt            <  2:896664/150‑1 (MQ=255)                                                                                                                                           ttCGCGGGCTTGCCCGCGAAGAGGCCA‑GACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCACGTc           >  1:1234080/1‑150 (MQ=255)                                                                                                                                           ttCGCGGGCAAGCCCGCAAAGAGGCCAGGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTcggtcggt                                                               >  1:1172462/1‑99 (MQ=255)                                                                                                                                           ttCGCGGGCAAGCCCGCAAAGAGGCCAGGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTcggtcggt                                                               <  2:1172462/99‑1 (MQ=255)                                                                                                                                                    ttGCCCGCGAAGAGGCCA‑GACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCACGTCATTCTCCTg  <  2:1261275/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                      CCTGAGCGAAGGCAGCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCA‑GACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCACGTCATTCTCCTG  >  NC_002947/284440‑284736

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