BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,408,101	G→A	5.8%	V162V (GTC→GTT)	PP_3021 ←	LysE family transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,408,101	0	G	A	5.8%	303.4 / 5.1	101	V162V (GTC→GTT)	PP_3021	LysE family transporter
Reads supporting (aligned to +/- strand): ref base G (46/49); new base A (3/3); total (49/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

Only 100 of 101 total aligned reads displayed.

ATGCGCACCTTGCAACAAGTCAGCCCCGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAGGCTGGCATAGAACAGCGCGTTCTTGGGGTTGAGGATGCCAGACAAGAAGCCCATGC  >  NC_002947/3407953‑3408249

aTGCGCACCTTGCAACAAGTCAGCCCCGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGa                                                                                                                                                     <  2:423888/150‑1 (MQ=255)

 tGCGCACCTTGCAACAAGTCAGCCCCGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGAc                                                                                                                                                    <  2:23283/150‑1 (MQ=255)

   cgcACCTTGCATCAAGTCAGCCCCGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCTGCTGGTCGCCG‑TGGCCATCGGCAACCa                                                                                                                                                  <  1:1225726/150‑1 (MQ=255)

      aCCTTGCAACAAGTCAGCCCCGGGCGAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCa                                                                                                                                               >  1:706794/1‑150 (MQ=255)

       ccTTGCAACAAGTCAGCCCCGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCaa                                                                                                                                              <  1:288167/150‑1 (MQ=255)

       ccTTGCAACAAGTCAGCCCCGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCACGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCaa                                                                                                                                              <  2:1281302/150‑1 (MQ=255)

            caacaaGTCAGCCCCGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTccc                                                                                                                                         >  2:33192/1‑150 (MQ=255)

             aacaaGTCAGCCCCGGGCGAGATGCAGCAACAGGGCCCTCACCAACAGGACCGGGATTACCCCCGGAAGGGGCGCCAGCCATGCCTTGTTGCTGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCC‑CCGGAGCAAGCAAGTCCCa                                                                                                                                        >  1:1077638/1‑150 (MQ=255)

             aacaaGTCAGCCCCGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCGCa                                                                                                                                        >  2:828622/1‑150 (MQ=255)

                aaGCCAGCCCCGGGCCAGATGCAGCAACAACGCCCCCGACAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCagag                                                                                                                                     <  1:271249/150‑1 (MQ=255)

                 aGTCAGCCCCGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGc                                                                                                                                    >  1:10284/1‑150 (MQ=255)

                  gTCAGCCCCGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCa                                                                                                                                   >  1:221330/1‑150 (MQ=255)

                     aGCCCCGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGa                                                                                                                                >  2:1497388/1‑150 (MQ=255)

                         ccGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGAt                                                                                                                            <  1:1677747/150‑1 (MQ=255)

                         ccGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGAt                                                                                                                            >  2:481314/1‑150 (MQ=255)

                          cGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATg                                                                                                                           >  2:666142/1‑150 (MQ=255)

                          cGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATg                                                                                                                           >  1:1168706/1‑150 (MQ=255)

                               cAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGAAGATGCTGaa                                                                                                                      >  1:1022051/1‑150 (MQ=255)

                                aGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCCGAGCAGGACGATGCTGAAc                                                                                                                     >  2:1759240/1‑150 (MQ=255)

                                  aTGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACAt                                                                                                                   >  1:1062826/1‑150 (MQ=255)

                                         aacaacGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAgc                                                                                                            <  2:1171913/150‑1 (MQ=255)

                                          acaacGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAgcg                                                                                                           >  1:1064933/1‑150 (MQ=255)

                                          acaacGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAgcg                                                                                                           >  1:1099664/1‑150 (MQ=255)

                                            accGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGcgcg                                                                                                         <  1:173179/148‑1 (MQ=255)

                                             acGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGt                                                                                                        <  2:656381/150‑1 (MQ=255)

                                                 ccccGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGc                                                                                                    >  2:686745/1‑150 (MQ=255)

                                                   ccGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCg                                                                                                  >  1:1517916/1‑150 (MQ=255)

                                                     gCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGtt                                                                                                <  1:864792/150‑1 (MQ=255)

                                                         aCATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTcc                                                                                            >  1:1464812/1‑150 (MQ=255)

                                                         aCATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTcc                                                                                            >  1:355803/1‑150 (MQ=255)

                                                           aTCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGGAGCTCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAg                                                                                          >  1:1572978/1‑150 (MQ=255)

                                                           aTCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCTCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGttt                                                                                               >  2:519663/1‑145 (MQ=255)

                                                           aTCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCTCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGttt                                                                                               <  1:519663/145‑1 (MQ=255)

                                                           aTCACCAGCATGACCCCCGAGGCACGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAg                                                                                          >  1:573118/1‑150 (MQ=255)

                                                              aCCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCg                                                                                       >  1:1030551/1‑150 (MQ=255)

                                                               ccAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCgg                                                                                      <  1:828622/150‑1 (MQ=255)

                                                                cAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTc                                                                                             >  2:1226739/1‑142 (MQ=255)

                                                                cAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTc                                                                                             <  1:1226739/142‑1 (MQ=255)

                                                                cAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCggg                                                                                     >  2:1643543/1‑150 (MQ=255)

                                                                cAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCggg                                                                                     >  1:285802/1‑150 (MQ=255)

                                                                cAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCggg                                                                                     >  1:283841/1‑150 (MQ=255)

                                                                cAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCggg                                                                                     <  2:446001/150‑1 (MQ=255)

                                                                cAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCggg                                                                                     <  2:1062826/150‑1 (MQ=255)

                                                                 aGCATGACCCCCGCGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGGTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCgggg                                                                                    <  2:655059/150‑1 (MQ=255)

                                                                 aGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCgggg                                                                                    <  1:666142/150‑1 (MQ=255)

                                                                     tGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGa                                                                                <  2:925956/150‑1 (MQ=255)

                                                                      gACCCCCGCGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGAc                                                                               >  2:1662436/1‑150 (MQ=255)

                                                                        cccccGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGc                                                                             <  1:686745/150‑1 (MQ=255)

                                                                        cccccGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGc                                                                             >  2:124248/1‑150 (MQ=255)

                                                                         ccccGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCt                                                                            >  1:1886743/1‑150 (MQ=255)

                                                                           ccGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTgg                                                                          >  2:1782216/1‑150 (MQ=255)

                                                                              aGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGa                                                                       <  2:1168706/150‑1 (MQ=255)

                                                                                 cccGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCTGCTGGTCGCCG‑TGGCCATCGGCAACCAGcgg                                                                                                                                              >  2:1720501/1‑74 (MQ=18)

                                                                                 cccGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCTGCTGGTCGCCG‑TGGCCATCGGCAACCAGcgg                                                                                                                                              <  1:1720501/76‑3 (MQ=18)

                                                                                  ccGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCaa                                                                                                                                              <  2:575883/75‑1 (MQ=255)

                                                                                  ccGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCaa                                                                                                                                              >  1:575883/1‑75 (MQ=255)

                                                                                   cGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATg                                                                  <  2:1749101/150‑1 (MQ=255)

                                                                                   cGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATg                                                                  >  2:1431183/1‑150 (MQ=255)

                                                                                       ccAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTgg                                                              <  1:1782216/150‑1 (MQ=255)

                                                                                         aGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGcc                                                            >  2:1440020/1‑150 (MQ=255)

                                                                                          gCCACACCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTcc                                                                                                                                          <  2:595439/71‑1 (MQ=255)

                                                                                          gCCACACCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTcc                                                                                                                                          >  1:595439/1‑71 (MQ=255)

                                                                                            cACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCt                                                                                                                         <  1:533946/86‑1 (MQ=255)

                                                                                            cACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCt                                                                                                                         >  2:533946/1‑86 (MQ=255)

                                                                                            cACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCggg                                                         >  1:200745/1‑150 (MQ=255)

                                                                                             aCGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAGGc                                                        <  2:285802/150‑1 (MQ=255)

                                                                                             aCGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAGGc                                                        <  2:283841/150‑1 (MQ=255)

                                                                                                 cAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAGGCTGGc                                                    <  2:355803/150‑1 (MQ=255)

                                                                                                  aGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAGGCTGGCa                                                   <  1:1497388/150‑1 (MQ=255)

                                                                                                    gCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAGGCTGGCATa                                                 <  2:221330/150‑1 (MQ=255)

                                                                                                      tGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAGGCTgacttgtt                                               <  2:1077638/150‑8 (MQ=255)

                                                                                                      tGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAGGCTGGCATAGa                                               >  2:645493/1‑150 (MQ=255)

                                                                                                                 cgccgcAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAGGCTGGCATAGAACAGCGCGTTc                                    <  1:124248/150‑1 (MQ=255)

                                                                                                                    cgcAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAGGCTGGCATAGAACAGCGCGTTCTTg                                 <  2:1548837/150‑1 (MQ=255)

                                                                                                                    cgcAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAGGCTGGCATAGAACAGCACGGTCTTg                                 >  2:583524/1‑150 (MQ=255)

                                                                                                                     gcAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAGGCTGGCATAGAACAGCGCGTTCTTgg                                <  2:10284/150‑1 (MQ=255)

                                                                                                                     gcAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAGGCTGGCATAGAACAGCGCGTTCTTgg                                >  2:516511/1‑150 (MQ=255)

                                                                                                                     gcAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAGGCTGGCATAGAACAGCGCGTTCTTgg                                <  2:1517916/150‑1 (MQ=255)

                                                                                                                      cAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAGGCTGGCATAGAACAGCGCGTTCTTggg                               <  1:481314/150‑1 (MQ=255)

                                                                                                                        gTACCCGCGGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGATGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTTGCGACCATGCTGGCCAGGCTGGCATAGAACAGCGCGTTCTTGGGGt                             <  1:1525056/150‑1 (MQ=255)

                                                                                                                         tACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAGGCTGGCATAGAACAGCGCGTTCTTGGGGtt                            <  1:1733001/150‑1 (MQ=255)

                                                                                                                           cccGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAGGCTGGCATAGAACAGCGCGTTCTTGGGGTTGa                          <  2:573118/150‑1 (MQ=255)

                                                                                                                            ttgCTGGTCGCCG‑TGGCCATCGGCAACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAGGCTGGCATAGaa                                              >  2:1784327/3‑129 (MQ=38)

                                                                                                                            ttgCTGGTCGCCG‑TGGCCATCGGCAACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAGGCTGGCATAGaa                                              <  1:1784327/127‑1 (MQ=38)

                                                                                                                            ccGCTGGTTGCCGCTGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGAAGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAGGCTGGCATAGAACAGCGCGTTCTTGGGGTTGAg                         <  2:1350182/150‑1 (MQ=255)

                                                                                                                            ccGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCa                                                                    >  2:1556307/1‑107 (MQ=255)

                                                                                                                            ccGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCa                                                                    <  1:1556307/107‑1 (MQ=255)

                                                                                                                            ccGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAGGCTGGCATAGAACAGCGCGTTCTTGGGGTTGAg                         <  2:1218911/150‑1 (MQ=255)

                                                                                                                            ccGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAGGCTGGCATAGAACAGCGCGTTCTTGGGGTTGAg                         >  1:91972/1‑150 (MQ=255)

                                                                                                                                 ggTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAGGCTGGCATAGAACAGCGCGTTCTTGGGGTTGAGGATGc                    <  1:33192/150‑1 (MQ=255)

                                                                                                                                   ttGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAGGCTGGCATAGAACAGCGCGTTCTTGGGGTTGAGGATGCCa                  <  1:1643543/150‑1 (MQ=255)

                                                                                                                                     gCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAGGCTGGCATAGAACAGCGCGTTCTTGGGGTTGAGGATGCCAGa                <  2:1572978/150‑1 (MQ=255)

                                                                                                                                      ccGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGtt                                                                                                >  1:577010/1‑69 (MQ=255)

                                                                                                                                      ccGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGtt                                                                                                <  2:577010/69‑1 (MQ=255)

                                                                                                                                       cGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAGGCTGGCATAGAACAGCGCGTTCTTGGGGTTGAGGATGCCAGACa              >  1:875585/1‑150 (MQ=255)

                                                                                                                                             ccA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAGGCTGGCCTCCACTTCCACGTTCTTGGGGTTGAGGATGCCAGACAAGAAGc        >  2:1894201/1‑150 (MQ=255)

                                                                                                                                              cA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAGGCTGGCATAGAACAGCGCGTTCTTGGGGTTGAGGATGCCAGACAAGAAGcc       <  2:1030551/150‑1 (MQ=255)

                                                                                                                                                  ggCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAGGCTGGCATAGAACAGCGCGTTCTTGGGGTTGAGGATGCCAGACAAGAAGCCCAt    >  2:1693894/1‑150 (MQ=255)

                                                                                                                                                   gCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAGGCTGGCATAGAACAGCGCGTTCTTGGGGTTGAGGATGCCAGACAAGAAGCCCATg   <  1:1684026/150‑1 (MQ=255)

                                                                                                                                                    cGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAGGCTGGCATAGAACAGCGCGTTCTTGAGTTTGAGGATGCCAGACAAGAAGCCCATGc  <  1:348566/150‑1 (MQ=255)

ATGCGCACCTTGCAACAAGTCAGCCCCGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAGGCTGGCATAGAACAGCGCGTTCTTGGGGTTGAGGATGCCAGACAAGAAGCCCATGC  >  NC_002947/3407953‑3408249

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