Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,625,537	T→A	11.6%	intergenic (‑50/+27)	PP_3195 ← / ← PP_3196	hypothetical protein/nodulation protein d1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,625,537	0	T	A	11.6%	169.2 / 10.7	67	intergenic (‑50/+27)	PP_3195/PP_3196	hypothetical protein/nodulation protein d1
Reads supporting (aligned to +/- strand):  ref base T (26/33);  new base A (4/4);  total (30/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 = 9.98e-01

CTTCTGCGGTCGTGGAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACAT‑CGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTCGCCGTCCATGCTGCGGTGCCACAGCATGA  >  NC_002947/3625393‑3625682                                                                                                                                                  |                                                                                                                                                   cTTCTGCGGTCGTGGAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAGACACCGCAGCg                                                                                                                                                 >  5:305233/1‑149 (MQ=255) cTTCTGCGGTCGTGGAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCgg                                                                                                                                                >  5:151824/1‑149 (MQ=255)    cTGCGGTCGTGGAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGGGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGgt                                                                                                                                               <  5:45729/147‑1 (MQ=255)    cTGCGGTCGTGGAGTCCGGCGCCCACTTGAGAGACAGGCAGGCAGGGCGCGGTTTGATGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCATTCGTGATCAGGGATGCGAGGTATCGGGTAACGCGATGCCCCAG‑CACCGCTGCGgt                                                                                                                                               <  2:326363/147‑1 (MQ=255)       cGGTCGTGGAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTg                                                                                                                                          <  7:19963/149‑1 (MQ=255)        ggTCGTGGAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAGACACCGCAGCGgtgctgg                                                                                                                                          <  7:2575/149‑5 (MQ=255)        ggTCGTGGAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAGACACCGCAGCGgtgctgg                                                                                                                                          <  6:305233/149‑5 (MQ=255)               gAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAGACACCGCAGCGGTG‑CTGGGGcatcg                                                                                                                                  <  1:174771/149‑5 (MQ=255)                  tCCGGCGCCCACTTGAGCGACAGGCAGGCAGGTCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCCGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGgcc                                                                                                                               <  6:130885/149‑1 (MQ=255)                    cGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGgccgc                                                                                                                             >  5:49943/1‑149 (MQ=255)                    cGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGgccgc                                                                                                                             >  6:122069/1‑149 (MQ=255)                          ccACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCt                                                                                                                         >  6:187479/1‑147 (MQ=255)                            aCTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAgg                                                                                                                     >  3:77467/1‑149 (MQ=255)                                     aCAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGgcgc                                                                                                            <  6:234840/149‑1 (MQ=255)                                        ggcaggcaggGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCa                                                                                                         <  8:44421/149‑1 (MQ=255)                                         gcgggcaggGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCGGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGACGCCTCTCAGGAGGTGgcgcg                                                                                                           >  3:37016/4‑146 (MQ=255)                                          caggcaggGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCaga                                                                                                       <  5:372146/149‑1 (MQ=255)                                           aggcaggGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCagag                                                                                                      <  5:42962/149‑1 (MQ=255)                                             gcaggGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAACCGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCa                                                                                                    >  4:141135/1‑149 (MQ=255)                                                gggCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCt                                                                                                 <  7:356486/149‑1 (MQ=255)                                                 ggCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTg                                                                                                >  2:164945/1‑149 (MQ=255)                                                  gcgcGGTTTGAATGCAAGGTACGCACGCAGGGCAGGGTGGGTCCAATCGTGATCGGGGATCCGAGGGATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTgg                                                                                               <  4:37016/149‑1 (MQ=255)                                                    gcgGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCt                                                                                             <  6:343453/149‑1 (MQ=255)                                                       gTTTGCAGGCAAGGTACGAAGTCACGGCAGGGTGGGTCCAATCCTGATCAGGGATTCGAGTTATCGAGATGCGCGATGCCCCAG‑CACCGCTGCTGTGTCTGCGGCGATGGCAGAAGCTCATGCGGTGCCGCACAGAGCATCTGGCtcag                                                                                          >  6:189657/1‑146 (MQ=255)                                                            aaGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTg                                                                                     >  6:52511/1‑149 (MQ=255)                                                              ggCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGgcgcgc                                                                                                          >  5:88034/1‑126 (MQ=255)                                                              ggCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGgcgcgc                                                                                                          <  6:88034/126‑1 (MQ=255)                                                              ggCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGgc                                                                                   >  7:17552/1‑149 (MQ=255)                                                              ggCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGgc                                                                                   >  8:45673/1‑149 (MQ=255)                                                              ggCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGAGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGATCAGGCGGTGGCTCGCATAGCATATGGCTGCGTGCTGgc                                                                                   >  8:108957/1‑149 (MQ=255)                                                                 aaGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGGTGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCa                                                                                >  1:74419/1‑149 (MQ=255)                                                                       cggacgcagggcaggGTGGGGCCAATCGTGATCCGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGaa                                                                           <  8:360853/143‑1 (MQ=255)                                                                           ggcagggcagGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACa                                                                      <  7:45673/149‑1 (MQ=255)                                                                                  cagggTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACt                                                                              >  1:292439/1‑134 (MQ=255)                                                                                  cagggTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACt                                                                              <  2:292439/134‑1 (MQ=255)                                                                                              aaTCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACAT‑c                                                                   <  4:233499/132‑1 (MQ=255)                                                                                              aaTCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACAT‑c                                                                   >  3:233499/1‑132 (MQ=255)                                                                                              aaTCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACAT‑CGCTGATGCGCTCGCGaa                                                  <  5:122069/149‑1 (MQ=255)                                                                                              aaTCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACAT‑CGCTGATGCGCTCGCGaa                                                  <  4:172484/149‑1 (MQ=255)                                                                                                tCGTGATCATTGATGCGAGGTATCTCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACAT‑CGCTGATGCGCTCGCGAAGc                                                <  6:9057/149‑1 (MQ=255)                                                                                                   tGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACAT‑CGCTGATGCGCTCGCGAAGCCAt                                             <  6:49943/149‑1 (MQ=255)                                                                                                     aTCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACAT‑CGCTGATGCGCTCGCGAAGCCATCt                                           >  2:301233/1‑149 (MQ=255)                                                                                                     aTCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACAT‑CGCTGATGCGCTCGCGAAGCCATCt                                           >  8:167881/1‑149 (MQ=255)                                                                                                                 aGGTAGCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCTATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACGGAATGACAT‑CGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTc                               <  6:282756/149‑1 (MQ=255)                                                                                                                    tATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACAT‑CGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTCGcc                            <  7:108957/149‑1 (MQ=255)                                                                                                                    tATCGCGAAACGCGATGCCCCAG‑CACAGCAGCGGTGTCTGGGGCGATTGCCGCCGCTCAGGCGGTGACGAGCAGAGCATCTGGATGCGTGCTGGCGCACTGAATGACAT‑AGCTGATGCGCTCGCGAAGCAATCTGTGCATGGGGTCGcc                            >  8:354376/1‑149 (MQ=255)                                                                                                                      tCGCGAAACGCGATGCCCCAG‑CACCACTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACAT‑CGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTCGCCGt                          >  7:317550/1‑149 (MQ=255)                                                                                                                        gcgAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACAT‑CGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTCGCCg                           <  1:346692/146‑1 (MQ=255)                                                                                                                        gcgAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACAT‑CGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTCGCCg                           >  2:346692/1‑146 (MQ=255)                                                                                                                        gcgAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACAT‑CGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTCGCCGTc                         >  6:15841/1‑148 (MQ=255)                                                                                                                            aaCGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACAT‑CGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTCGCCGTCCAtgc                    >  3:3458/1‑149 (MQ=255)                                                                                                                                cgATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACAT‑CGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTCGCCGTCCATGCTGCg                <  6:151824/149‑1 (MQ=255)                                                                                                                                   tGCGCCAG‑CTCCGCTGCGGTGTGTTGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGTGCATCTGGCTGCGTGCTGGCGCACTGAATGACAT‑CGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTCGCCGTCCATGCTGCGGTg             <  4:77467/149‑1 (MQ=255)                                                                                                                                   tGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACAT‑CGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTCGCCGTCCATGCTGCGGTg             <  3:141135/149‑1 (MQ=255)                                                                                                                                     ccccAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACAT‑CGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTCGCCGTCCATGCTGCGGTGc            <  1:293329/148‑1 (MQ=255)                                                                                                                                      cccagacACCGCAGCGGTG‑CTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACAT‑CGCTGATGCGCTCGCGAAGCCATCTCTGCATGGGGTCGCCGTCCATGCTGCGGTGCca          >  7:340661/7‑149 (MQ=255)                                                                                                                                      cccAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACAT‑CGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTCGCCGTCCATGCTGCGGTGCca          <  5:52511/149‑1 (MQ=255)                                                                                                                                       ccagacACCGCAGCGGTG‑CTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACAT‑CGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTCGCCGTCCATGCTGCGGTGCca          >  4:237001/6‑148 (MQ=255)                                                                                                                                       ccAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCCGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTCGCCGTCCATGCTGCGGTGCca          <  1:301233/149‑1 (MQ=255)                                                                                                                                       ccAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACAT‑CGc                                                                 <  1:109715/93‑1 (MQ=255)                                                                                                                                       ccAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACAT‑CGc                                                                 >  2:109715/1‑93 (MQ=255)                                                                                                                                         aG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACAT‑CGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTCGCCGTCCATGCTGCGGTGCcaca        <  5:156636/148‑1 (MQ=255)                                                                                                                                         aG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACAT‑AGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTCGCCGTCCATGCTGCGGTGCca          >  6:156636/1‑146 (MQ=255)                                                                                                                                          gacACCGCAGCGGTG‑CTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACAT‑CGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTCGCCGTCCATGCTGCGGTGCCACAGc      <  3:237001/147‑1 (MQ=255)                                                                                                                                             aCCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACAT‑CGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTCGCCGTCCATGCTGCGGTGCCACAGCAt    <  7:167881/149‑1 (MQ=255)                                                                                                                                              ccGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACAT‑CGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTCGCCGTCCATGCTGCGGTGCCACAGCAt    >  4:159481/1‑148 (MQ=255)                                                                                                                                               cGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACAT‑CGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTCGCCGTCCATGCTGCGGTGCCACAGCATGa  <  1:164945/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                   CTTCTGCGGTCGTGGAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACAT‑CGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTCGCCGTCCATGCTGCGGTGCCACAGCATGA  >  NC_002947/3625393‑3625682

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