Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,977,516	G→A	5.3%	L439L (TTG→TTA)	yifB →	ATP‑dependent chelatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,977,516	0	G	A	5.3%	231.2 / 2.5	75	L439L (TTG→TTA)	yifB	ATP‑dependent chelatase
Reads supporting (aligned to +/- strand):  ref base G (35/36);  new base A (2/2);  total (37/38)
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.29e-01

GAGCACTACATTGAACAACCAGCCTTGTGGTGAAACCAGTGCCGATGTCGCCGCCAGGGTGGCTGAAGCGCGGGAGCTGCAGCACAGGCGACAAGGCTGTGCGAATGCATTT‑CTCGACCTTGAAGGTTTAAGGCAT‑‑‑CACTGTGGGTTGGTAGCGGCGGACCAAGTCTGGCTGGAAGGGGCCTGCGAGCGGCTGACCTTGTCGCTACGGGCGGCGCACCGGTTGTTGAAGGTGGCGCGGACACTGGCGGACCTGGAGGGATGCGAGACGATTGGCCGGGCGCACCTGGCCGAGGCC  >  NC_002947/5977369‑5977663                                                                                                                                                        |                                                                                                                                                    gagCACTACATTGAACAACCAGCCTTGTGGTGAAACCAGTGCCGATGTCGCCGCCAGGGTGGCTGAAGCGCGGGAGCTGCAGCACAGGCGACAAGGCTGTGCGAATGCATTT‑CTCGACCTTGAAGGTTTAAGGCAT‑‑‑CACTGTGGGTTgg                                                                                                                                                    <  8:8120/149‑1 (MQ=255)    cACTACATTGAACAACCAGCCTTGTGGTGAAACCAGTGCCGATGTCGCCGCCAGGGTGGCTGAAGCGCGGGAGCTGCAGCACAGGCGACAAGGCTGTGCGAATGCATTT‑CTCGACCTTGAAGGTTTAAGGCAT‑‑‑CACTGTGGGTTGGTAg                                                                                                                                                 >  1:301573/1‑149 (MQ=255)      cTACATTGAACAACCAGCCTTGTGGTGAAACCAGTGCCGATGTCGCCGCCAGGGTGGCTGAAGCGCGGGAGCTGCAGCACAGGCGACAAGGCTGTGCGAATGCATTT‑CTCGACCTTGAAGGTTTAAGGCAT‑‑‑CACTGTGGGTTGGTAgcg                                                                                                                                               <  5:74203/149‑1 (MQ=255)          aTTGAACAACCAGCCTTGTGGTGAAACCAGTGCCGATGTCGCCGCCAGGGTGGCTGAAGCGCGGGAGCTGCAGCACAGGCGACAAGGCTGTGCGAATGCATTT‑CTCGACCTTGAAGGTTTAAGGCAT‑‑‑CACTGTGGGTTGGTAGcggcgg                                                                                                                                           <  4:7015/149‑1 (MQ=255)              aacaacCAGCCTTGTGGTGAAACCAGTGCCGATGTCGCCGCCAGGGTGGCTGAAGCGCGGGAGCTGCAGCACAGGCGACAAGGCTGTGCGAATGCATTT‑CTCGACCTTGAAGGTTTAAGGCAT‑‑‑CACTGTGGGTTGGTAGCGGCGGAcc                                                                                                                                        >  6:269181/1‑148 (MQ=255)                      gCCTTGTGGTGAAACCAGTGCCGATGTCGCCGCCAGGGTGGCTGAAGCGCGGGAGCTGCAGCACAGGCGACAAGGCTGTGCGAATGCATTT‑CTCGACCTTGAAGGTTTAAGGCAT‑‑‑CACTGTGGGTTGGTAGCGGCGGACCAAGTctggc                                                                                                                               >  3:32389/1‑149 (MQ=255)                              gtgAAACCAGTGCCGATGTCGCCGCCAGGGTGGCTGAAGCGCGGGAGCTGCAGCACAGGCGACAAGGCTGTGCGAATGCATTT‑CTCGACCTTGAAGGTTTAAGGCAT‑‑‑CACTGTGGGTTGGTAGCGGCGGACCAAGTCTGGCTGGAAggg                                                                                                                       >  1:77242/1‑149 (MQ=255)                               tgAAACCAGTGCCGATGTCGCCGCCAGGGTGGCTGAAGCGCGGGAGCTGCAGCACAGGCGACAAGGCTGTGCGAATGCATTT‑CTCGACCTTGAAGGTTTAAGGCAT‑‑‑CACTGTGGGTTGGTAGCGGCGGACCAAGTCTGGCTGGAAgggg                                                                                                                      <  8:85497/149‑1 (MQ=255)                                gAAACCAGTGCCGATGTCGCCGCCAGGGTGGCTGAAGCGCGGGAGCTGCAGCACAGGCGACAAGGCTGTGCGAATGCATTT‑CTCGACCTTGAAGGTTTAAGGCAT‑‑‑CACTGTGGGTTGGTAGCGGCGGACCAAGTCTGGCTGGAAGGGGc                                                                                                                     >  5:26542/1‑149 (MQ=255)                                 aaaCCAGTGCCGATGTCGCCGCCAGGGTGGCTGAAGCGCGGGAGCTGCAGCACAGGCGACAAGGCTGTGCGAATGCATTT‑CTCGACCTTGAAGGTTTAAGGCAT‑‑‑CACTGTGGGTTGGTAGCGGCGGACCAAGTCTGGCTGGTAGGGGc                                                                                                                     <  3:86912/148‑1 (MQ=255)                                     cAGTGCCGATGTCGCCGCCAGGGTGGCTGAAGCGCGGGAGCTGCAGCACAGGCGACAAGGCTGTGCGAATGCATTT‑CTCGACCTTGAAGGTTTAAGGCAT‑‑‑CACTGTGGGTTGGTAGCGGCGGACCAAGTCTGGCTGGAAGGGGCCTGc                                                                                                                 >  1:87847/1‑148 (MQ=255)                                          ccGATGTCGCCGCCAGGGTGGCTGAAGCGCGGGAGCTGCAGCACAGGCGACAAGGCTGTGCGAATGCATTT‑CTCGACCTTGAAGGTTTAAGGCAT‑‑‑CACTGTGGGTTGGTAGCGGCGGACCAAGTCTGGCTGGAAGGGGCCTGCGAGCgg                                                                                                           <  7:219639/149‑1 (MQ=255)                                          ccGATGTCGCCGCCAGGGTGGCTGAAGCGCGGGAGCTGCAGCACAGGCGACAAGGCTGTGCGAATGCATTT‑CTCGACCTTGAAGGTTTAAGGCAT‑‑‑CACTGTGGGTTGGTAGCGGCGGACCAAGTCTGGCTGGAAGGGGCCTGCGAGCgg                                                                                                           <  8:266029/149‑1 (MQ=255)                                              tGTCGCCGCCAGGGTGGCTGAAGCGCGGGAGCTGCAGCACAGGCGACAAGGCTGTGCGAATGCATTT‑CTCGACCTTGAAGGTTTAAGGCAT‑‑‑CACTGTGGGTTGGTAGCGGCGGACCAAGTCTGGCTGGAAGGGGCCTGCGAGCGGCTGa                                                                                                       <  7:145517/149‑1 (MQ=255)                                              tGTCGCCGCCAGGGTGGCTGAAGCGCGGGAGCTGCAGCACAGGCGACAAGGCTGTGCGAATGCATTT‑CTCGACCTTGAAGGTTTAAGGCAT‑‑‑CACTGTGGGTTGGTAGCGGCGGACCAAGTCTGGCTGGAAGGGGCCTGCGAGCGGCTGa                                                                                                       <  3:222108/149‑1 (MQ=255)                                                tCGCCGCCAGGGTGGCTGAAGCGCGGGAGCTGCAGCACAGGCGACAAGGCTGTGCGAATGCATTTGTTCGACCTTGAAGGTTTAAGGCAT‑‑‑CACTGTGGGTTGGTAGCGGCGGACCAAGTCTGGCTGGAAGGGGCCTGCGAGCGGCTGAc                                                                                                      >  1:96437/1‑149 (MQ=255)                                                tCGCCGCCAGGGTGGCTGAAGCGCGGGAGCTGCAGCACAGGCGACAAGGCTGTGCGAATGCATTT‑CTCGACCTTGAAGGTTTAAGGCAT‑‑‑CACTGTGGGTTGGTAGCGGCGGACCAAGTCTGGCTGGAAGGGGCCTGCGAGCGGCTGAcc                                                                                                     <  6:2301/149‑1 (MQ=255)                                                  gccgccAGGGTGGCTGAAGCGCGGGAGCTGCAGCAAAGGCGACAAGGATGTGCGAATGAATTT‑GTCGACCTTGAAGGTTTAAGGCAT‑‑‑CACTGTGGGTTGGTAGCGGCGGACCAAGTCTGGCTGGAAGGGGCGTGCGAGCGGCTGACCtt                                                                                                   >  4:208264/1‑149 (MQ=255)                                                     gccAGGGTGGCTGAAGCGCGGGAGCTGCAGCACAGGCGACAAGGCTGTGCGAATGCATTT‑CTCGACCTTGAAGGTTTAAGGCAT‑‑‑CACTGTGGGTTGGTAGCGGCGGACCAAGTCTGGCTGGAAGGGGCCTGCGAGCGGCTGACCTTGTc                                                                                                <  5:269181/149‑1 (MQ=255)                                                      ccAGGGTGGCTGAAGCGCGGGAGCTGCAGCACAGGCGACAAGGCTGTGCGAATGCATTT‑CTCGACCTTGAAGGTTTAAGGCAT‑‑‑CACTGTGGGTTGGTAGCGGCGGACCAAGTCTGGCTGGAAGGGGCCTGCGAGCGGCTGACCTTGTCg                                                                                               <  2:77242/149‑1 (MQ=255)                                                                  aaGCGCGGGAGCTGCAGCACAGGCGACAAGGCTGTGCGAATGCATTT‑CTCGACCTTGAAGGTTTAAGGCAT‑‑‑CACTGTGGGTTGGTAGCGGCGGACCAAGTCTGGCTGGAAGGGGCCTGCGAGCGGCTGACCTTGTCGCTACGggcgg                                                                                     <  2:210183/147‑1 (MQ=255)                                                                      gcgGGAGCTGCAGCACAGGCGACAAGGCTGTGCGAATGCATTT‑CTCGACCTTGAAGGTTTAAGGCATTAGCACTGTGGGTTGGTAGCGGCGGACCAAGTCTGGCTGGAAGGGGCCTGCGAGCGGCTGACCTTGTCGCTACGGgcggcg                                                                                   >  1:100271/1‑148 (MQ=255)                                                                          gAGCTGCAGCACAGGCGACAAGGCTGTGCGAATGCATTT‑CTCGACCTTGAAGGTTTAAGGCAT‑‑‑CACTGTGGGTTGGTAGCGGCGGACCAAGTCTGGCTGGAAGGGGCCTGCGAGCGGCTGACCTTGTCGCTACGGGCGGCGCACCGgtt                                                                           >  4:250210/1‑149 (MQ=255)                                                                           aGCTGCAGCACAGGCGACAAGGCTGTGCGAATGCATTT‑CTCGACCTTGAAGGTTTAAGGCAT‑‑‑CACTGTGGGTTGGTAGCGGCGGACCAAGTCTGGCTGGAAGGGGCCTGCGAGCGGCTGACCTTGTCGCTACGGGCGGCGCACCGgttg                                                                          <  6:61287/149‑1 (MQ=255)                                                                              tGCAGCACAGGCGACAAGGCTGTGCGAATGCATTT‑CTCGACCTTGAAGGTTTAAGGCAT‑‑‑CACTGTGGGTTGGt                                                                                                                                                   <  3:2421/73‑1 (MQ=255)                                                                              tGCAGCACAGGCGACAAGGCTGTGCGAATGCATTT‑CTCGACCTTGAAGGTTTAAGGCAT‑‑‑CACTGTGGGTTGGt                                                                                                                                                   >  4:2421/1‑73 (MQ=255)                                                                               gcagcaCAGGCGACAAGGCTGTGCGAATGCATTT‑CTCGACCTTGAAGGTTTAAGGCAT‑‑‑CACTGTGGGTTAGTAGCGGCGGACCAAGTCTGGCTGGAAGGGGCCTGCGAGCGGCTGACCTTGTc                                                                                                >  8:287456/1‑123 (MQ=255)                                                                               gcagcaCAGGCGACAAGGCTGTGCGAATGCATTT‑CTCGACCTTGAAGGTTTAAGGCAT‑‑‑CACTGTGGGTTAGTAGCGGCGGACCAAGTCTGGCTGGAAGGGGCCTGCGAGCGGCTGACCTTGTc                                                                                                <  7:287456/123‑1 (MQ=255)                                                                               gcagcaCAGGCGACAAGGCTGTGCGAATGCATTT‑CTCGACCTTGAAGGTTTAAGGCAT‑‑‑CACTGTGGGTTAGTAGCGGCGGACCAAGTCTGGCTGGAAGGGGCCTGCGAGCGGCTGACCTTGTc                                                                                                >  8:291735/1‑123 (MQ=255)                                                                               gcagcaCAGGCGACAAGGCTGTGCGAATGCATTT‑CTCGACCTTGAAGGTTTAAGGCAT‑‑‑CACTGTGGGTTAGTAGCGGCGGACCAAGTCTGGCTGGAAGGGGCCTGCGAGCGGCTGACCTTGTc                                                                                                <  7:291735/123‑1 (MQ=255)                                                                               gcagcaCAGGCGACAAGGCTGTGCGAATGCATTT‑CTAGACCTTGAAGGTTTAAGGCAT‑‑‑CACTGTGGGTTGGTAGCGGCGGACCAAGTCTGGCTGGAAGGGGCCTGCGAGCGGCTGACCTTGTCGCTACGGGCGGCGCACCGGttgttg                                                                       >  2:308281/1‑148 (MQ=255)                                                                                      aGGCGACAAGGCTGTGCGAATGCATTT‑CTCGACCTTGATGGTTTAAGGCAT‑‑‑CACTGTGGGTTGGTAGCGGCTGACCAAGTCTGTCTGGAAGGGGCCTGCGTGCGGCTGACCTTGTCGCTACGGGCGGCGCACCGGTTGTTGAAGGTGgc                                                               <  2:301573/149‑1 (MQ=255)                                                                                      aGGCGACAAGGCTGTGCGAATGCATTT‑CTCGACCTTGAAGGTTTAAGGCAT‑‑‑CACTGTGGGTTGGTAGCGGCGGACCAAGTCTGGCTGGAAGGGGCCTGCGAGCGGCTGACCTTGTCGCTACGGGCGGCGCACCGGTTGTTGAAGGTGgc                                                               >  4:75344/1‑149 (MQ=255)                                                                                        gCGACAAGGCTGTGCGAATGCATTT‑CTCGACCTTGAAGGTTTAAGGCAT‑‑‑CACTGTGGGTTGGTAGCGGCGGACCAAGTCTGGCTGGAAGGGGCCTGCGAGCGGCTGACCTTGTCGCTACGGGCGGCGCACCg                                                                              >  3:83413/1‑132 (MQ=255)                                                                                        gCGACAAGGCTGTGCGAATGCATTT‑CTCGACCTTGAAGGTTTAAGGCAT‑‑‑CACTGTGGGTTGGTAGCGGCGGACCAAGTCTGGCTGGAAGGGGCCTGCGAGCGGCTGACCTTGTCGCTACGGGCGGCGCACCg                                                                              <  4:8005/132‑1 (MQ=255)                                                                                        gCGACAAGGCTGTGCGAATGCATTT‑CTCGACCTTGAAGGTTTAAGGCAT‑‑‑CACTGTGGGTTGGTAGCGGCGGACCAAGTCTGGCTGGAAGGGGCCTGCGAGCGGCTGACCTTGTCGCTACGGGCGGCGCACCg                                                                              <  4:83413/132‑1 (MQ=255)                                                                                        gCGACAAGGCTGTGCGAATGCATTT‑CTCGACCTTGAAGGTTTAAGGCAT‑‑‑CACTGTGGGTTGGTAGCGGCGGACCAAGTCTGGCTGGAAGGGGCCTGCGAGCGGCTGACCTTGTCGCTACGGGCGGCGCACCg                                                                              >  3:8005/1‑132 (MQ=255)                                                                                        gCGACAAGGCTGTGCGAATGCATTT‑CTCGACCTTGAAGGTTTAAGGCAT‑‑‑CACTGTGGGTTGGTAGCGGCGGACCAAGTCTGGCTGGAAGGGGCCTGCGAGCGGCTGACCTTGTCGCTACGGGCGGCGCACCGGTTGTTGaa                                                                     >  7:167161/1‑141 (MQ=255)                                                                                        gCGACAAGGCTGTGCGAATGCATTT‑CTCGACCTTGAAGGTTTAAGGCAT‑‑‑CACTGTGGGTTGGTAGCGGCGGACCAAGTCTGGCTGGAAGGGGCCTGCGAGCGGCTGACCTTGTCGCTACGGGCGGCGCACCGGTTGTTGaa                                                                     <  8:167161/141‑1 (MQ=255)                                                                                        gCGACAAGGCTGTGCGAATGCATTT‑CTCGACCTTGAAGGTTTAAGGCAT‑‑‑CACTGTGGGTTGGTAGCGGCGGACCAAGTCTGGCTGGAAGGGGCCTGCGAGCGGCTGACCTTGTCGCTACGGGCGGCGCACCGGTTGTTGAAGGTGgcgc                                                             >  3:264115/1‑149 (MQ=255)                                                                                             aaGGCTGTGCGAATGCATTT‑CTCGACCTTGAAGGTTTAAGGCAT‑‑‑CACTGTGGGTTGGTAGCGGCGGACCAAGTCTGGCTGGAAGGGGCCTGCGAGCGGCTGACCTTGTCGCTACGGGCGGCGCACCGGTTGTTGAAGGTGGCGCGGaca                                                        >  5:328261/1‑149 (MQ=255)                                                                                             aaGGCTGTGCGAATGCAATT‑CTCGACCTTGAAGGTTTAAGGCAT‑‑‑CACTGTGGGTTGGTAGCGGCGGACCAAGTCTGGCTGGAAGGGGCCTGAGAGCGGCTGACCTTTTAGCTACGGGCGGCGCACCGGTTGTTGAAGGAGGCCCGGaca                                                        >  7:237414/1‑149 (MQ=255)                                                                                                  tgtgCGAATGCATTT‑CTCGACCTTGAAGGTTTAAGGCAT‑‑‑CACTGTGGGTTGGTAGCGGCGGACCAAGTCTGGCTGGAAGGGGCCTGCGAGCGGCTGACCTTGTCGCTACGGGCGGCGCACCGGTTGTTGAAGTTGGCGCGGACACTGGc                                                   >  8:230661/1‑149 (MQ=255)                                                                                                    tgCGAATGCATTT‑CTCGACCTTGAAGGTTTAAGGCAT‑‑‑CACTGTGGGTTGGTAGCGGCGGACCAAGTCTGGCTGGAAGGGGCCTGCGAGCGGCTGACCTTGTCGCTACGGGCGGCGCACCGGTTGTTGAAGGTGGCGCGGACACTGGc                                                   <  3:75344/147‑1 (MQ=255)                                                                                                       gAATGCATTT‑CTCGACCTTGAAGGTTTAAGGCAT‑‑‑CACTGTGGGTTGGTAGCGGCGGACCAAGTCTGGCTGGAAGGGGCCTGCGAGCGGCTGACCTTGTCGCTACGGGCGGCGCAGCGGTTGTTGAAGGTGGCGCGGACACTGGCGGAcc                                              >  8:166752/1‑149 (MQ=255)                                                                                                        aaTGCATTT‑CTCGACCTTGAAGGTTTAAGGCAT‑‑‑CACTGTGGGTTGGTAGCGGCGGACCAAGTCTGGCTGGAAGGGGCCTGCGAGCGGCTGACCTTGTCGCTACGGGCGGCGCACCGGTTGTTGAAGGTGGCGCGGACACTg                                                     <  7:147550/141‑1 (MQ=255)                                                                                                        aaTGCATTT‑CTCGACCTTGAAGGTTTAAGGCAT‑‑‑CACTGTGGGTTGGTAGCGGCGGACCAAGTCTGGCTGGAAGGGGCCTGCGAGCGGCTGACCTTGTCGCTACGGGCGGCGCACCGGTTGTTGAAGGTGGCGCGGACACTg                                                     >  8:147550/1‑141 (MQ=255)                                                                                                        aaTGCATTT‑CTCGACCTTGAAGGTTTAAGGCAT‑‑‑CACTGTGGGTTGGTAGCGGCGGACCAAGTCTGGCTGGAAGGGGCCTGCGAGCGGCTGACCTTGTCGCTACGGGCGGCGCACCGGTTGTTGAAGGTGGCGCGGACACTGGCGGACCt                                             <  4:264115/149‑1 (MQ=255)                                                                                                         aTGCATTT‑CTCGACCTTGAAGGTTTAAGGCAT‑‑‑CACTGTGGGTTGGTAGCGGCGGACCAAGTCTGGCTGGAAGGGGCCTGCGAGCGGCTGACCTTGTCGCTACGGGCGGCGCACCGGTTGTTGAAGGTGGCGCGGACACTGGCGGACCTg                                            >  5:291701/1‑149 (MQ=255)                                                                                                          tGCATTT‑CTCGACCTTGAAGGTTTAAGGCAT‑‑‑CACTGTGGGTTGGTAGCGGCGGACCAAGTCTGGCTGGAAGGGGCCTGCGAGCGGCTGACCTTGTCGCTACGGGCGGCGCACCGGTTGTTGAAGGTGGCGCGGACACTGGCGGACCTgg                                           <  4:32389/149‑1 (MQ=255)                                                                                                            cATTT‑CTCGACCTTGAAGGTTTAAGGCAT‑‑‑CACTGTGGGTTGGTAGCGGCGGACCAAGTCTGGCTGGAAGGGGCCTGCGAGCGGCTGACCTTGTCGCTGCGGGCGGCGCACCGGTTGTTGAAGGTGGCGCGGACACTGGCGGACCTGGAg                                         <  2:87847/149‑1 (MQ=255)                                                                                                             aTTT‑CTCGACCTTGAAGGTTTAAGGCAT‑‑‑CACTGTGGGTTGGTAGCGGCGGACCAAGTCTGGCTGGAAGGGGCCTGCGAGCGGCTGACCTTGTCGCTACGGGCGGCGCACCGGTTGTTGAAGGTGGcgcg                                                            <  6:284971/129‑1 (MQ=255)                                                                                                             aTTT‑CTCGACCTTGAAGGTTTAAGGCAT‑‑‑CACTGTGGGTTGGTAGCGGCGGACCAAGTCTGGCTGGAAGGGGCCTGCGAGCGGCTGACCTTGTCGCTACGGGCGGCGCACCGGTTGTTGAAGGTGGcgcg                                                            >  5:284971/1‑129 (MQ=255)                                                                                                             aTTT‑CTCGACCTTGAAGGTTTAAGGCAT‑‑‑CACTGTGGGTTGGTAGCGGCGGACCAAGTCTGGCTGGAAGGGGCCTGCGAGCGGCTGACCTTGTCGCTACGGGCGGCGCACCGGTTGTTGAAGGTGGCGCGGACACTGGCGGACCTGGAgg                                        <  3:250210/149‑1 (MQ=255)                                                                                                              ttt‑CTCGACCTTGAAGGTTTAAGGCAT‑‑‑CACTGTGGGTTGGTAGCGGCGGACCAAGTCTGGCTGGAAGGGGCCTGCGAGCGGCTGACCTTGTCGCTACGGGCGGCGCACCGGTTGTTGAAGGTGGCGCGGACACTGGCgg                                                 >  2:37748/1‑139 (MQ=255)                                                                                                              ttt‑CTCGACCTTGAAGGTTTAAGGCAT‑‑‑CACTGTGGGTTGGTAGCGGCGGACCAAGTCTGGCTGGAAGGGGCCTGCGAGCGGCTGACCTTGTCGCTACGGGCGGCGCACCGGTTGTTGAAGGTGGCGCGGACACTGGCgg                                                 <  1:37748/139‑1 (MQ=255)                                                                                                              ttt‑CTCGACCTTGAAGGTTTAAGGCAT‑‑‑CACTGTGGGTTGGTAGCGGCGGACCAAGTCTGGCTGGAAGGGGCCTGCGAGCGGCTGACCTTGTCGCTACGGGCGGCGCACCGGTTGTTGAAGGTGGCGCGGACACTGGCGGACCTGGAgg                                        >  7:23487/1‑148 (MQ=255)                                                                                                                    cGACCTTGAAGGTTTAAGGCAT‑‑‑CACTGTGGGTTGGTAGCGGCGGACCAAGTCTGGCTGGAAGGGGCCTGCGAGCGGCTGACCTTGTCGCTACGGGCGGCGCACCGGttgttg                                                                       >  7:217816/1‑112 (MQ=255)                                                                                                                    cGACCTTGAAGGTTTAAGGCAT‑‑‑CACTGTGGGTTGGTAGCGGCGGACCAAGTCTGGCTGGAAGGGGCCTGCGAGCGGCTGACCTTGTCGCTACGGGCGGCGCACCGGttgttg                                                                       <  8:217816/112‑1 (MQ=255)                                                                                                                      aCCTTGAAGGTTTAAGGCAT‑‑‑CACTGTGGGTTGGTAGCGGCGGACCAAGTCTGGCTGGAAGGGGCCTGCGAGCGGCTGACCTTGTCGCTACGGGCGGCGCACCGGTTGTTGAAGGTGGCGCGGACACTg                                                     >  3:128124/1‑128 (MQ=255)                                                                                                                      aCCTTGAAGGTTTAAGGCAT‑‑‑CACTGTGGGTTGGTAGCGGCGGACCAAGTCTGGCTGGAAGGGGCCTGCGAGCGGCTGACCTTGTCGCTACGGGCGGCGCACCGGTTGTTGAAGGTGGCGCGGACACTg                                                     <  4:128124/128‑1 (MQ=255)                                                                                                                      aCCTTGAAGGTTTAAGGCAT‑‑‑CACTGTGGGTTGGTAGCGGCGGACCAAGTCTGGCTGGAAGGGGCCTGCGAGCGGCTGACCTTGTCGCTACGGGCGGCGCACCGGTTGTTGAAGGTGGCGCGGACACTGGCGGACCTGGAGGGATGc                                   >  4:248219/1‑146 (MQ=255)                                                                                                                         ttGAAGGTTTAAGGCAT‑‑‑CACTGTGGGTTGGTAGCGGCGGACCAAGTCTGGCTGGAAGGGGCCTGCGAGCGGCTGCCCTTGTCGCTACGGGCGGCGCACCGGTTGTTGAAGGTGGCGCGGACACTGGCGGACCGGGAGGGATGCGAGAc                              >  8:67359/1‑148 (MQ=255)                                                                                                                           gAAGGTTTAAGGCAT‑‑‑CACTGTGGGTTGGTAGCGGCGGACCAAGTCTGGCTGGAAGGGGCCTGCGAGCGGCTGACCTTGTCGCTACGGGCGGCGCACCGGTTGTTGAAGGTGGCGCGGACACTGGCGGACCTGGAGGGATGCGAGACGAt                           >  7:14638/1‑149 (MQ=255)                                                                                                                               ttttATGGCAT‑‑‑GACTGTGGGTTGGTAGCGGCGGACCAAGGGTGGCTGGAAGGGGCCTGCGAGCGGGTGACCTTGTCGCTACGGGCGGCGCACCGGTTGTTGAAGGTGGCGCGGACACTGGCGGACCTGGAGGGATGCGAGACGATTgg                        <  8:237414/147‑1 (MQ=255)                                                                                                                                     ggCAT‑‑‑CACTGTGGGTTGGTAGCGGCGGACCAAGTCTGGCTGGAAGGGGCCTGCGAGCGGCTGACCTTGTCGCTACGGGCGGCGCACCGGTTGTTGAAGGTGGCGCGGACACTGGCGGACCTGGAGGGATGCGAGACGATTGGCCGGgcg                 >  3:148069/1‑149 (MQ=255)                                                                                                                                      gCAT‑‑‑CACTGTGGGTTGGTAGCGGCGGACCAAGTCTGGCTGGAAGGGGCCTGCGAGCGGCTGAcc                                                                                                     >  1:301115/1‑64 (MQ=255)                                                                                                                                      gCAT‑‑‑CACTGTGGGTTGGTAGCGGCGGACCAAGTCTGGCTGGAAGGGGCCTGCGAGCGGCTGAcc                                                                                                     <  2:301115/64‑1 (MQ=255)                                                                                                                                      gCAT‑‑‑CACTGTGGGTTGGTAGCGGCGGACCAAGTCTGGCTGGAAGGGGCCTGCGAGCGGCTGACCTTGTCGCTACGGGCGGCGCACCGGTTGTTGAAGGTGGCGCGGACACTGGCGGACCTGGAGGGATGCGAGACGATTGGCCGGgcgc                <  3:248219/149‑1 (MQ=255)                                                                                                                                             cACTGTGGGTTGGTCGCGGCGGCCCACGTCTGGCGGGAAGGGGCCTGCGAGCGGCTGACCTTGTCGCTACGGGCGGCGCACCGGTGGTTGACGGTGGCGCGGACACTGGCGGACCTGGAGGGATGCGAGACGATTGGCCGGGCGCAc              <  7:67359/147‑1 (MQ=255)                                                                                                                                              aCTGTGGGTTGGTAGCGGCGGACCAAGTCTGGCTGGAAGGGGCCTGCGAGCGGCTGACCTTGTCGCTACGGGCGGCGCACCGGTTGTTGAAGGTGGCGCGGACACTGGCGGACCTGGAGGGATGCGAGACGATTGGCCGGGCGCACCTg           <  2:96437/149‑1 (MQ=255)                                                                                                                                               cTGTGGGTTGGTAGCGGCGGACCAAGTCTGGCTGGAAGGGGCCTGCGAGCGGCTGACCTTGTCGCTACGGGCGGCGCACCGGTTGTTGAAGGTGGCGCGGACACTGGCGGACCTGGAGGGATGCGAGACGATTGGCCGGGCGCACCTgg          >  4:41168/1‑149 (MQ=255)                                                                                                                                                    ggTTGGTAGCGGCGGACCAAGTCTGGCTGGAAGGGGCCTGCGAGCGGCTGACCTTGTCGCTACGGGCGGCGCACCGGTTGTTGAAGGTGGCGCGGACACTGGCGGACCTGGAGGGATGCGAGACGATTGGCCGGGCGCACCTGGCCGa      <  2:260247/148‑1 (MQ=255)                                                                                                                                                    ggTTGGTAGCGGCGGACCAAGTCTGGCTGGAAGGGGCCTGCGAGCGGCTGACCTTGTCGCTACGGGCGGCGCACCGGTTGTTGAAGGTGGCGCGGACACTGGCGGACCTGGAGGGATGCGAGACGATTGGCCGGGCGCACCTGGCCGAg     <  7:166752/149‑1 (MQ=255)                                                                                                                                                       tGGTAGCGGCGGACCAAGTCTGGCTGGAAGGGGCCTGCGAGCGGCTGACCTTGTCGCTACGGGCGGCGCACCGGTTGTTGAAGGTGGCGCGGACACTGGCGGACCTGGAGGGATGCGAGACGATTGGCCGGGCGCACCTGGCCGAGGcc  >  8:220443/1‑149 (MQ=255)                                                                                                                                                       tGGTAGCGGCGGACCAAGTCTGGCTGGAAGGGGCCTGCGAGCGGCTGACCTTGTCGCTACGGGCGGCGCACCGGTTGTTGAAGGTGGCGCGGACACTGGCGGACCTGGAGGGATGCGAGACGATTGGCCGGGCGCACCTGGCCGAGGcc  <  6:26542/149‑1 (MQ=255)                                                                                                                                                        ggTAGCGGCGGACCAAGTCTGGCTGGAAGGGGCCTGCGAGCGGCTGACCTTGTCGCTACGGGCGGCGCACCGGTTGTTGAAGGTGGCGCGGACACTGGCGGACCTGGAGGGATGCGAGACGATTGGCCGGGCGCACCTGGCCGAGGcc  <  1:145097/148‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                                    GAGCACTACATTGAACAACCAGCCTTGTGGTGAAACCAGTGCCGATGTCGCCGCCAGGGTGGCTGAAGCGCGGGAGCTGCAGCACAGGCGACAAGGCTGTGCGAATGCATTT‑CTCGACCTTGAAGGTTTAAGGCAT‑‑‑CACTGTGGGTTGGTAGCGGCGGACCAAGTCTGGCTGGAAGGGGCCTGCGAGCGGCTGACCTTGTCGCTACGGGCGGCGCACCGGTTGTTGAAGGTGGCGCGGACACTGGCGGACCTGGAGGGATGCGAGACGATTGGCCGGGCGCACCTGGCCGAGGCC  >  NC_002947/5977369‑5977663

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