BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,892,398	C→G	5.9%	E59D (GAG→GAC)	PP_2544 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,892,398	0	C	G	5.9%	749.5 / 24.5	226	E59D (GAG→GAC)	PP_2544	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (101/110); new base G (2/12); total (103/123)
Fisher's exact test for biased strand distribution p-value = 2.35e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 235 total aligned reads displayed.

TGGTCACCATGGCCTTCATGCCTGACTACCAGATCCTGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTGGTTCATCCACGGTGAACCGAGCAGGGCTGCCGCAGCGATCAGGGCCAGGATGGCGAGGGCGATGTAGGTCTTCACGG  >  NC_002947/2892249‑2892543

   tCACCATGGCCTTCATGCCTGACTACCAGATCCTGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTcc                                                                                                                                                >  1:6354615/1‑150 (MQ=255)

   tCACCATGGCCTTCATGCCTGACTACCAGATCCTGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTcc                                                                                                                                                <  1:7591409/150‑1 (MQ=255)

     aCCATGGCCTTCATGCCTGACTACCAGATCCTGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCt                                                                                                                                              >  2:7264777/1‑150 (MQ=255)

      ccATGGCCTTCATGCCTGACTACCAGATCCTGGTCATTTCCACCGGTATTGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTcccgt                                                                                                                                             >  2:4104601/1‑148 (MQ=255)

      ccATGGCCTTCATGCCTGACTACCAGATCCTGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCtt                                                                                                                                             <  2:1909879/150‑1 (MQ=255)

        aTGGCCTTCATGCCTGACTACCTGATCCTGGTCATTTCCCCCGGTATCGCTGCGGGGTTCGTGGTGGCGTTGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCTGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTgg                                                                                                                                           <  1:4036725/150‑1 (MQ=255)

        aTGGCCTTCATGCCTGACTACCAGATCCTGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTgg                                                                                                                                           <  2:7295876/150‑1 (MQ=255)

         tGGCCTTCATGCCTGACTACCAGATCCTGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTggg                                                                                                                                          <  1:2541925/150‑1 (MQ=255)

           gCCTTCATGCCTGACTACCAGATCCTGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGa                                                                                                                                        >  2:6899027/1‑150 (MQ=255)

                cATGCCTGACTACCAGATCCTGGTCATTTCCCCCGGTATCGCGGCGGCGTTGGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGAACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCgg                                                                                                                                   <  1:3796937/150‑1 (MQ=255)

                cATGCCTGACTACCAGATCCTGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGGTCCCCAAGGGAGCCgg                                                                                                                                   <  2:479186/150‑1 (MQ=255)

                cATGCCTGACTACCAGATCCTGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCgg                                                                                                                                   >  2:5308828/1‑150 (MQ=255)

                cATGCCTGACTACCAGATCCTGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCgg                                                                                                                                   <  2:8135904/150‑1 (MQ=255)

                 aTGCCTGACTACCAGATCCTGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGc                                                                                                                                  <  2:4236306/150‑1 (MQ=255)

                 aTGCCTGACTACCAGATCCTGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGc                                                                                                                                  <  2:5648781/150‑1 (MQ=255)

                 aTGCCTGACTACCAGATCCTGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCGTGGGGACCCGGc                                                                                                                                  >  1:1915502/1‑150 (MQ=255)

                  tGCCTGACTACCAGATCCTGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGcc                                                                                                                                 <  1:914723/150‑1 (MQ=255)

                   gCCTGACTACCAGATCCTGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCg                                                                                                                                >  1:2996608/1‑150 (MQ=255)

                    ccTGACTACCAGATCCTGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGGTCCCCAAGGGAGCCGGCCGt                                                                                                                               <  1:3644377/150‑1 (MQ=255)

                     cTGACTACCAGATCCTGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGtt                                                                                                                              >  1:8915205/1‑150 (MQ=255)

                      tGACTACCAGATCCTGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGttt                                                                                                                             >  1:7862238/1‑150 (MQ=255)

                       gACTACCATATCCTGTTCATTTCCACGGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGGTCCCCAAGGGAGCCGGCCGTTTc                                                                                                                            <  2:7370320/150‑1 (MQ=255)

                       gACTACCAGATCCTGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTc                                                                                                                            <  2:1268223/150‑1 (MQ=255)

                        aCTACCAGATCCTGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTcg                                                                                                                           <  1:7196750/150‑1 (MQ=255)

                             cAGATCCTGGTCATTTCCACCGGTATCGCTGCGGCGGTCATGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCgttg                                                                                                                      >  2:777190/1‑150 (MQ=255)

                                 tCCTGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTgg                                                                                                                  <  2:5147772/150‑1 (MQ=255)

                                 tCCTGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTgg                                                                                                                  <  2:8732878/150‑1 (MQ=255)

                                 tCCTGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTgg                                                                                                                  <  1:5253861/150‑1 (MQ=255)

                                  ccTGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGGTGTTggg                                                                                                                 >  2:3943924/1‑150 (MQ=255)

                                   cTGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGc                                                                                                                >  2:2209255/1‑150 (MQ=255)

                                   cTGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGc                                                                                                                <  1:4887064/150‑1 (MQ=255)

                                   cTGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTGGTTGGGc                                                                                                                >  2:2804208/1‑150 (MQ=255)

                                     ggTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTGGGGTTGTTGGGCCt                                                                                                              >  2:3156878/1‑150 (MQ=255)

                                            tCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGCTCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCttt                                                                                                       <  2:8193038/150‑1 (MQ=255)

                                             ccACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCtttt                                                                                                      >  2:6048366/1‑150 (MQ=255)

                                             ccACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCGTGGGGACCCGGCCGTTTCGCGGTGTTGGGGCGGGGCtttt                                                                                                      >  1:5003712/1‑150 (MQ=255)

                                              tcggggTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGGTTCGCCGTGGTGGGCCTGGGCttttt                                                                                                     >  2:5720719/5‑150 (MQ=255)

                                                  ggTATCGCTGCGGCGGTCGTGGTGGCGTTGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCCGCACTACGCTTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGc                                                                                                 <  1:216353/150‑1 (MQ=255)

                                                  ggTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGc                                                                                                 <  1:1567871/150‑1 (MQ=255)

                                                   gTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCg                                                                                                >  2:3008086/1‑150 (MQ=255)

                                                     agcGCTGCGGCGGTCGTGGGGGCGCGGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTGTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCggg                                                                                              <  1:3178747/148‑1 (MQ=255)

                                                      tCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGc                                                                                             >  1:655679/1‑150 (MQ=255)

                                                        gCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCAt                                                                                           >  2:3310886/1‑150 (MQ=255)

                                                          tGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTggg                                                                                                                                          >  2:1690520/1‑101 (MQ=255)

                                                          tGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTggg                                                                                                                                          <  1:1690520/101‑1 (MQ=255)

                                                          tGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCAtgt                                                                                         <  1:369022/150‑1 (MQ=255)

                                                            cggcggTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCATTTTCTGCGGGCATGTgg                                                                                       <  2:8524614/150‑1 (MQ=255)

                                                              gcggTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCggg                                                                                              <  1:5587269/141‑1 (MQ=255)

                                                               cggTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGaaa                                                                                    <  2:655679/150‑1 (MQ=255)

                                                                    gtggtgGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTg                                                                               <  2:3426967/150‑1 (MQ=255)

                                                                     tggtggCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTgg                                                                                                            <  1:6962281/120‑1 (MQ=255)

                                                                     tggtggCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTgg                                                                              <  1:3348440/150‑1 (MQ=255)

                                                                      ggtggCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTTTTTGGCCTGGGGTTTTTATTCGGGCATGTGGAAACCCAGGt                                                                             >  1:6348200/1‑150 (MQ=255)

                                                                           cGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTGGTTCATc                                                                        <  1:1035729/150‑1 (MQ=255)

                                                                            gATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTGGTTCATcc                                                                       <  2:6742467/150‑1 (MQ=255)

                                                                            gATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGGTGTTGGGCCAGGGCTTTTTCTTCGGGCCTGTGGAAAAGCTGGTTAAtac                                                                       >  2:1156885/1‑148 (MQ=255)

                                                                                    tggtgCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTGGTTCATCCACGGTGaa                                                               <  1:2889469/150‑1 (MQ=255)

                                                                                     ggtgCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTGGTTCATCCACGGTGAAc                                                              <  1:3975627/150‑1 (MQ=255)

                                                                                      gtgCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGGTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTGGTTCATCCACGGTGAAcc                                                             >  2:6979044/1‑150 (MQ=255)

                                                                                       tgCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTGGTTCATCCACGGTGAACCg                                                            <  1:5821313/150‑1 (MQ=255)

                                                                                         cACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTGGTTCATCCACGGTGAACCGAg                                                          <  1:7440878/150‑1 (MQ=255)

                                                                                           ccAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTGGTTCATCCACGGTGAACCGAGCa                                                        >  2:5009883/1‑150 (MQ=255)

                                                                                            ccgATCCGATCTGGTAGCATGCTGCACTAAGCGTCAGTCGCATGTGTAAAGGGCGGGCTCCCTGGGGGACCCGGCATTTTGGCTTTGTTGGGCGTGGGATTTTTGTGCGGGCATTTTGAAACGCTGGTTCATCCACGGTGACCCGAGCAg                                                       <  2:6627205/148‑1 (MQ=255)

                                                                                                tCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTGGTTCATCCACGGTGAACCGAGCAGGGCt                                                   <  1:8753125/150‑1 (MQ=255)

                                                                                                        ggTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCATGGGGACCCGGCCGTGTGGGCTTGTTGGGCCGGGGGGTTTTGTTCGGGGAGGTGGGAAACCTGGTTAATCCACCCTGAACCCAGCCCGGCGGCCGCAGc                                           >  1:7916089/1‑150 (MQ=255)

                                                                                                         gTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTGGTTCATCCACGGTGAACCGAGCAGGGCTGCCGCAGCg                                          >  2:4543595/1‑150 (MQ=255)

                                                                                                           agcaagcaGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTGGTTCATCCACGGTGAACCGAGCAGGGCTGCCGCAGCGAt                                        <  2:4847709/150‑1 (MQ=255)

                                                                                                             caagcaGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTGGTTCATCCACGGTGAACCGAGCAGGGCTGCCGCAGCGATCa                                      <  1:6278857/150‑1 (MQ=255)

                                                                                                              aagcaGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCCTGGGGAAACGCTGGTTCATCCACGGTGAACCGCGCCGGGCTGCCGCcgcgcgcag                                     >  2:5446598/1‑145 (MQ=255)

                                                                                                              aagcaGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTGGTTCATCCACGGTGAACCGAGCAGGGCTGCCGCAGCGATCAg                                     <  2:359264/150‑1 (MQ=255)

                                                                                                               agcagcACTACGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCGTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCGTGGTCTGCGGGCCTGGGGACACGCTGGTTCATCCACGGTGAACCGAGCAGGGCTGCCGCAGCGATCAgg                                    >  2:5831917/1‑150 (MQ=255)

                                                                                                               agcagcACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTGGTTCATCCACGGTGAACCGAGCAGGGCTGCCGCAGCGATCAgg                                    <  2:4486261/150‑1 (MQ=255)

                                                                                                                  agcaCTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTGGTTCATCCACGGTGAACCGAGCAGGGCTGCCGCAGCGATCAGGGcc                                 <  2:5635282/150‑1 (MQ=255)

                                                                                                                   gcaCTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTGGTTCATCCACGGTGAACCGAGCAGGGCTGCCGCAGCGATCAGGGCCa                                <  1:6034321/150‑1 (MQ=255)

                                                                                                                      cTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTGGTTCATCCACGGTGAACCGAGCAGGGCTGCCGCAGCGATCAGGGCCAGGa                             <  1:4225721/150‑1 (MQ=255)

                                                                                                                        aaGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTGGTTCATCCACGGTGAACCGAGCAGGGCTGCCGCAGCGATCAGGGCCAGGATg                           <  1:9118957/150‑1 (MQ=255)

                                                                                                                        aaGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTGGTTCATCCACGGTGAACCGAGCAGGGCTGCCGCAGCGATCAGGGCCAGGATg                           <  1:532913/150‑1 (MQ=255)

                                                                                                                         aGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCGTGGGGACCCGGCCGGGTCGCGGTGTTGGGCCTGGGCGTGTTCGGGGGGGCAGTGGTAAAGCTGGTTGATACCCGGGGCACCGGGCCGGGCTGCCGCCGCGCTCAGGGGCACGATgg                          >  2:1565151/1‑150 (MQ=255)

                                                                                                                          gCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTGGTTCATCCACGGTGAACCGAGCAGGGCTGCCGCAGCGATCAGGGCCAGGATGGc                         <  2:5238648/150‑1 (MQ=255)

                                                                                                                          gCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTGGTTCATCCACGGTGAACCGAGCAGGGCTGCCGCAGCGATCAGGGCCAGGATGGc                         <  1:6146636/150‑1 (MQ=255)

                                                                                                                            gTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAGCGCTGGTTCATCCACGGTGAACCGAGCAGGGCTGCCGCAGCGATCAGGGCCAGGATGGCGa                       >  1:9101386/1‑150 (MQ=255)

                                                                                                                            gTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTGGTTCATCCACGGTGAACCGAGCCGGGCTGCCGCAGCGATCAGGGCCAGGTTggagt                       >  1:6843742/1‑147 (MQ=255)

                                                                                                                            gTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTGGTTCATCCACGGTGAACCGAGCAGGGCTGCCGCAGCGATCAGGGCCAGGATGGCGa                       <  1:1209335/150‑1 (MQ=255)

                                                                                                                            gTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTGGTTCATCCACGGTGAACCGAGCAGGGCTGCCGCAGCGATCAGGGCCAGGATGGCGa                       <  1:4453819/150‑1 (MQ=255)

                                                                                                                             tCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTGGTTCATCCACGGTGAACCGAGCAGGGCTGCCGCAGCGATCAGGGCCAGGATGGCGAg                      <  1:8377201/150‑1 (MQ=255)

                                                                                                                              cACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTGGTTCATCCACGGTGAACCGAGCAGGGCTGCCGCAGCGATCAGGGCCAGGATGGCGAgg                     <  2:2495072/150‑1 (MQ=255)

                                                                                                                              cACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTGGTTCATCCACGGTGAACCGAGCAGGGCTGCCGCAGCGATCAGGGCCAGGATGGCGAgg                     <  1:1412475/150‑1 (MQ=255)

                                                                                                                               aCTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTGGTTCATCCACGGTGAACCGAGCAGGGCTGCCGCAGCGATCAGGGCCAGGATGGCGAggg                    <  1:2520127/150‑1 (MQ=255)

                                                                                                                               aCTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTGGTTCATCCACGGTGAACCGAGCAGGGCTGCCGCAGCGATCAGGGCCAGGATGGCGAggg                    <  1:6740117/150‑1 (MQ=255)

                                                                                                                               aCTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTGGTTCATCCACGGTGAACCGAGCAGGGCTGCCGCAGCGATCAGGGCCAGGATGGCGAggg                    <  1:777190/150‑1 (MQ=255)

                                                                                                                                 tCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGGTTCGCGGTGGGGGGCCTGGGCGTTTTCTGCGGGCATGTGGAAACGCTGGTTCATCCACGGTGCACCGAGCAGGGCTGCCGCAGCGAGCAGGGCCAGGATGGCGAGGGCg                  >  2:6635011/1‑150 (MQ=255)

                                                                                                                                      tctGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTGGCGTTGTTTGGCCGGGGCGTTTTCTGCGGGCATGTGGTAACGCTGGTTCATCCACGGTGAACCGAGCAGGGCTGCCGCAGCGATCAGGGCCAGGATGGCGAGGGCGATGTa             >  2:1944621/1‑150 (MQ=255)

                                                                                                                                      tctGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTGGTTCATCCACGGTGAACCGAGCAGGGCTGCCGCAGCGATCAGGGCCAGGATGGCGAGGGCGATGTa             <  2:2366204/150‑1 (MQ=255)

                                                                                                                                       ctGGACACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTGGTTCATCCACGGTGAACCGAGCAGGGCTGCCGCAGCGATCAGGGCCAGGATGGCGAGGGCGATGTAg            <  1:5423510/150‑1 (MQ=255)

                                                                                                                                       ctGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTGGTTCATCCACGGTGAACCGAGCAGGGCTGCCGCAGCGATCAGGGCCAGGATGGCGAGGGCGATGTAg            <  2:3553193/150‑1 (MQ=255)

                                                                                                                                       ctGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTGGTTCATCCACGGTGAACCGAGCAGGGCTGCCGCAGCGATCAGGGCCAGGATGGCGAGGGCGATGTAg            <  1:8327846/150‑1 (MQ=255)

                                                                                                                                               ggccgggTCCCCAAGGGAGCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTGGTTCATCCACGGTGAACCGAGCAGGGCTGCCGCAGCGATCAGGGCCAGGATGGCGAGGGCGATGTAGGTCTTCAc    >  1:6781982/1‑150 (MQ=255)

                                                                                                                                               ggccggcTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTGGTTCATCCACGGTGAACCGAGCAGGGCTGCCGCAGCGATCAGGGCCAGGATGGCGAGGGCGATGTAGGTCTTCAc    <  2:324370/150‑1 (MQ=255)

                                                                                                                                               ggccggcTCCCTTGGGGACCCGGCCGTTTAGCGTTGTTTGGCCTGGGCTTTTTCTGCGGGCATGTGGAAAAGCTGGTTAATCCACGGTGAAACGAGCAGGGCTGCCGCAGCGATCAGGGCCAGGATGGCGAGGGCGATGTAGGTCTTGAc    >  1:8185931/1‑150 (MQ=255)

                                                                                                                                                 ccggcTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTGGTTCATCCACGGTGAACCGAGCAGGGCTGCCGCAGCGATCAGGGTCAGGATGGCGAGGGCGATGTAGGTCTTCACgg  <  1:8697638/150‑1 (MQ=255)

TGGTCACCATGGCCTTCATGCCTGACTACCAGATCCTGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTGGTTCATCCACGGTGAACCGAGCAGGGCTGCCGCAGCGATCAGGGCCAGGATGGCGAGGGCGATGTAGGTCTTCACGG  >  NC_002947/2892249‑2892543

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 11 ≤` `ATCG/ATCG` `< 25 ≤` `ATCG/ATCG` `< 37 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`