BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP009974	1,444,631	C→T	D26N (GAC→AAC)	RPPX_06380 ←	NUDIX hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009974	1,444,631	0	C	T	93.3%	846.6 / 28.0	255	D26N (GAC→AAC)	RPPX_06380	NUDIX hydrolase
Reads supporting (aligned to +/- strand): ref base C (11/6); new base T (129/109); total (140/115)
Fisher's exact test for biased strand distribution p-value = 4.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

Only 100 of 260 total aligned reads displayed.

TCCCAATACCCGGGGTACAGCGCCTTGCTCAGCGTACGCCGGTGTACACACAATTCACCGGCACTGTTGAACAGCAGAATGAACGTGCAGCGACCGATCAGGCCACGCTCGCGCAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTCGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGGGCGGCCTCGGTTGGGCTGATAGCCATCATCAACCCTGCGACAGCAGCTGACGCAGATCGATAACCGCAGCGTTGGCCCGGGAAATGTAGTTGGCCATG  >  CP009974/1444483‑1444778

   cAATACCCGGGGTACAGCGCCTTGCTCAGCGTACGCCGGTGTACACACAATTCACCGGCACTGTTGAACAGCAGAATGAACGTGCAGCGACCGATCAGGCCACGCTCGCGCAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTTGGttt                                                                                                                                                <  1:4117537/151‑1 (MQ=255)

       aCCCGGGGTACAGCGCCTTGCTCAGCGTACGCCGGTGTACACACAATTCACCGGCACTGTTGAACAGCAGAATGAACGTGCAGCGACCGATCAGGCCACGCTCGCGCAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTc                                                                                                                                            <  2:5170210/151‑1 (MQ=255)

        cccGGGGTACAGCGCCTTGCTCAGCGTACGCCGGTGTACACACAATTCACCGGCACTGTTGAACAGCAGAATGAACGTGCAGCGACCGATCAGGCCACGCTCGCGCAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTcc                                                                                                                                           >  2:4512807/1‑151 (MQ=255)

          cGGGGTACAGCGCCTTGCTCAGCGTACGCCGGTGTACACACAATTCACCGGCACTGTTGAACAGCAGAATGAACGTGCAGCGACCGATCAGGCCACGCTCGCGCAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCAc                                                                                                                                         >  2:2206628/1‑151 (MQ=255)

           ggggTACAGCGCCTTGCTCAGCGTACGCCGGTGTACACACAATTCACCGGCACTGTTGAACAGCAGAATGAACGTGCAGCGACCGATCAGGCCACGCTCGCGCAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCAcc                                                                                                                                        >  1:1476749/1‑151 (MQ=255)

              gTACAGCGCCTTGCTCAGCGTACGCCGGTGTACACACAATTCACCGGCACTGTTGAACAGCAGAATGAACGTGCAGCGACCGATCAGGCCACGCTCGCGCAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAg                                                                                                                                     <  2:3985090/151‑1 (MQ=255)

               tACAGCGCCTTGCTCAGCGTACGCCGGTGTACACACAATTCACCGGCACTGTTGAACAGCAGAATGAACGTGCAGCGACCGATCAGGCCACGCTCGCGCAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAgg                                                                                                                                    >  2:2107312/1‑151 (MQ=255)

                 cAGCGCCTTGCTCAGCGTACGCCGGTGTACACACAATTCACCGGCACTGTTGAACAGCAGAATGAACGTGCAGCGACCGATCAGGCCACGCTCGCGCAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTg                                                                                                                                  >  1:511006/1‑151 (MQ=255)

                  aGCGCCTTGCTCAGCGTACGCCGGTGTACACACAATTCACCGGCACTGTTGAACAGCAGAATGAACGTGCAGCGACCGATCAGGCCACGCTCGCGCAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTgc                                                                                                                                 >  2:3804040/1‑150 (MQ=255)

                   gcgcCTTGCTCAGCGTACGCCGGTGTACACACAATTCACCGGCACTGTTGAACAGCAGAATGAACGTGCAGCGACCGATCAGGCCACGCTCGCGCAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGAc                                                                                                                                >  2:1111475/1‑151 (MQ=255)

                             cAGCGTACGCCGGTGTACACACAATTCACCGGCACTGTTGAACAGCAGAATGAACGTGCAGCGACCGATCAGGCCACGCTCGCGCAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCgg                                                                                                                      >  1:2804513/1‑151 (MQ=255)

                             cAGCGTACGCCGGTGTACACACAATTCACCGGCACTGTTGAACAGCAGAATGAACGTGCAGCGACCGATCAGGCCACGCTCGCGCAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCgg                                                                                                                      <  1:2054335/151‑1 (MQ=255)

                                cGTACGCCGGTGTACACACAATTCACCGGCACTGTTGAACAGCAGAATGAACGTGCAGCGACCGATCAGGCCACGCTCGCGCAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCg                                                                                                                   <  2:678900/151‑1 (MQ=255)

                                cGTACGCCGGTGTACACACAATTCACCGGCACTGTTGAACAGCAGAATGAACGTGCAGCGACCGATCAGGCCACGCTCGCGCAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCg                                                                                                                   >  1:4230489/1‑151 (MQ=255)

                                 gTACGCCGGCGTACACACAATTCACCGGCCCTGTCGAACAGCAGAATGAACGTGCAGCGACCGATCCGGCCACGCTCGCGCAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTTGGCTTCGTCCACCCAGGTGCCCAGCCCGCCGTCgc                                                                                                                  >  2:5065808/1‑150 (MQ=255)

                                    cGCCGGTGTACACACAATTCACCGGCACTGTTGAACAGCAGAATGAACGTGCAGCGACCGATCAGGCCACGCTCGCGCAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAg                                                                                                               >  1:3962247/1‑151 (MQ=255)

                                      ccGGTGTACACACAATTCACCGGCACTGTTGAACAGCAGAATGAACGTGCAGCGACCGATCAGGCCACGCTCGCGCAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAgcc                                                                                                             >  2:1387346/1‑151 (MQ=255)

                                      ccGGTGTACACACAATTCACCGGCACTGTTGAACAGCAGAATGAACGTGCAGCGACCGATCAGGCCACGCTCGCGCAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAgcc                                                                                                             >  1:613914/1‑151 (MQ=255)

                                      ccGGTGTACACACAAGTCACCGGCACTGTTGAACAGCAGAATGAACGTGCAGCGACCGATCAGGCCACGCTCGCGCAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAgcc                                                                                                             <  2:1266014/151‑1 (MQ=255)

                                        ggTGTACACACAATTCACCGGCACTGTTGAACAGCAGAATGAACGTGCAGCGACCGATCAGGCCACGCTCGCGCAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAgccgc                                                                                                           >  1:703390/1‑151 (MQ=255)

                                          tgtACACACAATTCACCGGCACTGTTGAACAGCAGAATGAACGTGCAGCGACCGATCAGGCCACGCTCGCGCAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGccc                                                                                                         <  1:2831087/151‑1 (MQ=255)

                                             acacacAATTCACCGGCACTGTTGAACAGCAGAATGAACGTGCAGCGACCGATCAGGCCACGCTCGCGCAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGAt                                                                                                      >  1:2897986/1‑151 (MQ=255)

                                               acacaATTCACCGGCACTGTTGAACAGCAGAATGAACGTGCAGCGACCGATCAGGCCACGCTCGCGCAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATgg                                                                                                    <  1:3682076/151‑1 (MQ=255)

                                                cacaATTCACCGGCACTGTTGAACAGCAGAATGAACGTGCAGCGACCGATCAGGCCACGCTCGCGCAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATggg                                                                                                   <  1:4196689/151‑1 (MQ=255)

                                                cacaATTCACCGGCACTGTTGAACAGCAGAATGAACGTGCAGCGACCGATCAGGCCACGCTCGCGCAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATggg                                                                                                   >  2:1235128/1‑151 (MQ=255)

                                                 acaATTCACCGGCACTGTTGAACAGCAGAATGAACGTGCAGCGACCGATCAGGCCACGCTCGCGCAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGCTGggc                                                                                                  >  2:1979017/1‑151 (MQ=255)

                                                  caATTCACCGGCACTGTTGAACAGCAGAATGAACGTGCAGCGACCGATCAGGCCACGCTCGCGCAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGggcg                                                                                                 <  2:471336/151‑1 (MQ=255)

                                                  caATTCACCGGCACTGTTGAACAGCAGAATGAACGTGCAGCGACCGATCAGGCCACGCTCGCGCAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTCGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGggcg                                                                                                 >  2:2805258/1‑151 (MQ=255)

                                                   agttgagcGGCGCTGTGGGACAGCAGGATGAGCGGGCAGCGACCGGTCAGGCCACGCGCGCGCAGCTCGGCACGGGGCAAGGCGCCCAGTACCGGGTTGGTTTCGTCCACCCAGGTGACCAGATCGCGGTCGGAAGGCGCCCGATGggcgg                                                                                                <  2:1715638/144‑1 (MQ=255)

                                                     ttCACCGGCACTGTTGAACAGCAGAATGAACGTGCAGCGACCGATCAGGCCACGCTCGCGCAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTCGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGGGCGGcc                                                                                              >  1:2724621/1‑151 (MQ=255)

                                                      tCACCGGCACTGTTGAACAGCAGAATGAACGTGCAGCGACCGATCAGGCCACGCTCGCGCAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGGGCGGCCt                                                                                             <  1:891798/151‑1 (MQ=255)

                                                      tCACCGGCACTGTTGAACAGCAGAATGAACGTGCAGCGACCGATCAGGCCACGCTCGCGCAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGGGCGGCCt                                                                                             <  1:1147623/151‑1 (MQ=255)

                                                      ggggcGGCACTGTGGAACGGGGGAGTGAAGGTGCAGCGAGCGATCGGGCGACGCTCGCGCAGCTCGGGACGGGGGAGGGCGCGCAGTACCTGGTCGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGGGCGGCCt                                                                                             <  2:361531/147‑1 (MQ=255)

                                                       cACCGGCACTGTTGAACAGCAGAATGAACGTGCAGCGACCGATCAGTCCACGCTCGCGCAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGGGCGGCCTc                                                                                            <  2:3728365/151‑1 (MQ=255)

                                                       cACCGGCACTGTTGAACAGCAGAATGAACGTGCAGCGACCGATCAGGCCACGCTCGCGCAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGGGCGGCCTc                                                                                            <  1:3750064/151‑1 (MQ=255)

                                                       cACCGGCACTGTTGAACAGCAGAATGAACGTGCAGCGACCGATCAGGCCACGCTCGCGCAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGGGCGGCCTc                                                                                            <  1:530955/151‑1 (MQ=255)

                                                        aCCGGCACTGTTGAACAGCAGAATGAACGTGCAGCGACCGATCAGGCCACGCTCGCGCAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGGGCGGCCTCg                                                                                           >  1:1347233/1‑151 (MQ=255)

                                                         ccGGCACTGTTGAACAGCAGAATGAACGTGCAGCGACCGATCAGGCCACGCTCGCGCAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGGGCGGCCTCgg                                                                                          <  2:3715401/151‑1 (MQ=255)

                                                            gCACTGGTGAACAGCAGAATGAACGTGCAGCGACCGATCAGGCCACGCTCGCGCAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGGGCGGCCTCGGTTg                                                                                       <  1:4618436/151‑1 (MQ=255)

                                                             cACTGTTGAACAGCAGAATGAACGTGCAGCGACCGATCAGGCCACGCTCGCGCAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGGGCGGCCTCGGTTgg                                                                                      >  2:656792/1‑151 (MQ=255)

                                                              aCTGTTGAACAGCAGAATGAACGTGCAGCGACCGATCAGGCCACGCTCGCGCAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGGGCGGCCTCGGTTggg                                                                                     >  2:2105518/1‑151 (MQ=255)

                                                              aCTGTTGAACAGCAGAATGAACGTGCAGCGACCGATCAGGCCACGCTCGCGCAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTCGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGGGCGGCCTCGGTTggg                                                                                     >  1:2028072/1‑151 (MQ=255)

                                                                 gTTGAACAGCAGAATGAACGTGCAGCGACCGATCAGGCCACGCTCGCGCAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGGGCGGCCTCGGTTGGGCTg                                                                                  >  1:3806142/1‑151 (MQ=255)

                                                                  ttGAACAGCAGAATGAACGTGCAGCGACCGATCAGGCCACGCTCGCGCAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGGGCGGCCTCGGTTGGGCTGa                                                                                 >  1:872167/1‑151 (MQ=255)

                                                                    taaCAGCAGAATGAACGTGCAGCGACCGATCAGGCCACGCTCGCGCAGCTCGGCACGGCGCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCCCGCGGTCGGAAGCCGCCCGATGGGCGGCCTCGGTTGGGCTGATa                                                                               <  1:1875472/150‑1 (MQ=255)

                                                                     aaCAGCAGAATGAACGTGCAGCGACCGATCAGGCCACGCTCGCGCAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGGGCGGCCTCGGTTGGGCTGATAg                                                                              >  2:114095/1‑151 (MQ=255)

                                                                     aaCAGCAGAATGAACGTGCAGCGACCGATCAGGCCACGCTCGCGCAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGGGCGGCCTCGGTTGGGCTGATAg                                                                              <  1:1896625/151‑1 (MQ=255)

                                                                     aaCAGCAGAATGAACGTGCAGCGACCGATCAGGCCACGCTCGCGCAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGGGCGGCCTCGGTTGGGCTGATAg                                                                              >  2:113916/1‑151 (MQ=255)

                                                                      aCAGCAGAATGAACGTGCAGCGACCGATCAGGCCACGCTCGCGCAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGGGCGGCCTCGGTTGGGCTGATAGc                                                                             <  1:2123438/151‑1 (MQ=255)

                                                                      aCAGCAGAATGAACGTGCAGCGACCGATCAGGCCACGCTCGCGCAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGGGCGGCCTCGGTTGGGCTGATAGc                                                                             >  2:1509606/1‑151 (MQ=255)

                                                                        agcagAATGAACGTGCAGCGACCGATCAGGCCACGCTCGCGCAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGGGCGGCCTCGGTCGGGCTGATAGCca                                                                           >  2:5131838/1‑151 (MQ=255)

                                                                             aaTGAACGTGCAGCGACCGATCAGGCCACGCTCGCGCAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGGGCGGCCTCGGTTGGGCTGATAGCCatcatc                                                                      <  1:3868245/151‑1 (MQ=255)

                                                                             aaTGAACGTGCAGCGACCGATCAGGCCACGCTCGCGCAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGGGCGGCCTCGGTTGGGCTGATAGCCatcatc                                                                      >  1:1678305/1‑151 (MQ=255)

                                                                             aaTGAACGTGCAGCGACCGATCAGGCCACGCTCGCGCAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGGGCGGCCTCGGTTGGGCTGATAGCCatcatc                                                                      <  1:1773115/151‑1 (MQ=255)

                                                                                   cGTGCAGCGACCGATCAGGCCACGCTCGCGCAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGGGCGGCCTCGGTTGGGCTGATAGCCATCATCAACCCt                                                                <  1:5041123/151‑1 (MQ=255)

                                                                                     ggCAGCGACCGATCAGGCCACGCTCGCGCAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGGGCGGCCTCGGTTGGGCTGATAGCCATCATCAACCCTGc                                                              <  1:459659/150‑1 (MQ=255)

                                                                                       cAGCGACCGATCAGGCCACGCTCGCGCAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGGGCGGCCTCGGTTGGGCTGATAGCCATCATCAACCCTGCGa                                                            >  1:3272900/1‑151 (MQ=255)

                                                                                         gCGACCGATCAGGCCACGCTCGCGCAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGGGCGGCCTCGGTTGGGCTGATCGCCATCATCAACCCTGCGAca                                                          >  2:3957522/1‑151 (MQ=255)

                                                                                               gATCAGGCCACGCTCGCGCAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGGGCGGCCTCGGTTGGGCTGATAGCCATCATCAACCCTGCGACAGCAGCt                                                    >  1:2610482/1‑151 (MQ=255)

                                                                                               gATCAGGCCACGCTCGCGCAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGGGCGGCCTCGGTTGGGCTGATAGCCATCATCAACCCTCCGACagcagcc                                                    >  2:3299053/1‑150 (MQ=255)

                                                                                                aTCAGGCCACGCTCGCGCAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGGGCGGCCTCGGTTGGGCTGATAGCCATCATCAACCCTGCGACAGCAGCTg                                                   >  1:3398690/1‑151 (MQ=255)

                                                                                                 tCAGGCCACGCTCGCGCAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGGGCGGCCTCGGTTGGGCTGATAGCCATCATCAACCCTGCGACAGCAGCTGa                                                  <  1:3427467/151‑1 (MQ=255)

                                                                                                  cAGGCCACGCTCGCGCAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTCGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGGGCGGCCTCGGTTGGGCTGATAGCCATCATCAACCCTGCGACAGCAGCTGAc                                                 >  1:4823591/1‑151 (MQ=255)

                                                                                                   aGGCCACGCTCGCGCGGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTCGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGGGCGGCCTCGGTTGGGCTGATAGCCATCCTCAACCCTCCGCCAGCAGCTGACg                                                >  2:576112/1‑151 (MQ=255)

                                                                                                   aGGCCACGCTCGCGCAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTCGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGGGCGGCCTCGGTTGGGCTGATAGCCATCATCAACCCTGCGACAGCAGCTGACg                                                <  2:2492436/151‑1 (MQ=255)

                                                                                                     gCCACGCTCGCGCAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGGGCGGCCTCGGTTGGGCTGATAGCCATCATCAACCCTGCGACAGCAGCTGACGCa                                              <  1:1011257/151‑1 (MQ=255)

                                                                                                      ccACGCTCGCGCAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGGGCGGCCTCGGTTGGGCTGATAGCCATCATCAACCCTGCGACAGCAGCTGACGCAg                                             <  1:3597528/151‑1 (MQ=255)

                                                                                                            tCGCGCAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGGGCGGCCTCGGTTGGGCTGATAGCCATCATCAACCCTGCGACCGCAGCTGACGCAgcccgcc                                       >  1:4304012/1‑145 (MQ=255)

                                                                                                            tCGCGCAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGGGCGGCCTCGGTTGGGCTGATAGCCATCATCAACCCTGCGACAGCAGCTGACGCAGATCGAt                                       <  1:2386316/151‑1 (MQ=255)

                                                                                                             cgcgCAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGCTGCGCGGCCTCGGTTGGGCTGATAGCCATCATCAACCCGGCAACAGCAGCTGACGCAGATCGATa                                      >  2:2925589/1‑151 (MQ=255)

                                                                                                             cgcgCAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGGGCGGCCTCGGTTGGGCTGATAGCCATCATCAACCCTGCGACAGCAGCTGACGCAGATCGATa                                      >  2:1781282/1‑151 (MQ=255)

                                                                                                              gcgcAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGGGCGGCCTCGGTTGGGCTGATAGCCATCATCAACCCTGCGACAGCAGCTGACGCAGATCGATaa                                     <  1:1227658/151‑1 (MQ=255)

                                                                                                               cgcAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGGGCGGCCTCGGTTGGGCTGATAGCCATCATCAACCCTGCGACAGCAGCTGACGCAGATCGATAAc                                    >  1:2299972/1‑151 (MQ=255)

                                                                                                               cgcAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTCGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGGGCGGCCTCGGTTGGGCTGATAGCCATCATCAACCCTg                                                               >  1:3335678/1‑124 (MQ=255)

                                                                                                               cgcAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTCGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGGGCGGCCTCGGTTGGGCTGATAGCCATCATCAACCCTg                                                               <  2:3335678/124‑1 (MQ=255)

                                                                                                                gcAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGGGCGGCCTCGGTTGGGCTGATAGCCATCATCAACCCTGCGACAGCAGCTGACGCAGATCGATAAcc                                   >  2:594245/1‑151 (MQ=255)

                                                                                                                 cAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGGGCGGCCTCGGTTGGGCTGATAGCCATCATCAACCCTGCGACAGCAGCTGACGCAGATCGATAACCg                                  <  2:2051828/151‑1 (MQ=255)

                                                                                                                 cAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGGGCGGCCTCGGTTGGGCTGATAGCCATCATCAACCCTGCGACAGCAGCTGACGCAGATCGATAACCg                                  <  2:4182824/151‑1 (MQ=255)

                                                                                                                  aGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGGGCGGCCTCGGTTGGGCTGATAGCCATCATCAACCCTGCGACAGCAGCTGACGCAGATCGATAACCGc                                 >  2:3157653/1‑151 (MQ=255)

                                                                                                                    cTCGGCACGGGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGGGCGGCCTCGGTTGGGCTGATAGCCATCATCAACCCTGCGACAGCAGCTGACGCAGATCGATAACCGCAg                               >  2:840517/1‑151 (MQ=255)

                                                                                                                    cTCGGCACGAGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGGGCGGCCTCGGTTGGGCTGATAGCCATCATCAACCCTGCGACAGCAGCTGACGCAGATCGATAACCGCAg                               >  1:3130489/1‑151 (MQ=255)

                                                                                                                      cGGCACGGGGCAAGTCGCCGTGTCCCTGGTTGGTTGCGTCGACCCAGGGGACCAGCTCGCGGTCGGAAGCCGCCCGATGGGCGGCCTCGGTTGGGCTGATAGCCATCATCAACCCTGCGACAGCAGCTGACGCAGATCGATAACCGCAGCg                             <  1:583714/151‑1 (MQ=255)

                                                                                                                         cACGGTGCAAGGCGCCCAGTACCTGGTCGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGGGCGGCCTCGGTTGGGCTGATAGCCATCATCAACCCTGCGACAGCAGCTGACGCAGATCGATAACCGCAGCGTTg                          <  2:268269/151‑1 (MQ=255)

                                                                                                                         cACGGGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGGGCGGCCTCGGTTGGGCTGATAGCCATCATCAACCCTGCGACAGCAGCTGACGCAGATCGATAACCGCAGCGTTg                          <  1:2918433/151‑1 (MQ=255)

                                                                                                                          aCGGGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGGGCGGCCTCGGTTGGGCTGATAGCCATCATCAACCCTGCGACAGCAGCTGACGCAGATCGATAACCGCAGCGTTgg                         >  1:3828366/1‑151 (MQ=255)

                                                                                                                           cGGGGCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGGGCGGCCTCGGTTGGGCTGATAGCCATCATCAACCCTGCGACAGCAGCTGACGCAGATCGATAACCGCAGCGTTGGc                        >  1:4069119/1‑151 (MQ=255)

                                                                                                                            ggggCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGGGCGGCCTCGGTTGGGCTGATAGCCATCATCAACCCTGCGACAGCAGCTGACGCAGATCGATAACCGCAGCGTTCGcc                       >  2:3635710/1‑151 (MQ=255)

                                                                                                                             gggCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGGGCGGCCTCGGTTGGGCTGATAGCCATCATCAACCCTGCGACAGCAGCTGACGCAGATCGATAACCGCAGCGTTGGccc                      <  1:2193484/151‑1 (MQ=255)

                                                                                                                              ggCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGGGCGGCCTCGGTTGGGCTGATAGCCATCATCAACCCTGCGACAGCAGCTGACGCAGATCGATAACCGCAGCGTTGGCCCg                     <  1:1702791/151‑1 (MQ=255)

                                                                                                                              ggCAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGGGCGGCCTCGGTTGGGCTGATAGCCATCATCAACCCTGCGACAGCAGCTGACGCAGATCGATAACCGCAGCGTTGGCCCg                     >  1:1307576/1‑151 (MQ=255)

                                                                                                                                cAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGGGCGGCCTCGGTTGGGCTGATAGCCATCATCAACCCTGCGACAGCAGCTGACGCAGATCGATAACCGCAGCGTTGGCCCggg                   <  2:491354/151‑1 (MQ=255)

                                                                                                                                cAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGGGCGGCCTCGGTTGGGCTGATAGCCATCATCAACCCTGCGACAGCAGCTGACGCAGATCGATAACCGCAGCGTTGGCCCggg                   >  2:3046867/1‑151 (MQ=255)

                                                                                                                                cAAGGCGCCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGGGCGGCCTCGGTTGGGCTGATAGCCATCATCAACCCTGCGACAGCAGCTGACGCAGATCGATAACCGCAGCGTTGGCCCggg                   <  2:3526572/151‑1 (MQ=255)

                                                                                                                                 atagtgcCCAGTATCTTGTTGGTTTCGTACACCCAGGTGACCAGCTTGCGGACGGTAGCCGCCCGATGGGCGGCCTTGGTTGGGCTGATAGTCATCATTAACTCTGCGACAGCAGCTGACTCAGTTTGATAATCGCTGAGTTGGACAgggt                  >  2:1742415/6‑150 (MQ=255)

                                                                                                                                     cgcCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGGGCGGCCTCGGTTGGGCTGATAGCCATCATCAACCCTGCGACAGCAGCTGACGCAGATCGATAACCGCAGCGTTGGCCCGGGAAATg              >  2:4436394/1‑151 (MQ=255)

                                                                                                                                     cgcCCAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGGGCGGCCTCGGTTGGGCTGATAGCCATCATCAACCCTGCGACAGCAGCTGACGCAGATCGATAACCGCAGCGTTGGCCCGGGAAATg              <  1:1082255/151‑1 (MQ=255)

                                                                                                                                        ccAGTACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGGGCGGCCTCGGTTGGGCTGATAGCCATCATCAACCCTGCGACAGCAGCTGACGCAGATCGATAACCGCAGCGTTGGCCCGGGAAATGTAg           >  1:3502364/1‑151 (MQ=255)

                                                                                                                                            tACCTGGTTGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGGGCGGCCTCGGTTGGGCTGATAGCCATCATCAACCCTGCGACAGCAGCTGACGCAGATCGATAACCGCAGCGTTGGCCCGGGAAATGTAGTTgg       <  1:711597/151‑1 (MQ=255)

                                                                                                                                                tggttggtTTCGTCCACCTAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGGGGGGCCTCGGTTGGGCTGATAGCCATCATCAACCCTGCGACAGCAGCTGACGCAGATCGATAACCGCAGCGTTGGCCCGGGAAATGTAGTTGGCCAt   <  2:1574936/151‑1 (MQ=255)

                                                                                                                                                 ggttggttTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGGGCGGCCTCGGTTGGGCTGATAGCCATCATCAACCCTGCGACAGCAGCTGACGCAGATCGATAACCGCAGCGTTGGCCCGGGAAATGTAGTTGGCCATg  >  2:2430811/1‑151 (MQ=255)

TCCCAATACCCGGGGTACAGCGCCTTGCTCAGCGTACGCCGGTGTACACACAATTCACCGGCACTGTTGAACAGCAGAATGAACGTGCAGCGACCGATCAGGCCACGCTCGCGCAGCTCGGCACGGGGCAAGGCGCCCAGTACCTGGTCGGTTTCGTCCACCCAGGTGACCAGCTCGCGGTCGGAAGCCGCCCGATGGGCGGCCTCGGTTGGGCTGATAGCCATCATCAACCCTGCGACAGCAGCTGACGCAGATCGATAACCGCAGCGTTGGCCCGGGAAATGTAGTTGGCCATG  >  CP009974/1444483‑1444778

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