BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,095,940	G→A	9.4%	R15C (CGT→TGT)	rhaA ←	L‑rhamnose isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,095,940	0	G	A	9.4%	274.9 / 12.7	95	R15C (CGT→TGT)	rhaA	L‑rhamnose isomerase
Reads supporting (aligned to +/- strand): ref base G (42/43); new base A (6/3); total (48/47)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

Only 100 of 101 total aligned reads displayed.

ATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTAAAGGTGGTCAGATTCGCGGTGGTGCTGAC  >  NC_000913/4095791‑4096089

aTTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACg                                                                                                                                                       <  1:1041530/150‑1 (MQ=255)

    ccTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGt                                                                                                                                                   >  1:991719/1‑150 (MQ=255)

       gtgGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTc                                                                                                                                                <  1:137371/150‑1 (MQ=255)

         ggCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGc                                                                                                                                              >  2:1456898/1‑150 (MQ=255)

          gCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCt                                                                                                                                             <  1:306046/150‑1 (MQ=255)

           ccTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTa                                                                                                                                            >  2:1331808/1‑150 (MQ=255)

           ccTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACAGCGTTTCGcgg                                                                                                                                            <  1:821493/150‑3 (MQ=255)

           ccTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACAGCGTTTCGcgg                                                                                                                                            >  2:821493/1‑148 (MQ=255)

             tGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACg                                                                                                                                                       <  1:639139/137‑1 (MQ=255)

             tGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACg                                                                                                                                                       >  2:639139/1‑137 (MQ=255)

              gAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCt                                                                                                                                         >  1:1323694/1‑150 (MQ=255)

              gAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACAGCGTTTCGcggcgg                                                                                                                                         <  2:319304/150‑6 (MQ=255)

              gAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACAGCGTTTCGcggcgg                                                                                                                                         <  2:700461/150‑6 (MQ=255)

                aTCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTcc                                                                                                                                       >  1:829946/1‑150 (MQ=255)

                  ccccccGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCa                                                                                                                                     >  1:888418/1‑150 (MQ=255)

                   cccccGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAg                                                                                                                                    >  2:1375365/1‑150 (MQ=255)

                    ccccGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAgg                                                                                                                                   <  2:308770/150‑1 (MQ=255)

                       cGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACg                                                                                                                                                       <  1:378378/127‑1 (MQ=255)

                       cGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACg                                                                                                                                                       >  2:378378/1‑127 (MQ=255)

                            aGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTcc                                                                                                                           >  2:622211/1‑150 (MQ=255)

                              cGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAg                                                                                                                         <  1:1601896/150‑1 (MQ=255)

                               gAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGt                                                                                                                        >  2:443899/1‑150 (MQ=255)

                                aaCCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGtt                                                                                                                       <  1:1008095/150‑1 (MQ=255)

                                 tgtttCCGGTTTTGAAAACCCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTg                                                                                                                      >  1:1062330/4‑150 (MQ=255)

                                 tgtttCCGGTTTTGAAAACCCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTg                                                                                                                      >  1:1024876/4‑150 (MQ=255)

                                  gtttCCGGTTTTGAAAACCCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGa                                                                                                                     >  1:837361/3‑150 (MQ=255)

                                   cTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAg                                                                                                                    >  1:1109530/1‑150 (MQ=255)

                                        gggTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTc                                                                                                               <  1:232527/150‑1 (MQ=255)

                                         ggTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCa                                                                                                              >  1:603463/1‑150 (MQ=255)

                                              tCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATgcgc                                                                                                         >  1:1253943/1‑150 (MQ=255)

                                               cAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTNGNCATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATgcgcn                                                                                                        <  2:726345/150‑2 (MQ=255)

                                                aaaaCCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCaa                                                                                                       >  2:1424431/1‑150 (MQ=255)

                                                aaaaCCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCAAGTTGAGTGGTCATGCGCaa                                                                                                       >  2:743626/1‑150 (MQ=255)

                                                   aCCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGc                                                                                                    <  1:1477616/150‑1 (MQ=255)

                                                      ggAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCGCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCGcc                                                                                                 <  2:1468048/150‑1 (MQ=255)

                                                        aaaCATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCtt                                                                                               <  1:873067/150‑1 (MQ=255)

                                                          aCATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTtg                                                                                             <  1:622211/150‑1 (MQ=255)

                                                          aCATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTtg                                                                                             <  2:1029828/150‑1 (MQ=255)

                                                          aCATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTtg                                                                                             <  1:1331808/150‑1 (MQ=255)

                                                           catcatCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTtgt                                                                                            >  1:1513661/1‑150 (MQ=255)

                                                           catcatCGCCCTGCCAGCAGTACATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTtgt                                                                                            >  2:861962/1‑150 (MQ=255)

                                                              catcGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGTGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTtgtctg                                                                                         >  2:971218/1‑150 (MQ=255)

                                                              catcGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTtgtctg                                                                                         >  2:1392150/1‑150 (MQ=255)

                                                              catcGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTtgtctg                                                                                         >  2:310164/1‑150 (MQ=255)

                                                              catcGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTtgtctg                                                                                         <  1:1456898/150‑1 (MQ=255)

                                                                   cccTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCgtg                                                                                    <  2:991719/150‑1 (MQ=255)

                                                                    ccTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCgtgt                                                                                   <  2:888418/150‑1 (MQ=255)

                                                                      tGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTag                                                                                 >  2:510555/1‑150 (MQ=255)

                                                                       gCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTaga                                                                                <  2:947516/150‑1 (MQ=255)

                                                                           gcagTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGa                                                                            >  2:863292/1‑150 (MQ=255)

                                                                                gCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTg                                                                       >  1:665791/1‑150 (MQ=255)

                                                                                  aTTGGAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGATCCTTTGTCTGTCGTGTAGAGTGAATCTgcg                                                                     >  2:864308/1‑150 (MQ=255)

                                                                                   ttGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTgcgc                                                                    <  2:1591498/150‑1 (MQ=255)

                                                                                    tGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGcc                                                                   <  2:711884/150‑1 (MQ=255)

                                                                                     gAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCca                                                                  >  2:417210/1‑150 (MQ=255)

                                                                                            gTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTg                                                           <  1:997420/150‑1 (MQ=255)

                                                                                             tAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTgg                                                          >  2:1576500/1‑150 (MQ=255)

                                                                                              aaaCGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTggg                                                         <  1:596120/150‑1 (MQ=255)

                                                                                                 cGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCaa                                                      >  1:73368/1‑150 (MQ=255)

                                                                                                     cAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTc                                                  <  2:648422/150‑1 (MQ=255)

                                                                                                     cAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTc                                                  <  2:29360/150‑1 (MQ=255)

                                                                                                      aaGTTGGCGCAGTGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCa                                                 <  1:971218/150‑1 (MQ=255)

                                                                                                      aaGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCa                                                 <  1:310164/150‑1 (MQ=255)

                                                                                                       aGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGt                                                                              >  2:1502378/1‑120 (MQ=255)

                                                                                                       aGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGt                                                                              <  1:1502378/120‑1 (MQ=255)

                                                                                                       aGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCAc                                                <  2:829946/150‑1 (MQ=255)

                                                                                                       aGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCAc                                                <  1:417210/150‑1 (MQ=255)

                                                                                                               cAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGa                                        <  1:743626/150‑1 (MQ=255)

                                                                                                                 gcgcCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGAtt                                      >  1:592957/1‑150 (MQ=255)

                                                                                                                 gcgcCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGAtt                                      <  2:1105936/150‑1 (MQ=255)

                                                                                                                    cctcctCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGTCCTGTTCCCGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTTCGCAACATAGTGGGCAATTTCACTGTCAGGATTAgg                                   <  2:30757/150‑1 (MQ=255)

                                                                                                                    cctcctCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAgg                                   >  2:221798/1‑150 (MQ=255)

                                                                                                                     ctcctcGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAggg                                  <  2:874068/150‑1 (MQ=255)

                                                                                                                     ctcctcGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAggg                                  >  1:160031/1‑150 (MQ=255)

                                                                                                                      tcctcGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAgggg                                 <  2:1109530/150‑1 (MQ=255)

                                                                                                                      tcctcGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAgggg                                 <  1:1564595/150‑1 (MQ=255)

                                                                                                                      tcctcGACATAAATCCCCACCGCCGCGAAACTCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAgggg                                 <  2:1426928/150‑1 (MQ=255)

                                                                                                                        ctcGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTa                               <  2:73368/150‑1 (MQ=255)

                                                                                                                        ctcGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTa                               <  1:1392150/150‑1 (MQ=255)

                                                                                                                        ctcGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTa                               >  1:1303207/1‑150 (MQ=255)

                                                                                                                        ctcGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTa                               >  2:705452/1‑150 (MQ=255)

                                                                                                                           gACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAACGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTAAAg                            >  1:865114/1‑150 (MQ=255)

                                                                                                                             cATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTAAAggt                          <  1:443899/150‑1 (MQ=255)

                                                                                                                               tCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTAAAggtgg                        >  1:871873/1‑150 (MQ=255)

                                                                                                                                  aTCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTAAAGGTGGTCa                     >  1:960878/1‑150 (MQ=255)

                                                                                                                                     cccACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTAAAGGTGGTCAGAt                  <  2:1323694/150‑1 (MQ=255)

                                                                                                                                      ccACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTAAAGGTGGTCAGAtt                 >  2:210261/1‑150 (MQ=255)

                                                                                                                                       cACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTAAAGGTGGTCAGATTc                >  2:1211929/1‑150 (MQ=255)

                                                                                                                                          gctagcgAAACAGCGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCATATAGTGGGCAATTTCACTGTCAGGATTAGGGGTAAAGGTGGTCAGATTcgcg             >  2:1172556/5‑150 (MQ=255)

                                                                                                                                          gctagcgAAACAGCGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTAAAGGTGGTCAGATTcgcg             >  1:157740/5‑150 (MQ=255)

                                                                                                                                           ctagcgAAACAGCGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTAAAGGTGGTCAGATTCGCgg            >  2:954256/4‑150 (MQ=255)

                                                                                                                                            tagcgAAACAGCGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTAAAGGTGGTCAGATTCGCggt           >  1:1228420/3‑150 (MQ=255)

                                                                                                                                             cgcgAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTAAAGGTGGTCAGATTCGCggtg          <  1:705452/150‑1 (MQ=255)

                                                                                                                                             agcgAAACAGCGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTAAAGGTGGTCAGATTCGCggtg          >  1:327216/2‑150 (MQ=255)

                                                                                                                                              gcgAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTAAAGGTGGTCAGATTCGCggtgg         >  2:208352/1‑150 (MQ=255)

                                                                                                                                                 aaaCGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTAAAGGTGGTCAGATTCGCGGTGGTGc      <  1:1375365/150‑1 (MQ=255)

                                                                                                                                                  aaCGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTAAAGGTGGTCAGATTCGCGGTGGTGCt     <  2:603463/150‑1 (MQ=255)

                                                                                                                                                    cGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTAAAGGTGGTCAGATTCGCGGTGGTGCTga   <  2:1024876/150‑1 (MQ=255)

                                                                                                                                                    cGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGCCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTAAAGGTGGTCAGATTCGCGGTGGTGCTga   <  1:1550182/150‑1 (MQ=255)

                                                                                                                                                     gCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTAAAGGTGGTCAGATTCGCGGTGGTGCTgac  >  1:879870/1‑150 (MQ=255)

ATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTAAAGGTGGTCAGATTCGCGGTGGTGCTGAC  >  NC_000913/4095791‑4096089

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