BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,095,842	A→C	8.2%	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	8.2%	276.0 / 16.4	98	G47G (GGT→GGG)	rhaA	L‑rhamnose isomerase
Reads supporting (aligned to +/- strand): ref base A (43/47); new base C (5/3); total (48/50)
Fisher's exact test for biased strand distribution p-value = 4.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.01e-01

Only 100 of 102 total aligned reads displayed.

TGGCATGTAAATTAAGCCGTTTCGGCCCCGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGC  >  NC_000913/4095694‑4095991

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:399384/150‑1 (MQ=255)

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

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

                         ccccGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCcagca                                                                                                                             <  2:1237495/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:993368/150‑1 (MQ=255)

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

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

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

                                                   ccAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGTTTTGAAAACCCGGa                                                                                                                                                  >  1:1235892/1‑103 (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                                                                                  <  2:401699/150‑1 (MQ=255)

                                                                    tCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGcctcc                                                                                  <  1:773282/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                                                                                  <  2:1500100/138‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

                                                                                           cGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACcgccgc                                                           >  1:711884/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                                          <  1:821493/150‑3 (MQ=255)

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

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

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

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

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

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

                                                                                                               gAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACAGCGTTTCGcggcgg                                       <  2: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                                                     >  2:378378/1‑127 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                                                                                                                cAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGaaa                                                       >  2:1009935/1‑101 (MQ=255)

                                                                                                                                                cAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGaaa                                                       <  1:1009935/101‑1 (MQ=255)

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

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

                                                                                                                                                    aCCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGc  <  1:1477616/150‑1 (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` `< 32 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`