Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,772,157	A→C	5.0%	I269S (ATC→AGC)	betT‑I ←	choline transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,772,157	0	A	C	5.0%	253.6 / 4.2	80	I269S (ATC→AGC)	betT‑I	choline transporter
Reads supporting (aligned to +/- strand):  ref base A (41/35);  new base C (2/2);  total (43/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.01e-01

GCCGGCTTCGCCTTCATACACATAGAGGTCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGAACAGGTATTCCAGGCCGGACGATACTT  >  NC_002947/5772009‑5772304                                                                                                                                                     |                                                                                                                                                    gCCGGCTTCGCCTTCATACACATAGAGGTCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCatga                                                                                                                                                     >  8:63500/1‑149 (MQ=255)  cctgCTTCGCCTTCATACACATAGAGGTCGAAGGTTTTCAGCACGATGCCGTTCATGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCAATAAGATCAACAGGCCACTGAACAGCatgat                                                                                                                                                    <  3:49690/146‑1 (MQ=255)     gCTTCGCCTTCATACACATAGAGGTCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGtt                                                                                                                                                 >  1:335086/1‑149 (MQ=255)       ttCGCCTTCATACACATAGAGGTCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGtt                                                                                                                                                 >  6:378146/1‑147 (MQ=255)       ttCGCCTTCATACACATAGAGGTCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGtt                                                                                                                                                 >  6:76355/1‑147 (MQ=255)         cGCCTTCATACACATAGAGGTCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAg                                                                                                                                              <  1:314683/148‑1 (MQ=255)         cGCCTTCATACACATAGAGGTCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAg                                                                                                                                              <  5:378146/148‑1 (MQ=255)          gCCATCATACACATAGAGGTCGAAGGTTTTCAGCACGATGCCGTTCAGCTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTATCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCAGGATGTTCAgga                                                                                                                                            >  4:68077/1‑148 (MQ=255)           ccTTCATACACATAGAGGTCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGtt                                                                                                                                           <  2:116017/149‑1 (MQ=255)            cTTCATACACATAGAGGTCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTg                                                                                                                                          >  4:2854/1‑149 (MQ=255)               cATACACATAGAGGTCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTgg                                                                                                                                         <  4:41578/147‑1 (MQ=255)               cATACACATAGAGGGCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACATCATGCTGTTCAgtggcc                                                                                                                                         <  2:205523/147‑6 (MQ=255)                    acaTAGAGGTCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACATCATGCTGTTCAgt                                                                                                                                             <  4:67328/138‑2 (MQ=255)                    acaTAGAGGTCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACATCATGCTGTTCAgt                                                                                                                                             >  3:67328/1‑137 (MQ=255)                    acaTAGAGGTCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAggcg                                                                                                                                  >  5:166271/1‑149 (MQ=255)                      agagagGTCGAAGGTTTTCAGCACGATGCCGTTCCGGTAGTCGCCAACGTTCTGCACGAAGCGGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGGACGAAGATCAACAGGCCACTGAACAGCATGATGGTCAGGTTGGACAggcggc                                                                                                                                >  6:181799/3‑149 (MQ=255)                      aTAGAGGTCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAggcggc                                                                                                                                >  3:264612/1‑149 (MQ=255)                      aTAGAGGTCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAggcggc                                                                                                                                >  1:324585/1‑149 (MQ=255)                        agagGTCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGAAGCACGAAGATAATCAGGCCACTGAACAGCATGATGTTCAGGTTGGAAAGgaggcg                                                                                                                               >  4:313489/1‑148 (MQ=255)                               gAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGt                                                                                                                       >  8:165035/1‑149 (MQ=255)                                aaGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGtt                                                                                                                      <  2:324585/149‑1 (MQ=255)                                  ggTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGATGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTtct                                                                                                                    <  3:313489/149‑1 (MQ=255)                                   gTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAggcggc                                                                                                                                >  2:333637/1‑136 (MQ=255)                                   gTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAggcggc                                                                                                                                <  1:333637/136‑1 (MQ=255)                                    ttttCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTtctc                                                                                                                   >  7:168818/1‑148 (MQ=255)                                         aGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGc                                                                                                              <  3:2854/148‑1 (MQ=255)                                                    gTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCgg                                                                                                                       >  8:157423/1‑127 (MQ=255)                                                    gTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGt                                                                                                                       <  7:157423/128‑1 (MQ=255)                                                    gTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTg                                                                                                  >  7:317709/1‑149 (MQ=255)                                                          gTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCTGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCa                                                                                                                                               <  5:350842/98‑1 (MQ=255)                                                          gTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCAAGCAGCACGAAGATCAACAGGCCACTGAACAGCATTATGTTCa                                                                                                                                               >  6:350842/1‑98 (MQ=255)                                                           tAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACg                                                                                             >  3:324499/1‑147 (MQ=255)                                                           tAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACg                                                                                             <  4:324499/147‑1 (MQ=255)                                                             gTCGCCAACGTTCTTCACGAAGCCGTTGAGCAGGTGCAGTGTTTCGCCGCCCAGCAGCACGAAGATCAACATGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTa                                                                                          <  3:68077/148‑1 (MQ=255)                                                              tCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGc                                                                                        >  7:195354/1‑149 (MQ=255)                                                                gCCAACGTTCTGCACGATGCCGTTGAGCAGGTGCAGGGTTTCGGCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGc                                                                                                 <  5:230584/138‑1 (MQ=255)                                                                gCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGcc                                                                                       <  2:335086/148‑1 (MQ=255)                                                                gCCAACGTTCTGCACGAAGCCGATGAGCACGTGCACGGTTTCGCCGCCCAGAAGCACGAACAACAACAGTCCACTGAACACCATGAAGTTCAGGTTGGACAGGCGTCGCATGACGTAATCAACGCCCGATACCGCTGc                                                                                                 >  6:230584/1‑138 (MQ=255)                                                                   aaCGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGc                                                                                   >  6:283156/1‑149 (MQ=255)                                                                   aaCGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGc                                                                                   >  7:123528/1‑149 (MQ=255)                                                                      gTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCgg                                                                                <  4:264612/149‑1 (MQ=255)                                                                        tCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGc                                                                               >  8:339222/1‑148 (MQ=255)                                                                            cACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCAt                                                                            <  8:317709/147‑1 (MQ=255)                                                                            cACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGa                                                                          >  8:77823/1‑149 (MQ=255)                                                                             aCGAAGCCGTTGAGCAGGTGCTGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGAc                                                                         >  1:116431/1‑149 (MQ=255)                                                                             aCGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGAc                                                                         <  7:63500/149‑1 (MQ=255)                                                                               gAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGAc                                                                         >  5:141771/1‑147 (MQ=255)                                                                                aaGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAg                                                                      <  5:76355/149‑1 (MQ=255)                                                                                 aGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAgg                                                                     <  8:123528/149‑1 (MQ=255)                                                                                    cGTTGAGCATGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATc                                                                  <  6:203849/149‑1 (MQ=255)                                                                                          gCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACTGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCaccacc                                                            <  5:319146/149‑1 (MQ=255)                                                                                          gCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGcc                                                                              <  8:254973/131‑1 (MQ=255)                                                                                           cAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCAccacca                                                           >  3:84023/1‑149 (MQ=255)                                                                                           cAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCAGCCATGACCAGGATCAccacca                                                           >  4:131401/1‑149 (MQ=255)                                                                                             ggTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCAGCCATGACCAGGATCACCACcagc                                                         <  3:131401/149‑1 (MQ=255)                                                                                               tGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACcagcag                                                       <  8:168818/149‑1 (MQ=255)                                                                                                  aGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCCCGCCCGATACCGCAGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCCGCAGGtt                                                    >  3:85474/1‑149 (MQ=255)                                                                                                  aGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGCTTGGACAGGCGGCGGATGCCGATCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCCCCACCAGCAGGtt                                                    >  2:192423/1‑149 (MQ=255)                                                                                                      tttCACCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACAAGGATCACCACCAGCAGGTTGGTc                                                >  7:166672/1‑149 (MQ=255)                                                                                                           ccgccCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCt                                           >  7:113017/1‑149 (MQ=255)                                                                                                                cagcagCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGt                                      <  4:84023/149‑1 (MQ=255)                                                                                                                   cagcaCGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCa                                   <  7:165035/149‑1 (MQ=255)                                                                                                                       aagatgATCAACATGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGcc                               <  2:89629/144‑1 (MQ=255)                                                                                                                       aCGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGcc                               >  5:261435/1‑149 (MQ=255)                                                                                                                        cGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCAGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGcc                               >  1:200751/1‑148 (MQ=255)                                                                                                                         gAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCg                              >  1:134162/1‑148 (MQ=255)                                                                                                                         gAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCg                              >  3:355333/1‑148 (MQ=255)                                                                                                                          aaGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGaa                            <  5:283156/149‑1 (MQ=255)                                                                                                                           aGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGAAc                           <  7:296022/149‑1 (MQ=255)                                                                                                                            gATCAACAGGCCCCTGAACAGCATGATGAGCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGAACa                          <  6:166271/149‑1 (MQ=255)                                                                                                                             aTCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGAACAg                         <  4:355333/149‑1 (MQ=255)                                                                                                                              tcatcaGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCGACGGTAGCCACGCCGGCCATGACCAGGATCACCAGCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGAACAgg                        <  1:192423/149‑1 (MQ=255)                                                                                                                                aaCAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTg                                                   <  4:321966/120‑1 (MQ=255)                                                                                                                                aaCAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTg                                                   >  3:321966/1‑120 (MQ=255)                                                                                                                                   aGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGTTTTGTCTTGCTGTGGTCCATGCCGAACAGGTATTc                   <  2:116431/149‑1 (MQ=255)                                                                                                                                   aGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGAACAGGTATTc                   >  3:89696/1‑149 (MQ=255)                                                                                                                                       cACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCGACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGAAGAGGTATTCCAgg               >  8:916/1‑149 (MQ=255)                                                                                                                                       cACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGAACAGGTATTCCAgg               >  7:126449/1‑149 (MQ=255)                                                                                                                                       accTGAACATCATGCTGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGAACAGGTATTCCAgg               >  1:95422/3‑149 (MQ=255)                                                                                                                                          tGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGAACAGGTATTCCAGGCCg            >  1:293608/1‑149 (MQ=255)                                                                                                                                             aCAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGAACAGGTATTCCAGGCCGGAc         <  8:195354/149‑1 (MQ=255)                                                                                                                                                 cATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGAACAGGTATTCCAGGCCGGACGAt      >  3:326735/1‑148 (MQ=255)                                                                                                                                                  atgatgTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGAACAGGTATTCCAGGCCGGACGATAc    >  4:190138/1‑149 (MQ=255)                                                                                                                                                    gatgTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGAACAGGTATTCCAGGCCGGACGATACtt  >  3:33507/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                    GCCGGCTTCGCCTTCATACACATAGAGGTCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGAACAGGTATTCCAGGCCGGACGATACTT  >  NC_002947/5772009‑5772304

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