BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,546,717	T→A	6.0%	I327F (ATC→TTC)	dnaB ←	replicative DNA helicase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,546,717	0	T	A	6.0%	1355.5 / 50.6	413	I327F (ATC→TTC)	dnaB	replicative DNA helicase
Reads supporting (aligned to +/- strand): ref base T (205/183); new base A (14/11); total (219/194)
Fisher's exact test for biased strand distribution p-value = 8.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.04e-01

Only 100 of 419 total aligned reads displayed.

CAGCTGCGAAAGGGCAATGACCGGGCAGTTGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAGTCGTCGTCGTCCAGCTGGCC  >  NC_002947/5546568‑5546864

cAGCTGCGAAAGGGCAATGACCGGGCAGTTGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGAt                                                                                                                                                     >  1:4700212/1‑150 (MQ=255)

   cTGCGAAAGGGCAATGACCGGGCAGTTGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTc                                                                                                                                                  >  1:7139325/1‑150 (MQ=255)

     gCGAAAGGGCAATGACCGGGCAGTTGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGc                                                                                                                                                >  2:5902644/1‑150 (MQ=255)

     gCGAAAGGGCAATGACCGGGCAGTTGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGc                                                                                                                                                <  1:2148084/150‑1 (MQ=255)

      cGAAAGGGCAATGACCGGGCAGTTGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGcc                                                                                                                                               >  1:4972444/1‑150 (MQ=255)

      cGAAAGGGCAATGACCGGGCAGTTGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGAAATCGcc                                                                                                                                               <  2:7473443/150‑1 (MQ=255)

        aaaGGGCAATGACCGGGCAGTTGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCAt                                                                                                                                             >  2:7422414/1‑150 (MQ=255)

        aaaGGGCAATGACCGGGCAGTTGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGAAATCGCCAt                                                                                                                                             <  1:2335547/150‑1 (MQ=255)

          agtgCAATGACCGGGCAGTTGAATTCTTTGTCCAGGGCCTTGAGCGATCGGGAAATCTCGTAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCTTGATCATGGCGAAATCGCCATga                                                                                                                                           <  2:2307538/147‑2 (MQ=255)

             gCAATGACCGGGCAGTTGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGAAATCGCCATGATCa                                                                                                                                        <  1:7547700/150‑1 (MQ=255)

              cAATGACCGGGCAGTTGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCAc                                                                                                                                       >  2:5505731/1‑150 (MQ=255)

              cAATGACCGGGCAGTTGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCAc                                                                                                                                       >  2:8367920/1‑150 (MQ=255)

               aaTGACCGGGCAGTTGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACg                                                                                                                                      >  2:1530485/1‑150 (MQ=255)

               aaTGACCGGGCAGTTGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACg                                                                                                                                      >  2:5623500/1‑150 (MQ=255)

               aaTGACCGGGCAGTTGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGAGCATGGCGAAATCGCCATGATCatg                                                                                                                                      <  2:2248073/150‑3 (MQ=255)

                       ggCAGTTGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAggc                                                                                                                              >  2:5730876/1‑150 (MQ=255)

                        gCAGTTGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAggcg                                                                                                                             <  1:7813608/150‑1 (MQ=255)

                            ttGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGcggcgg                                                                                                                         <  2:5340356/150‑1 (MQ=255)

                            ttGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGcggcgg                                                                                                                         <  2:407498/150‑1 (MQ=255)

                            ttGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGcggcgg                                                                                                                         <  2:1730077/150‑1 (MQ=255)

                            ttGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGcggcgg                                                                                                                         >  1:8036154/1‑150 (MQ=255)

                            ttGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGcggcgg                                                                                                                         >  1:6225941/1‑150 (MQ=255)

                              gAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGt                                                                                                                       >  2:2587487/1‑150 (MQ=255)

                              gAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCGCACGGGCCAGGCGGCGGGt                                                                                                                       >  2:6347998/1‑150 (MQ=255)

                                aTTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTAc                                                                                                                     >  1:5882256/1‑150 (MQ=255)

                                  tCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTAcgc                                                                                                                   <  2:2405125/150‑1 (MQ=255)

                                    tttGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTAcgcgc                                                                                                                 <  1:5794426/150‑1 (MQ=255)

                                      ttgCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTAcgcgcgc                                                                                                               <  2:2819139/148‑1 (MQ=255)

                                             ggCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTc                                                                                                        >  1:8466585/1‑150 (MQ=255)

                                              gCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTcc                                                                                                       <  2:5936523/150‑1 (MQ=255)

                                                cTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGa                                                                                                     <  2:1862948/150‑1 (MQ=255)

                                                cTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGa                                                                                                     <  1:3122450/150‑1 (MQ=255)

                                                 ttGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAg                                                                                                    <  2:8597646/150‑1 (MQ=255)

                                                    agcgagcgGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAgggg                                                                                                 >  2:4492770/1‑150 (MQ=255)

                                                      cgagcgGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCt                                                                                               <  1:2831884/150‑1 (MQ=255)

                                                       gagcgGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTg                                                                                              <  2:7015427/150‑1 (MQ=255)

                                                        agctGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGa                                                                                             <  2:1106595/150‑1 (MQ=255)

                                                        agcgGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGa                                                                                             >  2:6699474/1‑150 (MQ=255)

                                                          cgGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGTATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATa                                                                                           >  2:1988545/1‑150 (MQ=255)

                                                          cgGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATa                                                                                           >  1:2648704/1‑150 (MQ=255)

                                                           gggAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATAc                                                                                          >  2:292747/1‑150 (MQ=255)

                                                             gAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATAccc                                                                                        <  1:3830569/150‑1 (MQ=255)

                                                                aTCTCGGACATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACccgcc                                                                                     <  2:7291142/150‑1 (MQ=255)

                                                                   tcGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTa                                                                                  >  2:468032/1‑150 (MQ=255)

                                                                   tcGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTa                                                                                  >  2:3711740/1‑150 (MQ=255)

                                                                     ggAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTAtc                                                                                >  2:4259964/1‑150 (MQ=255)

                                                                      gAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTAtcg                                                                               >  1:4135610/1‑150 (MQ=255)

                                                                         aTTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTAtcgtcg                                                                            <  1:4291393/150‑1 (MQ=255)

                                                                          tttCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGa                                                                           >  1:6074092/1‑150 (MQ=255)

                                                                              gTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGaa                                                                       >  1:7325799/1‑150 (MQ=255)

                                                                              gTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGaa                                                                       <  2:791335/150‑1 (MQ=255)

                                                                                 ggTACGGTTGTCACCGGTCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACaa                                                                    >  2:6140208/1‑150 (MQ=255)

                                                                                 ggTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACaa                                                                    <  1:7513166/150‑1 (MQ=255)

                                                                                 ggTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACaa                                                                    >  1:4010358/1‑150 (MQ=255)

                                                                                   tACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAAtt                                                                  >  1:1592686/1‑150 (MQ=255)

                                                                                     cGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAAtttt                                                                <  2:2685293/150‑1 (MQ=255)

                                                                                      ggTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTc                                                               <  2:5219120/150‑1 (MQ=255)

                                                                                       gTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCg                                                              <  2:8139208/150‑1 (MQ=255)

                                                                                       gTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCg                                                              <  1:6838502/150‑1 (MQ=255)

                                                                                        ttGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCgg                                                             >  1:8248112/1‑150 (MQ=255)

                                                                                          gTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTc                                                           <  1:2034366/150‑1 (MQ=255)

                                                                                            cACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGt                                                         >  1:1736553/1‑150 (MQ=255)

                                                                                             aCCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGtt                                                        >  2:7722933/1‑150 (MQ=255)

                                                                                                    gACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAgg                                                 >  2:8444440/1‑150 (MQ=255)

                                                                                                      ccccGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGtt                                               >  2:3411334/1‑150 (MQ=255)

                                                                                                       cccGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCCTGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTc                                              >  1:226020/1‑150 (MQ=255)

                                                                                                       cccGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTc                                              >  2:2206966/1‑150 (MQ=255)

                                                                                                        ccGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCa                                             >  1:5505138/1‑150 (MQ=255)

                                                                                                         cGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCAc                                            >  1:1460434/1‑150 (MQ=255)

                                                                                                         cGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCAc                                            >  1:6285495/1‑150 (MQ=255)

                                                                                                              tCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCg                                       >  2:3527004/1‑150 (MQ=255)

                                                                                                              tCTGCATCAGCTGCAGGGAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCg                                       <  1:6373669/150‑1 (MQ=255)

                                                                                                               cTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGa                                      <  2:7309077/150‑1 (MQ=255)

                                                                                                                tGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAg                                     <  2:6010043/150‑1 (MQ=255)

                                                                                                                tGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAg                                     >  2:4347878/1‑150 (MQ=255)

                                                                                                                tGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAg                                     <  1:6642715/150‑1 (MQ=255)

                                                                                                                tGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAg                                     <  2:2428515/150‑1 (MQ=255)

                                                                                                                tGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAg                                     <  1:3160859/150‑1 (MQ=255)

                                                                                                                    tCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATGCCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCa                                 >  1:538066/1‑150 (MQ=255)

                                                                                                                    tCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCa                                 >  2:6820763/1‑150 (MQ=255)

                                                                                                                         tGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGgcgc                            >  1:1803678/1‑150 (MQ=255)

                                                                                                                         tGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGgcgc                            >  2:209860/1‑150 (MQ=255)

                                                                                                                              gTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCa                       <  2:3611191/150‑1 (MQ=255)

                                                                                                                              gTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCa                       <  1:5926898/150‑1 (MQ=255)

                                                                                                                               tAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAg                      >  2:1027380/1‑150 (MQ=255)

                                                                                                                                 gTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAgtc                    <  1:4352450/150‑1 (MQ=255)

                                                                                                                                    gACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAgtcgtc                 >  1:8634164/1‑150 (MQ=255)

                                                                                                                                    gACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAgtcgtc                 >  2:3851165/1‑150 (MQ=255)

                                                                                                                                    gACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAgtcgtc                 <  1:5613641/150‑1 (MQ=255)

                                                                                                                                    gACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAgtcgtc                 >  1:4225705/1‑150 (MQ=255)

                                                                                                                                    gACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAgtcgtc                 <  1:3454727/150‑1 (MQ=255)

                                                                                                                                    gACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAgtcgtc                 >  1:1086180/1‑150 (MQ=255)

                                                                                                                                      ccgtGAGCATGGCGAAATCGCCATGATCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAgtcgtcgt               >  1:2973115/4‑150 (MQ=255)

                                                                                                                                      ccATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAgtcgtcgt               <  1:1229706/150‑1 (MQ=255)

                                                                                                                                      ccATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAgtcgtcgt               >  1:7069137/1‑150 (MQ=255)

                                                                                                                                            tCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAGTCGTCGTCGTCCa         <  2:1398459/150‑1 (MQ=255)

                                                                                                                                            gcATGGCGAAATCGCCATGATCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAGTCGTCGTCGTCCa         >  2:5024419/2‑150 (MQ=255)

                                                                                                                                               tGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAGTCGTCGTCGTCCAGCt      >  1:7204840/1‑150 (MQ=255)

                                                                                                                                                  cGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAGTCGTCGTCGTCCAGCTGGc   <  1:5140461/150‑1 (MQ=255)

                                                                                                                                                   gATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAGTCGTCGTCGTCCAGCTGGcc  <  2:4864585/150‑1 (MQ=255)

CAGCTGCGAAAGGGCAATGACCGGGCAGTTGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAGTCGTCGTCGTCCAGCTGGCC  >  NC_002947/5546568‑5546864

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