BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	618,358	A→G	5.8%	S71P (TCG→CCG)	PP_0532 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	618,358	0	A	G	5.8%	303.5 / 5.9	103	S71P (TCG→CCG)	PP_0532	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (46/51); new base G (4/2); total (50/53)
Fisher's exact test for biased strand distribution p-value = 4.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.12e-01

Only 100 of 105 total aligned reads displayed.

GCCGGCCCTGTCGATGCGCAGTTCGGTCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTA‑GTGAC‑GCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCTGTAACACCGACGTGTCCATCTTCTGGTCATGCTCGGGCTGG  >  NC_002947/618212‑618506

gccggccCTGTCGATGCGCAGTTCGGTCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATg                                                                                                                                                      <  8:42956/149‑1 (MQ=255)

gccggccCTGTCGATGCGCAGTTCGGGCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATg                                                                                                                                                      <  6:382620/149‑1 (MQ=255)

    gccCTGTCGATGCGCAGTTCGGTCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCATCGGGGTACCCGGTGCTgt                                                                                                                                                  <  6:63799/149‑2 (MQ=255)

          tCGATGCGCAGTTCGGTCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCt                                                                                                                                              <  2:47671/147‑1 (MQ=255)

          tCGATGCGCAGTTCGGTCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTg                                                                                                                                             <  2:164425/148‑1 (MQ=255)

            gATGCGCAGTTCGGTCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGa                                                                                                                                           <  1:93745/148‑1 (MQ=255)

            gATGCGCAGTTCGGTCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGAc                                                                                                                                          >  6:393055/1‑149 (MQ=255)

            gATGCGCAGTTCGGTCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGAc                                                                                                                                          >  5:385499/1‑149 (MQ=255)

             aTGCGCAGTTCGGTCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGTGTACCCCGATGCTGGCGCTGGAc                                                                                                                                          >  3:233600/1‑148 (MQ=255)

                 gcAGTTCGGTCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTc                                                                                                                                     >  7:81364/1‑149 (MQ=255)

                 gcAGTTCGGTCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTc                                                                                                                                     <  6:293621/149‑1 (MQ=255)

                    gTTCGGTCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGc                                                                                                                                  >  3:167868/1‑149 (MQ=255)

                     ttCGGTCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCt                                                                                                                                 >  1:276732/1‑149 (MQ=255)

                     ttCGGTCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCt                                                                                                                                 >  4:24164/1‑149 (MQ=255)

                        ggTCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGc                                                                                                                              >  6:183460/1‑149 (MQ=255)

                         gTCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAAGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGc                                                                                                                              <  1:348048/149‑1 (MQ=255)

                          tCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCa                                                                                                                            <  6:27559/149‑1 (MQ=255)

                            cGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAAc                                                                                                                          >  6:70660/1‑149 (MQ=255)

                             ggTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAAcc                                                                                                                         >  6:242848/1‑149 (MQ=255)

                              gTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGG‑ATCATGGGAATCGAAGGGACCACACCCCCGGGTACCCCCATGAGGGCGATGGACGATTCGCCTTGCCAAcc                                                                                                                         >  1:142348/1‑149 (MQ=255)

                              gTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAAccc                                                                                                                        <  2:320070/149‑1 (MQ=255)

                                   ccACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAgg                                                                                                                   >  8:157809/1‑149 (MQ=255)

                                   ccACCGCCACCAGCAGGAGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAgg                                                                                                                   <  4:397712/149‑1 (MQ=255)

                                     aCCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCg                                                                                                                 >  8:402601/1‑149 (MQ=255)

                                        gCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGTGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTa                                                                                                              <  4:233600/149‑1 (MQ=255)

                                        gCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTa                                                                                                              <  3:24164/149‑1 (MQ=255)

                                        gCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTa                                                                                                              <  7:225663/149‑1 (MQ=255)

                                         ccaccaGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTA‑g                                                                                                            <  2:338822/149‑1 (MQ=255)

                                         ccaccaGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTA‑g                                                                                                            >  5:326569/1‑149 (MQ=255)

                                          caccaGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTA‑Gt                                                                                                           <  5:393055/149‑1 (MQ=255)

                                             cagcagGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTA‑GTGAc                                                                                                        <  7:402601/149‑1 (MQ=255)

                                             cagcagGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTA‑GTGAc                                                                                                        >  8:125830/1‑149 (MQ=255)

                                             cagcagGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCCC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTA‑GTGAc                                                                                                        >  2:216416/1‑149 (MQ=255)

                                                cagGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGGACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTA‑GTGAC‑g                                                                                                      >  3:269550/1‑147 (MQ=255)

                                                 agGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTA‑GTGAC‑GCAg                                                                                                   <  6:298758/149‑1 (MQ=255)

                                                   gggCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTA‑GTGAC‑GCAGGc                                                                                                 >  3:292022/1‑149 (MQ=255)

                                                        ccTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTA‑GTGAC‑GCAGGCCATCg                                                                                            <  5:183460/149‑1 (MQ=255)

                                                        ccTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTA‑GTGAC‑GCAGGCCATCg                                                                                            <  3:130011/149‑1 (MQ=255)

                                                          tGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCt                                                                                                                                              <  3:156951/99‑1 (MQ=255)

                                                          tGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCt                                                                                                                                              >  4:156951/1‑99 (MQ=255)

                                                              ccACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTA‑g                                                                                                            >  8:374048/1‑128 (MQ=255)

                                                              ccACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTA‑g                                                                                                            <  7:374048/128‑1 (MQ=255)

                                                              ccACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTA‑GTGAC‑GCAGGCCATCGGCCGAg                                                                                      <  4:269550/149‑1 (MQ=255)

                                                              ccACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTA‑GTGAC‑GCAGGCCATCGGCCGAg                                                                                      >  2:13355/1‑149 (MQ=255)

                                                              ccACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTA‑GTGAC‑GCAGGCCATCGGCCGAg                                                                                      >  6:338841/1‑149 (MQ=255)

                                                                 cGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTA‑GTGAC‑GCAGGCCATCGGCCGAGTcc                                                                                   <  8:170627/149‑1 (MQ=255)

                                                                    cGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTACGTGAC‑GCAGGCCATCGGCCGAGTCCAg                                                                                 <  1:13355/149‑1 (MQ=255)

                                                                      aGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTA‑GTGAC‑GCAGGCCATCGGCCGAGTCCAGATc                                                                              <  4:351316/149‑1 (MQ=255)

                                                                        ccTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTA‑GTGAC‑GCAGGCCATCGGCCGAGTCCAGATCaa                                                                            <  7:157809/149‑1 (MQ=255)

                                                                         cTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTA‑GTGAC‑GCAGGCCATCGGCCGAGTCCAGATCAAg                                                                           >  4:185490/1‑149 (MQ=255)

                                                                           gTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTA‑GTGAC‑GCAGGCCATCGGCCGAGTCCAGATCAAGCt                                                                         >  4:40258/1‑149 (MQ=255)

                                                                              cAGTTGCTGCTGATCCAGCGTCCCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTA‑GTGACAGCAGGCCATCGGCCGAGTCCAGATCAAGCTgg                                                                       <  1:52547/149‑1 (MQ=255)

                                                                              cAGTTGCTGCTGATCCAGCGTCCCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTA‑GTGAC‑GCAGGCCATCGGCCGAGTCCAGATCAAGCTggg                                                                      <  4:167868/149‑1 (MQ=255)

                                                                              cAGTTGCTGCTGATCCAGCGTCCCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTA‑GTGAC‑GCAGGCCATCGGCCGAGTCCAGATCAAGCTgg                                                                       <  4:292022/148‑1 (MQ=255)

                                                                              cAGTTGCTGCTGATCCAGCGTCCCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTA‑GTGAC‑GCAGGCCATCGGCCGAGTCCAGATCAAGCTgg                                                                       <  6:161985/148‑1 (MQ=255)

                                                                               aGTTGCTGCTGATCCAGCGTCCCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTA‑GTGAC‑GCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGt                                                                     >  1:262302/1‑149 (MQ=255)

                                                                                    ctgctgATCCAGCGTCCCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTA‑GTGAC‑GCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAg                                                                 <  7:125830/148‑1 (MQ=255)

                                                                                    ctgctgATCCAGCGTCCCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTA‑GTGAC‑GCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAg                                                                 <  5:392758/148‑1 (MQ=255)

                                                                                    ctgctgATCCAGCGTCCCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTA‑GTGAC‑GCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAg                                                                 <  2:142348/148‑1 (MQ=255)

                                                                                      gctgATCCAGCGTCCCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCgg                                                                                                                                   <  2:313116/82‑1 (MQ=255)

                                                                                      gctgATCCAGCGTCCCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCgg                                                                                                                                   >  1:313116/1‑82 (MQ=255)

                                                                                      gctgATCCAGCGTCCCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTA‑GTGAC‑GCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGa                                                              >  4:112877/1‑149 (MQ=255)

                                                                                          aTCCAGCGTCCCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTA‑GTGAC‑GCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAgg                                                          <  3:345730/149‑1 (MQ=255)

                                                                                          aTCCAGCGTCCCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTA‑GTGAC‑GCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAg                                                           <  3:40258/148‑1 (MQ=255)

                                                                                           tCCAGCGTCCCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTA‑GTGAC‑GCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGc                                                         >  2:359444/1‑149 (MQ=255)

                                                                                           tCCAGCGTCCCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGATGGCGCTGGACGATTCGGCTTGCCAACCCACCGGCGGTA‑GTGAC‑GCAGGCCTTCGGCCGAGTCCAGAGCCAGCTGGGTGAAGAGATCGGc                                                         >  4:40520/1‑149 (MQ=255)

                                                                                           tCCAGCGTACCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTA‑GTGAC‑GCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGc                                                         >  8:138628/1‑149 (MQ=255)

                                                                                            ccAGCGTCCCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCTGCTAGCCAACCCACAGGCGGTA‑GTGAC‑GCAGGCCATCGGCCGAGTCCACATCAAGCTGGGTGAAGTGATAGGc                                                         >  2:233038/1‑148 (MQ=255)

                                                                                             cAGCGTCCCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTA‑GTGAC‑GCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGcc                                                        <  8:353887/148‑1 (MQ=255)

                                                                                             cAGCGTCCCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTA‑GTGAC‑GCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGc                                                         >  1:21524/1‑147 (MQ=255)

                                                                                             cAGCGTCACCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTA‑GTGAC‑GCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCa                                                       <  5:352719/149‑1 (MQ=255)

                                                                                               gCGTCCCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTA‑GTGAC‑GCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAgg                                                     >  7:330002/1‑149 (MQ=255)

                                                                                                 gTCCCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTA‑GTGAC‑GCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTc                                                   <  2:276732/149‑1 (MQ=255)

                                                                                                     ccAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTA‑GTGAC‑GCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTgg                                               >  3:298261/1‑149 (MQ=255)

                                                                                                     ccAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTA‑GTGAC‑GCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTgg                                               >  5:345712/1‑149 (MQ=255)

                                                                                                         cgc‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTA‑GTGAC‑GCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCt                                           <  3:185490/149‑1 (MQ=255)

                                                                                                          gc‑agcaTTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTA‑GTGAC‑GCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCTg                                          >  2:35294/1‑149 (MQ=255)

                                                                                                                  tCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTA‑GTGAC‑GCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCTGTAACAcc                                   >  6:221865/1‑149 (MQ=255)

                                                                                                                    ccATCG‑AGCATCCACAGCACCGGGTACCCCGATGGTGGCGCTGGACGATTCGGGTTGCCAACCCACAGGCGGTA‑GTGAC‑GCAGGCCATCGGGCGAGTCCAGATCAAGCTGGGGGAAGTGATAGGCCAGGTCCTGGCGCTGGAACACCGa                                 >  2:259683/1‑149 (MQ=255)

                                                                                                                    ccATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTA‑GTGAC‑GCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCt                                           >  4:118658/1‑139 (MQ=255)

                                                                                                                    ccATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTA‑GTGAC‑GCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCt                                           <  3:118658/139‑1 (MQ=255)

                                                                                                                    ccATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTA‑GTGAC‑GCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCTGTAACACCGa                                 >  4:118244/1‑149 (MQ=255)

                                                                                                                         g‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTA‑GTGAC‑GCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCTGTAACACCGACgtgt                            >  3:277814/1‑149 (MQ=255)

                                                                                                                         g‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTA‑GTGAC‑GCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCTGTAACACCGACgtg                             <  1:197011/148‑1 (MQ=255)

                                                                                                                               tCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTA‑GTGAC‑GCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCTGTAACACCGACGTGTCCAtc                       <  7:138628/149‑1 (MQ=255)

                                                                                                                               tCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTA‑GTGAC‑GCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCTGTAACACCGACGTGTCCAtc                       <  2:215940/149‑1 (MQ=255)

                                                                                                                                  ccaGCATCGGGGTACCCGGTGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTA‑GTGAC‑GCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCTGTAACACCGACGTGTCCAtcttc                    >  4:144189/2‑149 (MQ=255)

                                                                                                                                  ccaGCATCGGGGTACCCGGTGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTA‑GTGAC‑GCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCTGTAACACCGACGTGTCCAtcttc                    >  3:413089/2‑149 (MQ=255)

                                                                                                                                  acaGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTA‑GTGAC‑GCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCTGTAACACCGACGTGTCCAtcttc                    >  5:383590/1‑149 (MQ=255)

                                                                                                                                  acaGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTA‑GTGAC‑GCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCTGTAACACCGACGTGTCCAtcttc                    <  6:383590/149‑1 (MQ=255)

                                                                                                                                   caGCATCGGGGTACCCGGTGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTA‑GTGAC‑GCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCTGTAACACCGACGTGTCCAtcttct                   >  5:294259/1‑149 (MQ=255)

                                                                                                                                   caGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTA‑GTGAC‑GCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCTGTAACACCGACGTGTCCAtcttct                   <  5:15039/149‑1 (MQ=255)

                                                                                                                                   agcatcggggTA‑CCCGGTGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTA‑GTGAC‑GCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCTGTAACACCGACGTGTCCATCTTCTg                  <  3:144189/142‑1 (MQ=255)

                                                                                                                                    gcatcggggTA‑CCCGGTGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTA‑GTGAC‑GCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCTGTAACACCGACGTGTCCATCTTCTgg                 >  5:108704/7‑149 (MQ=255)

                                                                                                                                    aGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTA‑GTGAC‑GCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCTGTAACACCGACGTGTCCATCTTCTg                  <  2:347871/149‑1 (MQ=255)

                                                                                                                                       aCCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTA‑GTGAC‑GCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCTGTAACACCGACGTGTCCATCTTCTGGt                >  3:311776/1‑148 (MQ=255)

                                                                                                                                       aCCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTA‑GTGAC‑GCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCTGTAACACCGACGTGTCCATCTTCTGGTc               <  6:345712/149‑1 (MQ=255)

                                                                                                                                           ggTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTA‑GTGAC‑GCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCTGTAACACCGACGTGTCCATCTTCTGGTCATGc           <  4:10544/149‑1 (MQ=255)

                                                                                                                                                cccGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGAA‑GTGAC‑GCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTCATAGGCCAGGTCCTGGCGCTGTAACCCCGACGTGTCCATCTTCTGGTCATGCTCggg      >  8:300323/1‑149 (MQ=255)

                                                                                                                                                    aTGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTA‑GTGAC‑GCAGGCCATCGGCCGAGTCGAGATCAAGCTGGGTGAAGTGATAGGCGAGGTCCTGGCGCTGTAACACCGACGTGTCCATCTTCTGGTCATGCTCGGGCTgg  <  1:259683/149‑1 (MQ=255)

GCCGGCCCTGTCGATGCGCAGTTCGGTCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGC‑AGCATTCCCATCG‑AGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTA‑GTGAC‑GCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCTGTAACACCGACGTGTCCATCTTCTGGTCATGCTCGGGCTGG  >  NC_002947/618212‑618506

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