Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	284,586	T→A	17.9%	intergenic (+33/+30)	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,586	0	T	A	17.9%	117.9 / 17.2	61	intergenic (+33/+30)	betT‑II/atsK	choline/carnitine/betaine transporter/alpha‑ketoglutarate‑dependent sulfonate dioxygenase
Reads supporting (aligned to +/- strand):  ref base T (27/21);  new base A (5/6);  total (34/27)
Fisher's exact test for biased strand distribution p-value = 7.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

AGCGAAGGCAGCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAA‑GGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCACGTCATTCTCC  >  NC_002947/284443‑284733                                                                                                                                                |                                                                                                                                                     aGCGAAGGCAGCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGccc                                                                                                                                                 <  3:63555/149‑1 (MQ=255) aGCGAAGGCAGCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGccc                                                                                                                                                 <  3:147894/149‑1 (MQ=255)     aaGGCAGCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCAAGCCcgcg                                                                                                                                              >  8:116336/1‑148 (MQ=255)      aGGCAGCCAGAACTATGACCTGGTGGGCTACAGCAATGAGCAGATCATCAACGACATCCTCGACCAGTCCTAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCAAGCCCGCGaa                                                                                                                                            <  7:209449/149‑1 (MQ=255)      aGGCAGCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCAAGCCCGCGaa                                                                                                                                            <  8:234689/149‑1 (MQ=255)       ggCAGCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGT‑CTGGCCTCTTCGCGGGCAAGCCCGCGAaga                                                                                                                                          <  6:273706/149‑1 (MQ=255)         cAGCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCAAGCCCGCGAagag                                                                                                                                         <  7:116336/149‑1 (MQ=255)              aGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAag                                                                                                                                           <  1:19130/142‑1 (MQ=255)              aGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAag                                                                                                                                           >  2:19130/1‑142 (MQ=255)                   tATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAgg                                                                                                                               >  8:302950/1‑149 (MQ=255)                   tATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAgg                                                                                                                               <  6:57036/149‑1 (MQ=255)                            gtggGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCTGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACATGCAACCGc                                                                                                                        <  2:29416/147‑1 (MQ=255)                             tggGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGc                                                                                                                        >  2:263937/1‑146 (MQ=255)                                gCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCAAGCCCGCGAaga                                                                                                                                          <  1:198464/125‑1 (MQ=255)                                gCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCAGTCCTGGCCTCTTCGCGGGCAAGCCCGCGAaga                                                                                                                                          >  2:198464/1‑125 (MQ=255)                                      gCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAA‑GGTGCa                                                                                                           >  2:52247/1‑149 (MQ=255)                                       cAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAA‑GGTGCAt                                                                                                          >  5:179995/1‑149 (MQ=255)                                          ggAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCAGACAGGCAACCGCAAGAAA‑AGTGCATccc                                                                                                       >  2:317343/1‑149 (MQ=255)                                           gAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAA‑GGTGCATcccc                                                                                                      <  1:263937/149‑1 (MQ=255)                                           gAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTAGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAA‑GGTGCGTccc                                                                                                       >  2:137252/1‑148 (MQ=255)                                             gCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAA‑GGTGCATcccc                                                                                                      >  5:280117/1‑147 (MQ=255)                                                     tcaACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAA‑GGTGCATCCCCAAGGATGGTc                                                                                            >  2:42259/1‑149 (MQ=255)                                                            cATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAA‑GGTGCATCCCCAAGGATGGTCGCCCGGt                                                                                     >  6:102990/1‑149 (MQ=255)                                                                     ccAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAA‑GGTGCATCCCCAAGGATGGTCGCCCGGTGCATCAccc                                                                            >  3:262533/1‑149 (MQ=255)                                                                      cAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCGTGCCCGCGCAGAGGCCAGACAGGAAGCCGCTAGAAA‑GGTGAATCCCCAACGATGGTCGCCCGGTGCATCAAccg                                                                           >  8:125192/1‑149 (MQ=255)                                                                            gAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGCCCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAAAGGCCAGCCAGGGAACCGCTTGAAA‑GGTGGATCCCCAAGCATGGGCGGCCGGTGGCTCTCccgccgcgt                                                                     >  8:262796/1‑146 (MQ=37)                                                                                 gCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAA‑GGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTtcggt                                                                 <  3:130698/148‑1 (MQ=255)                                                                                 gCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAA‑GGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTtcggt                                                                 >  4:130698/1‑148 (MQ=255)                                                                                  cacaTGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAA‑GGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTtcggtcg                                                               <  4:262533/149‑1 (MQ=255)                                                                                   acaTGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAA‑GGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTtcggtcgg                                                              <  1:100096/149‑1 (MQ=255)                                                                                    caggCAGTCCCTGCCCCTGGCCCGCTGCCCATGAGCACGTGCCCTGCCCCCTGCGCGGGCTGGCCGGCTAAGAGGCCAGACAGGCAACCGCTAGAAA‑GGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTcggtcggt                                                             <  1:317343/146‑1 (MQ=255)                                                                                    caTGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAA‑GGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTcggtcggt                                                             <  5:144072/149‑1 (MQ=255)                                                                                    caTGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAA‑GGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTcggtcggt                                                             <  6:280117/149‑1 (MQ=255)                                                                                       gCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAA‑GGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTaa                                                           >  6:115769/1‑148 (MQ=255)                                                                                         aGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAA‑GGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAatc                                                         >  3:66562/1‑148 (MQ=255)                                                                                             ccTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGTGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAA‑GGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCa                                                    >  3:113289/1‑149 (MQ=255)                                                                                              cTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAA‑GGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCAc                                                   >  5:109073/1‑149 (MQ=255)                                                                                              cTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAA‑GGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCAc                                                   >  3:312101/1‑149 (MQ=255)                                                                                               tGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAA‑GGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCAcc                                                  >  5:272242/1‑149 (MQ=255)                                                                                                      ggTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAA‑GGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAAt                                           >  2:187244/1‑149 (MQ=255)                                                                                                          cGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAA‑GGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTg                                       <  7:302950/149‑1 (MQ=255)                                                                                                           gCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAA‑GGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTgg                                      >  5:84918/1‑149 (MQ=255)                                                                                                                ccAAGCTCTCCTGTCGTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAA‑GTTGCATCCCCAAGGATGGTCTCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCaaaa                                            <  5:267878/138‑1 (MQ=255)                                                                                                                ccAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGAAAGCGAAGAGGCCAGACAGGCAACCGCTAGAAA‑GGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACACGCCGGTTCGGTAGGTAATCATCCACCGCaaaa                                            >  6:267878/1‑138 (MQ=255)                                                                                                                 cAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAA‑GGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCAcg                                >  4:110152/1‑149 (MQ=255)                                                                                                                    gCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAA‑GGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCgg                             <  1:187244/149‑1 (MQ=255)                                                                                                                        tcctgtcctgGCCTCTTCGTGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAA‑GGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGtt                           <  4:113289/147‑1 (MQ=255)                                                                                                                         cctgtcctgGCCTCTTCGCGGGCTTGCCCCCGAAGAGGCCAGACAGGCAACCGCTAGAAA‑GGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCCGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCAcgccggagtc                        >  2:135795/1‑142 (MQ=255)                                                                                                                            gtcctgGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAA‑GGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTccc                      <  7:238232/148‑1 (MQ=255)                                                                                                                            gtcctgGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAA‑GGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTccc                      >  8:238232/1‑148 (MQ=255)                                                                                                                            gtcctgGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAA‑GGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTcc                       >  2:193427/1‑147 (MQ=255)                                                                                                                                     tctTCGCGGGCAAGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAA‑GGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCAc            >  7:89597/1‑149 (MQ=255)                                                                                                                                       ttCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAA‑GGTGCATCCCCAAGGATGGTCGCGCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGATCGCCAc            <  4:66562/147‑1 (MQ=255)                                                                                                                                       ttCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAA‑GGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCAc            <  5:115769/147‑1 (MQ=255)                                                                                                                                       ttCGCGGGCAAGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCACg           >  1:76506/1‑149 (MQ=255)                                                                                                                                       ttCGCGGGCAAGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAA‑GGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAAc                 >  6:8610/1‑142 (MQ=255)                                                                                                                                       ttCGCGGGCAAGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAA‑GGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAAc                 <  5:8610/142‑1 (MQ=255)                                                                                                                                         cgcgGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAA‑GGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCACGTc         <  6:109073/148‑1 (MQ=255)                                                                                                                                           ctggCTTGCCCGCGAAGAGGCCAGACAGGCAAACGCTAGAAA‑GGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCTGTGCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCACGTCAtt      <  7:125192/147‑1 (MQ=255)                                                                                                                                            gggCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAA‑GGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAGAATGCTGGGTCACGCGGTTGTCCCAGAACGCCACGTCATtc     >  5:276590/1‑149 (MQ=255)                                                                                                                                               cTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAA‑GGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCACGTCATTCTcc  >  8:55937/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                                     AGCGAAGGCAGCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAA‑GGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCACGTCATTCTCC  >  NC_002947/284443‑284733

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑