BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,970,426	A→C	6.5%	W42G (TGG→GGG)	PP_5723 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,970,426	0	A	C	6.5%	511.2 / 9.0	163	W42G (TGG→GGG)	PP_5723	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (40/109); new base C (11/2); total (51/112)
Fisher's exact test for biased strand distribution p-value = 5.96e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 167 total aligned reads displayed.

ATTGGCCAGTCGCCAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATT‑TGT‑GCCTTGATGCTGATATTGCCGGGGCCGCTGT‑GCAGCCCAAT‑CGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTGGCGGCCATTTCTCCGTTTCAGTGGACAGCCCCAT  >  NC_002947/5970278‑5970575

aTTGGCCAGTCGCCAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCaa                                                                                                                                                          <  2:4361922/150‑1 (MQ=255)

 ttGGCCAGTCGCCAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAg                                                                                                                                                         <  1:6697335/150‑1 (MQ=255)

 ttGGCCAGTCGCCAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAg                                                                                                                                                         <  1:6697280/150‑1 (MQ=255)

 ttGGCCAGTCGCCAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTGGAAAACTGCACGCGATCTTTATGGGAGCGGCCTTGGGGGGCGGGTGGGGGGGGGAGCGCCCCCCCAg                                                                                                                                                         >  2:4007800/1‑150 (MQ=255)

     ccAGTCGCCAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGAtt‑t                                                                                                                                                    <  2:1741253/150‑1 (MQ=255)

     ccAGTCGCCAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGAtt‑t                                                                                                                                                    <  1:4006419/150‑1 (MQ=255)

     ccAGTCGCCAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGAtt‑t                                                                                                                                                    <  2:5917650/150‑1 (MQ=255)

     ccAGTCGCCAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCGTGCGTGACGATTTGGGTGCGCAGCGCCCCCCCAGaagt                                                                                                                                                     >  1:7754742/1‑147 (MQ=255)

      cAGTCGCCAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGGCGGGGGTGGGGGGCGCGGCGCCCCCCCAgggttg                                                                                                                                                    >  2:6027490/1‑145 (MQ=255)

         tCGCCAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATT‑TGT‑Gc                                                                                                                                               <  1:5299278/150‑1 (MQ=255)

           gCCAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATT‑TGT‑GCCt                                                                                                                                             <  1:1645256/150‑1 (MQ=255)

            ccAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAACTGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATT‑TGT‑GCCtt                                                                                                                                            <  2:7674150/150‑1 (MQ=255)

            ccAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATT‑TGT‑GCCtt                                                                                                                                            <  1:6262341/150‑1 (MQ=255)

            ccAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATT‑TGT‑GCCtt                                                                                                                                            <  2:2886250/150‑1 (MQ=255)

            ccAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATT‑TGT‑GCCtt                                                                                                                                            <  1:8199968/150‑1 (MQ=255)

            ccAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCGTGGGGGGGGGGTTGGGGGGGGGGCGCCCCCCAACAAT‑tgt‑ggggg                                                                                                                                            >  2:2502323/1‑146 (MQ=255)

               gCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGGGGCGCGGGAGGGGGGCGGAGCGCCCCCCCAGAGT‑ttt‑ggcgtttt                                                                                                                                         >  2:5751070/1‑143 (MQ=255)

                cACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATT‑TGT‑GCCTTGATg                                                                                                                                        <  2:9043242/150‑1 (MQ=255)

                   gggAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATT‑TGT‑GCCTTGATGCTg                                                                                                                                     <  2:440458/150‑1 (MQ=255)

                    ggAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATT‑TGT‑GCCTTGATGCTGa                                                                                                                                    <  2:1599910/150‑1 (MQ=255)

                    ggAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATT‑TGT‑GCCTTGATGCTGa                                                                                                                                    <  2:3952319/150‑1 (MQ=255)

                     gAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGGGGCGGGGTTTTGGTGGGGAGCGCCCCCCACCAAT‑TTT‑GGCATGATGCTggt                                                                                                                                   >  2:3462704/1‑148 (MQ=255)

                     gAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGGCGCGTTTGGGGGGCGAAGCGCCCCCCCACCTG‑TGT‑GCCGTGATGCTGat                                                                                                                                   >  2:7477772/1‑150 (MQ=255)

                       aaCATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATT‑TGT‑GCCTTGATGCTGatat                                                                                                                                 <  1:4113554/150‑1 (MQ=255)

                          aTGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATT‑TGT‑GCCTTGATGCTGATATTGc                                                                                                                              <  1:8623826/150‑1 (MQ=255)

                           tGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATT‑TGT‑GCCTTGATGCTGATATTGcc                                                                                                                             <  1:3234955/150‑1 (MQ=255)

                             cacaATTTCGGTGCGCAGTAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATT‑TGT‑GCCTTGATGCTGATATTGCCgg                                                                                                                           <  2:4798352/150‑1 (MQ=255)

                             cacaATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATT‑TGT‑GCCTTGATGCTGATATTGCCgg                                                                                                                           <  1:1395543/150‑1 (MQ=255)

                             cacaATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATT‑TGT‑GCCTTGATGCTGATATTGCCgg                                                                                                                           <  1:8730044/150‑1 (MQ=255)

                              acaATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATT‑TGT‑GCCTTGATGCTGATATTGCCggg                                                                                                                          <  1:5641411/150‑1 (MQ=255)

                               caATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATT‑TGT‑GCCTTGATGCTGATATTGCCgggg                                                                                                                         <  2:9115211/150‑1 (MQ=255)

                                aaTTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATT‑TGT‑GCCTTGATGCTGATATTGCCGGGGc                                                                                                                        <  2:5439568/150‑1 (MQ=255)

                                  tttCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATT‑TGT‑GCCTTGATGCTGATATTGCCGGGGCCg                                                                                                                      <  1:1044604/150‑1 (MQ=255)

                                    tCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATT‑TGT‑GCCTTGATGCTGATATTGCCGGGGCCGCt                                                                                                                    <  2:420166/150‑1 (MQ=255)

                                    tCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATT‑TGT‑GCCTTGATGCTGATATTGCCGGGGCCGCt                                                                                                                    <  1:9019623/150‑1 (MQ=255)

                                     cGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATT‑TGT‑GCCTTGATGCTGATATTGCCGGGGCCGCtg                                                                                                                   <  2:7527879/150‑1 (MQ=255)

                                     cGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATT‑TGT‑GCCTTGATGCTGATATTGCCGGGGCCGCtg                                                                                                                   <  1:2372537/150‑1 (MQ=255)

                                         gcgcAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATT‑TGT‑GCCTTGATGCTGATATTGCCGGGGCCGCTGT‑GCa                                                                                                              <  1:8110614/150‑1 (MQ=255)

                                           gcagcaGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATT‑TGT‑GCCTTGATGCTGATATTGCCGGGGCCGCTGT‑GCAGc                                                                                                            <  1:9119125/150‑1 (MQ=255)

                                               gagaaatggccgccagtgtagcacgcacaaacctggcgcctgagtgcacattggcctttgtgtGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATT‑TGT‑GCCTTGATGCTGATATTGCCGGGGCCGCTGT‑GCAGCCCaa                                                                                                        >  2:1067401‑M1/64‑150 (MQ=255)

                                               gagaaatggccgccagtgtagcacgcacaaacctggcgcctgagtgcacattggcctttgtgtGCGGCCTTGCGTCGCGATTGCGGTGCGAAGCGCCCCCCAAGATT‑TGT‑GCCTTGATGCTGATATTGCCGGGGCCGCTGT‑GCAGCCCaa                                                                                                        >  2:1067361‑M1/64‑150 (MQ=255)

                                               cagcGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATT‑TGT‑GCCTTGATGCTGATATTGCCGGGGCCGCTGT‑GCAGCCCaa                                                                                                        <  1:7611895/150‑1 (MQ=255)

                                                agcGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATT‑TGT‑GCCTTGATGCTGATATTGCCGGGGCCGCTGT‑GCAGCCCAAt                                                                                                       <  2:4901596/150‑1 (MQ=255)

                                                agcGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATT‑TGT‑GCCTTGATGCTGATATTGCCGGGGCCGCTGT‑GCAGCCCAAt                                                                                                       <  2:5802605/150‑1 (MQ=255)

                                                agaaatggccgccagtgtagcacgcacaaacctggcgcctgagtgcacattggcctttgtgtGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATT‑TGT‑GCCTTGATGCTGATATTGCCGGGGCCGCTGT‑GCAGCCCAAt                                                                                                       >  2:6490043‑M1/63‑150 (MQ=255)

                                                   ggggCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATT‑TGT‑GCCTTGATGCTGATATTGCCGGGGCCGCTGT‑GCAGCCCAAT‑cgc                                                                                                   <  1:5324215/150‑1 (MQ=255)

                                                   ggggCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATT‑TGT‑GCCTTGATGCTGATATTGCCGGGGCCGCTGT‑GCAGCCCAAT‑cgc                                                                                                   <  2:2595527/150‑1 (MQ=255)

                                                     tggccgccagtgtagcacgcacaaacctggcgcctgagtgcacattggcctttgtgtGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATT‑TGT‑GCCTTGATGCTGATATTGCCGGGGCCGCTGT‑GCAGCCCAAT‑CGCGa                                                                                                 >  1:3232133‑M1/58‑150 (MQ=255)

                                                     ggCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATT‑TGT‑GCCTTGATGCTGATATTGCCGGGGCCGCTGT‑GCAGCCCAAT‑CGCGa                                                                                                 <  2:6273573/150‑1 (MQ=255)

                                                     ggCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATT‑TGT‑GCCTTGATGCTGATATTGCCGGGGCCGCTGT‑GCAGCCCAAT‑CGCGa                                                                                                 <  2:619418/150‑1 (MQ=255)

                                                     ggCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATT‑TGT‑GCCTTGATGCTGATATTGCCGGGGCCGCTGT‑GCAGCCCAAT‑CGCGa                                                                                                 <  1:8996746/150‑1 (MQ=255)

                                                     ggCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATT‑TGT‑GCCTTGATGCTGATATTGCCGGGGCCGCTGT‑GCAGCCCAAT‑CGCGa                                                                                                 <  1:7851767/150‑1 (MQ=255)

                                                     ggCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATT‑TGT‑GCCTTGATGCTGATATTGCCGGGGCCGCTGT‑GCAGCCCAAT‑CGCGa                                                                                                 <  2:7185001/150‑1 (MQ=255)

                                                      ggccgccagtgtagcacgcacaaacctggcgcctgagtgcacattggcctttgtgtGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATT‑TGT‑GCCTTGATGCTGATATTGCCGGGGCCGCTGT‑GCAGCCCAAT‑CGCGAc                                                                                                >  2:8440210‑M1/57‑150 (MQ=255)

                                                        aaGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATT‑TGT‑GCCTTGATGCTGATATTGCCGGGGCCGCTGT‑GCAGCCCAAT‑CGCGACGc                                                                                              <  2:1747968/150‑1 (MQ=255)

                                                          gCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATT‑TGT‑GCCTTGATGCTGATATTGCCGGGGCCGCTGT‑GCAGCCCAAT‑CGCGACGCaa                                                                                            <  2:8159421/150‑1 (MQ=255)

                                                               gcgAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGGCGCGAGTGGGGTGCGCAGCCCCCCCCCACAAT‑TTT‑TCCGTGATGCTGGTATTGCCCGGGCCCCGGG‑GCAGCCCAAA‑AGCGCCGCAACGCCg                                                                                       >  2:8957847/1‑150 (MQ=255)

                                                                 gAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCGTTCGTGGCGGGTTGGGTGCGGAGCGCCCCCCCAGATT‑TGT‑TGCGTGTTGGTGATATTTGCGGGGGCGCGGT‑GGGGCCCCAA‑AGAGACGCCACGCCGca                                                                                     >  2:1422905/1‑150 (MQ=255)

                                                                            aaCAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGGGGGGCGGTTTGGGGGGGGAGCGCCCCCCAAGATT‑TTT‑GGCGTTATGGTGTTAATGCCGGGGGCGCGGG‑GCAGCCCACT‑AGCGACGCAACGCCGCACACACAGGccc                                                                          >  2:5199059/1‑149 (MQ=255)

                                                                                   tcaaaaCGGCCCGCGATCTTGGGGGGTGCGGCCTGGCGGCGCGTGGGGGGGGCGACCCCCCCCCCCAG‑TT‑TGTGGCCTTCTTGCTCATGTGGCCGGGGCGGCTGG‑GCCGCCCACC‑CGCCCCGCAAGGCCGCACACAAAGGCCAATGTGCa                                                                   <  2:5903101/148‑1 (MQ=255)

                                                                                               gctttttttGTGGGTGGGGCTTTGGGTCGGTTTGGGGGTGCGACGCGCCCCCCATTTTT‑TGT‑GCCTTTTTGCT‑TTTTTGGCGGGGGCGGGGTGGCCGCCCAAT‑CGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAg                                                       <  2:3753742/144‑1 (MQ=255)

                                                                                               gcgATCTTTGTGGGAGCGGCCATGCGTGGCCGTTGGGGGGCGGAGCCCCCCCCCAGATT‑TGT‑TGCATGGTGGTTGTTTTGGCGGGGCCGCGGG‑GCAGCCCCAG‑CGCGGCCCCAGGCCCCCCCCCCAGGCCCAGGTGGCCGCAGGCGcccg                                                       >  2:5262659/1‑148 (MQ=255)

                                                                                                     tttGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATT‑TGT‑GCCTTGATGCTGATATTGCCGGGGCCg                                                                                                                      <  2:6842107/83‑1 (MQ=255)

                                                                                                          gggAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATT‑TGT‑GCCTTGATGCTGATATTGCCGGGGCCGCt                                                                                                                    >  1:6478969/1‑80 (MQ=255)

                                                                                                           gctgCGGCCTTGCGTCGGTTGTGGGGGGCTACGC‑CCCCCCAATTTTGTGT‑GCGTTTTGGCTGATTTGGGGGGGGCCGCGGT‑GCACCCCAAG‑CGCGCCCCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTgtgcgtgc                                           <  2:5488253/147‑1 (MQ=255)

                                                                                                                        acaaatcTTGGGGGGCGCTTCGCACCCCAATATT‑TGT‑GCCTTGATGCTGATATTGCCGGGGCCGCTGT‑GCAGCCCAAT‑CGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTGGCGGcc                              <  2:7660591‑M1/143‑1 (MQ=255)

                                                                                                                        GTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATT‑TGT‑GCCTTGATGCTGATATTGCCGGGGCCGCTGT‑GCAGCCCAAT‑CGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTGGCGGcc                              <  1:8450004‑M1/150‑1 (MQ=255)

                                                                                                                        GTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATT‑TGT‑GCCTTGATGCTGATATTGCCGGGGCCGCTGT‑GCAGCCCAAT‑CGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTGGCGGcc                              <  1:377231‑M1/150‑1 (MQ=255)

                                                                                                                         caaatcTTGGGGTGCGAAGCGCCCCCCAAGATT‑TGT‑GCCTTGATGCTGATATTGCCGGGGCCGCTGT‑GCAGCCCAAT‑CGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTGGCGGCCa                             <  2:1443955‑M1/144‑1 (MQ=255)

                                                                                                                         TCGCGATTGGGGTGCGAAGCGCCCCCCAAGATT‑TGT‑GCCTTGATGCTGATATTGCCGGGGCCGCTGT‑GCAGCCCAAT‑CGCGACGCAAGGCCGCACACAATGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTGGCGGCCa                             <  2:3569599‑M1/150‑1 (MQ=255)

                                                                                                                         TCGCGATTGGGGTGCGAAGCGCCCCCAAAGATT‑TGT‑GCCTTGATGCTGATATTGCCGGGGCCGCTGT‑GCAGCCCAAT‑CGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTGGCGGCCa                             <  2:6794758‑M1/150‑1 (MQ=255)

                                                                                                                          aaatcTTGGGGGGCGCTTCGCACCCCAATATT‑TGT‑GCCTTGATGCTGATATTGCCGGGGCCGCTGT‑GCAGCCCAAT‑CGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTGGCGGCCAt                            <  2:4092114‑M1/145‑1 (MQ=255)

                                                                                                                            atcTTGGGGGGCGCTTCGCACCCCAAGATT‑TGT‑GCCTTGATGCTGATATTGCCGGGGCCGCTGT‑GCAGCCCAAT‑CGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTGGCGGCCAttt                          <  1:8976381‑M1/147‑1 (MQ=255)

                                                                                                                            atcTTGGGGGGCGCTTCGCACCCCAAGATT‑TGT‑GCCTTGATGCTGATATTGCCGGGGCCGCTGT‑GCAGCCCAAT‑CGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTGGCGGCCAttt                          <  1:9157706‑M1/147‑1 (MQ=255)

                                                                                                                             GATTGGGGTGCGAAGCGCCCCCCAAGATT‑TGT‑GCCTTGATGCTGATATTGCCGGGGCCGCTGT‑GCAGCCCAAT‑CGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTGGCGGCCATTtc                         <  1:4988944‑M1/150‑1 (MQ=255)

                                                                                                                              ATTGGGGTGCGAAGCGCCCCCCAAGATT‑TGT‑GCCTTGATGCTGATATTGCCGGGGCCGCTGT‑GCAGCCCAAT‑CGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTGGCGGCCATTtct                        <  2:5896732‑M1/150‑1 (MQ=255)

                                                                                                                                tGGGGGGCGACGCGCCCCCAACTTTT‑TGT‑GCCTTTTTGCTTTTTGGGGCGGGGCGGCTGT‑GCCCCCCACG‑CGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTGGCGGCCATTTCTcc                      <  2:3515063/150‑1 (MQ=255)

                                                                                                                                  GGGTGCGAAGCGCCCCCCAAGATT‑TGT‑GCCTTGATGCTGATATTGCCGGGGCCGCTGT‑GCAGCCCAAT‑CGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTa                                         >  2:4955999‑M1/1‑129 (MQ=255)

                                                                                                                                  GGGTGCGAAGCGCCCCCCAAGATT‑TGT‑GCCTTGATGCTGATATTGCCGGGGCCGCTGT‑GCAGCCCAAT‑CGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTGGCGGCCATTTCTCCGt                    <  1:7510585‑M1/150‑1 (MQ=255)

                                                                                                                                  GGGTGCGAAGCGCCCCCCAAGATT‑TGT‑GCCTTGATGCTGATATTGCCGGGGCCGCTGT‑GCAGCCCAAT‑CGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTGGCGGCCATTTCTCCGt                    <  1:1123844‑M1/150‑1 (MQ=255)

                                                                                                                                   GGTGCGAAGCGCCCCCCAAGATT‑TGT‑GCCTTGATGCTGATATTGCCGGGGCCGCTGT‑GCAGCCCAAT‑CGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGgc                                                            >  1:966086‑M1/1‑109 (MQ=255)

                                                                                                                                   GGTGCGAAGCGCCCCCCAAGATT‑TGT‑GCCTTGATGCTGATATTGCCGGGGCCGCTGT‑GCAGCCCAAT‑CGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGgc                                                            >  1:5535799‑M1/1‑109 (MQ=255)

                                                                                                                                   GGTGCGAAGCGCCCCCCAAGATT‑TGT‑GCCTTGATGCTGATATTGCCGGGGCCGCTGT‑GCAGCCCAAT‑CGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTGGCGACCATTTCTCCgtg                   <  1:6388326‑M1/150‑2 (MQ=255)

                                                                                                                                      gcgccgcgcCCCCCATTTTT‑TGG‑GCGTTGTTGCTGTTTTTGCCGGGGCCGGGGC‑GCA‑CCCAATGCGCGCCGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTGGCGGCCATTTCTCCGTTTCa                <  2:1246737/145‑1 (MQ=255)

                                                                                                                                      gCGACGCCCCCCCCACTTTT‑GGT‑GCCTTGAGGCTTTTTGTGCCGGGGCCGCTGC‑GCAGCCAAAT‑CGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTGGCGGCCATTTCTCCGTTTCa                <  2:8253174/150‑1 (MQ=255)

                                                                                                                                        gccgcgcCCCCCCCGTTT‑TGG‑GCGTTGGTGGTGTTTGGGCGGGGGCGGCGGG‑GCCCCCCACG‑CGCGCCCCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTGGCGGCCATTTCTCCGTTTCAGt              <  2:1997041/147‑1 (MQ=255)

                                                                                                                                        gAAGCGCCCCCCAAGATT‑TGT‑GCCTTGATGCTGATATTGCCGGGGCCGCTGT‑GCAGCCCAAT‑CGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTGGCGGCCATTTCTCCGTTTCAGt              <  2:7754742/150‑1 (MQ=255)

                                                                                                                                          aGCGCCCCCCAAGATT‑TGT‑GCCTTGATGCTGATATTGCCGGGGCCGCTGT‑GCAGCCCAAT‑CGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTGGCGGCCATTTCTCCGTTTCAGTgg            >  1:2223774/1‑150 (MQ=255)

                                                                                                                                           gcgcCCCCCAAGATT‑TGT‑GCCTTGATGCTGATATTGCCGGGGCCGCTGT‑GCAGCCCAAT‑CGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTGGCGGCCATTTCTCCGTTTCAGTGGa           >  1:5700161/1‑150 (MQ=255)

                                                                                                                                               cccccAAGATT‑TGT‑GCCTTGATGCTGATATTGCCGGGGCCGCTGT‑GCAGCCCAAT‑CGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTGGCGGCCATTTCTCCGTTTCAGTGGACAGc       >  1:3009031/1‑150 (MQ=255)

                                                                                                                                                 cccAAGATT‑TGT‑GCCTTGATGCTGATATTGCCGGGGCCGCTGT‑GCAGCCCAAT‑CGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTGGCGGCCATTTCTCCGTTTCAGTGGACAGccc     >  2:3800146/1‑150 (MQ=255)

                                                                                                                                                  ccAAGATT‑TGT‑GCCTTGATGCTGATATTGCCGGGGCCGCTGT‑GCAGCCCAAT‑CGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTGGCGGCCATTTCTCCGTTTCAGTGGACAGcccc    >  1:1623476/1‑150 (MQ=255)

                                                                                                                                                  ccAAGATT‑TGT‑GCCTTGATGCTGATATTGCCGGGGCCGCTGT‑GCAGCCCAAT‑CGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTGGCGGCCATTTCTCCGTTTCAGTGGACAGcccc    >  2:7947618/1‑150 (MQ=255)

                                                                                                                                                  ccAAGATT‑TGT‑GCCTTGATGCTGATATTGCCGGGGCCGCTGT‑GCAGCCCAAT‑CGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTGGCGGCCATTTCTCCGTTTCAGTGGACAGcccc    >  1:2621194/1‑150 (MQ=255)

                                                                                                                                                  ccAAGATT‑TGT‑GCCTTGATGCTGATATTGCCGGGGCCGCTGT‑GCAGCCCAAT‑CGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTGGCGGCCATTTCTCCGTTTCAGTGGACAGcccc    >  1:6064937/1‑150 (MQ=255)

                                                                                                                                                   cAAGATT‑TGT‑GCCTTGATGCTGATATTGCCGGGGCCGCTGT‑GCAGCCCAAT‑CGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTGGCGGCCATTTCTCCGTTTCAGTGGACAGCCCCa   >  2:7383734/1‑150 (MQ=255)

                                                                                                                                                   cAAGATT‑TGT‑GCCTTGATGCTGATATTGCCGGGGCCGCTGT‑GCAGCCCAAT‑CGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTGGCGGCCATTTCTCCGTTTCAGTGGACAGCCCCa   >  1:8130960/1‑150 (MQ=255)

                                                                                                                                                   cAAGATT‑TGT‑GCCTTGATGCTGATATTGCCGGGGCCGCTGT‑GCAGCCCAAT‑CGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTGGCGGCCATTTCTCCGTTTCAGTGGACAGCCCCa   >  1:4836267/1‑150 (MQ=255)

                                                                                                                                                   cAAGATT‑TGT‑GCCTTGATGCTGATATTGCCGGGGCCGCTGT‑GCAGCCCAAT‑CGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTGGCGGCCATTTCTCCGTTTCAGTGGACAGCCCCa   >  1:4452997/1‑150 (MQ=255)

                                                                                                                                                    aaGATT‑TGT‑GCCTTGATGCTGATATTGCCGGGGCCGCTGT‑GCAGCCCAAT‑CGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTGGCGGCCATTTCTCCGTTTCAGTGGACAGCCCCAt  >  1:3641920/1‑150 (MQ=255)

ATTGGCCAGTCGCCAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATT‑TGT‑GCCTTGATGCTGATATTGCCGGGGCCGCTGT‑GCAGCCCAAT‑CGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTGGCGGCCATTTCTCCGTTTCAGTGGACAGCCCCAT  >  NC_002947/5970278‑5970575

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