BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,095,965	T→G	5.6%	E6D (GAA→GAC)	rhaA ←	L‑rhamnose isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,095,965	0	T	G	5.6%	1709.2 / 9.6	615	E6D (GAA→GAC)	rhaA	L‑rhamnose isomerase
Reads supporting (aligned to +/- strand): ref base T (267/296); new base G (2/50); total (269/346)
Fisher's exact test for biased strand distribution p-value = 2.21e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 625 total aligned reads displayed.

GTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTAAAGGTGGTCAGATTCGCGGTGGTGCTGACGACCTGACGGAAATCATCCACATTG  >  NC_000913/4095816‑4096114

  cAGCGAACCGTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTc                                                                                                                                                     <  1:5743725/150‑1 (MQ=255)

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

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

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

     cGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAg                                                                                                                                                  >  2:4367223/1‑150 (MQ=255)

      gAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGt                                                                                                                                                 <  2:610076/150‑1 (MQ=255)

         ccTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGa                                                                                                                                              <  2:5839690/150‑1 (MQ=255)

         ccTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGa                                                                                                                                              <  2:2400001/150‑1 (MQ=255)

           ttCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGt                                                                                                                                            >  1:8094058/1‑150 (MQ=255)

             ccGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTgg                                                                                                                                          >  1:2626886/1‑150 (MQ=255)

                ggTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGGTGAGTGGTCa                                                                                                                                       >  2:1057757/1‑150 (MQ=255)

                 gTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCAt                                                                                                                                      >  1:8816059/1‑150 (MQ=255)

                  ttttCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATg                                                                                                                                     <  1:4757078/150‑1 (MQ=255)

                  ttttCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTAGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATg                                                                                                                                     >  2:8652992/1‑150 (MQ=255)

                   tttCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCTCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATgc                                                                                                                                    <  1:4437754/150‑1 (MQ=255)

                        aaaCCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCaaa                                                                                                                               >  1:6802517/1‑150 (MQ=255)

                        aaaCCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCaaa                                                                                                                               <  1:6502736/150‑1 (MQ=255)

                        aaaCCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCGGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCaaa                                                                                                                               >  1:1018025/1‑150 (MQ=255)

                         aaCCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAg                                                                                                                              <  1:888192/150‑1 (MQ=255)

                         aaCCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTCGAGTGGTCATGCGCAAAg                                                                                                                              >  1:9162429/1‑150 (MQ=255)

                            cGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTc                                                                                                                           >  1:8204628/1‑150 (MQ=255)

                            cGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTc                                                                                                                           >  1:2976332/1‑150 (MQ=255)

                             ggAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTcc                                                                                                                          <  1:6133387/150‑1 (MQ=255)

                               aaaCATCATCGCCCTGCCAGCAGTGCATCGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGGTCCAGTTGAGTGGTCATGCGCAAAGCTCCtt                                                                                                                        <  1:2684932/150‑1 (MQ=255)

                                   atcatcGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTtgtc                                                                                                                    >  1:3410530/1‑150 (MQ=255)

                                    tcatcGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTtgtct                                                                                                                   <  1:6426088/150‑1 (MQ=255)

                                      atcGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTtgtctgt                                                                                                                 >  1:7902110/1‑150 (MQ=255)

                                      atcGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTtgtctgt                                                                                                                 >  2:1466172/1‑150 (MQ=255)

                                       tcGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTtgtctgtc                                                                                                                >  1:3061120/1‑150 (MQ=255)

                                       tcGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTtgtctgtc                                                                                                                <  1:2209730/150‑1 (MQ=255)

                                         gCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCgt                                                                                                              >  1:5525688/1‑150 (MQ=255)

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

                                           ccTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCgtgt                                                                                                            >  2:6356441/1‑150 (MQ=255)

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

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

                                              gCCAGCAGTGCATTGAAACGGGTAAACGATCACGTTGGCGCAGCGCCTCCTCGACATCACTCCCCACCGCCGCGAAACGCTGTTTCGCTAGCGCCCAGGCCGGGTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTaga                                                                                                         <  2:6991724/150‑1 (MQ=255)

                                              gCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCACTCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGGTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTaga                                                                                                         <  2:7706596/150‑1 (MQ=255)

                                              gCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGGTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTaga                                                                                                         <  1:7605718/150‑1 (MQ=255)

                                               ccAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGTGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTagag                                                                                                        >  1:7192676/1‑150 (MQ=255)

                                               ccAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTagag                                                                                                        >  2:888514/1‑150 (MQ=255)

                                                        cATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTgc                                                                                               >  1:8511250/1‑150 (MQ=255)

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

                                                              aaCGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCcaca                                                                                         >  2:3248126/1‑150 (MQ=255)

                                                                 gggTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAg                                                                                      >  1:1553510/1‑150 (MQ=255)

                                                                  cctAAACGATCAAGTAGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCAGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGGGTAGAGTGAATCTGCGCCACAGAGt                                                                                     >  2:4192565/3‑150 (MQ=255)

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

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

                                                                   gTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGGTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTg                                                                                    <  2:7368543/150‑1 (MQ=255)

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

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

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

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

                                                                               gTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACt                                                                        <  2:1702449/150‑1 (MQ=255)

                                                                                    cgcAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAg                                                                   <  1:4658970/150‑1 (MQ=255)

                                                                                    cgcAGCGCCTCCTCGACATCAATCCCCACCGACGCGAAACGCTGTTTCGCTAGCTACCAGGCCTGTTCCTGTTGAGTGGTCATGCGCAAATCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAg                                                                   >  1:2770962/1‑150 (MQ=255)

                                                                                         cgcCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTa                                                              <  1:9214039/150‑1 (MQ=255)

                                                                                         cgcCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTa                                                              <  1:488426/150‑1 (MQ=255)

                                                                                           cctcctCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAgg                                                            <  2:2937676/150‑1 (MQ=255)

                                                                                           cctcctCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAgg                                                            <  2:210498/150‑1 (MQ=255)

                                                                                            ctcctcGACATCAATCCCCACCGCCGCGAACCGCTGTTTCGCTAGCTCCCAGGCCTGGGCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAggg                                                           <  2:7599969/150‑1 (MQ=255)

                                                                                             tcctcGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTc                                                                                                                                        >  1:7733947/1‑72 (MQ=255)

                                                                                              cctcGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGt                                                         >  1:4004769/1‑150 (MQ=255)

                                                                                              cctcGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGGGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGt                                                         >  1:6287835/1‑150 (MQ=255)

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

                                                                                                tcGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTaa                                                       >  2:5221635/1‑150 (MQ=255)

                                                                                                 cGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTaaa                                                      >  2:7848176/1‑150 (MQ=255)

                                                                                                 cGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTaaa                                                      <  1:3896956/150‑1 (MQ=255)

                                                                                                  gACATCACTCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTAAAg                                                     <  1:2424095/150‑1 (MQ=255)

                                                                                                  gACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTAAAg                                                     <  1:7338897/150‑1 (MQ=255)

                                                                                                     aTCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTAAAggtg                                                  <  1:8476527/150‑1 (MQ=255)

                                                                                                     aTCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTAAAggtg                                                  >  1:4115025/1‑150 (MQ=255)

                                                                                                      tCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTAAAggtgg                                                 <  2:3873779/150‑1 (MQ=255)

                                                                                                       cAATCCCCCCCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTAAAggtggt                                                <  2:1894304/150‑1 (MQ=255)

                                                                                                       cAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTAAAggtggt                                                >  2:3711481/1‑150 (MQ=255)

                                                                                                       cAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTGCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTAAAggtggt                                                <  2:3457640/150‑1 (MQ=255)

                                                                                                       cAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGGTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTAAAggtggt                                                <  1:1224953/150‑1 (MQ=255)

                                                                                                        aaTCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTAAAGGTGGTc                                               <  1:5907603/150‑1 (MQ=255)

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

                                                                                                              cACCGCCGCGACACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTAAAGGTGGTCAGATTc                                         <  2:7123913/150‑1 (MQ=255)

                                                                                                                ccgccgCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTAAAGGTGGTCAGATTcgc                                       >  1:6864371/1‑150 (MQ=255)

                                                                                                                ccgccgCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGGTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTAAAGGTGGTCAGATTcgc                                       <  2:9148464/150‑1 (MQ=255)

                                                                                                                   ccgcGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTAAAGGTGGTCAGATTCGCggt                                    <  2:4841914/150‑1 (MQ=255)

                                                                                                                   ccgcGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTAAAGGTGGTCAGATTCGCggt                                    >  2:276513/1‑150 (MQ=255)

                                                                                                                   ccgcGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTAAAGGTGGTCAGATTCGCggt                                    <  2:6647769/150‑1 (MQ=255)

                                                                                                                   ccgcGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTAAAGGTGGTCAGATTCGCggt                                    >  2:2507165/1‑150 (MQ=255)

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

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

                                                                                                                     gcgAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGGTGAGTGGTCCTGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGAGTAGGGGTAAAGGTGGTCAGATTCGCggtgg                                  >  2:6277520/1‑150 (MQ=255)

                                                                                                                       gAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTAAAGGTGGTCAGATTCGCGgtggtg                                <  2:6006653/150‑1 (MQ=255)

                                                                                                                       gAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTAAAGGTGGTCAGATTCGCGgtggtg                                <  2:4539416/150‑1 (MQ=255)

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

                                                                                                                           cGCTGTTTCGCTAGCTCCCAGGCCTGGTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTAAAGGTGGTCAGATTCGCGGTGGTGCTga                            <  1:3227016/150‑1 (MQ=255)

                                                                                                                                    gCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTAAAGGTGGTCAGATTCGCGGTGGTGCTGACGACCTGAc                   <  1:5147020/150‑1 (MQ=255)

                                                                                                                                        gCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTAAAGGTGGTCAGATTCGCGGTGGTGCTGACGACCTGACGGaa               <  2:3978234/150‑1 (MQ=255)

                                                                                                                                          tCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTAAAGGTGGTCAGATTCGCGGTGGTGCTGACGACCTGACGGAAat             >  1:2029949/1‑150 (MQ=255)

                                                                                                                                          tCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTAAAGGTGGTCAGATTCGCGGTGGTGCTGACGACCTGACGGAAat             >  1:1697501/1‑150 (MQ=255)

                                                                                                                                           cccAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTAAAGGTGGTCAGATTCGCGGTGGTGCTGACGACCTGACGGAAatc            <  2:4762048/150‑1 (MQ=255)

                                                                                                                                            ccAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTAAAGGTGGTCAGATTCGCGGTGGTGCTGACGACCTGACGGAAatca           <  1:9003419/150‑1 (MQ=255)

                                                                                                                                             cAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTAAAGGTGGTCAGATTCTCGGTGGTGCTGACGACCTGACGGAAatcat          <  2:5344621/150‑1 (MQ=255)

                                                                                                                                                gCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTAAAGGTGGTCAGATTCGCGGTGGTGCTGACGACCTGACGGAAATCATCca       >  2:8392957/1‑150 (MQ=255)

GTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTAAAGGTGGTCAGATTCGCGGTGGTGCTGACGACCTGACGGAAATCATCCACATTG  >  NC_000913/4095816‑4096114

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