Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	5,674,751	+G	intergenic (+113/+27)	rhlE‑II → / ← yceI	ATP‑dependent RNA helicase/polyisoprene‑binding acidic stress response factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,674,750	1	.	G	98.6%	230.5 / ‑4.8	70	intergenic (+112/+28)	rhlE‑II/yceI	ATP‑dependent RNA helicase/polyisoprene‑binding acidic stress response factor
Reads supporting (aligned to +/- strand):  ref base . (0/1);  major base G (32/36);  minor base A (0/1);  total (32/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 = 8.14e-01

CGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGAT‑A‑‑GGGCCGCAAA‑GC‑GCGCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGC  >  NC_002947/5674602‑5674890                                                                                                                                                          |                                                                                                                                             cgGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGAT‑A‑‑GGGCCGCAAA‑gc                                                                                                                                               <  6:363571/149‑1 (MQ=255)                  tgctgACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGAT‑A‑‑GGGCCGCAAA‑GCGGCCCCCTCGGtttt                                                                                                                                >  6:104684/1‑147 (MQ=255)                    ctgACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGAT‑A‑‑GGGCCGCAAA‑GCGGCCCCCTCGGTTTTCAGc                                                                                                                            <  1:56666/149‑1 (MQ=255)                       aCCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGAT‑A‑‑GGGCCGCAAA‑GCGGCCCCCTCGGTTTTCAGCTg                                                                                                                          <  3:145477/148‑1 (MQ=255)                        cccGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGAT‑A‑‑GGGCCGCAAA‑GCGGCCCCCTCGGTTTTCAGCTGGc                                                                                                                        <  6:240156/149‑1 (MQ=255)                         ccGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGAT‑A‑‑GGGCCGCAAA‑GCGGCCCCCTCGGTTTTCAGCTgg                                                                                                                         >  4:218075/1‑147 (MQ=255)                             gtgggCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGAT‑A‑‑GGGCCGCAAA‑GCGGCCCCCTCGGTTTTCAGCTGGCAGaaa                                                                                                                   <  2:56724/147‑1 (MQ=255)                              aGGGCTGTTAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGAT‑A‑‑GGGCCGCAAA‑GCGGCCCCCTCGGTTTTCAGCTGGCAGAAAt                                                                                                                  <  2:370105/149‑1 (MQ=255)                              aGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGAT‑A‑‑GGGCCGCAAA‑GCGGCCCCCTCGGTTTTCAGCTGGCAGAAAt                                                                                                                  <  4:243584/149‑1 (MQ=255)                                         aaGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGAT‑A‑‑GGGCCGCAAA‑GCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCtt                                                                                                       <  5:184004/149‑1 (MQ=255)                                          aGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCTGCTCCTGCAGAGAATTTCTTGATCTCTACAGGTGCGGCCTTGTGTCGTGAT‑A‑‑GGGCCGCAAA‑GCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTc                                                                                                      <  6:254380/149‑1 (MQ=38)                                            tACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGAT‑A‑‑GGGCCGCAAA‑GCGGCCCCCGCGGGTTTCAGCTGGCAGAAATTAGTTCTGCTTca                                                                                                     >  8:86556/1‑148 (MQ=37)                                            tACCGAGGGGGCTGCGTCGCAGCCCTATCCCGACACCAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGAT‑A‑‑GGGCCGCAAA‑GCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTcac                                                                                                    <  8:38320/149‑1 (MQ=38)                                                  gggggCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGAT‑A‑‑GGGCCGCAAA‑GCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTc                                                                                              <  6:214576/149‑1 (MQ=255)                                                   ggggCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGAT‑A‑‑GGGCCGCAAA‑GCGGCCCCCTCGGTTTTCAGCt                                                                                                                           <  7:148732/119‑1 (MQ=255)                                                   ggggCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGAT‑A‑‑GGGCCGCAAA‑GCGGCCCCCTCGGTTTTCAGCt                                                                                                                           >  8:148732/1‑119 (MQ=255)                                                   ggggCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGAT‑A‑‑GGGCCGCAAA‑GCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCt                                                                                                >  2:371051/1‑146 (MQ=255)                                                     ggCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGAT‑A‑‑GGGCCGCAAA‑GCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGa                                                                                            >  6:29993/1‑148 (MQ=255)                                                      gctgctTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGAT‑A‑‑GGGCCGCAAA‑GCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGaca                                                                                          <  3:8844/149‑1 (MQ=255)                                                          ctTCGCAGCCCTATCACGACACAAGACCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGAT‑A‑‑GGGCCGCAAA‑GCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGa                                                                                      >  8:31101/1‑149 (MQ=255)                                                           ttCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGAT‑A‑‑GGGCCGCAAA‑GCGGCCCCCTCGGTTTTGAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGa                                                                                      >  7:140956/1‑148 (MQ=255)                                                              gCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGAT‑A‑‑GGGCCGCAAA‑GCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATAtc                                                                                  <  2:43374/149‑1 (MQ=255)                                                                  cccTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGAT‑A‑‑GGGCCGCAAA‑GCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATAtccagg                                                                              >  6:206185/1‑149 (MQ=255)                                                                       tCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGAT‑A‑‑GGGCCGCAAA‑GCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATAtccaggtccag                                                                         <  1:248943/149‑1 (MQ=255)                                                                       tCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGAT‑A‑‑GGGCCGCAAA‑GCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTAGACACTGATAtccaggtcaag                                                                         >  5:171802/1‑149 (MQ=38)                                                                         aCGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGAT‑A‑‑GGGCCGCAAA‑GCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATAtccaggtccagg                                                                        <  3:218895/148‑1 (MQ=255)                                                                         aCGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGAT‑A‑‑GGGCCGCAAA‑GCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATAtccaggtccagg                                                                        <  1:371051/148‑1 (MQ=255)                                                                         aCGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGAT‑A‑‑GGGCCGCAAA‑GCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAggg                                                                       >  3:74491/1‑149 (MQ=255)                                                                         aCGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGAT‑A‑‑GGGCCGCAAA‑GCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAggg                                                                       >  4:211637/1‑149 (MQ=255)                                                                          cGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGGCTTGTGTCGTGAT‑A‑‑GGGCCGCAAA‑GCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGt                                                                      <  5:86297/149‑1 (MQ=255)                                                                           gACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGAT‑A‑‑GGGCCGCAAA‑GCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTc                                                                     >  1:288076/1‑149 (MQ=255)                                                                                aaGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGAT‑A‑‑GGGCCGCAAA‑GCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGc                                                                  <  1:166720/147‑1 (MQ=255)                                                                                aaGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGAT‑A‑‑GGGCCGCAAA‑GCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGc                                                                  >  2:166720/1‑147 (MQ=255)                                                                                    ccGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGAT‑A‑‑GGGCCGCAAA‑GCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTa                                                            <  3:218075/149‑1 (MQ=255)                                                                                      gCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGAT‑A‑‑GGGCCGCAAA‑GCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGc                                                          <  4:74491/149‑1 (MQ=255)                                                                                           tGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGAT‑A‑‑GGGCCGCAAA‑GCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCg                                                     >  7:200926/1‑149 (MQ=255)                                                                                             cAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGAT‑A‑‑GGGCCGCAAA‑GCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCggg                                                   >  4:73875/1‑149 (MQ=255)                                                                                             cAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGAT‑A‑‑GGGCCGCAAA‑GCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCgg                                                    <  3:1834/148‑1 (MQ=255)                                                                                             cAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGAT‑A‑‑GGGCCGCAAA‑GCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCgg                                                    <  7:31101/148‑1 (MQ=255)                                                                                              agagAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGAT‑A‑‑GGGCCGCAAA‑GCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCTAAGTAGCACACGGGc                                                  >  2:307840/1‑149 (MQ=39)                                                                                              agagAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGAT‑A‑‑GGGCCGCAAA‑GCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGc                                                  >  8:16690/1‑149 (MQ=255)                                                                                                   atttatttATCTCTACAGCAGCGGGCTTGTGTGGTGATCATTGGGCCGCAAA‑GCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCt                                                <  1:312180/149‑1 (MQ=255)                                                                                                   aTTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGAT‑A‑‑GGGCCGCAAA‑GCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTg                                             <  7:16690/149‑1 (MQ=255)                                                                                                    tttCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGAT‑A‑‑GGGCCGCAAA‑GCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTg                                             <  2:308044/148‑1 (MQ=255)                                                                                                     aagaaattctcTGCAGGAGCGGCCTTGTGTCGTGAT‑A‑‑GGGCTGCGAA‑GCAGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGAt                                           <  2:325036/142‑1 (MQ=25)                                                                                                                  cAGGAGCGGCCTTGTGTCGTGAT‑A‑‑GGGCCGCAAA‑GCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTc                                  <  8:80041/145‑1 (MQ=255)                                                                                                                  cAGGAGCGGCCTTGTGTCGTGAT‑A‑‑GGGCCGCAAA‑GCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTc                                  >  7:80041/1‑145 (MQ=255)                                                                                                                    ggAGCGGCCTTGTGTCGTGAT‑A‑‑GGGCCGCAAA‑GCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCgttcag                            <  5:104684/149‑1 (MQ=255)                                                                                                                     gAGCGGCCTTGTGTCGTGAT‑A‑‑GGGCCGCAAA‑GCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCgttcagg                           >  5:137030/1‑149 (MQ=255)                                                                                                                     gAGCGGCCTTGTGTCGTGAT‑A‑‑GGGCCGCAAA‑GCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCgttcagg                           >  6:257331/1‑149 (MQ=255)                                                                                                                         ggCCTTGTGTCGTGAT‑A‑‑GGGCCGCAAA‑GCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCgttcaggttca                       <  5:126981/149‑1 (MQ=255)                                                                                                                           ccTTGTTTCGTGAT‑A‑‑GGGCCGCAAATGG‑GCCCCCTCGGTTTTCAGCTGGCATAAATTTCTTCTTCTTCACACCTTCGTCACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGc                     <  3:14969/149‑1 (MQ=255)                                                                                                                           ccTTGTGTCGTGAT‑A‑‑TGGCCGCAAA‑GCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGc                     <  8:140956/149‑1 (MQ=255)                                                                                                                               gtgtCGTGAT‑A‑‑GGGCCGCAAA‑GCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCgtgg                 >  3:181788/1‑149 (MQ=255)                                                                                                                               gtgtCGTGAT‑A‑‑GGGCCGCAAA‑GCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCgtg                  <  8:200926/148‑1 (MQ=255)                                                                                                                                  tCGTGAT‑A‑‑GGGCCGCAAA‑GCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCgtggtgg              >  1:71963/1‑149 (MQ=255)                                                                                                                                  tCGTGAT‑A‑‑GGGCCGCAAA‑GCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCgtggtgg              >  7:258177/1‑149 (MQ=255)                                                                                                                                    gTGAT‑A‑‑GGGCCGCAAA‑GCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTATCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCgtggtggtg            <  4:317359/149‑1 (MQ=255)                                                                                                                                    gTGAT‑A‑‑GGGCCGCAAA‑GCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCgtggtggtg            >  3:317359/1‑149 (MQ=255)                                                                                                                                     tGAT‑A‑‑GGGCCGCAAA‑GCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGtggtggtgg           >  3:87326/1‑149 (MQ=255)                                                                                                                                          a‑‑GGGCCGCAAA‑GCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGtt       >  5:302561/1‑149 (MQ=255)                                                                                                                                          a‑‑GGGCCGCAAA‑GCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGACCTTTGATGCCGAAGTCgttcaggtt                         <  2:290942/131‑1 (MQ=255)                                                                                                                                          a‑‑GGGCCGCAAA‑GCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGACCTTTGATGCCGAAGTCgttcaggtt                         >  1:290942/1‑131 (MQ=255)                                                                                                                                               gCCGCAAA‑GCGGCCCCCTCGGTTTTCAGCTGGCAGATATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGaa    <  3:225151/149‑1 (MQ=255)                                                                                                                                                ccGCAAA‑GCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCgtggtg               >  6:252459/1‑137 (MQ=255)                                                                                                                                                ccGCAAA‑GCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCgtggtg               <  5:252459/137‑1 (MQ=255)                                                                                                                                                ccGCAAA‑GCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGa     >  7:129839/1‑147 (MQ=255)                                                                                                                                                ccGCAAA‑GCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGa     <  8:129839/147‑1 (MQ=255)                                                                                                                                                 cGCAAA‑GCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATg                                          <  8:250941/109‑1 (MQ=255)                                                                                                                                                   cAAA‑GCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGc  >  4:94035/1‑147 (MQ=255)                                                                                                                                                   cAAA‑GCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGc  >  2:143345/1‑147 (MQ=255)                                                                                                                                                          |                                                                                                                                             CGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGAT‑A‑‑GGGCCGCAAA‑GC‑GCGCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGC  >  NC_002947/5674602‑5674890

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