Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,001,315	T→C	5.6%	Y329C (TAC→TGC)	amaA ←	L‑pipecolate oxidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,001,315	0	T	C	5.6%	209.1 / 2.3	71	Y329C (TAC→TGC)	amaA	L‑pipecolate oxidase
Reads supporting (aligned to +/- strand):  ref base T (32/35);  new base C (2/2);  total (34/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.45e-01

GTGCCGGTCGTGTTACCCCGCCCGTTGTACCCAGTCACCGCCACCAACCCAGGCGCCGGCTCGAACAGGCGCATCAGGTGGTCA‑GGGGTGAAGTCGATGCAGCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCGGGTAATAGTGGCTCTGGATGCGATCGGCCCAGCTGCGCACGAACCAGGCTGGCTTGTTGTCGACCCGGCCCAGGCTGCCGAGCAGCAGGCGGCCCTGATCGTCGCGGCGGATGCTGCTGAGCACGGTGCGGGTATCCCACGAACCTTGG  >  NC_002947/6001171‑6001461                                                                                                                                                  |                                                                                                                                                   gTGCCGGTCGTGTTACCCCGCCCGTTGTACCCAGTCACCGCCACCAACCCAGGCGCCGGCTCGAACAGGCGCATCAGGTGGTCA‑GGGGTGAAGTCGATGCAGCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCGGGTAATa                                                                                                                                                <  8:238088/149‑1 (MQ=255)     cgttCGTGTTACCCAGCCAGTTGTACCCAGTCACCGCCACCAACCCAGGCGCCGGCTCGAACAGGCGCATCAGGTTGTCA‑GGGGTGAAGTCGATGCAGCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCGGGTAATAGTgg                                                                                                                                            <  8:99003/146‑1 (MQ=255)             ttACCCCGCCCGTTGTACCCAGTCACCGCCACCAACCCAGGCGCCGGCTCGAACAGGCGCATCAGGTGGTCA‑GGGGTGAAGTCGATGCAGCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCGGGTAATAGTGGCTCTGGa                                                                                                                                     <  7:223856/148‑1 (MQ=255)             ttACCCCGCCCGTTGTACCCAGTCACCGCCACCAACCCAGGCGCCGGCTCGAACAGGCGCATCAGGTGGTCA‑GGGGTGAAGTCGATGCAGCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCGGGTAATAGTGGCTCTGGAt                                                                                                                                    <  2:187635/149‑1 (MQ=255)              tACCCCGCCCGTTGTACCCAGTCACCGCCACCAACCCAGGCGCCGGCTCGAACAGGCGCATCAGGTGGTCA‑GGGGTGAAGTCGATGCAGCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCGGGTAATAGTGGCTCTGGATg                                                                                                                                   >  3:220143/1‑149 (MQ=255)                 cccgcccgTTGTACCCAGTCACCGCCACCAACCCAGGCGCCGGCTCGAACAGGCGCATCAGGTGGTCA‑GGGGTGAAGTCGATGCAGCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCGGGTAATAGTGGCTCTGGATg                                                                                                                                   >  1:113157/1‑146 (MQ=255)                 cccgcccgTTGTACCCAGTCACCGCCACCAACCCAGGCGCCGGCTCGAACAGGCGCATCAGGTGGTCA‑GGGGTGAAGTCGATGCAGCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCGGGTAATAGTGGCTCTGGATg                                                                                                                                   <  2:113157/146‑1 (MQ=255)                      ccgTTGTACCCAGTCACCGCCACCAACCCAGGCGCCGGCTCGAACAGGCGCATCAGGTGGTCA‑GGGGTGAAGTCGATGCAGCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCGGGTAATAGTGGCTCTGGATGCGATCGGc                                                                                                                           <  2:187516/149‑1 (MQ=255)                         ttGTACCCAGTCACCGCCACCAACCCAGGCGCCGGCTCGAACAGGCGCATCAGGTGGTCA‑GGGGTGAAGTCGATGCAGCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCGGGTAATAGTGGCTCTGGATGCGATCGGCCCa                                                                                                                        <  4:220143/149‑1 (MQ=255)                              cccAGTCACCGCCACCAACCCAGGCGCCGGCTCGAACAGGCGCATCAGGTGGTCA‑GGGGTGAAGTCGATGCAGCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCGGGTAATAGTGGCTCTGGATGCGATCGGCCCAGCTgc                                                                                                                   <  7:33646/149‑1 (MQ=255)                               ccAGTCACCGCCACCAACCCAGGCGCCGGCTCGAACAGGCGCATCAGGTGGTCA‑GGGGTGAAGTCGATGCAGCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCGGGTAATAGTGGCTCTGGATGCGATCGGCCCAGCTgcg                                                                                                                  >  3:270760/1‑149 (MQ=255)                               ccAGTCACCGCCACCAACCCAGGCGCCGGCTCCAACAGGCGCATCAGGTGGTCA‑GGGGTGAAGTCGATGCAGCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCGGGTAATAGTGGCTCTGGATGCGATCGGCCCAGCTgcg                                                                                                                  >  1:35812/1‑149 (MQ=255)                                 aGTCACCGCCACCAACCCAGGCGCCGGCTCGAACAGGCGCATCAGGTGGTCA‑GGGGTGAAGTCGATGCAGCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCGGGTAATAGTGGCTCTGGATGCGATCGGCCCAGCTGCGCa                                                                                                                <  2:144835/149‑1 (MQ=255)                                      ccGCCACCAACCCAGGCGCCGGCTCGAACAGGCGCATCAGGTGGTCA‑GGGGTGAAGTCGATGCAGCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCGGGTAATAGTGGCTCTGGATGCGATCGGCCCAGCTGCGCACGAAc                                                                                                           >  2:181422/1‑149 (MQ=255)                                      ccGCCACCAACCCAGGCGCCGGCTCGAACAGGCGCATCAGGTGGTCA‑GGGGTGAAGTCGATGCAGCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCGGGTAATAGTGGCTCTGGATGCGATCGGCCCAGCTGCGCACGAAc                                                                                                           >  6:125030/1‑149 (MQ=255)                                       cGCCACCAACCCAGGCGCCGGCTCGAACAGGCGCATCAGGTGGTCA‑GGGGTGAAGTCGATGCAGCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCGGGTAATAGTGGCTCTGGATGCGATCGGCCCAGCTGCGCACGAAc                                                                                                           >  2:269434/1‑148 (MQ=255)                                       cGCCACCAACCCAGGCGCCGGCTCGAACAGGCGCATCAGGTGGTCA‑GGGGTGAAGTCGATGCAGCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCGGGTAATAGTGGCTCTGGATGCGATCGGCCCAGCTGCGCACGAAc                                                                                                           >  7:139318/1‑148 (MQ=255)                                           accaACCCAGGCGCCGGCTCGAACAGGCGCATCAGGTGGTCA‑GGGGTGAAGTCGATGCAGCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCGGGTAATAGTGGCTCTGGATGCGATCGGCCCAGCTGCGCACGAACCAggc                                                                                                      <  6:45955/149‑1 (MQ=255)                                            ccaACCCAGGCGCCGGCTCGAACAGGCGCATCAGGTGGTCA‑GGGGTGAAGTCGATGCAGCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCGGGTAATAGTGGCTCTGGATGCGATCGGCCCAGCTGCGCACGAACCAggct                                                                                                     <  5:255332/149‑1 (MQ=255)                                             caACCCAGGCGCCGGCTCGAACAGGCGCATCAGGTGGTCA‑GGGGTGAAGTCGATGCAGCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCGGGTAATAGTGGCTCTGGATGCGATCGGCCCAGCTGCGCACGAACCAggctg                                                                                                    >  3:205366/1‑149 (MQ=255)                                               cctgatgcgcCTGTTCGAGCCGGCGCATCAGGTGGTCA‑GGGGTGAAGTCGATGCAGCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCGGGTAATAGTGGCTCTGGATGCGATCGGCCCAGCTGCGCa                                                                                                                >  8:14220/7‑135 (MQ=255)                                               aCCCAGGCGCCGGCTCGAACAGGCGCATCAGGTGGTCA‑GGGGTGAAGTCGATGCAGCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCGGGTAATAGAGGCTCTGGATGCGATCGGCCCAGCTGAGCACGATCCAggctgga                                                                                                  >  7:125089/1‑148 (MQ=255)                                                    ggCGCCGGCTCGACCAGGCGCATCAGGTGGTCA‑GGGGTGAAGTCGATGCAGCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCGGGTAATAGTGGCTCTGGATGCGATCGGCCCAGCTGCGCACGAACCAGGCTGGCttgtt                                                                                             <  4:270760/149‑1 (MQ=255)                                                    ggCGCCGGCTCGAACAGGCGCATCAGGTGGTCA‑GGGGTGAAGTCGATGCAGCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCGGGTAATAGTGGCTCTGGATGCGATCGGCCCAGCTGCGCACGAACCAGGCTGGCttgtt                                                                                             <  7:58808/149‑1 (MQ=255)                                                      cgcCGGCTCGAACAGGCGCATCAGGTGGTCA‑GGGGTGAAGTCGATGCAGCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCGGGTAATAGTGGCTCTGGATGCGATCGGCCCAGCTGCGCACGAACCAGGCTGGCTtgttgt                                                                                           <  1:181422/149‑1 (MQ=255)                                                       gcCGGCTCGAACAGGCGCATCAGGTGGTCAGGGGGTGAAGTCGATGCAGCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCGGGTAATAGTGGCTCTGGATGCGATCGGCCCAGCTGCGCACGAACCAGGCTGGCtt                                                                                                >  1:40233/1‑144 (MQ=255)                                                       gcCGGCTCGAACAGGCGCATCAGGTGGTCA‑GGGGTGAAGTCGATGCAGCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCGGGTAATAGTGGCTCTGGATGCGATCGGCCCAGCTGCGCACGAACCAGGCTGGCtt                                                                                                <  2:40233/143‑1 (MQ=255)                                                         cGGCTCGAACAGGCGCATCAGGTGGTCA‑GGGGTGAAGTCGATGCAGCCGGTCCGGTGCATTTCCCACTCGACCTTGCCCAGCTCCGGGTAATAGTGGCTCTGGATGCGATCGGCCCAGCTGCGCACGAACCAGGCTGGCTTGTTGTCgc                                                                                        >  1:235160/1‑148 (MQ=255)                                                               gAACAGGCGCATCAGGTGGTCA‑GGGGTGAAGTCGATGCAGCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCGGGTAATAGTGGCTCTGGATGCGATCGGCCCAGCTGCGCACGAACCAGGCTGGCTTGTTGTCGACCCgg                                                                                   <  1:269434/148‑1 (MQ=255)                                                               gAACAGGCGCATCAGGTGGTCA‑GGGGTGAAGTCGATGCAGCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCGGGTAATAGTGGCTCTGGATGCGATCGGCCCAGCTGCGCACGAACCAGGCTGGCTTGTTGTCGACCCgg                                                                                   <  5:125030/148‑1 (MQ=255)                                                               gAACAGGCGCATCAGGTGGTCA‑GGGGTGAAGTCGATGCAGCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCGGGTAATAGTGGCTCTGGATGCGATCGGCCCAGCTGCGCACGAACCAGGCTGGCTTGTTGTCGACCCg                                                                                    <  6:116912/147‑1 (MQ=255)                                                                    ggCGCATCAGGTGGTCA‑GGGGTGAAGTCGATGCAGCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCGGGTAATAGTGGCTCTGGATGCGAGCGGCCCAGCTGCGCACGAACCAGGCTGGCTTGTTGTCGACCCGGCCCAgg                                                                             <  4:234315/149‑1 (MQ=255)                                                                      cgaaTCAGGTGGTCA‑GGCGTGAAGTCGATGCAGCCGGTCCTGTGCATTTCCCACTCGACCTTGCCCAGCTCCGGGTAATAGTGGCTCTGGATGCGATCGGCCCAGCTGCGCACGAACCAGGCTGGCTTGTTGTCGACCCGGCCCATGCt                                                                           <  4:207690/146‑1 (MQ=255)                                                                       gcATCAGGTGGTCA‑GGGGTGAAGTCGATGCAGCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCGGGTAATAGTGGCTCTGGATGCGATCGGCCCAGCTGCGCACGAACCAGGCTGGCTTGTTGTCGACCCGGCCCAGGCTg                                                                          <  4:205366/149‑1 (MQ=255)                                                                       gcATCAGGTGGTCA‑GGGGTGAAGTCGATGCAGCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCGGGTAATAGTGGCTCTGGATGCGATCGGCCCAGCTGCGCACGAACCAGGCTGGCTTGTTGTCGACCCGGCCCAGGCTg                                                                          >  7:177203/1‑149 (MQ=255)                                                                            aGGTGGTCA‑GGGGTTAAGTCGATGCAGCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCGGGTAATAGTGGCTCTGGATGCGATCGGCCCAGCTGCGCACGAACCAGGCTGGCTTGTTGTCGACCCGGCCCAGGCTGCCGa                                                                      <  2:162735/148‑1 (MQ=255)                                                                                  tCA‑GGGGTGAAGTTGATGCAGCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCGGGTAATAGTGGCTCTGGATGCGATCGGCCCAGCTGCGCACGAACCAggct                                                                                                     <  8:24810/111‑1 (MQ=255)                                                                                  tCA‑GGGGTGAAGTTGATGCAGCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCGGGTAATAGTGGCTCTGGATGCGATCGGCCCAGCTGCGCACGAACCAggct                                                                                                     >  7:24810/1‑111 (MQ=255)                                                                                       gggTGAAGTCGATGCAGCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCGGGTAATAGTGGCTCTGGATGCGATCGGCCCAGCTGCGCACGAACCAGGCTGGCTTGTTGTCGACCCGGCCCAGGCTGCCGAGCAGCAggcgg                                                           >  4:68555/1‑149 (MQ=255)                                                                                       gggTGAAGTCGATGCAGCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCGGGTAATAGTGGCTCTGGATGCGATCGGCCCAGCTGCGCACGAACCAGGCTGGCTTGTTGTCGACCCGGCCCAGGCTGCCGAGCAGCAggcgg                                                           <  8:177203/149‑1 (MQ=255)                                                                                         gTGAAGTCGATGCAGCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCGGGTAATAGTGGCTCTGGATGCGATCGGCCCAGCTGCGCACGAACCAGGCTGGCTTGTTGTCGACCCGGCCCAGGCTGCCGAGCAGCAGGCGGcc                                                         <  8:139318/149‑1 (MQ=255)                                                                                              gTCGATGCAGCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCGGGTAATAGTGGCTCTGGATGCGATCGGCCCAGCTGCGCACGAACCAGGCTGGCTTGTTGTCGACCCGGCCCAGGCTGCCGAGCAGCAGGCGGCCCTGAt                                                    >  7:246362/1‑149 (MQ=255)                                                                                                 gATGCAGCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCGGGTAATAGTGGCTCTGGATGCGATCGGCCCAGCTGCGCACGAACCAGGCTGGCTTGTTGTCGACCCGGCCCAGGCTGCCGAGCAGCAGGCGGCCCTGAtcgt                                                 >  6:9554/1‑149 (MQ=255)                                                                                                  aTGCAGCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCGGGTAATAGTGGCTCTGGATGCGATCGGCCCAGCTGCGCACGAACCAGGCTGGCTTGTTGTCGACCCGGCCCAGGCTGCCGAGCAGCAGGCGGCCCTGAtcgtc                                                <  5:109257/149‑1 (MQ=255)                                                                                                      aGCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCGGGTAATAGTGGCTCTGGATGCGATCGGCCCAGCTGCGCACGAACCAGGCTGGCTTGTTGTCGACCCGGCCCAGGCTGCCGAGCAGCAggcggc                                                          >  1:281883/1‑135 (MQ=255)                                                                                                      aGCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCGGGTAATAGTGGCTCTGGATGCGATCGGCCCAGCTGCGCACGAACCAGGCTGGCTTGTTGTCGACCCGGCCCAGGCTGCCGAGCAGCAggcggc                                                          <  2:281883/135‑1 (MQ=255)                                                                                                       gCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCGGGCAATAGTGGCTCTGGATGCGATCGGCCCAGCTGCGCACGAACCAGGCTGGCTTGTTGTCGACCTGGCCCAGGCTGCCGAGCAGCAGGCGGCCCTGATCGTCgcggc                                           >  3:281433/1‑149 (MQ=255)                                                                                                       gCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCGGGCAATAGTGGCTCTGGATGCGATCGGCCCAGCTGCGCACGAACCAGGCTGGCTTGTTGTCGACCTGGCCCAGGCTGCCGAGCAGCAGGCGGACCTGATCGTCgcggc                                           >  3:285964/1‑149 (MQ=255)                                                                                                         cGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCGGGTAATAGTGGCTCTGGATGCGATCGGCCCAGCTGCGCACGAACCAGGCTGGCTTGTTGTCGACCCGGCCCAGGCTGCCg                                                                       >  1:258354/1‑119 (MQ=255)                                                                                                         cGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCGGGTAATAGTGGCTCTGGATGCGATCGGCCCAGCTGCGCACGAACCAGGCTGGCTTGTTGTCGACCCGGCCCAGGCTGCCg                                                                       <  2:258354/119‑1 (MQ=255)                                                                                                            tCCATTGCATTTCCCACTCGACCTTGCCCATCTCCGGGTAATAGTGGCTCTGGATGCGATCGGCCCAGCTGCGCACGAACCAGGCTGGCTTGTTGTCGACCCGGCCCAGGCTGCCGAGCAGCAGGCGGCCCTGATCGTCGCGGCGGAtg                                      <  8:125089/149‑1 (MQ=255)                                                                                                                gTGCATTTCCCACTCGACCTTGCCCAGCTCCGGGTAATAGTGGCTCTGGATGCGATCGGCCCAGCTGCGCACGAACCAGGCTGGCTTGTTGTCGACCCGGCCCAGGCTGCCGAGCAGCAGGCGGCCCTGAt                                                    >  5:251517/1‑131 (MQ=255)                                                                                                                gTGCATTTCCCACTCGACCTTGCCCAGCTCCGGGTAATAGTGGCTCTGGATGCGATCGGCCCAGCTGCGCACGAACCAGGCTGGCTTGTTGTCGACCCGGCCCAGGCTGCCGAGCAGCAGGCGGCCCTGAt                                                    >  3:284151/1‑131 (MQ=255)                                                                                                                 tGCATTTCCCACTCGACCTTGCCCAGCTCCGGGTAATAGTGGCTCTGGATGCGATCGGCCCAGCTGCGCACGAACCAGGCTGGCTTGTTGTCGACCCGGCCCAGGCTGCCGAGCAGCAGGCGGCCCTGATCGTCGCGGCGGATgctgct                                 >  6:144065/1‑149 (MQ=255)                                                                                                                    aTTTCCCACTCGACCTTGCCCAGCTCCGGGTAATAGTGGCTCTGGATGCGATCGGCCCAGCTGCGCACGAACCAGGCTGGCTTGTTGTCGACCCGGCCCAGGCTGCCGAGCAGCAGGCGGCCCTGATCGTCGCGGCGGATGCTGCTGAg                              <  2:35812/149‑1 (MQ=255)                                                                                                                     tttCCCACTCGACCTTGCCCAGCTCCGGGTAATAGTGGCTCTGGATGCGATCGGCCCAGCTGCGCACGAACCAGGCTGGCTTGTTGTCGACCCGGCCCAGGCTGCCGAGCAGCAGGCGGCCCTGATCGTCGCGGCGGATGCTGCTGAg                              >  1:211973/1‑148 (MQ=255)                                                                                                                      ttCCCACTCGACCTTGCCCAGCTCCGGGTAATAGTGGCTCTGGATGCGATCGGCCCAGCTGCGCACGAACCAGGCTGGCTTGTTGTCGACCCGGCCCAGGCTGCCGAGCAGCAGGCGGCCCTGATCGTCGCGGCGGATGCTGCTGAGCa                            >  6:170064/1‑149 (MQ=255)                                                                                                                       tCCCACTCGACCTTGCCCAGCTCCGGGTAATAGTGGCTCTGGATGCGATCGGCCCAGCTGCGCACGAACCAGGCTGGCTTGTTGTCGACCCGGCCCAGGCTGCCGAGCAGCAGGCGGCCCTGATCGTCGCGGCGGATGCTGCTGAGCAc                           >  7:78186/1‑149 (MQ=255)                                                                                                                        cccACTCGACCTTGCCCAGCTCCGGGTAATAGTGGCTCTGGATGCGATCGGCCCAGCTGCGCACGAACCAGGCTGGCTTGTTGTCGACCCGGCCCAGGCTGCCGAGCAGCAGGCGGCCCTGATCGTCGCGGCGGATGCTGCTGAGCACg                          >  4:260174/1‑149 (MQ=255)                                                                                                                          cACTCGACCTTGCCCAGCTCCGGGTAATAGTGGCTCTGGATGCGATCGGCCCAGCTGCGCACGAACCAGGCTGGCTTGTTGTCGACCCGGCCCAGGCTGCCGAGCAGCAGGCGGCCCTGATCGTCGCGGCGGATGCTGCTGAGCACGGt                        >  1:58479/1‑149 (MQ=255)                                                                                                                               gACCTTGCCCAGCTCCGGGTAATAGTGGCTCTGGATGCGATCGGCCCAGCTGCGCACGAACCAGGCTGGCTTGTTGTCGACCCGGCCCAGGCTGCCGAGCAGCAGGCGGCCCTGATCGTCGCGGCGGATGCTGCTGAGCACGGTGCggg                   >  5:46745/1‑149 (MQ=255)                                                                                                                               gACCTTGCCCAGCTCCGGGTAATAGTGGCTCTGGATGCGATCGGCCCAGCTGCGCACGAACCAGGCTGGCTTGTTGTCGACCCGGCCCAGGCTGCCGAGCAGCAGGCGGCCCTGATCGTCGCGGCGGATGCTGCTGAGCACGGTGCgg                    <  5:144065/148‑1 (MQ=255)                                                                                                                               gACCTTGCCCAGCTCCGGGTAATAGTGGCTCTGGATGCGATCGGCCCAGCTGCGCACGAACCAGGCTGGCTTGTTGTCGACCCGGCCCAGGCTGCCGAGCAGCAGGCGGCCCTGATCGTCGCGGCGGATGCTGCTGAGCACGGTGCgg                    <  3:260174/148‑1 (MQ=255)                                                                                                                                  cTTGCCCAGCTCCGGGTAATAGTGGCTCTGGATGCGATCGGCCCAGCTGCGCACGAACCAGGCTGGCTTGTTGTCGACCCGGCCCAGGCTGCCGAGCAGCAGGCGGCCCTGATCGTCGCGGCGGATGCTGCTGAGCACGGTGCGGGTAt                >  7:250304/1‑149 (MQ=255)                                                                                                                                       ccAGCTCCGGGTAATAGTGGCTCTGGATGCGATCGGCCCAGCTGCGCACGAACCAGGCTGGCTTGTTGTCGACCCGGCCCAGGCTGCCGAGCAGCAGGCGGCCCTGATCGTCGCGGCGGATGCTGCTGAGCACGGTGCGGGTATCCCAc           <  5:170064/149‑1 (MQ=255)                                                                                                                                         aGCTCCGGGCAATAGTGGCTCTGGATGCGATCGGCCCAGCTGCGCACGAACCAGGCTGGCTTGTTGTCGACCTGGCCCAGGCTGCCGAGCAGCAGGCGGCCCTGATCGTCGCGGCGGATGCTGCTGAGCACGGTGCGGGTATCCCACGa         <  4:281433/149‑1 (MQ=255)                                                                                                                                         aGCTCCGGGCAATAGTGGCGCTGGATGCGATCGGCCCAGCTGCGCACGAACCAGGCTGGCTTGTTGTCGACCTGGCCCAGGCTGCCGAGCAGCAGGCGGCCCTGATCGTCGCGGCGGATGATGCTGAGCACGGTGCGGGTATCCCACGa         <  4:285964/149‑1 (MQ=255)                                                                                                                                              cGGGTAATAGTGGCTCTGGATGCGATCGGCCCAGCTGCGCACGAACCAGGCTGGCTTGTTGTCGACCCGGCCCAGGCTGCCGAGCAGCAGGCGGCCCTGATCGTCGCGGCGGATGCTGCTGAGCACGGTGCGGGTATCCCACGAACAtt    >  2:17659/1‑149 (MQ=255)                                                                                                                                               gggTAATAGTGGCTCTGGATGCGATCGGCCCAGCTGCGCACGAACCAGGCTGGCTTGTTGTCGACCCGGCCCAGGCTGCCGAGCAGCAGGCGGCCCTGATCGTCGCGGCGGATGCTGCTGAGCACGGTGCGGGTATCCCACGAACCTTg   <  8:246362/149‑1 (MQ=255)                                                                                                                                               gggTAATAGTGGCTCTGGATGCGATCGGCCCAGCTGCGCACGAACCAGGCTGGCTTGTTGTAGACCCGGCCCAGGCTCCagag                                                                     >  5:164555/1‑83 (MQ=255)                                                                                                                                                ggTAATAGTGGCTCTGGATGCGATCGGCCCAGCTGCGCACGAACCAGGCTGGCTTGTTGTCGACCCGGCCCAGGCTGCCGAGCAGCAGGCGGCCCTGATCGTCGCGGCGGATGCTGCTGAGCACGGTGCGGGTATCCCACGAACCTTgg  <  3:68555/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                   GTGCCGGTCGTGTTACCCCGCCCGTTGTACCCAGTCACCGCCACCAACCCAGGCGCCGGCTCGAACAGGCGCATCAGGTGGTCA‑GGGGTGAAGTCGATGCAGCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCGGGTAATAGTGGCTCTGGATGCGATCGGCCCAGCTGCGCACGAACCAGGCTGGCTTGTTGTCGACCCGGCCCAGGCTGCCGAGCAGCAGGCGGCCCTGATCGTCGCGGCGGATGCTGCTGAGCACGGTGCGGGTATCCCACGAACCTTGG  >  NC_002947/6001171‑6001461

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