Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,786,023	T→A	8.1%	intergenic (+39/‑147)	rpsB → / → tsf	30S ribosomal protein S2/elongation factor Ts

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,786,023	0	T	A	8.1%	200.6 / 6.0	71	intergenic (+39/‑147)	rpsB/tsf	30S ribosomal protein S2/elongation factor Ts
Reads supporting (aligned to +/- strand):  ref base T (33/32);  new base A (2/4);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TACATGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAA‑AGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAA‑CGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCAAGAGGATTT  >  NC_002947/1785883‑1786165                                                                                                                                              |                                                                                                                                                tACATGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAA‑AGGGGGCTTTGCCCCCt                                                                                                                                         >  7:299587/1‑149 (MQ=255)  acctGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAA‑AGGGGGCTTTGCCCCCtt                                                                                                                                        <  7:210808/146‑1 (MQ=255)  aCATGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAA‑AGGGGGCTTTGcccccgt                                                                                                                                        >  2:242728/1‑147 (MQ=255)   cATGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAA‑AGGGGGCTTTGCCCCCtt                                                                                                                                        <  4:91054/148‑1 (MQ=255)     tGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAA‑AGGGGGCTTTGCCCCCttttt                                                                                                                                     <  7:270954/149‑1 (MQ=255)     tGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAA‑AGGGGGCTTTGCCCCCttttt                                                                                                                                     <  4:195496/149‑1 (MQ=255)     tGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAA‑AGGGGGCAAAGCCCCCttttt                                                                                                                                     <  8:169585/149‑1 (MQ=255)      gACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAA‑AGGGGGCTTTGCCCCCTTTTTg                                                                                                                                    >  8:76567/1‑149 (MQ=255)        cTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCAGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAA‑AGGGGGCTTTGCCCCCTTTTTGc                                                                                                                                   >  5:306162/1‑148 (MQ=255)         ttCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAA‑AGGGGGCAAAGCCCCCTTTTTGc                                                                                                                                   >  3:218074/1‑147 (MQ=255)         ttCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAA‑AGGGGGCAAAGCCCCCTTTTTGc                                                                                                                                   <  4:218074/147‑1 (MQ=255)           cGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAA‑AGGGGGCAAAGCCCCCTTTTTGcgt                                                                                                                                 <  8:183049/147‑3 (MQ=255)             aTGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAA‑AGGGGGCTTTGCCCCCTTTTTGCCACCtt                                                                                                                             <  7:76567/149‑1 (MQ=255)                 cTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGCCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAA‑AGGGGGCTTTGCCCCCTTTTTGCCACCTTGaaa                                                                                                                         >  5:135539/1‑149 (MQ=255)                      gCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAA‑AGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCctg                                                                                                                    >  8:243823/1‑149 (MQ=255)                      gCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAA‑AGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCctg                                                                                                                    >  5:6446/1‑149 (MQ=255)                      gCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAA‑AGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCctg                                                                                                                    <  1:267950/149‑1 (MQ=255)                                gcgGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAA‑AGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCAt                                                                                                          >  2:63703/1‑149 (MQ=255)                                 cgGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAA‑AGGGGGCTTTGCCCCCTTTTTGCCACCt                                                                                                                              <  5:321910/128‑1 (MQ=255)                                 cgGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAA‑AGGGGGCTTTGCCCCCTTTTTGCCACCt                                                                                                                              >  6:321910/1‑128 (MQ=255)                                 cgGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAA‑AGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATc                                                                                                         <  5:175376/149‑1 (MQ=255)                                        aacaacGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAA‑AGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCg                                                                                                  >  5:88266/1‑149 (MQ=255)                                          caacGTTGCCGCCGGAACCGAAGTTTACGCTCAAGACGCGGCTGCACCTGCTGCTGAGTAAAAGCCGCTAGCGTCAAATTGGCACGCAAAA‑AGGGGCCATTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGca                                                                                                <  1:81963/149‑1 (MQ=255)                                               ttGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAA‑AGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATg                                                                                           <  6:224953/149‑1 (MQ=255)                                               ttGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTAAGTAATAGGCGCGAGCGGATACTTGGCCCGCAAAA‑AGGGGGGTTTGCCCCCTTTTTGCCCCCATGAAATACGGCGGTAATCATAGGCCCCGcatcatt                                                                                           >  1:247660/1‑148 (MQ=255)                                                         cACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAA‑AGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCctgctg                                                                                  >  4:72569/1‑148 (MQ=255)                                                          aCCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAA‑AGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCtga                                                                                 >  3:199642/1‑148 (MQ=255)                                                             gAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAA‑AGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCtgaaat                                                                              >  2:151203/1‑148 (MQ=255)                                                                 tttACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAA‑AGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCGGCATCATGGGCCTGCTGAAATGaaaa                                                                         <  6:6446/149‑1 (MQ=255)                                                                      gCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAA‑AGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATc                                                                                                                       <  2:57697/98‑1 (MQ=255)                                                                      gCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAA‑AGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATc                                                                                                                       >  1:57697/1‑98 (MQ=255)                                                                        tGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAA‑AGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGa                                                                  >  7:283682/1‑149 (MQ=255)                                                                        tGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAA‑AGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGa                                                                  >  2:54861/1‑149 (MQ=255)                                                                                       aCCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAA‑AGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCa                                                            <  4:187967/140‑1 (MQ=255)                                                                                       aCCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAA‑AGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCa                                                            >  3:187967/1‑140 (MQ=255)                                                                                       aCCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAA‑AGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAA‑c                                                  >  5:301264/1‑149 (MQ=255)                                                                                         ctgctgctgAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAGAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCACCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGaag                                                   >  1:68732/1‑147 (MQ=255)                                                                                         ctgctgctgAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAA‑AGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAAtt                                                       <  8:201846/143‑1 (MQ=255)                                                                                         ctgctgctgAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAA‑AGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAAtt                                                       >  7:201846/1‑143 (MQ=255)                                                                                         ctgctgctgAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAA‑AGGGGGCTTTGCCCCCTTTTTGACACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCATGATGAAATGAAAACAGCGGATTTGCAGAATTTAA‑CGc                                                >  4:4916/1‑149 (MQ=255)                                                                                         ctgctgctgAGTAATACGCTCTAGCGTCTACTTGGCACGCAAAA‑AGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCATGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTg                                                      <  2:68732/144‑1 (MQ=255)                                                                                           gctgctgAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAA‑AGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAA‑CGccc                                              >  2:235977/1‑149 (MQ=255)                                                                                           gctgctgAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAA‑AGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAA‑CGccc                                              >  2:240410/1‑149 (MQ=255)                                                                                            ctgctgAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAA‑AGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAA‑CGCCCg                                             >  6:142626/1‑149 (MQ=255)                                                                                               ctgAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAA‑AGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAA‑CGCCCGTGa                                          <  1:54861/149‑1 (MQ=255)                                                                                                 gAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAA‑AGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAA‑CGCCCGTGAtt                                        >  5:63385/1‑149 (MQ=255)                                                                                                     aaTAGGCGGTAGCGTCTACTTGGCACGCAAAA‑AGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAA‑CGCCCGTGATTGAc                                     <  3:4916/148‑1 (MQ=255)                                                                                                      aTAGGCGCTAGCGTCTACTTGGCACGCAAAA‑AGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAA‑CGCCCGTGATTGACgg                                   <  3:72569/149‑1 (MQ=255)                                                                                                      aTAGGCGCTAGCGTCTACTTGGCACGCAAAA‑AGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAA‑CGCCCGTGATTGACgg                                   <  4:199642/149‑1 (MQ=255)                                                                                                         ggCGCTAGCGTCTACTTGGCACGCAAAA‑AGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAA‑CGCCCGTGATTGACGGGt                                 <  7:243823/148‑1 (MQ=255)                                                                                                          gcgcTAGCGTCTACTTGGCACGCAAAA‑AGGGGGCTTTGCCCCCTTTTTGACACCTTTAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAATACAGCGGATTTGCAGAATTGAA‑CGCCCGTGATTGACGGGTg                                <  6:306162/148‑1 (MQ=255)                                                                                                              tAGCGTCTACTTGGCACGCAAAA‑AGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAA‑CGCCCGTGATTGACGGGTGGAATg                           >  8:27252/1‑149 (MQ=255)                                                                                                                gCGTCTACTTGGCACGCAAAA‑AGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAA‑CGCCCGTGATTGACGGGTGGAATg                           <  8:283682/147‑1 (MQ=255)                                                                                                                gCGTCTACTTGGCACGCAAAA‑AGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAA‑CGCCCGTGATTGACGGGTGGAATg                           <  1:132909/147‑1 (MQ=255)                                                                                                                gCGTCTACTTGGCACGCAAAA‑AGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAA‑CGCCCGTGATTGACGGGTGGAATGGt                         <  8:299587/149‑1 (MQ=255)                                                                                                                             aCGCAAAA‑AGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAA‑CGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTc             >  8:87948/1‑148 (MQ=255)                                                                                                                             aCGCAAAA‑AGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAA‑CGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTc             >  2:302972/1‑148 (MQ=255)                                                                                                                              cGCAAAA‑AGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCTCCGCATCATGGGCCTGCTGAAATGAAAAGAGCGGATTTGCAGAATTGAA‑CGCCCGTGATTGACGGGTGGTATGGTTGAAAAACTTTCCa           <  6:135539/149‑1 (MQ=255)                                                                                                                                 aaaa‑aGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAA‑CGCCCGt                                            <  2:40507/113‑1 (MQ=255)                                                                                                                                 aaaa‑aGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAA‑CGCCCGt                                            >  1:40507/1‑113 (MQ=255)                                                                                                                                 aaaa‑aGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAA‑CGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCaa          >  7:219712/1‑147 (MQ=255)                                                                                                                                 aaaa‑aGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAA‑CGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCAaga        <  7:27252/149‑1 (MQ=255)                                                                                                                                 aaaa‑aGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAA‑CGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCAaga        <  2:247660/149‑1 (MQ=255)                                                                                                                                 aaaa‑aGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAA‑CGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCAaga        <  1:240410/149‑1 (MQ=255)                                                                                                                                 aaaa‑aGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAA‑CGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCAaga        <  1:235977/149‑1 (MQ=255)                                                                                                                                 aaaa‑aGGGGGCAAAGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAA‑CGCCCGTGATTGACGGGTGGAATGGTTGaaaaa                  <  7:294101/139‑1 (MQ=255)                                                                                                                                 aaaa‑aGGGGGCAAAGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAA‑CGCCCGTGATTGACGGGTGGAATGGTTGaaaaa                  >  8:294101/1‑139 (MQ=255)                                                                                                                                        ggggCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCGGAATTGAA‑CGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCAAGAGGAttt  >  4:262861/1‑149 (MQ=255)                                                                                                                                         gggCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAATCGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCAAGAGGAttt  <  1:242728/149‑1 (MQ=255)                                                                                                                                         gggCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAA‑CGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCa           <  2:312490/139‑1 (MQ=255)                                                                                                                                         gggCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAA‑CGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCa           >  1:312490/1‑139 (MQ=255)                                                                                                                                              |                                                                                                                                                TACATGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAA‑AGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAA‑CGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCAAGAGGATTT  >  NC_002947/1785883‑1786165

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