BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,095,833	G→T	8.0%	N50K (AAC→AAA)	rhaA ←	L‑rhamnose isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,095,833	0	G	T	8.0%	328.3 / 8.4	98	N50K (AAC→AAA)	rhaA	L‑rhamnose isomerase
Reads supporting (aligned to +/- strand): ref base G (44/46); new base T (5/3); total (49/49)
Fisher's exact test for biased strand distribution p-value = 7.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 101 total aligned reads displayed.

TCCAGATAGATGGCATGTAAATTAAGCCGTTTCGGCCCCGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCA  >  NC_000913/4095684‑4095981

tCCAGATAGATGGCATGTAAATTAAGCCGTTTCGGCCCCGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCggg                                                                                                                                                      <  2:407713/150‑1 (MQ=255)

 ccAGATAGATGGCATGTAAATTAAGCCGTTTCGGCCCCGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGt                                                                                                                                                     >  1:399384/1‑150 (MQ=255)

  cAGATAGATGGCATGTAAATTAAGCCGTTTCGGCCCCGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGtt                                                                                                                                                    >  2:1471383/1‑150 (MQ=255)

  cAGATAGATGGCATGTAAATTAAGCCGTTTCGGCCCCGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGtt                                                                                                                                                    >  1:477207/1‑150 (MQ=255)

  cAGATAGATGGCATGTAAATTAAGCCGTTTCGGCCCCGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGtt                                                                                                                                                    >  2:773282/1‑150 (MQ=255)

   agatagatGGCATGTAAATTAAGCCGTTTCGGCCCCGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGttt                                                                                                                                                   <  1:1380339/150‑1 (MQ=255)

    gatagatGGCATGTAAATTAAGCCGTTTCGGCCCCGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGtttt                                                                                                                                                  <  2:769870/150‑1 (MQ=255)

          tGGCATGTAAATTAAGCCGTTTCGGCCCCGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAAc                                                                                                                                            >  1:700461/1‑150 (MQ=255)

          tGGCATGTAAATTAAGCCGTTTCGGCCCCGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAAc                                                                                                                                            >  2:137371/1‑150 (MQ=255)

             cATGTAAATTAAGCCGTTTCGGCCCCGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCgg                                                                                                                                         >  1:1029828/1‑150 (MQ=255)

             cATGTAAATTAAGCCGTTTCGGCCCCGGAAGCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCgg                                                                                                                                         <  2:867769/150‑1 (MQ=255)

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

               tGTAAATTAAGCCGTTTCGGCCCCGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGaa                                                                                                                                       >  2:306046/1‑150 (MQ=255)

                gTAAATTAAGCGGTTTCGGCCCCGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACTTAGCTCACTGGCATTAGGCGCTTTTCCCGGATAATTGCCTGTGGCCTGAATCCCCGCTTTCAGCGAACCTTCCGGTTTTGAAAACCCGGaag                                                                                                                                      <  2:771673/150‑2 (MQ=255)

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

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

                gTAAATTAAGCCGTTTCGGCCCCGGAATCAGCCGCATAGCCTCTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGaaa                                                                                                                                      >  1:1064681/1‑150 (MQ=255)

                 tAAATTAAGCCGTTTCGGCCCCGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAAc                                                                                                                                     >  1:1591498/1‑150 (MQ=255)

                    aTTAAGCCGTTTCGGCCCCGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAAcatc                                                                                                                                  <  1:1501635/150‑1 (MQ=255)

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

                         gCCGTTTCGGCCCCGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGc                                                                                                                             <  2:399384/150‑1 (MQ=255)

                              ttCGGCCCCGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGc                                                                                                                        >  2:79621/1‑150 (MQ=255)

                                   ccccGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCcagca                                                                                                                   <  2:1237495/150‑1 (MQ=255)

                                     ccGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGt                                                                                                                 <  1:907055/150‑1 (MQ=255)

                                      cGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTg                                                                                                                >  1:29360/1‑150 (MQ=255)

                                      cGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTg                                                                                                                >  1:648422/1‑150 (MQ=255)

                                         aaTCAGCCGCATTGCCTGTTCCAGTTCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCAt                                                                                                             <  2:492629/150‑1 (MQ=255)

                                             aGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGaa                                                                                                         <  2:1235395/150‑1 (MQ=255)

                                              gCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGaaa                                                                                                        >  1:947516/1‑150 (MQ=255)

                                                cGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACg                                                                                                      <  2:477207/150‑1 (MQ=255)

                                                   aTAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGt                                                                                                   <  1:359261/150‑1 (MQ=255)

                                                   aTAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGt                                                                                                   <  1:930018/150‑1 (MQ=255)

                                                      gCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCTTCGCCCTGCCAGCAGTGCATTGAAACGGGTaaa                                                                                                >  2:1564595/1‑150 (MQ=255)

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

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

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

                                                            tCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCa                                                                                          <  2:907731/150‑1 (MQ=255)

                                                            tCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCa                                                                                          <  2:993368/150‑1 (MQ=255)

                                                             ccAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGTTTTGAAAACCCGGa                                                                                                                                        <  2:1235892/103‑1 (MQ=255)

                                                             ccAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGTTTTGAAAACCCGGa                                                                                                                                        >  1:1235892/1‑103 (MQ=255)

                                                             ccAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGTTTTGAAAACCCGGa                                                                                                                                        <  2:1241655/103‑1 (MQ=255)

                                                             ccAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGTTTTGAAAACCCGGa                                                                                                                                        >  1:1241655/1‑103 (MQ=255)

                                                             ccAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCaa                                                                                         <  1:821967/150‑1 (MQ=255)

                                                              cAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAg                                                                                        <  2:714717/150‑1 (MQ=255)

                                                                    ggCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGgcg                                                                                  >  1:589588/1‑150 (MQ=255)

                                                                      cACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCa                                                                                >  1:1426928/1‑150 (MQ=255)

                                                                              tCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGcctcc                                                                        <  1:773282/150‑1 (MQ=255)

                                                                              tCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGcctcc                                                                        <  2:401699/150‑1 (MQ=255)

                                                                               cACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGgc                                                                                   >  2:306714/1‑138 (MQ=255)

                                                                               cACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGgc                                                                                   <  1:306714/138‑1 (MQ=255)

                                                                               cACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGcctcct                                                                       >  1:726345/1‑150 (MQ=255)

                                                                               cACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGcctcct                                                                       >  1:1105936/1‑150 (MQ=255)

                                                                                 cTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCg                                                                     <  2:68894/150‑1 (MQ=255)

                                                                                   ggCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGAc                                                                   <  2:589588/150‑1 (MQ=255)

                                                                                    gCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACa                                                                  >  2:1477616/1‑150 (MQ=255)

                                                                                     cATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACAt                                                                 >  2:232527/1‑150 (MQ=255)

                                                                                      aTTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTACATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATc                                                                <  1:44619/150‑1 (MQ=255)

                                                                                       ttACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCa                                                               >  2:958265/1‑150 (MQ=255)

                                                                                          cgcgCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGcctcc                                                                        >  1:1500100/1‑138 (MQ=255)

                                                                                          cgcgCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGcctcc                                                                        <  2:1500100/138‑1 (MQ=255)

                                                                                           gcgcTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATcc                                                           <  1:1471383/150‑1 (MQ=255)

                                                                                             gcTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATcccc                                                         <  2:294852/150‑1 (MQ=255)

                                                                                                ttGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCAcc                                                      <  1:1535149/150‑1 (MQ=255)

                                                                                                   cccGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGgc                                                                                   >  2:223188/1‑118 (MQ=255)

                                                                                                   cccGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGgc                                                                                   <  1:223188/118‑1 (MQ=255)

                                                                                                   cccGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCAccgcc                                                   >  2:1550182/1‑150 (MQ=255)

                                                                                                     cGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACcgccgc                                                 >  1:711884/1‑150 (MQ=255)

                                                                                                     cGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACcgccgc                                                 >  1:1532112/1‑150 (MQ=255)

                                                                                                        aTAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGaa                                              <  1:1174548/150‑1 (MQ=255)

                                                                                                           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                                >  2:821493/1‑148 (MQ=255)

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

                                                                                                                       cTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATc                                                                >  1:71153/1‑117 (MQ=255)

                                                                                                                       cTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATc                                                                <  2:71153/117‑1 (MQ=255)

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

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

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

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

                                                                                                                         gAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACAGCGTTTCGcggcgg                             <  2:319304/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)

TCCAGATAGATGGCATGTAAATTAAGCCGTTTCGGCCCCGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCA  >  NC_000913/4095684‑4095981

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