BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,709,090	T→C	5.4%	L377P (CTG→CCG) ‡	PP_5009 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,709,090	0	T	C	5.4%	586.8 / 13.7	183	L377P (CTG→CCG) ‡	PP_5009	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (78/95); new base C (7/3); total (85/98)
Fisher's exact test for biased strand distribution p-value = 1.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.99e-01

Only 100 of 189 total aligned reads displayed.

CTTCGACCTGCCTTCGCAGGAACGAACTGCGCCCGCTGTCAGCCCACCGCGATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCATTTCCCAAGCCATCATCGGCGCCGGCACACAAGGGCAGT  >  NC_002947/5708941‑5709239

 ttcGACCTGCCTTCGCAGGAACGAACTGCGCCCGCTGTCAGCCCACCGCGATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTg                                                                                                                                                      >  2:2184368/1‑150 (MQ=255)

 ttcGACCTGCCTTCGCAGGAACGAACTGCGCCCGCTGTCAGCCCACCGCGATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTg                                                                                                                                                      >  2:2985958/1‑150 (MQ=255)

 ttcGACCTGCCTTCGCAGGAACGAACTGCGCCCGCTGTCAGCCCACCGCGATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTg                                                                                                                                                      >  2:1026058/1‑150 (MQ=255)

 ttcGACCTGCCTTCGCAGGAACGAACTGCGCCCGCTGTCAGCCCACCGCGATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTg                                                                                                                                                      >  2:5229211/1‑150 (MQ=255)

  tcGACCTGCCTTCGCAGGAACGAACTGCGCCCGCTGTCAGCCCACCGCGATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTgg                                                                                                                                                     <  2:4912148/150‑1 (MQ=255)

    gACCTGCCTTCGCAGGAACGAACTGCGCCCGCTGTCAGCCCACCGCGATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGc                                                                                                                                                   >  1:13467/1‑150 (MQ=255)

       cTGCCTTCGCAGGAACGAACTGCGCCCGCTGTCAGCCCACCGCGATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGAt                                                                                                                                                >  2:3472219/1‑150 (MQ=255)

         gCCTTCGCAGGAACGAACTGCGCCCGCTGTCAGCCCACCGCGATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCATACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGt                                                                                                                                              <  2:1456764/150‑1 (MQ=255)

          ccTTCGCAGGAACGAACTGCGCCCGCTGTCAGCCCACCGCGATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCCCAGCGATGtg                                                                                                                                             <  2:607212/150‑2 (MQ=255)

            ttCGCAGGAACGAACTGCGCCCGCTGTCAGCCCACCGCGATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGa                                                                                                                                           <  2:2718997/150‑1 (MQ=255)

             tCGCAGGAACGAACTGCGCCCGCTGTCAGCCCACCGCGATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGaa                                                                                                                                          <  1:4663171/150‑1 (MQ=255)

             tCGCAGGAACGAACTGCGCCCGCTGTCAGCCCACCGCGATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGaa                                                                                                                                          <  1:4739243/150‑1 (MQ=255)

             tCGCAGGAACGAACTGCGCCCGCTGTCAGCCCACCGCGATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGaa                                                                                                                                          <  2:2317577/150‑1 (MQ=255)

              cGCAGGAACGAACTGCGCCCGCTGTCAGCCCACCGCGATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAg                                                                                                                                         <  2:1176579/150‑1 (MQ=255)

                   gaacgaacTGCGCCCGCTGTCAGCCCACCGCGATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATc                                                                                                                                    >  2:2474892/1‑150 (MQ=255)

                    aacgaacTGCGCCCGCTGTCAGCCCACCGCGATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATcc                                                                                                                                   <  2:616110/150‑1 (MQ=255)

                      cgaacTGCGCCCGCTGTCAGCCCACCGCGATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCGCACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGt                                                                                                                                 <  1:815562/150‑1 (MQ=255)

                       gaacTGCGCCCGCTGTCAGCCCACCGCGATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTg                                                                                                                                >  2:4057265/1‑150 (MQ=255)

                          cTGCGCCCGCTGTCAGCCCACCGCGATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGAtt                                                                                                                             <  1:3472219/150‑1 (MQ=255)

                           tGCGCCCGCTGTCAGCCCACCGCGATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTa                                                                                                                            <  1:203737/150‑1 (MQ=255)

                              gcCCGCTGTCAGCCCACCGCGATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTAcgc                                                                                                                         <  1:2930379/150‑1 (MQ=255)

                              gcCCGCAGTCAGCCCACCGCGAGGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTAcgc                                                                                                                         <  2:5330055/150‑1 (MQ=255)

                                ccGCTGTCAGCCCACCGCGATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTAcgccg                                                                                                                       >  2:1610526/1‑150 (MQ=255)

                                ccGCTGTCAGCCCACCGCGATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTAcgccg                                                                                                                       >  1:1204740/1‑150 (MQ=255)

                                   cTGTCAGCCCACCGCGATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTAcgccgctc                                                                                                                    <  1:4614968/150‑1 (MQ=255)

                                     gTCAGCCCACCGCGATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTAcgccgctcgc                                                                                                                  >  1:4177068/1‑150 (MQ=255)

                                      tCAGCCCACCGCGATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTAcgccgctcgcc                                                                                                                 >  1:5225567/1‑150 (MQ=255)

                                          cccACCGCGATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCtt                                                                                                             <  1:2989599/150‑1 (MQ=255)

                                             aCCGCGATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCg                                                                                                          >  1:857680/1‑150 (MQ=255)

                                             aCCGCGATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCg                                                                                                          >  1:4176006/1‑150 (MQ=255)

                                              ccGCGATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGt                                                                                                         <  1:695545/150‑1 (MQ=255)

                                              ccGCGATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGt                                                                                                         <  1:1553432/150‑1 (MQ=255)

                                              ccGCGATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGt                                                                                                         >  1:4614361/1‑150 (MQ=255)

                                               cgcgATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTg                                                                                                        <  1:1020546/150‑1 (MQ=255)

                                                gcgATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTga                                                                                                       <  1:1009466/150‑1 (MQ=255)

                                                 cgATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTgag                                                                                                      >  2:1130784/1‑150 (MQ=255)

                                                  gATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTgaga                                                                                                     >  2:2119118/1‑150 (MQ=255)

                                                     ggCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGAccc                                                                                                  <  1:3428089/150‑1 (MQ=255)

                                                       ctgctCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTg                                                                                                >  1:3380727/1‑150 (MQ=255)

                                                       ctgctCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTg                                                                                                >  1:278254/1‑150 (MQ=255)

                                                        tgctCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGc                                                                                               >  1:2700360/1‑150 (MQ=255)

                                                         gctCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATcc                                                                                                                                   >  1:4298806/1‑113 (MQ=255)

                                                         gctCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCa                                                                                              <  1:4760576/150‑1 (MQ=255)

                                                         gctCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCa                                                                                              <  1:5367046/150‑1 (MQ=255)

                                                           tCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGa                                                                                            >  2:3576820/1‑150 (MQ=255)

                                                             aaTGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGAcc                                                                                          >  1:1774468/1‑150 (MQ=255)

                                                              aTGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCa                                                                                         <  1:3283434/150‑1 (MQ=255)

                                                              aTGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCa                                                                                         <  2:1003735/150‑1 (MQ=255)

                                                                  cAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTa                                                                                     >  1:1103688/1‑150 (MQ=255)

                                                                     cgcTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATcc                                                                                  <  1:210358/150‑1 (MQ=255)

                                                                      gcTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCa                                                                                 <  1:2119118/150‑1 (MQ=255)

                                                                      gcTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCGTGCGTGAGACCCTGCAGACCAGCTATCCa                                                                                 <  2:3337205/150‑1 (MQ=255)

                                                                            gACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAg                                                                           >  2:4353247/1‑150 (MQ=255)

                                                                             aCCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAgg                                                                          <  2:2436796/150‑1 (MQ=255)

                                                                               cGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCa                                                                        <  2:947557/150‑1 (MQ=255)

                                                                               cGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCa                                                                        <  1:3778546/150‑1 (MQ=255)

                                                                                  tACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTa                                                                     <  2:2166488/150‑1 (MQ=255)

                                                                                       cTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGa                                                                <  1:2184368/150‑1 (MQ=255)

                                                                                             cacGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGaa                                                          <  1:142478/150‑1 (MQ=255)

                                                                                                   caggaatgATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGAtgt                                                    <  2:2892740/143‑1 (MQ=255)

                                                                                                    cTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGAtgtg                                                   <  2:3687673/150‑1 (MQ=255)

                                                                                                    cTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGAtgtg                                                   <  2:4176006/150‑1 (MQ=255)

                                                                                                      gCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGca                                                 >  2:3939584/1‑150 (MQ=255)

                                                                                                       cTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGcac                                                >  2:321117/1‑150 (MQ=255)

                                                                                                        tCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGcaca                                               >  2:5386841/1‑150 (MQ=255)

                                                                                                           aTCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATcc                                            <  2:4616750/150‑1 (MQ=255)

                                                                                                           aTCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATcc                                            <  2:1774468/150‑1 (MQ=255)

                                                                                                              tGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCa                                         >  2:5259633/1‑150 (MQ=255)

                                                                                                              tGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCa                                         >  1:3644938/1‑150 (MQ=255)

                                                                                                               ggcggcGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCAt                                        >  1:664313/1‑150 (MQ=255)

                                                                                                               ggcggcGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCAt                                        <  2:1271292/150‑1 (MQ=255)

                                                                                                                   gcgATCCAAGTTCTGTAGAATGTCTGCAAATCTCTGGTCTTTGTCGAAGATATCCGTTATTACGAAGCTCGACTCTTGATTTAGACCCTGCAGAAAATCTATCCAGACAAGTCTGTACACGACCCTAACGATGTGCACATCCGCAtttaa                                    <  2:2803198/150‑3 (MQ=255)

                                                                                                                    cgAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCATTTccc                                   <  1:5259633/150‑1 (MQ=255)

                                                                                                                     gAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCATTTCCCa                                  <  1:321117/150‑1 (MQ=255)

                                                                                                                      aaCCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCATTTCCCaa                                 <  1:1638143/150‑1 (MQ=255)

                                                                                                                      aaCCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCATTTCCCaa                                 <  1:1881292/150‑1 (MQ=255)

                                                                                                                      aaCCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCATTTCCCaa                                 >  1:3186742/1‑150 (MQ=255)

                                                                                                                      aaCCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCATTTCCCaa                                 <  1:2305397/150‑1 (MQ=255)

                                                                                                                      aaCCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCATTTCCCaa                                 <  1:2474892/150‑1 (MQ=255)

                                                                                                                       aCCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCATTTCCCAAg                                >  2:4768709/1‑150 (MQ=255)

                                                                                                                          aaGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCATTTCCCAAGCca                             >  2:2065967/1‑150 (MQ=255)

                                                                                                                             gtcggtcgAAGGGCTGCACATAGCTGGGCGATGTCGAAGATATCCGTGATTACTCCGCTCGCCGCTTGCGTGTGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCATTTCACAATCcatca                          >  2:3878127/1‑150 (MQ=255)

                                                                                                                               cggtcgAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCATTTCCCAAGCGatcatc                        <  2:4177068/150‑1 (MQ=255)

                                                                                                                                ggtcgAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCATTTCCCAAGCCATCATCg                       >  2:4944967/1‑150 (MQ=255)

                                                                                                                                    gAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCATTTCCCAAGCcat                            >  1:1690025/1‑141 (MQ=255)

                                                                                                                                    gAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCATTTCCCAAGCcat                            <  2:1690025/141‑1 (MQ=255)

                                                                                                                                    gAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCATTTCCCAAGCCATCATCGgcgc                   >  2:312362/1‑150 (MQ=255)

                                                                                                                                     aaGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCATTTCCCAAGCCATCATCGGCGcc                  >  1:3412425/1‑150 (MQ=255)

                                                                                                                                      aGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGTACGATGTGCACATCCGCATTTCCCAAGCCATCATCGGCGCCg                 <  1:3697276/150‑1 (MQ=255)

                                                                                                                                          ttcgacaTCGCCCAGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCATTTCCCaa                                 >  1:5175336/5‑130 (MQ=255)

                                                                                                                                          ttcgacaTCGCCCAGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCATTTCCCAAGCCATCATCGGCGCCGGcac             >  2:3512817/5‑150 (MQ=255)

                                                                                                                                          ttcgacaTCGCCCAGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCATTTCCCAAGCCATCATCGGCGCCGGcac             >  2:3417383/5‑150 (MQ=255)

                                                                                                                                          ttcgacaTCGCCCAGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCATTTCCCAAGCCATCATCGGCGCCGGcac             >  1:2724521/5‑150 (MQ=255)

                                                                                                                                                 tCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCATTTCCCAAGCCATCATCGGCGCCGGCACACAAggg      <  2:799510/150‑1 (MQ=255)

                                                                                                                                                   gCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCATTTCCCAAGCCATCATCGGCGCCGGcacac           <  2:160747/143‑1 (MQ=255)

                                                                                                                                                   gCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCATTTCCCAAGCCATCATCGGCGCCGGCACACAAGGGCa    <  2:4644973/150‑1 (MQ=255)

                                                                                                                                                   gCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCATTTCCCAAGCCATCATCGGCGCCGGCACACAAGGGCa    <  2:90230/150‑1 (MQ=255)

                                                                                                                                                     tGGGCGATGTCGAAGATATCCGTGTTTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCATTTCCCAAGCCATCATCGGCGCCGGCACACAAGGGCAGt  <  2:1735235/150‑1 (MQ=255)

                                                                                                                                                     tGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCATTTCCCAAGCCATCATCGGCGCCGGCACACAAGGGCAGt  <  1:2964491/150‑1 (MQ=255)

                                                                                                                                                     tGGGCGATGTCGAAGATATCCGTGATTACGCCGCACGCCGCTTGCGTGAGACCCTCCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACCATGTGCCCATCCGCATTTCCCAAGCCATCAACGCCGCCGCCACACAAGGGCAga  >  2:5190462/1‑149 (MQ=255)

CTTCGACCTGCCTTCGCAGGAACGAACTGCGCCCGCTGTCAGCCCACCGCGATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCATTTCCCAAGCCATCATCGGCGCCGGCACACAAGGGCAGT  >  NC_002947/5708941‑5709239

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