Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	284,587	T→A	6.7%	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,587	0	T	A	6.7%	203.5 / 4.2	75	intergenic (+34/+29)	betT‑II/atsK	choline/carnitine/betaine transporter/alpha‑ketoglutarate‑dependent sulfonate dioxygenase
Reads supporting (aligned to +/- strand):  ref base T (32/38);  new base A (3/2);  total (35/40)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.88e-01

CCTGAGCGAAGGCAGCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAA‑GAGGCCAGACA‑GGCAACCGCTAGAAAGGTGCAT‑CCCCAAGGAT‑GGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCG‑‑‑C‑AAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCACGTCATTCTCCT  >  NC_002947/284439‑284734                                                                                                                                                     |                                                                                                                                                            ccTGAGCGAAGGCAGCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCtt                                                                                                                                                             >  2:242068/1‑149 (MQ=255)    gAGCGAAGGCAGCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGcc                                                                                                                                                          >  8:49568/1‑149 (MQ=255)     aGCGAAGGCAGCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGc                                                                                                                                                           <  6:329023/147‑1 (MQ=255)      gCGAAGGCAGCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGccc                                                                                                                                                         <  8:325877/148‑1 (MQ=255)        gAAGGCAGCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCcgc                                                                                                                                                       >  7:18783/1‑148 (MQ=255)          aGGCAGCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCTTCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGaa                                                                                                                                                    <  8:69793/149‑1 (MQ=255)              aGCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCACGGGCTTGCCCGCGAA‑GAgg                                                                                                                                               <  2:153619/149‑1 (MQ=255)               gCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAA‑GAGGc                                                                                                                                              >  5:154584/1‑149 (MQ=255)                ccAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGa                                                                                                                                                     >  1:172717/1‑142 (MQ=255)                ccAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGa                                                                                                                                                     <  2:172717/142‑1 (MQ=255)                 cAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAA‑GAGGcc                                                                                                                                             >  8:57494/1‑148 (MQ=255)                 cAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAA‑GAGGCCa                                                                                                                                            <  6:368093/149‑1 (MQ=255)                  aGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGaa                                                                                                                                                    >  4:176030/1‑141 (MQ=255)                  aGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGaa                                                                                                                                                    <  3:176030/141‑1 (MQ=255)                         tGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAA‑GAGGCCAGACA‑GGc                                                                                                                                    <  3:117721/148‑1 (MQ=255)                           aCCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAA‑GAGGCCAGACA‑GGCAAc                                                                                                                                 <  2:120274/149‑1 (MQ=255)                           aCCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGT‑CTGGCCTCTTCGCGGGCAAGCCCGCGaa                                                                                                                                                    >  7:88808/1‑131 (MQ=255)                           aCCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGT‑CTGGCCTCTTCGCGGGCAAGCCCGCGaa                                                                                                                                                    <  8:88808/131‑1 (MQ=255)                                gtggGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAA‑GAGGCCAGACA‑GGCAACCGCTa                                                                                                                            <  3:259304/149‑1 (MQ=255)                                    gCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGCCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAA‑GAGGCCAGACA‑GGCAACCGCTAGaaa                                                                                                                        <  7:85542/149‑1 (MQ=255)                                     cTACAGCAAGGAGCAGATCATCAACGCCCTCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAA‑GAGGCCAGACA‑GGCAACCGCTAGAAAg                                                                                                                       <  7:168980/149‑1 (MQ=255)                                      tACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGCA‑GAGGCCAGACA‑GGCAACCGCTAGAAAgg                                                                                                                      >  7:325880/1‑149 (MQ=255)                                       aCAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAA‑GAGGCCAGACA‑GGCAACCGCTAGAAAGGt                                                                                                                     >  7:172948/1‑149 (MQ=255)                                         aGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCAAGCCCGCGAA‑GAGGCCa                                                                                                                                            <  1:346507/125‑1 (MQ=255)                                         aGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCAAGCCCGCGAA‑GAGGCCa                                                                                                                                            >  2:346507/1‑125 (MQ=255)                                           cAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAA‑GAGGCCAGACA‑GGCAACCGCTAGAAAGGTGCAt                                                                                                                 >  2:285784/1‑149 (MQ=255)                                             aGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAa‑gag                                                                                                                                                >  3:123827/1‑117 (MQ=255)                                             aGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAa‑gag                                                                                                                                                <  4:123827/117‑1 (MQ=255)                                              ggAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAATGAGGCCAGACA‑GGCAACCGCTAGAAAGGTGCAT‑cc                                                                                                              <  1:242068/149‑1 (MQ=255)                                                 gCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAA‑GAGGCCAGACA‑GGCAACCGCTAGAAAGGTGCAT‑CCCCaa                                                                                                          <  8:172948/149‑1 (MQ=255)                                                 gCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAA‑GAGGCCAGACA‑GGCAACCGCTAGAAAGGTGCAT‑CCCCaa                                                                                                          <  5:349967/149‑1 (MQ=255)                                                    gATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCACATAGCCGGCGTTCCCGGCAA‑TGGGCCAGAAA‑CGCAACCGGGATCAAGGTGCCT‑ACCCAAgg                                                                                                        >  2:70182/1‑148 (MQ=255)                                                        atcaACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAA‑GAGGCCAGACA‑GGCAACCGCTAGAAAGGTGCAT‑CCCCAAGGAT‑GGt                                                                                                  >  7:296906/1‑149 (MQ=255)                                                                     tCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAA‑GAGGCCAGACA‑GGCAACCGCTAGAAAGGTGCAT‑CCCCAAGGAT‑GGTCGCCCGGTGCAt                                                                                      >  6:167989/1‑148 (MQ=255)                                                                     tCGACAAGTAAGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCAGTACTGGCCTCTTAGCGGGCTTGCCCGCGAA‑GAGGCCAGACA‑GGCACCCGCAAGAAAGGTGAAT‑CCCCATGGCT‑GGTCGCCCGGTTCAt                                                                                      >  6:375628/1‑148 (MQ=255)                                                                       gACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAA‑GAGGCCAGACA‑GGCAACCGCTAGAAAGGTGCAT‑CCCCAAGGAT‑GGTCGCCCGGTGCATCa                                                                                    <  1:70182/148‑1 (MQ=255)                                                                          cAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAA‑GAGGCCAGACA‑GGCAACCGCTAGAAAGGTGCAT‑CCCCAAGGAT‑GGTCGCCCGGTg                                                                                         >  6:358966/1‑140 (MQ=255)                                                                          cAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAA‑GAGGCCAGACA‑GGCAACCGCTAGAAAGGTGCAT‑CCCCAAGGAT‑GGTCGCCCGGTg                                                                                         <  5:358966/140‑1 (MQ=255)                                                                             tACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAA‑GAGGCCAGACA‑GGCAACCGCTAGAAAGGTGCAT‑CCCCAAGGAT‑GGTCGCCCGGTGCATCACccgccg                                                                             >  2:255881/1‑149 (MQ=255)                                                                                       acaTGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTGTTCGCGGGCTTGCCCGCGAT‑GAGGCCAGACA‑GGCAACCGCTAGAAAGGTGCAT‑CCCCAAGGAT‑GGTCGCCCGGTGCATCACCCGCCGGTtcggtcgg                                                                   <  3:270981/149‑1 (MQ=255)                                                                                         aTGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAA‑GAGGCCAGACA‑GGCAACCGCTAGAAAGGTGCAT‑CCCCAAGGAT‑GGTCGCCCGGTGCATCACCCGCCGGTTcggtcggt                                                                  >  4:221039/1‑148 (MQ=255)                                                                                         aTGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAA‑GAGGCCAGACA‑GGCAACCGCTAGAAAGGTGCAT‑CCCCAAGGAT‑GGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTa                                                                 <  8:317297/149‑1 (MQ=255)                                                                                            cAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAA‑GAGGCCAGACA‑GGCAACCGCTAGAAAGGTGCAT‑CCCCAAGGAT‑GGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAat                                                               >  8:333301/1‑148 (MQ=255)                                                                                                  cTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAA‑GAGGCCAGACA‑GGCAACCGCTAGAAAGGTGCAT‑CCCCAAGGAT‑GGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCAc                                                        <  7:333301/149‑1 (MQ=255)                                                                                                   tGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGAGAA‑GAGGCCAGACA‑GGCAACCGCTAGAAAGGTGCAT‑CCCCAAGGAT‑GGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCAcc                                                       >  6:367038/1‑149 (MQ=255)                                                                                                     cACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAA‑GAGGCCAGACA‑GGCAACCGCTAGAAAGGTGCATCCCCCAAGGAT‑GGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCg                                                      <  1:285784/149‑1 (MQ=255)                                                                                                     cACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAA‑GAGGCCAGACA‑GGCAACCGCTAGAAAGGTGCAT‑CCCCAAGGAT‑GGTCGCCCGGTGCATCACCCGCCGGTtcggt                                                                      <  4:21693/132‑1 (MQ=255)                                                                                                     cACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAA‑GAGGCCAGACA‑GGCAACCGCTAGAAAGGTGCAT‑CCCCAAGGAT‑GGTCGCCCGGTGCATCACCCGCCGGTtcggt                                                                      >  3:21693/1‑132 (MQ=255)                                                                                                        cTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAA‑GAGGCCAGACA‑GGCAACCGCTAGAAAGGTGCAT‑CCCCAAGGAT‑GGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCG‑‑‑C‑aa                                               >  3:128698/1‑148 (MQ=255)                                                                                                        cTGGGCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAA‑GAGGCCAGACA‑GGCAACCGCTAGAAAGGTGCAT‑CCCCAAGGAT‑GGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAatc                                                              <  7:140778/137‑1 (MQ=255)                                                                                                           gTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAA‑GAGGCCAGACA‑GGCAACCGCTAGAATGGTGCAT‑CCCCAAGGAT‑GGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCG‑‑‑C‑AAAATg                                           >  1:118289/1‑149 (MQ=255)                                                                                                              cGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAA‑GAGGCCAGACA‑GGCAACCGCTAGAAAGGTGCAT‑CCCCAAGGAT‑GGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCG‑‑‑C‑AAAATGCTg                                        <  6:154584/149‑1 (MQ=255)                                                                                                               gCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAA‑GAGGCCAGACA‑GGCAACCGCTAGAAAGGTGCAT‑CCCCAAGGAT‑GGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCG‑‑‑C‑aaa                                              <  6:305477/142‑1 (MQ=255)                                                                                                               gCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAA‑GAGGCCAGACA‑GGCAACCGCTAGAAAGGTGCAT‑CCCCAAGGAT‑GGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCG‑‑‑C‑aaa                                              >  5:305477/1‑142 (MQ=255)                                                                                                               gCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAA‑GAGGCCAGACA‑GGCAACCGCTAGAAAGGTGCAT‑CCCCAAGGAT‑GGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCG‑‑‑C‑AAAATGCTgg                                       <  7:57494/149‑1 (MQ=255)                                                                                                               gCTGACCAAGCTCTCCTGTCCGGGCCTGTTGGCGGGCTTGCCCTCGAA‑GAGGCCAGACA‑GGCAACCGCTAGAAAGGTGCAT‑CCCCAAGGAT‑GGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCG‑‑‑C‑AAAATGCTgg                                       <  8:18783/149‑1 (MQ=255)                                                                                                                     cAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAA‑GAGGCCAGACA‑GGCAACCGCTAGAAAGGTGCAT‑CCCCAAGGATGGGTCGCCCGGTGCATCACCCGCCGGTtcggtcgg                                                                   <  1:45460/120‑1 (MQ=255)                                                                                                                     cAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAA‑GAGGCCAGACA‑GGCAACCGCTAGAAAGGTGCAT‑CCCCAAGGAT‑GGTCGCCCGGTGCATCACCCGCCGGTtcggtcgg                                                                   >  2:45460/1‑119 (MQ=255)                                                                                                                     cAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAA‑GAGGCCAGACA‑GGCAACCGCTAGAAAGGTGCAT‑CCCCAAGGAT‑GGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCG‑‑‑C‑AAAATGCTGGGTCAcg                                 >  3:104066/1‑149 (MQ=255)                                                                                                                      aaGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAA‑GAGGCCAGACA‑GGCAACCGCTAGAAAGGTGCAT‑CCCCAAGGAT‑GGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCG‑‑‑C‑AAAATGCTGGGTCAcgc                                <  4:104066/149‑1 (MQ=255)                                                                                                                              ctttcctgGCCTCTGCCCGGCCTTGCCCGCGAA‑TAGGCCAGACA‑GGCAACCGCTAGAAAGGTGCAT‑CCCCAATGAT‑GGTCGCCCGGTGCATCACCCGCCGGTTTGGTCGTTAATCATCCACCG‑‑‑C‑AAAATGCTGTGTCACGCGGTTGTcc                        <  8:232646/146‑1 (MQ=255)                                                                                                                              ctgtcctgGCCTCTTCGCGGGCTTGCCCGCGAA‑GAGGCCAGACA‑GGCAACCGCTAGAAAGGTGCAT‑CCCCAAGGAT‑GGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCG‑‑‑C‑AAAATGCTg                                        >  2:132827/1‑133 (MQ=255)                                                                                                                              ctgtcctgGCCTCTTCGCGGGCTTGCCCGCGAA‑GAGGCCAGACA‑GGCAACCGCTAGAAAGGTGCAT‑CCCCAAGGAT‑GGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCG‑‑‑C‑AAAATGCTg                                        <  1:132827/133‑1 (MQ=255)                                                                                                                              ctgtcctgGCCTCTTCGCGGGCTTGCCCGCGAA‑GAGGCCAGACA‑GGCAACCGCTAGAAAGGTGCAT‑CCCCAAGGAT‑GGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCG‑‑‑C‑AAAATGCTGGGTCACGCGGTTGTcc                        <  1:255881/149‑1 (MQ=255)                                                                                                                                gtcctgGCCTCTTCGCGGGCTTGCCCGCGAA‑GAGGCCAGACAGGGCAACCGCTAGAAAGGTGCAT‑CCCCAAGGAT‑GGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCG‑‑‑C‑AAAATGCTGGGTCACGCGGTTGTccc                       >  1:256102/1‑149 (MQ=255)                                                                                                                                gtcctgGCCTCTTCGCGGGCTTGCCCGCGAA‑GAGGCCAGACA‑GGCAACCGCTAGAAAGGTGCAT‑CCCCAAGGAT‑GGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCG‑‑‑C‑AAAATGCTgg                                       <  7:347269/132‑1 (MQ=255)                                                                                                                                gtcctgGCCTCTTCGCGGGCTTGCCCGCGAA‑GAGGCCAGACA‑GGCAACCGCTAGAAAGGTGCAT‑CCCCAAGGAT‑GGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCG‑‑‑C‑AAAATGCTgg                                       >  8:347269/1‑132 (MQ=255)                                                                                                                                    tgGCCTCTTCGCGGGCTTGCCCGCGAA‑GAGGCCAGACA‑GGCAACCGCTAGAAAGGTGCAT‑CCCCAAGGAT‑GGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCG‑‑‑C‑AAAATGCTGGGTCACGCGGTTGTCCCAGAAc                  >  3:326498/1‑149 (MQ=255)                                                                                                                                      gCCTCTTCGCGGGCTTGCCCGCGAA‑GAGGCCAGACA‑GGCAACCGCTAGAAAGGTGCAT‑CCCCAAGGAT‑GGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCG‑‑‑C‑AAAATGCTGGGTCACGCGGTTGTCCCAGAACGc                >  2:321835/1‑149 (MQ=255)                                                                                                                                        ctctTCGCGGGCTTGCCCGCGAA‑GAGGCCAGACA‑GGCAACCGCTAGAAAGGTGCAT‑CCCCAAGGAT‑GGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCG‑‑‑C‑AAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCa              >  1:233100/1‑149 (MQ=255)                                                                                                                                        ctctTCGCGGGCAAGCCCGCGAA‑GAGGCCAGACA‑GGCAACCGCTAGAAAGGTGCAT‑CCCCAAGGAT‑GGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCG‑‑‑C‑AAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCa              >  2:239989/1‑149 (MQ=255)                                                                                                                                             cgcgGGCTTGCCCGCGAA‑GAGGCCAGACA‑GGCAACCGCTAGAAAGGTGCAT‑CCCCAAGGAT‑GGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCG‑‑‑C‑AAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCACGTc          <  3:221039/148‑1 (MQ=255)                                                                                                                                              gcgGGCTTGCCCGCGAA‑GAGGCCAGACA‑GGCAACCGCTAGAAAGGTGCAT‑CCCCAAGGAT‑GGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCG‑‑‑C‑AAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCACGTCa         <  7:49568/148‑1 (MQ=255)                                                                                                                                               cgGGCTTGCCCGCGAA‑GAGGCCAGACA‑GGCAACCGCTAGAAAGGTGCAT‑CCCCAAGGAT‑GGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCG‑‑‑C‑AAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCACGTCAtt       >  4:240889/1‑149 (MQ=255)                                                                                                                                               cgGGCGTGCCCGCGAA‑GAGGCCAGACA‑GGCAACCGCTAGAAAGGTGCAT‑CCCCAAGGAT‑GGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCG‑‑‑C‑AAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCACGTCAtt       <  8:296906/149‑1 (MQ=255)                                                                                                                                                gggCTTGCCCGGGAA‑GAGGCCAGACA‑GGCAACCGCTTGAAAGTTGCAT‑CCCCAAGCAT‑GGTCGCCCGGTGCATAACCCGCCGGTTCGGTCGGTAATCATCCACCG‑‑‑C‑AAAATGCTGGGTCACGCGTTTGTCCCAGAACGCCACGTCATtc      <  5:375628/149‑1 (MQ=255)                                                                                                                                                    ttGCCCGCGAA‑GAGGCCAGACA‑GGCAACCGCTAGAAAGGTGCAT‑CCCCAAGGAT‑GGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCATGGGGTCGAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCACGTCAtt       >  1:308151/1‑148 (MQ=255)                                                                                                                                                    ttGCCCGCGAA‑GAGGCCAGACA‑GGCAACCGCTAGAAAGGTGCAT‑CCCCAAGGAT‑GGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCG‑‑‑C‑AAAATGCTGGGTCACGCGGTTGTcc                        >  8:264447/1‑127 (MQ=255)                                                                                                                                                    ttGCCCGCGAA‑GAGGCCAGACA‑GGCAACCGCTAGAAAGGTGCAT‑CCCCAAGGAT‑GGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCG‑‑‑C‑AAAATGCTGGGTCACGCGGTTGTcc                        <  7:264447/127‑1 (MQ=255)                                                                                                                                                    ttGCCCGCGAA‑GAGGCCAGACA‑GGCAACCGCTAGAAAGGTGCAT‑CCCCAAGGAT‑GGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCG‑‑‑C‑AAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCACGTCATTCTCCt  >  7:322334/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                            CCTGAGCGAAGGCAGCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAA‑GAGGCCAGACA‑GGCAACCGCTAGAAAGGTGCAT‑CCCCAAGGAT‑GGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCG‑‑‑C‑AAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCACGTCATTCTCCT  >  NC_002947/284439‑284734

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