Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,305,576	G→T	6.1%	A559D (GCC→GAC)	recC ←	multifunctional ATP‑dependent DNA helicase/dsDNA/ssDNA exonuclease/ssDNA endonuclease recombinase subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,305,576	0	G	T	6.1%	209.6 / 3.6	66	A559D (GCC→GAC)	recC	multifunctional ATP‑dependent DNA helicase/dsDNA/ssDNA exonuclease/ssDNA endonuclease recombinase subunit
Reads supporting (aligned to +/- strand):  ref base G (29/33);  new base T (2/2);  total (31/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.19e-01

ATGCAGTCGCTCGCCCCACTGGCTGACGGTGGCAGGCTCGGACAAGGCCTGGCAGGCGACTTCGAGCGCATCCAGCAACGCGACCAGCGGGCCGATCAGTGCCGCGTCCAGGCCACCGATTTCATCGTAGGGTTCAATGCCATCACAGGCCTCACCCACCCCCACGGCATAACCCAACAGCATGCGGCGCAG‑GCCGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGCGCGGCGTCAAGGCCCCAACGGATGCCAGCGCCCTCGATCCAGCGGTG  >  NC_002947/5305428‑5305724                                                                                                                                                     |                                                                                                                                                      aTGCAGTCGCTCGCCCCACTGGCTGACGGTGGCAGGCTCGGACAAGGCCTGGCAGGCGACTTCGAGCGCATCCAGCAACGCGACCAGCGGGCCGATCAGTGCCGCGTCCAGGCCACCGATTTCATCGTAGGGTTCAATGCCATCACAgg                                                                                                                                                       >  6:342141/1‑149 (MQ=255)     aGTCGCTCGCCCCACTGGCTGACGGTGGCAGGCTCGGACAAGGCCTGGCAGGCGACTTCGAGCGCATCCAGCAACGCGACCAGCGGGCCGATCAGTGCCGCGTCCAGGCCACCGATTTCATCGTAGGGTTCAATGCCATCACAGGCCt                                                                                                                                                    >  2:90458/1‑148 (MQ=255)      gTCGCTCGCCCCACTGGCTGACGGTGGCAGGCTCGGACAAGGCCTGGCAGGCGACTTCGAGCGCATCCAGCAACGCGACCAGCGGGCCGATCAGTGCCGCGTCCAGGCCACCGATTTCATCGTAGGGTTCAATGCCATCACAGGCCTca                                                                                                                                                  <  4:40765/149‑1 (MQ=255)      gTCGCTCGCCCCACTGGCTGACGGTGGCAGGCTCGGACAAGGCCTGGCAGGCGACTTCGAGCGCATCCAGCAACGCGACCAGCGGGCCGATCAGTGCCGCGTCCAGGCCACCGATTTCATCGTAGGGTTCAATGCCATCACAGGCCTca                                                                                                                                                  >  3:286429/1‑149 (MQ=255)      gTCGCTCGCCCCACTGGCGGACGGTGGCAGGCTCGGACAAGGCCTGGCAGGCGACTTCGAGCGCATCCAGCAACGCGACCGGCGGGCCGATCAGTGCCGCGTCCAGGCCACCGATTTCATCGTAGGGTTCAATGCCATCACAGGCCTca                                                                                                                                                  <  4:180276/149‑1 (MQ=255)       tcgctcgcCCCACTGGCTGACGGTGGCAGGCTCGGACAAGGCCTGGCAGGCGACTTCGAGCGCATCCAGCAACGCGACCAGCGGGCCGATCAGTGCCGCGTCCAGGCCACCGATTTCATCGTAGGGTTCAATGCCATCACAGGCCTcac                                                                                                                                                 <  7:145524/149‑1 (MQ=255)       tcgctcgcCCCACTGGCTGACGGTGGCAGGCTCGGACAAGGCCTGGCAGGCGACTTCGAGCGCATCCAGCAACGCGACCAGCGGGCCGATCAGTGCCGCGTCCAGGCCACCGATTTCATCGTAGGGTTCAATGCCATCACAGGCCTcac                                                                                                                                                 >  8:145524/1‑149 (MQ=255)              ccccACTGGCTGACGGTGGCAGGCTCGGACAAGGCCTGGCAGGCGACTTCGAGCGCATCCAGCAACGCGACCAGCGGGCCGATCAGTGCCGCGTCCAGGCCACCGATTTCATCGTAGGGTTCAATGCCATCACAGGCCTCACCCAcccc                                                                                                                                          >  1:340950/1‑149 (MQ=255)                  aCTGGCTGACGGTGGCAGGCTCGGACAAGGCCTGGCAGGCGACTTCGAGCGCATCCAGCAACGCGACCAGCGGGCCGATCAGTGCCGCGTCCAGGCCACCGATTTCATCGTAGGGTTCAATGCCATCACAGGCCTCACCCACCCCCACg                                                                                                                                      >  2:74205/1‑149 (MQ=255)                    tGGCTGACGGTGGCAGGCTCGGACAAGGCCTGGCAGGCGACTTCGAGCGCATCCAGCAACGCGACCAGCGGGCCGATCAGTGCCGCGTCCAGGCCACCGATTTCATCGTAGGGTTCAATGCCATCACAGGCCTCACCCACCCCCACgg                                                                                                                                     >  8:313039/1‑148 (MQ=255)                      gCTGACGGTGGCAGGCTCGGACAAGGCCTGGCAGGCGACTTCGAGCGCATCCAGCAACGCGACCAGCGGGCCGATCAGTGCCGCGTCCAGGCCACCGATTTCATCGTAGGGTTCAATGCCATCACAGGCCTCACCCACCCCCACGGc                                                                                                                                    <  1:90458/147‑1 (MQ=255)                          aCGGTGGCAGGCTCGGACAAGGCCTGGCAGGCGACTTCGAGCGCATCCAGCAACGCGACCAGCGGGCCGATCAGTGCCGCGTCCAGGCCACCGATTTCATCGTAGGGTTCAATGCCATCACAGGCCTCACCCACCCCCACGGCATaa                                                                                                                                <  4:172812/147‑1 (MQ=255)                             gTGGCAGGCTCGGACAAGGGCTGGCAGGCGACTTCGAGCTCATCCAGCAACGCGACCAGCGGGCCGATCCGTGCCGCGTCCAGGCCACCGATTTCATCGTAGGGTTCAATGCCATCACAGGCCTCACCCACCCCCACGGCATAACCCaa                                                                                                                           <  6:149468/149‑1 (MQ=255)                             gTGGCAGGCTCGGACAAGGCCTGGCAGGCGACTTCGAGCGCATCCAGCAACGCGACCAGCGGGCCGATCAGTGCCGCGTCCAGGCCACCGATTTCATCGTAGGGTTCAATGCCATCACAGGCCTCACCCACCCCCACGGCATAACCCaa                                                                                                                           <  7:313039/149‑1 (MQ=255)                                  aGGCTCGGACAAGGCCTGGCAGGCGACTTCGAGCGCATCCAGCAACGCGACCAGCGGGCCGATCAGTGCCGCGTCCAGGCCACCGATTTCATCGTAGGGTTCAATGCCATCACAGGCCTCACCCACCCCCACGGCATAACCCAACAGCa                                                                                                                      <  5:94006/149‑1 (MQ=255)                                       cGGACAAGGCCTGGCAGGCGACTTCGAGCGCATCCAGCAACGCGACCAGCGGGCCGATCAGTGCCGCGTCCAGGCCACCGATTTCATCGTAGGGTTCAATGCCATCACAGGCCTCACCCCCCCCCACGGCATAACCCAACAGCATgcgg                                                                                                                 >  3:61379/1‑149 (MQ=255)                                        ggACAAGGCCTGGCAGGCGACTTCGAGCGCATCCAGCAACGCGACCAGCGGGCCGATCAGTGCCGCGTCCAGGCCACCTATTTCATCGTAGGGTTCAATGCCATCACAGGCCTCACCCACCCCCACGGCATAACCCAACAGCATgcggc                                                                                                                >  5:245501/1‑149 (MQ=255)                                          aCAAGGCCTGGCAGGCGACTTCGAGCGCATCCAGCAACGCGACCAGCGGGCCGATCAGTGCCGCGTCCAGGCCACCGATTTCATCGTAGGGTTCAATGCCATCACAGGCCTCACCCACCCCCACGGCATAACCCAACAGCATGCGgcgc                                                                                                              <  2:114711/149‑1 (MQ=255)                                            atggCCTGGCAGGGGAATTAGAGCTCATCCAGCAACGCGACCAGCTTGCCGATCAGTGCCGCGTCCAGGCCACCGATTTCATCGTAGGGTTCAATGCCATCACAGGCCTCACCCACCCCCACGGCATAACCCAACAGCATGCGGCGCa                                                                                                             <  7:142453/146‑1 (MQ=255)                                            aaGGCCTGGCAGGCGACTTCGAGCGCATCCAGCAACGCGACCAGCGGGCCGATCAGTGCCGCGTCCAGGCCACCGATTTCATCGTAGGGTTCAATGCCATCACAGGCCTCACCCACCCCCACGGCATAACCCAACAGCATGCGGCGCAg                                                                                                            <  2:53866/149‑1 (MQ=255)                                               gCCTGGCAGGCGACTTCGAGCGCATCCAGCAACGCGACCAGCGGGCCGATCAGTGCCGCGTCCAGGCCACCGATTTCATCGTAGGGTTCAATGCCATCACAGGCCTCACCCACCCCCACGGCATAACCCAACAGCATGCGGCGCAG‑Gcc                                                                                                        <  5:342141/149‑1 (MQ=255)                                               gCCTGGCAGGCGACTTCGAGCGCATCCAGCAACGCGACCAGCGGGCCGATCAGTGCCGCGTCCAGGCCACCGATTTCATCGTAGGGATCAATGCCATCACAGGCCTCACCCACCCCCACGGCATAACCCAACAGCATGCGGCGCAG‑Gcc                                                                                                        >  2:88491/1‑149 (MQ=255)                                                 cTGGCAGGCGACTTCGAGCGCATCCAGCAACGCGACCAGCGGGCCGATCAGTGCCGCGTCCAGGCCACCGATTTCATCGTAGGGTTCAATGCCATCACAGGCCTCACCCACCCCCACGGCATAACCCAACAGCATGCGGCGCAG‑GCCGa                                                                                                      <  6:248532/149‑1 (MQ=255)                                                    gCAGGCGACTTCGAGCGCATCCAGCAACGCGACCAGCGGGCCGATCAGTGCCGCGTCCAGGCCACCGATTTCATCGTAGGGTTCAATGCCATCACAGGCCTCACCCACCCCCACGGCATAACCCAACAGCATGCGGCGCAG‑GCCGaaa                                                                                                    <  8:265360/148‑1 (MQ=255)                                                     cAGGCGACTTCGAGCGCATCCAGCAACGCGACCAGCGGGCCGATCAGTGCCGCGTCCAGGCCACCGATTTCATCGTAGGGTTCAATGCCATCACAGGCCTCACCCACCCCCACGGCATAACCCAACAGCATGCGGCGCAG‑GCCGAAACg                                                                                                  <  4:74002/149‑1 (MQ=255)                                                             ttCGAGCGCATCCAGCAACGCGACCAGCGGGCCGATCAGTGCCGCGTCCAGGCCACCGATTTCATCGTAGGGTTCAATGCCATCACAGGCCTCACCCACCCCCACGGCATAACCCAACAGCATGCGGCGCAg                                                                                                            >  3:28066/1‑132 (MQ=255)                                                             ttCGAGCGCATCCAGCAACGCGACCAGCGGGCCGATCAGTGCCGCGTCCAGGCCACCGATTTCATCGTAGGGTTCAATGCCATCACAGGCCTCACCCACCCCCACGGCATAACCCAACAGCATGCGGCGCAg                                                                                                            <  4:28066/132‑1 (MQ=255)                                                               cGAGCGCATCCAGCAACGCGACCAGCGGGCCGATCAGTGCCGCGTCCAGGCCACCGATTTCATCGTAGGGTTCAATGCCATCACAGGCCTCACCCACCCCCACGGCATAACCCAACAGCATGCGGCGCAG‑GCCGAAACGCCAActgtt                                                                                         <  1:74205/148‑1 (MQ=255)                                                                gAGCGCATCCAGCAACGCGACCAGCGGGCCGATCAGTGCCGCGTCCAGGCCACCGATTTCATCGTAGGGTTCAATGCCATCACAGGCCTCACCCACCCCCACGGCATAACCCAACAGCATGCGGCGCAG‑GCCGAAACGCCAActgttct                                                                                       >  6:16312/1‑149 (MQ=255)                                                                  gcgcATCCAGCAACGCGACCAGCGGGCCGATCAGTGCCGCGTCCAGGCCACCGATTTCATCGTAGGGTTCAATGCCATCCCAGGCCTCACCCCCCCCCACGGCAGAACCCAACAGCATGCGGCGCAG‑GCCGAAACGCCAActgttctgt                                                                                     >  6:114837/1‑149 (MQ=255)                                                                  gcgcATCCAGCAACGCGACCAGCGGGCCGATCAGTGCCGCGTCCAGGCCACCGATTTCATCGTAGGGTTCAATGCCATCACAGGCCTCACCCACCCCCACGGCATAACCCAACAGCATGCGGCGCAG‑GCCGAAACGCCAActgttctgt                                                                                     >  6:37502/1‑149 (MQ=255)                                                                   cgcATCCAGCAACGCGACCAGCGGGCCGATCAGTGCCGCGTCCAGGCCACCGATTTCATCGTAGGGATCAATGCCATCACAGGCCTCACCCACCCCCACGGCATAACCCAACAGCATGCGGCGCAG‑GCCGAAACGCCAActgttctgtt                                                                                    <  1:88491/149‑1 (MQ=255)                                                                      aTCCAGCAACGCGACCAGCGGGCCGATCAGTGCCGAGTCCAGGCCACCGATTTCATCGTAGGGTTCAATGCCATCACAGGCCTCACCCACCCCCACGGCATAACCCAACAGCATGCGGCGCAG‑GCCGAAACGCCAACtgttctgttc                                                                                   >  1:50765/1‑147 (MQ=255)                                                                         cAGCAACGCGACCAGCGGGCCGATCAGTGCCGCGTCCAGGCCACCCATTTCATCGTAGGGTTCAATGCCATCACACGCCTCACCCACCCCCACGGCATAACCCAACACCAAGCGGCGCAG‑CACGAAACGCCAACTGTTCTGTTCGAgac                                                                              >  4:200261/1‑147 (MQ=255)                                                                          aGCAACGCGACCAGCGGGCCGATCAGTGCCGCGTCCAGGCCACCGATTTCATCGTAGGGTTCAATGCCATCACAGGCCTCACCCACCCCCACGGCATAACCCAACAGCATGCGGCGCAG‑GCCGAAACGCCAACTGTTCTGTTCGAGccc                                                                             >  6:254779/1‑149 (MQ=255)                                                                              aCGCGACCAGCGGGCCGATCAGTGCCGCGTCCAGGCCACCGATTTCATCGTAGGGTTCAATGCCATCACAGGCCTCACCCACCCCCACGGCATAACCCAACAGCATGCGGCGCAG‑GCCGAAACGCCAACTGTTCTGTTCGAGCCcggcg                                                                         <  2:84090/149‑1 (MQ=255)                                                                               cgcgACCAGCGGGCCGATCAGTGCCGCGTCCAGGCCACCGATTTCATCGTAGGGTTCAATGCCATCACAGTCCTCACCCACCCCCACGGCATAACCCAACAGCATGCGGCGCAG‑GCCGAAACGCCAActgtt                                                                                         >  2:324226/1‑132 (MQ=255)                                                                               cgcgACCAGCGGGCCGATCAGTGCCGCGTCCAGGCCACCGATTTCATCGTAGGGTTCAATGCCATCACAGTCCTCACCCACCCCCACGGCATAACCCAACAGCATGCGGCGCAG‑GCCGAAACGCCAActgtt                                                                                         >  2:329300/1‑132 (MQ=255)                                                                               cgcgACCAGCGGGCCGATCAGTGCCGCGTCCAGGCCACCGATTTCATCGTAGGGTTCAATGCCATCACAGTCCTCACCCACCCCCACGGCATAACCCAACAGCATGCGGCGCAG‑GCCGAAACGCCAActgtt                                                                                         <  1:329300/132‑1 (MQ=255)                                                                               cgcgACCAGCGGGCCGATCAGTGCCGCGTCCAGGCCACCGATTTCATCGTAGGGTTCAATGCCATCACAGTCCTCACCCACCCCCACGGCATAACCCAACAGCATGCGGCGCAG‑GCCGAAACGCCAActgtt                                                                                         <  1:324226/132‑1 (MQ=255)                                                                               cgcgACCAGCGGGCCGATCAGTGCCGCGTCCAGGCCACCGATTTCATCGTAGGGTTCAATGCCATCACAGGCCTCACCCACCCCCACGGCATAACCCAACAGCATGCGGCGCAG‑GCCGAAACGCCAACTGTTCTGTTCGAGCCcggcgg                                                                        >  1:74056/1‑149 (MQ=255)                                                                                  gACCAGCGGGCCGATCAGTGCCGCGTCCAGGCCACCGATTTCATCGTAGGGTTCAATGCCATCACAGGCCTCACCCACCCCCACGGCATAACCCAACAGCATg                                                                                                                    <  6:314249/103‑1 (MQ=255)                                                                                  gACCAGCGGGCCGATCAGTGCCGCGTCCAGGCCACCGATTTCATCGTAGGGTTCAATGCCATCACAGGCCTCACCCACCCCCACGGCATAACCCAACAGCATg                                                                                                                    >  5:314249/1‑103 (MQ=255)                                                                                  gACCAGCGGGCCGATCAGTGCCGCGTCCAGGCCACCGATTTCATCGTAGGGTTCAATGCCATCACAGGCCTCACCCACCCCCACGGCATAACCCAACAGCATGCGGCGCAG‑GCCGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGGCa                                                                     <  7:26218/149‑1 (MQ=255)                                                                                     cAGCGGGCCGATCAGTGCCGCGTCCAGGCCACCGATTTCATGGTAGGGTTCAATGCCATCACAGGCCTCACCCACCCCCACGGCATAACCCAACAGCATGCGGCGCAG‑GCCGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGc                                                                  <  5:295069/149‑1 (MQ=255)                                                                                        cGGGCCGATCAGTGCCGCGTCCAGGCCACCGATTTCATCGTAGGGTTCAATGCCATCACAGGCCTCACCCACCCCCACGGCATAACCCACCAGCATGCGGCGCAGTGCCGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGccc                                                                <  1:67091/149‑1 (MQ=255)                                                                                          ggCCGATCAGTGCCGCGTCCAGGCCACCGATTTCATCGTAGGGTTCAATGCCATCACAGGCCTCACCCACCCCCACGGCATAACCCAACAGCATGCGGCGCAG‑GCCGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAgg                                                             <  4:61379/149‑1 (MQ=255)                                                                                          ggCCGATCAGTGCCGCGTCCAGGCCACCGATTTCATCGTAGGGTTCAATGCCATCACAGGCCTCACCCACCCCCACGGCATAACCCAACAGCATGCGGCGCAG‑GCCGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAgg                                                             <  2:74056/149‑1 (MQ=255)                                                                                           gCCGATCAGTGCCGCGTCCAGGCCACCGATTTCATCGTAGGGTTCAATGCCATCACAGGCCTCACCCACCCCCACGGCATAACCCAACAGCATGCGGCGCAG‑GCCGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGc                                                            <  2:340950/149‑1 (MQ=255)                                                                                                tCAGTGCCGCGTCCAGGCCACCGATTTCATCGTAGGGTTCAATGCCATCACAGGCCTCACCCACCCCCACGGCATAACCCAACAGCATGCGGCGCAG‑GCCGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGcc                                                       <  6:245501/149‑1 (MQ=255)                                                                                                              aGGCCACCGATTTCATCGTAGGGTTCAATGCCATCACAGGCCTCACCCACCCCCACGGCATAACCCAACAGCATGCGGCGCAG‑GCCGAAAAGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGCGCGGCGt                                          >  5:194379/1‑148 (MQ=255)                                                                                                                gCCACCGATTTCATCGTAGGGTTCAATGCCATCACAGGCCTCACCCACCCCCACGGCATAACCCAACAGCATGCGGCGCAG‑GCCGAAAAGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGCGCGGCGTCaa                                       <  6:194379/149‑1 (MQ=255)                                                                                                                   aCCGATTTCATCGTAGGGTTCAATGCCATCACAGGCCTCACCCACCCCCACGGCATAACCCAACAGCATGCGGCGCAG‑GCCGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGCGCGGCGTCAAGGc                                    >  3:52255/1‑149 (MQ=255)                                                                                                                    ccGATTTCATCGTAGGGTTCAATGCCATCACAGGCCTCACCCACCCCCACGGCATAACCCAACAGCATGCGGCGCAG‑GCCGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGCGCGGCGTCAAGGcc                                   >  7:32819/1‑149 (MQ=255)                                                                                                                     cGATTTCATCGTAGGGTTCAATGCCATCACAGGCCTCACCCACCCCCACGGCATAACCCAACAGCATGCGGCGCAG‑GCCGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGCGCGGCGTCAAGGccc                                  >  8:108360/1‑149 (MQ=255)                                                                                                                     cGATTTCATCGTAGGGTTCAATGCCATCACAGGCCTCACCCACCCCCACGGCATAACCCAACAGCATGCGGCGCAG‑GCCGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGCGCGGCGTCAAGGccc                                  <  4:52255/149‑1 (MQ=255)                                                                                                                        tttCATCGTAGGGTTCAATGCCATCACAGGCCTCACCCACCCCCACGGCATAACCCAACAGCATGCGGCGCAG‑GCCGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGCGCGGCGTCAAGGCCCCaa                               >  1:164142/1‑149 (MQ=255)                                                                                                                             tCGTAGGGTTCAATGCCATCACAGGCCTCACCCACCCCCTCGGCATAACCCAACAGCATGCGGCGCAG‑GCCGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGCGCGGCGTCAAGGCCCCAACGGAt                          <  3:200261/149‑1 (MQ=255)                                                                                                                               gTAGGGTTCAATGCCATCACAGGCCTCACCCACCCCCACGGCATAACCCAACAGCATGCGGCGCAG‑GCCGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGCGCGGCGTCAAGGCCCCAACGGATGc                        <  7:108360/149‑1 (MQ=255)                                                                                                                                tAGGGTTCAATGCCATCACAGGCCTCACCCACCCCCACGGCATAACCCAACAGCATGCGGCGCAG‑GCCGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGCGCGGCGTCAAGGCCCCAACGGATGc                        >  4:38430/1‑148 (MQ=255)                                                                                                                                tAGGGTTCAATGCCATCACAGGCCTCACCCACCCCCACGGCATAACCCAACAGCATGCGGCGCAG‑GCCGAAACGCCAAATGTTCTGATCGACCCCGGCCGGCAAGCCAACGCTCGCCCGCTGCGCGGCGTCAAGGCCCCAACGGAAGc                        >  3:134534/1‑148 (MQ=255)                                                                                                                                           gCCATCACAGGCCTCACCCACCCCCACGGCATAACCCAACAGCATGCGGCGCAG‑GCCGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGCGCGGCGTCAAGGCCCCAACGGATGCCAgcg                   >  7:66697/1‑142 (MQ=255)                                                                                                                                           gCCATCACAGGCCTCACCCACCCCCACGGCATAACCCAACAGCATGCGGCGCAG‑GCCGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGCGCGGCGTCAAGGCCCCAACGGATGCCAgcg                   <  8:66697/142‑1 (MQ=255)                                                                                                                                           gCCATCACAGGCCTCACCCACCCCCACGGCATAACCCAACAGCATGCGGCGCAG‑GCCGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGCGCGGCGTCAAGGCCCCAACGGATGCCAGCGCCCTCGa            >  3:188519/1‑149 (MQ=255)                                                                                                                                              aTCACAGGCCTCCCCCACCCCCACGGCATAACCCAACAGCATGCGGCGCAG‑GCCGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGCGCGGCGTCAAGGCCCCAACGGATGCCAGCGCCCTCGATc          >  3:199760/1‑148 (MQ=255)                                                                                                                                               tCACAGGCCTCACCCACCCCCACGGCATAACCCAACAGCATGCGGCGCAG‑GCCGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGCGCGGCGTCAAGGCCCCAACGGATGCCAGCGCCCTCGATCCa        <  5:37502/149‑1 (MQ=255)                                                                                                                                                  caGGCCTCACCCACCCCCACGGCATAACCCAACAGCATGCGGCGCAG‑GCCGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGCGCGGCGTCAAGGCCCCAACGGATGCCAGCGCCCTCGATCCAGCg     >  5:7476/1‑149 (MQ=255)                                                                                                                                                     gCCTCACCCACCCCCACGGCATAACCCAACAGCATGCGGCGCAG‑GCCGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGCGCGGCGTCAAGGCCCCAACGGATGCCAGCGCCCTCGATCCAGCGGTg  <  6:7476/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                      ATGCAGTCGCTCGCCCCACTGGCTGACGGTGGCAGGCTCGGACAAGGCCTGGCAGGCGACTTCGAGCGCATCCAGCAACGCGACCAGCGGGCCGATCAGTGCCGCGTCCAGGCCACCGATTTCATCGTAGGGTTCAATGCCATCACAGGCCTCACCCACCCCCACGGCATAACCCAACAGCATGCGGCGCAG‑GCCGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGCGCGGCGTCAAGGCCCCAACGGATGCCAGCGCCCTCGATCCAGCGGTG  >  NC_002947/5305428‑5305724

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