Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,772,152	G→T	6.1%	L271M (CTG→ATG)	betT‑I ←	choline transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,772,152	0	G	T	6.1%	210.8 / 3.5	66	L271M (CTG→ATG)	betT‑I	choline transporter
Reads supporting (aligned to +/- strand):  ref base G (30/32);  new base T (2/2);  total (32/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 = 1.93e-01

ACTTGCCGGCTTCGCCTTCATACACATAGAGGTCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGAACAGGTATTCCAGGCCGGACG  >  NC_002947/5772005‑5772298                                                                                                                                                    |                                                                                                                                                   aCTTGCCGGCTTCGCCTTCATACACATAGAGGTCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAg                                                                                                                                                    >  4:23244/1‑148 (MQ=255)   ttGCCGGCTTCGCCTTCATACACATAGAGGTCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCat                                                                                                                                                 >  5:155386/1‑149 (MQ=255)     gCCGGCTTCGCCTTCATACACATATAGGTCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCatg                                                                                                                                                <  4:162233/148‑1 (MQ=255)     gCCGGCTTCGCCTTCATACACATAGAGGTCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCat                                                                                                                                                 <  3:235910/147‑1 (MQ=255)       cGGCTTCGCCTTCATACACATAGAGGTCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCatgatg                                                                                                                                             >  8:142450/1‑149 (MQ=255)         gCTTCGCCTTCATACACATAGAGGTCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGtt                                                                                                                                           >  2:165906/1‑149 (MQ=255)              gCCTTCATACACATAGAGGTCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAgg                                                                                                                                       <  2:26537/148‑1 (MQ=255)               ccTTCATACACATAGAGGTCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTc                                                                                                                                          >  3:226614/1‑144 (MQ=255)               ccTTCATACACATAGAGGTCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTc                                                                                                                                          <  4:226614/144‑1 (MQ=255)                   cATCCACATCGAGGTCGACGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCTTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGa                                                                                                                                  <  6:130130/148‑1 (MQ=255)                     tACACATAGAGGTCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGAGGTTCAGGTTGGACAg                                                                                                                               >  4:174887/1‑149 (MQ=255)                      acacATAGAGGTCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAgg                                                                                                                              >  8:11773/1‑149 (MQ=255)                                gTCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGc                                                                                                                    >  7:131240/1‑149 (MQ=255)                                   gAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGt                                                                                                                 >  2:72696/1‑149 (MQ=255)                                   gAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGt                                                                                                                 >  3:286701/1‑149 (MQ=255)                                   gAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGACGt                                                                                                                 >  6:67095/1‑149 (MQ=255)                                   gAAGGTTTTCAGCACGATGCCGTTCAAGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGt                                                                                                                 >  1:9676/1‑149 (MQ=255)                                       gTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGa                                                                                                                                  >  5:46170/1‑128 (MQ=255)                                       gTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGa                                                                                                                                  <  6:46170/128‑1 (MQ=255)                                          ttCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCAc                                                                                                          >  5:247633/1‑149 (MQ=255)                                          ttCAGCACGATGCCGTTCAGGTAGTCGCCAACCTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCa                                                                                                           <  1:122650/148‑1 (MQ=255)                                           tCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCAc                                                                                                          <  5:128494/148‑1 (MQ=255)                                             aGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGcc                                                                                                       >  3:161002/1‑149 (MQ=255)                                             aGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGcc                                                                                                       <  3:292795/149‑1 (MQ=255)                                              gCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGccc                                                                                                      <  7:11773/149‑1 (MQ=255)                                               cACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACATCATGCTGTTCAgtg                                                                                                                                      >  1:37214/1‑114 (MQ=255)                                               cACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACATCATGCTGTTCAgtg                                                                                                                                      >  1:41945/1‑114 (MQ=255)                                               cACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACATCATGCTGTTCAgtg                                                                                                                                      <  2:37214/116‑3 (MQ=255)                                               cACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACATCATGCTGTTCAgtg                                                                                                                                      <  2:41945/116‑3 (MQ=255)                                                  gATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATa                                                                                                  >  7:73840/1‑149 (MQ=255)                                                            aGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCa                                                                                         <  7:142450/148‑1 (MQ=255)                                                              gTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACgg                                                                                      <  4:286701/149‑1 (MQ=255)                                                              gTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACgg                                                                                      >  5:144041/1‑149 (MQ=255)                                                               tAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACgg                                                                                      >  7:187408/1‑148 (MQ=255)                                                               tAGTCGCCAAAGTTCTGCACGATGCCGTTGAGCAGGTGCAGGGTTTCGACGCCCTGCAGCACGTAGATAAACAGGCCACTGATCAGCATGATGTTCAGGTTGGACAGGCGTCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGt                                                                                     <  2:9676/149‑1 (MQ=255)                                                                  tCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGc                                                                                  >  5:188245/1‑149 (MQ=255)                                                                    gCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCGGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCa                                                                                >  2:325776/1‑149 (MQ=255)                                                                    gCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCa                                                                                >  8:210642/1‑149 (MQ=255)                                                                          gTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTa                                                                                    <  2:24221/139‑1 (MQ=255)                                                                          gTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCgg                                                                          >  3:235447/1‑149 (MQ=255)                                                                          gTTCTGCACGAAGCCGTTGAGCAGGTGCACGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTa                                                                                    >  1:24221/1‑139 (MQ=255)                                                                                   gAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCa                                                                 <  3:174887/149‑1 (MQ=255)                                                                                          ttGAGCAGGTGCAGGGTTTCGCCGCCCAGCATCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATcac                                                          <  7:43227/149‑1 (MQ=255)                                                                                          ttGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATcac                                                          <  4:235447/149‑1 (MQ=255)                                                                                          ttGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCAACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATcac                                                          <  6:155386/149‑1 (MQ=255)                                                                                              gCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATcaccac                                                       >  4:58714/1‑148 (MQ=255)                                                                                                aGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACcag                                                    >  8:204682/1‑149 (MQ=255)                                                                                                   tGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGACATGACCAGGATCACCACcagcag                                                 >  7:292768/1‑149 (MQ=255)                                                                                                      aGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGtt                                              <  6:188245/149‑1 (MQ=255)                                                                                                        ggTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTgg                                            <  3:23244/149‑1 (MQ=255)                                                                                                          tttCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTc                                          >  2:168698/1‑149 (MQ=255)                                                                                                          tttCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTc                                          >  8:185466/1‑149 (MQ=255)                                                                                                              gcggccCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTg                                       <  5:67095/145‑1 (MQ=255)                                                                                                              gccgccCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTg                                       <  1:72696/148‑1 (MQ=255)                                                                                                              gcagccCAGCAGGAGGAAGATCAACAGGCCACTGAACAGGATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATTCCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTTGTCTTg                                       <  8:292768/145‑1 (MQ=255)                                                                                                               ccgccCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCt                                     <  6:247633/149‑1 (MQ=255)                                                                                                                     agcagcACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTc                               >  1:163057/1‑149 (MQ=255)                                                                                                                     agcagcACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTc                               <  3:58714/149‑1 (MQ=255)                                                                                                                           aCGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGcc                         <  8:73840/149‑1 (MQ=255)                                                                                                                              aaGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGc                          <  6:304888/145‑1 (MQ=255)                                                                                                                              aaGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGc                          >  5:304888/1‑145 (MQ=255)                                                                                                                              aaGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGaa                      <  4:230637/149‑1 (MQ=255)                                                                                                                              aaGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGaa                      <  1:168698/149‑1 (MQ=255)                                                                                                                                gATCAACAGGCCACTGAACAGCATTATGTTCAGGTTGGTCAGGCGTCGTATGCGGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGAACa                    <  2:163057/149‑1 (MQ=255)                                                                                                                                 aTCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGAACAg                   >  2:165608/1‑149 (MQ=255)                                                                                                                                    aaCAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGAACAGGTa                <  5:124772/149‑1 (MQ=255)                                                                                                                                                  cAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGAACAGGAATTCCAGGCCGGACg  <  8:131240/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                   ACTTGCCGGCTTCGCCTTCATACACATAGAGGTCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGAACAGGTATTCCAGGCCGGACG  >  NC_002947/5772005‑5772298

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