BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	1,070,248	+GA	intergenic (+48/+94)	aroP‑I → / ← PP_0928	aromatic amino acid transport protein/K+‑dependent Na+/Ca+ exchanger‑like protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,070,246	1	.	G	89.0%	380.3 / 54.5	163	intergenic (+46/+96)	aroP‑I/PP_0928	aromatic amino acid transport protein/K+‑dependent Na+/Ca+ exchanger‑like protein
Reads supporting (aligned to +/- strand): ref base . (8/10); new base G (81/64); total (89/74)
Fisher's exact test for biased strand distribution p-value = 4.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.78e-01
*	NC_002947	1,070,246	2	.	A	88.9%	351.8 / 48.1	163	intergenic (+46/+96)	aroP‑I/PP_0928	aromatic amino acid transport protein/K+‑dependent Na+/Ca+ exchanger‑like protein
Reads supporting (aligned to +/- strand): ref base . (8/10); new base A (80/64); total (89/74)
Fisher's exact test for biased strand distribution p-value = 4.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.06e-01

Only 100 of 167 total aligned reads displayed.

GTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGT‑‑GAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCC  >  NC_002947/1070123‑1070358

||

gTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCtgt                                                                                                                      <  1:866061/124‑1 (MQ=255)

gTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCtgt                                                                                                                      <  1:866060/124‑1 (MQ=255)

 tAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCtgt                                                                                                                      >  2:866277/1‑123 (MQ=255)

   aTGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGT‑‑Ga                                                                                                                  <  2:866062/123‑1 (MQ=255)

    tGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCtgt                                                                                                                      >  1:866464/1‑120 (MQ=255)

    tGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGT‑‑Gag                                                                                                                 >  2:866311/1‑122 (MQ=255)

     gCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGT‑‑Ga                                                                                                                  >  2:866489/1‑121 (MQ=255)

      cTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTCTGCCGTTAGGTTTTCGGCGCCTGT‑‑Gagat                                                                                                               <  1:866064/123‑4 (MQ=255)

        gATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGT‑‑Gagatcga                                                                                                            >  1:866259/1‑118 (MQ=255)

        gATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGT‑‑Gagatcga                                                                                                            >  1:866263/1‑118 (MQ=255)

        gATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGT‑‑Gagatc                                                                                                              <  1:866067/122‑5 (MQ=255)

        gATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGACCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGT‑‑Gagatcga                                                                                                            >  2:866172/1‑118 (MQ=255)

         aTCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGT‑‑Gagatcga                                                                                                            >  2:866365/1‑117 (MQ=255)

         aTCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGT‑‑Gagatcga                                                                                                            <  1:866070/123‑7 (MQ=255)

         aTCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGT‑‑Gagatcg                                                                                                             <  2:866068/122‑6 (MQ=255)

             cGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGcgcc                                                                                                       <  1:866073/124‑1 (MQ=255)

                cATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGcgccgc                                                                                                     >  2:866528/1‑123 (MQ=255)

                 aTGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCcgcgc                                                                                                   <  1:866075/124‑1 (MQ=255)

                   ggATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAgcgc                                                                                                        <  2:866076/117‑1 (MQ=255)

                   ggATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCCGCgcgg                                                                                                 >  1:866341/1‑124 (MQ=255)

                    gATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCCGCgcggg                                                                                                <  1:866077/124‑1 (MQ=255)

                        gctgtgctgtGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCCGCgcgggcgg                                                                                             <  1:866080/123‑1 (MQ=255)

                        gctgtgctgtGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCCGCgcgggcgg                                                                                             <  2:866081/123‑1 (MQ=255)

                         ctgtgctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGTGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCCGCGCGggcggc                                                                                            <  2:866084/123‑1 (MQ=38)

                         ctgtgctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCCGCGCGGgcggcg                                                                                           >  1:866274/1‑124 (MQ=255)

                          tgtgctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCCGCGCGGgcggcg                                                                                           >  2:866302/1‑123 (MQ=255)

                          tgtgctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCCGCGCGGgcggcg                                                                                           >  2:866294/1‑123 (MQ=255)

                            tgctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGT‑‑GATAGCGAGCGCCGCGCGGGCGGCGCt                                                                                         <  2:866086/123‑1 (MQ=255)

                              ctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCCGCGCGggcggc                                                                                            <  2:866088/118‑1 (MQ=255)

                              ctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTc                                                                                        >  1:866228/1‑122 (MQ=255)

                               tgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGAt                                                                                     >  2:866247/1‑124 (MQ=255)

                                 tgggtggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGAtct                                                                                   <  2:866091/124‑1 (MQ=255)

                                 tgggtggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGAtct                                                                                   <  2:866092/124‑1 (MQ=255)

                                 tgggtggCTTGCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGAtct                                                                                   <  1:866090/124‑1 (MQ=38)

                                  gggtggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGCGAT‑‑CGAGCGCCGCGCGGGCGGCGCTc                                                                                        >  2:866359/1‑118 (MQ=255)

                                    gtggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCCGCgcgg                                                                                                 >  2:866275/1‑107 (MQ=255)

                                    gtggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGAtctc                                                                                  >  1:866313/1‑122 (MQ=255)

                                     tggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTcaca                                                                               <  2:866095/124‑1 (MQ=255)

                                     tggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTcaca                                                                               >  1:866233/1‑124 (MQ=255)

                                       gCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAgg                                                                             <  1:866096/124‑1 (MQ=255)

                                              aGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGAt                                                                         >  1:866421/1‑121 (MQ=255)

                                              aGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGa                                                                       >  2:866344/1‑123 (MQ=255)

                                              aGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGACGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGaa                                                                      >  1:866195/1‑124 (MQ=37)

                                               gCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGaaa                                                                     <  1:866098/124‑1 (MQ=255)

                                                   aagaagGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGaa                                                                      <  2:866102/119‑1 (MQ=255)

                                                   aagaagGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGaa                                                                      >  1:866102/1‑119 (MQ=255)

                                                   aagaagGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGcc                                                                  <  2:866103/123‑1 (MQ=255)

                                                     gaagGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCt                                                                 >  1:866317/1‑122 (MQ=255)

                                                       agGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTa                                                                >  1:866557/1‑121 (MQ=255)

                                                          cTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCa                                                          <  2:866105/124‑1 (MQ=255)

                                                           tCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAgcgc                                                                                                        >  1:866438/1‑77 (MQ=38)

                                                             gCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCCGCGCGGGAGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATg                                                        <  1:866107/123‑1 (MQ=38)

                                                             gCCAGGCCAAGGGCTGATCGGGTTTGTTCGGAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAgg                                                            <  2:866108/119‑1 (MQ=37)

                                                              ccAGTCCCAGGTGTGATCGGATTTGTTCGTAAGGGATATGCCGTAAGGTTCTCGGAGCATGTGAGAT‑‑CGAGTGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATg                                                        >  1:866330/1‑122 (MQ=21)

                                                              ccAGGCCAAGGGCTGATCGGCTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGc                                                       >  1:866368/1‑123 (MQ=37)

                                                                  gCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCa                                                          >  1:866353/1‑116 (MQ=255)

                                                                  gCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTc                                                    <  2:866112/122‑1 (MQ=255)

                                                                  gCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTg                                                  <  1:866113/124‑1 (MQ=255)

                                                                  gCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTg                                                  >  1:866410/1‑124 (MQ=255)

                                                                  gCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTg                                                  <  2:866110/124‑1 (MQ=255)

                                                                       gggCTGTTCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGc                                                       <  2:866115/114‑1 (MQ=37)

                                                                        ggCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCt                                                   <  2:866117/117‑1 (MQ=255)

                                                                           tGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCt                                          >  1:866354/1‑123 (MQ=255)

                                                                              tCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTggcg                                               <  2:866119/115‑1 (MQ=255)

                                                                                gggTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTggc                                                <  1:866120/112‑1 (MQ=255)

                                                                                gggTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGAcg                                      >  2:866297/1‑122 (MQ=255)

                                                                                 ggTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGt                                   <  1:866122/124‑1 (MQ=255)

                                                                                  gTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGAc                                       >  2:866461/1‑119 (MQ=255)

                                                                                   tttGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAgcgc                                                                                                        <  2:866124/53‑1 (MQ=21)

                                                                                   tttGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTggcgg                                              >  1:866494/1‑111 (MQ=255)

                                                                                     tGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTccc                               <  2:866125/124‑1 (MQ=255)

                                                                                      gTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGc                                                                   >  2:866237/1‑87 (MQ=21)

                                                                                        tCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCtt                            >  1:866208/1‑124 (MQ=255)

                                                                                        tCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCt                             <  2:866127/123‑1 (MQ=255)

                                                                                          gTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGCcgcg                                    >  1:866531/1‑114 (MQ=37)

                                                                                          gTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGAc                                       <  1:866128/111‑1 (MQ=255)

                                                                                           tAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGAtct                                                                                   >  2:866319/1‑66 (MQ=17)

                                                                                           tAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTc                           >  1:866312/1‑122 (MQ=255)

                                                                                           tAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCa                         >  2:866405/1‑124 (MQ=255)

                                                                                            aaGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCa                         >  2:866457/1‑123 (MQ=255)

                                                                                              ggggTATGCCGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGc                      >  2:866402/1‑124 (MQ=255)

                                                                                               gggTATGCCGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGcc                     >  1:866270/1‑124 (MQ=255)

                                                                                               gggTATGCCGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGcc                     <  1:866131/124‑1 (MQ=255)

                                                                                                ggTATGCCGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCcgcg                                                                                                    >  1:866395/1‑44 (MQ=14)

                                                                                                 gTATGCCGTTAGGGTTTCGGCGCCTGTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTcccc                              <  2:866132/113‑1 (MQ=37)

                                                                                                  tATGCCGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGAcgc                                     >  2:866447/1‑105 (MQ=255)

                                                                                                  tATGCCGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCt                  >  2:866323/1‑124 (MQ=255)

                                                                                                  tATGCCGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCt                  >  1:866375/1‑124 (MQ=255)

                                                                                                    tGCCGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCCGCGCGGGCGCCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCa                >  1:866293/1‑124 (MQ=38)

                                                                                                     gCCGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTc                 >  2:866361/1‑122 (MQ=255)

                                                                                                      ccGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATg                                                        <  1:866133/82‑1 (MQ=21)

                                                                                                      ccGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGccc                    >  1:866431/1‑118 (MQ=255)

                                                                                                      ccGTTAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCCGCGCGGACGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCTTTCCATGCCCCt                  >  2:866271/1‑120 (MQ=37)

                                                                                                         ttAGGTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAtgc            >  1:866135/1‑123 (MQ=255)

                                                                                                             gTTTTCGGCGCCTGTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAtgctg          <  2:866135/121‑1 (MQ=255)

                                                                                                              ttttCGGCGCCTGTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGAcgcg                                    >  1:866553/1‑94 (MQ=25)

                                                                                                                ttCGGCGCCTGTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCg    >  1:866334/1‑124 (MQ=255)

                                                                                                                ttCGGCGCCTGTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCg    <  1:866137/124‑1 (MQ=255)

                                                                                                                 tCGGCGCCTGTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTcccc                              <  1:866138/97‑1 (MQ=25)

                                                                                                                   ggCGCCTGTGAGAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGcc  <  1:866140/123‑1 (MQ=255)

||

GTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGT‑‑GAT‑‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCC  >  NC_002947/1070123‑1070358

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 16 ≤` `ATCG/ATCG` `< 26 ≤` `ATCG/ATCG` `< 33 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`