BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,388,044	T→A	5.1%	D149V (GAC→GTC)	PP_1209 ←	cold‑shock protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,388,044	0	T	A	5.1%	295.2 / 4.4	99	D149V (GAC→GTC)	PP_1209	cold‑shock protein
Reads supporting (aligned to +/- strand): ref base T (45/49); new base A (3/2); total (48/51)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.64e-01

Only 100 of 102 total aligned reads displayed.

CAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCAT‑TTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGCCAGATGCA  >  NC_002947/1387897‑1388182

cAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGtc                                                                                                                                            >  7:16012/1‑149 (MQ=255)

cAACGGCGGTTTACTGCAACCACGTCTTCGCCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCGGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGtc                                                                                                                                            <  8:60551/149‑1 (MQ=255)

  aCGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGtctc                                                                                                                                          >  8:355818/1‑149 (MQ=255)

    ggcggTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCgg                                                                                                                                        <  5:227972/149‑1 (MQ=255)

    ggcggTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCgg                                                                                                                                        <  7:355818/149‑1 (MQ=255)

     gcggTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCggg                                                                                                                                       <  4:124171/149‑1 (MQ=255)

      cggTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGga                                                                                                                                      <  1:389216/149‑1 (MQ=255)

       ggTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGgag                                                                                                                                     <  8:124606/149‑1 (MQ=255)

        gTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCTCCCGAGACTCGGgcga                                                                                                                                    <  6:393827/149‑1 (MQ=255)

        gTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGgaga                                                                                                                                    <  5:101566/149‑1 (MQ=255)

         tttACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGAt                                                                                                                                   <  2:29827/149‑1 (MQ=255)

            aCTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGACATCGCCCGAGTCTCGGGAGATGa                                                                                                                                 <  1:159775/148‑1 (MQ=255)

                cAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAAcc                                                                                                                             >  7:42199/1‑148 (MQ=255)

                cAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAAc                                                                                                                              <  4:183095/147‑1 (MQ=255)

                   ccACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGTACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAg                                                                                                                         >  7:301171/1‑149 (MQ=255)

                     aCGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGcc                                                                                                                       >  8:124296/1‑149 (MQ=255)

                     aCGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGcc                                                                                                                       >  1:342282/1‑149 (MQ=255)

                        tCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCttt                                                                                                                    <  3:132513/149‑1 (MQ=255)

                         cTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTc                                                                                                                   >  8:52732/1‑149 (MQ=255)

                           tCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCttt                                                                                                                    <  1:357689/146‑1 (MQ=255)

                           tCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCttt                                                                                                                    >  2:357689/1‑146 (MQ=255)

                            cGGCCTGCAGGCCGTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTc                                                                                                                   <  8:16012/146‑1 (MQ=255)

                            cGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTc                                                                                                                   <  8:43164/146‑1 (MQ=255)

                             ggCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGAc                                                                                                                <  7:124296/148‑1 (MQ=255)

                             ggCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACg                                                                                                               >  4:17826/1‑149 (MQ=255)

                              gCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACgt                                                                                                              >  8:317433/1‑149 (MQ=255)

                                 tGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGtt                                                                                                           >  4:351768/1‑149 (MQ=255)

                                 tGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGtt                                                                                                           >  4:198006/1‑149 (MQ=255)

                                   cAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGa                                                                                                         <  8:42199/149‑1 (MQ=255)

                                         cTTGTCTCGGTGCATCACCGAAAACTCCACGCGCTCGCCTTCGCCCAGCATGCGGTGGCCCTCCCCGCGGATGGCCCGAAAGTGCACGAAGATGGCATCGCCCGAGTCTCGGga                                                                                                                                      <  6:175328/114‑1 (MQ=255)

                                         cTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCAt‑t                                                                                                  >  6:88381/1‑149 (MQ=255)

                                            gTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCAt‑tt                                                                                                 <  3:310867/147‑1 (MQ=255)

                                            gTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCAT‑TTGa                                                                                               >  6:334132/1‑149 (MQ=255)

                                             tCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCAT‑TTGa                                                                                               <  3:198006/148‑1 (MQ=255)

                                             tCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCAT‑TTGa                                                                                               <  4:352818/148‑1 (MQ=255)

                                             tCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCAT‑TTGAc                                                                                              <  3:17826/149‑1 (MQ=255)

                                               gcgGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCAT‑TTGACgg                                                                                            >  7:143660/1‑149 (MQ=255)

                                                cgGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCAT‑TTGACGGt                                                                                           >  5:124814/1‑149 (MQ=255)

                                                   tGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGTACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCAT‑TTGACGGTGcc                                                                                        <  8:301171/149‑1 (MQ=255)

                                                     cATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCAT‑ttgacgttgacg                                                                                       >  7:338567/1‑148 (MQ=255)

                                                      aTCAACTAAAAGTCCACGGGCTTGCGTTCGAGCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCTTCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCAT‑TTGACGGTGCCGGt                                                                                     <  8:338567/149‑1 (MQ=255)

                                                       tCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCAT‑TTGACGGTGCCGGTa                                                                                    >  3:369771/1‑149 (MQ=255)

                                                          ccGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGtt                                                                                                           <  8:306269/124‑1 (MQ=255)

                                                          ccGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGtt                                                                                                           >  7:306269/1‑124 (MQ=255)

                                                           cGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCAT‑TTGACGGTGCCGGTATCAc                                                                                <  5:88381/149‑1 (MQ=255)

                                                             aaaaCTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCAT‑TTGACGGTGCCGGTATCACgg                                                                              >  5:164961/1‑149 (MQ=255)

                                                              aaaCTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCAT‑TTGACGGTGCCGGTATCACg                                                                               <  6:236922/147‑1 (MQ=255)

                                                                  tCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCAT‑TTGACGGTGCCGGTATCACGGTTGcc                                                                         >  2:3721/1‑149 (MQ=255)

                                                                   ccACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCAGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCAT‑TTGACGGTGCCGGTATCACGGTTGCCg                                                                        >  5:62346/1‑149 (MQ=255)

                                                                    cACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATCTTGACGGTGCCGGTATCACGGTTGcc                                                                         <  1:3721/148‑1 (MQ=255)

                                                                    cACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCAT‑TTGACGGTGCCGGTATCACGGTTGcc                                                                         <  6:164961/147‑1 (MQ=255)

                                                                    cACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCAT‑TTGACGGTGCCGGTATCACGGTTGCCgg                                                                       >  8:198600/1‑149 (MQ=255)

                                                                     aCGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCAT‑TTGACGGTGCCGGTATCACGGTTGCCgg                                                                       >  1:353791/1‑148 (MQ=255)

                                                                          cTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACg                                                                                                               <  6:299581/104‑1 (MQ=255)

                                                                          cTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACg                                                                                                               >  5:299581/1‑104 (MQ=255)

                                                                            ggCCTTCGACCAGAATGCGGTGGCCCTCACCGCGTATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGtt                                                                                                           <  1:100608/106‑1 (MQ=255)

                                                                            ggCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGtt                                                                                                           >  2:100608/1‑106 (MQ=255)

                                                                             gCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCAT‑TTGACGGTGCCGGTATCACGGTTGCCGGCATCTtgt                                                               <  7:198600/148‑1 (MQ=255)

                                                                             gCCTTCGACCAGAATGCGGTGGCCCTCAACGCGGATGGCCCGAAAGTGCACGATGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGACTTTAGACGTGTTGAACCAT‑TTGACGGTGCAGGAATCACGGTAGCAGGAATCAtgt                                                               >  5:265601/1‑148 (MQ=255)

                                                                               cTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCAT‑TTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGcc                                                            >  2:45891/1‑149 (MQ=255)

                                                                               cTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCAT‑TTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGc                                                             >  5:160501/1‑148 (MQ=255)

                                                                                  cGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCAT‑TTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCAc                                                          <  7:140935/148‑1 (MQ=255)

                                                                                   gACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCAT‑TTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACgt                                                        >  7:216647/1‑149 (MQ=255)

                                                                                    aCCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCTGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCAT‑TTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACgtg                                                       <  6:62346/149‑1 (MQ=255)

                                                                                      cAGAATGCGGTGGCCCTCACCGCGTATGGCGCGAAAGTGCACGTAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCGTTTCGACGTGTTGAACCAT‑TTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCGACGTGtt                                                     <  6:265601/149‑1 (MQ=255)

                                                                                       aGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCAT‑TTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACgtgtg                                                     >  8:293091/1‑147 (MQ=255)

                                                                                        gAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCAT‑TTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGtt                                                     >  3:359718/1‑147 (MQ=255)

                                                                                        gAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCAT‑TTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGa                                                   <  6:160501/149‑1 (MQ=255)

                                                                                         aaTGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGATCCAT‑TTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAg                                                  >  1:261979/1‑149 (MQ=255)

                                                                                         aaTGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCAT‑TTGACGGTGCCGGTATCACGGTTGCCGGCAGCTTGTGCCACGTGTTGAg                                                  >  4:125141/1‑149 (MQ=255)

                                                                                           tGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCAT‑TTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTg                                                <  2:342282/149‑1 (MQ=255)

                                                                                                    ccTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCAT‑TTGACGGTGCCGGTATCACGGTTGCCGGCATCTtgt                                                               >  4:122895/1‑125 (MQ=255)

                                                                                                    ccTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCAT‑TTGACGGTGCCGGTATCACGGTTGCCGGCATCTtgt                                                               <  3:122895/125‑1 (MQ=255)

                                                                                                          cgcgGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCAT‑TTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTgcgc                                 >  6:108402/1‑149 (MQ=255)

                                                                                                                   cccGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCAT‑TTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCgg                        <  7:52732/149‑1 (MQ=255)

                                                                                                                    ccGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCGGTTCGACGTGTTGAAGCAT‑TTGACTGTGCCGGTATCACGGTTGCAGGCATCTTGTGCCAAGTGTTGATTGCCGCGTGCCTTGCGAGTGCTGCgg                        >  4:321712/1‑148 (MQ=255)

                                                                                                                    ccGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCAT‑TTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCggg                       >  5:27830/1‑149 (MQ=255)

                                                                                                                    ccGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCAT‑TTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCggg                       >  7:12208/1‑149 (MQ=255)

                                                                                                                       aaaGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCAT‑TTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGGTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCaa                    <  8:216647/149‑1 (MQ=255)

                                                                                                                             cACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCAT‑TTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCa              >  2:290243/1‑149 (MQ=255)

                                                                                                                              aCGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCAT‑TTGACGGTGCCGGTATCACGGTTGCCg                                                                        >  6:320/1‑90 (MQ=255)

                                                                                                                              aCGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCAT‑TTGACGGTGCCGGTATCACGGTTGCCg                                                                        <  5:320/90‑1 (MQ=255)

                                                                                                                              aCGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCAT‑TTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCAc             <  5:334132/149‑1 (MQ=255)

                                                                                                                              aCGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCAT‑TTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCAc             >  4:140613/1‑149 (MQ=255)

                                                                                                                              aCGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCAT‑TTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCAc             >  2:256413/1‑149 (MQ=255)

                                                                                                                               cGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCAT‑TTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGGCGCGTGCCTTGCGCGTGCTGCGGGCAAAGTCCACg            <  2:353791/149‑1 (MQ=255)

                                                                                                                                tttggtttCATCTCCCGAGACTCGGGCGATGAAACCAAAGCCTTTCGACGTGTTGAACCAT‑TTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACgg           >  1:210305/8‑149 (MQ=255)

                                                                                                                                gAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCAT‑TTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACgg           >  1:192810/1‑149 (MQ=255)

                                                                                                                                gAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCAT‑TTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGCGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACgg           <  4:369771/149‑1 (MQ=255)

                                                                                                                                 aaGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCAT‑TTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGc          <  7:293091/149‑1 (MQ=255)

                                                                                                                                   gATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCAT‑TTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTTCTGCGGGCAAAGCCCACGGCCa        <  8:143660/149‑1 (MQ=255)

                                                                                                                                   gATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCAT‑TTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGc          >  2:186441/1‑147 (MQ=255)

                                                                                                                                         aTCTCCCGAGACTCGGGCGATGAAACCAAAGCCTTTCGACGTGTTGAACCAT‑TTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGCCAGATGca  >  6:67982/1‑149 (MQ=255)

                                                                                                                                         aTCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCAT‑TTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCcgcg                                           >  5:14706/1‑108 (MQ=255)

                                                                                                                                         aTCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCAT‑TTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCcgcg                                           <  6:14706/108‑1 (MQ=255)

                                                                                                                                         aTCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCAT‑TTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGCCAGATGc   <  1:45891/148‑1 (MQ=255)

                                                                                                                                          tctcccGAGACTCGGGCGATGAAACCAAAGCCTTTCGACGTGTTGAACCAT‑TTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGt                                          <  8:69990/105‑1 (MQ=255)

                                                                                                                                          tctcccGAGACTCGGGCGATGAAACCAAAGCCTTTCGACGTGTTGAACCAT‑TTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGt                                          >  7:69990/4‑108 (MQ=255)

                                                                                                                                           cGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCAT‑TTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGc                 <  6:58955/132‑1 (MQ=255)

                                                                                                                                           cGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCAT‑TTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGc                 >  5:58955/1‑132 (MQ=255)

CAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCAT‑TTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGCCAGATGCA  >  NC_002947/1387897‑1388182

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