BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,095,842	A→C	5.9%	G47G (GGT→GGG)	rhaA ←	L‑rhamnose isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,095,842	0	A	C	5.9%	272.6 / 5.8	110	G47G (GGT→GGG)	rhaA	L‑rhamnose isomerase
Reads supporting (aligned to +/- strand): ref base A (45/57); new base C (2/5); total (48/62)
Fisher's exact test for biased strand distribution p-value = 6.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.12e-01

Only 100 of 114 total aligned reads displayed.

TGGCATGTAAATTAAGCCGTTTCGGCCCCGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGC  >  NC_000913/4095694‑4095991

tGGCATGTAAATTAAGCCGTTTCGGCCCCGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAAc                                                                                                                                                      <  1:1741380/150‑1 (MQ=255)

   cATGTAAATTAAGCCGTTTCGGCCCCGGAATCAGCCGCATAGCCGGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCgg                                                                                                                                                   <  2:1495586/150‑1 (MQ=255)

    aTGTAAATTAAGCCGTTTCGGCCCCGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGa                                                                                                                                                  >  1:833654/1‑150 (MQ=255)

     tGTAAATTAAGCCGTTTCGGCCCCGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGaa                                                                                                                                                 <  1:1300991/150‑1 (MQ=255)

      ttAAATTAAGCCGTTTCGGCCCCGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTTCGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGTTTTGAAAACCCGGaag                                                                                                                                                <  2:756711/149‑2 (MQ=255)

      gTAAATTAAGCCGTTTCGGCCCCGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGaaa                                                                                                                                                <  2:1781969/150‑1 (MQ=255)

       tAAATTAAGCCGTTTCGGCCCCGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGTTTTGAAAACCCGGaagg                                                                                                                                               <  2:133091/150‑3 (MQ=255)

       tAAATTAAGCCGTTTCGGCCCCGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGTTTTGAAAACCCGGaagg                                                                                                                                               <  2:1732549/150‑3 (MQ=255)

       tAAATTAAGCCGTTTCGGCCCCGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGTTTTGAAAACCCGGaagg                                                                                                                                               <  1:1480708/150‑3 (MQ=255)

       tAAATTAAGCCGTATCGGCCCCGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGTTTTGAAAACCCGGaagg                                                                                                                                               <  2:127709/150‑3 (MQ=255)

        aaaTTAAGCCGTTTCGGCCCCGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGATCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAAca                                                                                                                                              <  1:1201956/150‑1 (MQ=255)

         aaTTAAGCCGTTTCGGCCCCGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAAcat                                                                                                                                             <  1:555378/150‑1 (MQ=255)

          aTTAAGCCGTTTCGGCCCCGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGTATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAAcatc                                                                                                                                            <  1:1619840/150‑1 (MQ=255)

            tAAGCCGTTTCGGCCCCGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAAcatcat                                                                                                                                          <  1:227458/150‑1 (MQ=255)

               gCCGTTTCGGCCCCGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGc                                                                                                                                       >  1:17262/1‑150 (MQ=255)

                        gCCCCGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGGTCACAGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTATTCAAACCCGGAAACATCATCGCCCTGCcagc                                                                                                                              >  2:1530190/1‑150 (MQ=255)

                        gCCCCGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCcagc                                                                                                                              >  1:231830/1‑150 (MQ=255)

                         ccccGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAACCCGGAAACATCATCGCCCTGCcagca                                                                                                                             >  2:1488044/1‑150 (MQ=255)

                          cccGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCcagcag                                                                                                                            >  1:1515012/1‑150 (MQ=255)

                           ccGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGt                                                                                                                           >  1:1630068/1‑150 (MQ=255)

                           ccGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGt                                                                                                                           >  2:1034374/1‑150 (MQ=255)

                             ggAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGc                                                                                                                         <  1:941064/150‑1 (MQ=255)

                                 tCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTg                                                                                                                     >  2:1029370/1‑150 (MQ=255)

                                 tCAGCCGAATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGGCAGCGAACCTTCCGTGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTg                                                                                                                     <  2:301591/150‑1 (MQ=255)

                                   aGCCGCATAGCCTGTTCCAGATCTGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGaa                                                                                                                   >  2:1709516/1‑150 (MQ=255)

                                        cATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACggg                                                                                                              <  1:1196068/150‑1 (MQ=255)

                                           aGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAATACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTaa                                                                                                           >  1:240286/1‑150 (MQ=255)

                                            gCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTaaa                                                                                                          >  2:264701/1‑150 (MQ=255)

                                             ccTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAAc                                                                                                         <  2:1638668/150‑1 (MQ=255)

                                              cTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACg                                                                                                        <  1:545576/150‑1 (MQ=255)

                                                gTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGAt                                                                                                      >  2:945661/1‑150 (MQ=255)

                                                gTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGAt                                                                                                      <  1:1057498/150‑1 (MQ=255)

                                                gTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGAt                                                                                                      >  2:491613/1‑150 (MQ=255)

                                                 ttCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAGCCGGAAACATCATCGCCCTGCCAGCAGTGCATTGCAACGGGTAAACGATc                                                                                                     >  1:14050/1‑150 (MQ=255)

                                                 ttCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATc                                                                                                     <  1:1105776/150‑1 (MQ=255)

                                                  tCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGTCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATcc                                                                                                    >  2:1744726/1‑149 (MQ=255)

                                                  tCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCa                                                                                                    >  1:804544/1‑150 (MQ=255)

                                                   ccAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCaa                                                                                                   >  2:1263892/1‑150 (MQ=255)

                                                    cAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAg                                                                                                  >  2:493789/1‑150 (MQ=255)

                                                     aGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGt                                                                                                 >  2:1704949/1‑150 (MQ=255)

                                                     aGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGt                                                                                                 <  2:1565242/150‑1 (MQ=255)

                                                     aGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCAACGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGt                                                                                                 >  2:1681800/1‑150 (MQ=255)

                                                      gATCGGCACGTAGCTCACTGGCATTTCGCGCTTTGCCCGGATAATTGCCTGTGGCCTGACTCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGtt                                                                                                >  1:256465/1‑150 (MQ=255)

                                                      gATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGtt                                                                                                >  1:506185/1‑150 (MQ=255)

                                                        tCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTgg                                                                                              <  2:1040186/150‑1 (MQ=255)

                                                             aCGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAg                                                                                         >  1:1721871/1‑150 (MQ=255)

                                                               gTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAgcg                                                                                       >  1:1508090/1‑150 (MQ=255)

                                                                   cTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGcctc                                                                                   >  1:758331/1‑150 (MQ=255)

                                                                    tCACTGGCATTACGCGCTTTGCCCGGATAATTGACTGTGGCCTGAATCACACCGCTCAGCGAACCTACCGAGTTTTCAAAACAGGAAACATCATCGCCATGCCAGCAGTGAATTGAAACGGGTAAACGATCAAGTTTTCGCAGCGcctcc                                                                                  >  1:30708/1‑150 (MQ=255)

                                                                     cACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCAGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGcctcct                                                                                 >  2:383519/1‑150 (MQ=255)

                                                                       cTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCg                                                                               >  1:1710073/1‑150 (MQ=255)

                                                                        tGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGa                                                                              <  1:1034374/150‑1 (MQ=255)

                                                                         ggCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGAGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGAc                                                                             <  2:1258614/150‑1 (MQ=255)

                                                                         ggCATTACGCCCTTTGCCCGGTTAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGAc                                                                             <  2:1101972/150‑1 (MQ=255)

                                                                            aTTACGCGCTTTGCCCGGATAATTGCCTGTCGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCACCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATc                                                                          >  1:797505/1‑150 (MQ=255)

                                                                              tACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCaa                                                                        <  2:796448/150‑1 (MQ=255)

                                                                                cgcgCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATc                                                                      <  1:1247977/150‑1 (MQ=255)

                                                                                            ggATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGATACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCTGCGCCTCCTCGACATCTCTCCCCGCCGCcgcg                                                          >  2:1230631/1‑150 (MQ=255)

                                                                                                 aTTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGCAACGGGTAAACGATCAAGTTGGCGCAGCGCATCCTCGCCATCAATCCCCACCGCCGCGAAAcc                                                     >  2:1395405/1‑149 (MQ=255)

                                                                                                  ttGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGc                                                    <  2:592589/150‑1 (MQ=255)

                                                                                                   tGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCt                                                   <  1:800058/150‑1 (MQ=255)

                                                                                                    gCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTg                                                  <  1:1278235/150‑1 (MQ=255)

                                                                                                        gtgGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTc                                              >  2:461227/1‑150 (MQ=255)

                                                                                                         tgGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCg                                             <  1:1538371/150‑1 (MQ=255)

                                                                                                         gggCCTTAATACCACCGGTCAGCTAACCTGCCGGGTTTTCAAAACCGGAAAAATCATCGCCCTGCAAGCAGTTCAGTGAAACTGGTAAACGATCAATTTGGCGCAGCGACTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCg                                             <  1:1402446/149‑1 (MQ=255)

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

                                                                                                           gCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGTGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCt                                           <  2:1250586/150‑1 (MQ=255)

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

                                                                                                            ccTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTa                                          <  2:1575559/150‑1 (MQ=255)

                                                                                                             cTGAATCCCCCGGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACAGCGTTTCGcggc                                         <  1:512224/150‑4 (MQ=255)

                                                                                                              tGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGc                                        <  1:1358861/150‑1 (MQ=255)

                                                                                                              tGAATCCCCCCGGGCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGc                                        <  1:111857/150‑1 (MQ=255)

                                                                                                               gAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCt                                       <  1:572282/150‑1 (MQ=255)

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

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

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

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

                                                                                                                       ccGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCCGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCGCCCAGGcc                               >  1:1132503/1‑150 (MQ=255)

                                                                                                                       ccGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGcc                               <  1:899979/150‑1 (MQ=255)

                                                                                                                          gTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCCTCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGt                            <  2:286296/150‑1 (MQ=255)

                                                                                                                            cAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTc                          <  1:773375/150‑1 (MQ=255)

                                                                                                                             agagAACCTTCCGTGGGTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAAGGGGTAAACGATCAAGTTGGCGCAGCTCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTcc                         <  2:865725/147‑1 (MQ=255)

                                                                                                                             aGCGAACCTACCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTcc                         >  1:497938/1‑150 (MQ=255)

                                                                                                                              gCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCa                        >  2:47761/1‑150 (MQ=255)

                                                                                                                                    cTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAg                  <  1:1007805/150‑1 (MQ=255)

                                                                                                                                    cTTCAGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAAGGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAg                  <  1:332976/150‑1 (MQ=255)

                                                                                                                                        cGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAACCGCTGTATCGCTTGATCCCAGGCCTGTTCCAGTTGAGTGGt              >  2:102245/1‑150 (MQ=255)

                                                                                                                                           gTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGATACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCAt           >  2:1458014/1‑150 (MQ=255)

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

                                                                                                                                            ttttCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGATACGCTGTTTCGCTAGCTCACAGGCCTGTTCCAGTTGAGTGGTCATg          >  1:931718/1‑150 (MQ=255)

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

                                                                                                                                             tttCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATgc         <  2:1319459/150‑1 (MQ=255)

                                                                                                                                              ttCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATgcg        <  2:385656/150‑1 (MQ=255)

                                                                                                                                              ttCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATgcg        <  2:239637/150‑1 (MQ=255)

                                                                                                                                              ttCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATgcg        <  1:491613/150‑1 (MQ=255)

                                                                                                                                                cAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCa      <  2:1721871/150‑1 (MQ=255)

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

                                                                                                                                                 aaaaCCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATACCCACCGCCGCGAAACGCTGTATCGCTAGCTCCCAGGCCTGATCCAGTTGAGTGGTCATGCGCaa     >  2:759973/1‑150 (MQ=255)

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

                                                                                                                                                    aCCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGCCATCAATCCCCACCGCCGCGAACCGCTGTTTCGCTAGCTCCCCGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGc  >  1:61176/1‑150 (MQ=255)

TGGCATGTAAATTAAGCCGTTTCGGCCCCGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGC  >  NC_000913/4095694‑4095991

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