BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,095,942	T→C	6.9%	Q14R (CAG→CGG) ‡	rhaA ←	L‑rhamnose isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,095,942	0	T	C	6.9%	476.6 / 18.3	161	Q14R (CAG→CGG) ‡	rhaA	L‑rhamnose isomerase
Reads supporting (aligned to +/- strand): ref base T (70/79); new base C (5/6); total (76/85)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.17e-01

Only 100 of 163 total aligned reads displayed.

TGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTAAAGGTGGTCAGATTCGCGGTGGTGCTGACGA  >  NC_000913/4095793‑4096091

tGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCt                                                                                                                                                       >  2:2183438/1‑150 (MQ=255)

 gCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTg                                                                                                                                                      >  1:2209061/1‑150 (MQ=255)

    tgtgGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGttt                                                                                                                                                   >  2:1214231/1‑150 (MQ=255)

      tgGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCg                                                                                                                                                 >  2:1508442/1‑150 (MQ=255)

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

      tgGCCGGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCg                                                                                                                                                 <  2:958846/150‑1 (MQ=255)

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

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

         gcTGAATCCCCGCGGTCAGCGAACCTTCCGGTTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGATACGGGTAAAGGTTCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACAGCGTTTCGcgg                                                                                                                                              <  1:1052343/149‑3 (MQ=255)

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

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

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

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

         ccTGAATCCCCCCGGTCAGCGAACCTACCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACAGCGTTTCGcgg                                                                                                                                              <  1:576983/150‑3 (MQ=255)

          cTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAg                                                                                                                                             >  2:2233310/1‑150 (MQ=255)

          cTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACAGCGTTTCGcggc                                                                                                                                             <  2:766645/150‑4 (MQ=255)

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

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

              aTCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTcc                                                                                                                                         <  2:2049911/150‑1 (MQ=255)

               tCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTccc                                                                                                                                        <  2:330689/150‑1 (MQ=255)

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

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

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

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

                   cccGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGc                                                                                                                                    >  1:1538278/1‑150 (MQ=255)

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

                     cGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTGTCGCTAGCTCCCAGGCCt                                                                                                                                  <  1:981252/150‑1 (MQ=255)

                      ggTCAGCGAACCTTCCGGGTTTTCAAAAACGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAAAGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTCTTTCGCTAGCTCCCAGGCCTg                                                                                                                                 >  2:159958/1‑150 (MQ=255)

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

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

                             gatgtttCCGGTTTTGAAAACCCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGt                                                                                                                          >  2:1910167/6‑150 (MQ=255)

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

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

                               aCCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTg                                                                                                                        >  2:2075840/1‑150 (MQ=255)

                                gtttCCGGTTTTGAAAACCCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGa                                                                                                                       >  2:1729833/3‑150 (MQ=255)

                                ccTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGa                                                                                                                       >  2:67947/1‑150 (MQ=255)

                                     cGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGt                                                                                                                  <  1:2204497/150‑1 (MQ=255)

                                        gTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCAt                                                                                                               <  2:280903/150‑1 (MQ=255)

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

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

                                           ttCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCACGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATgcg                                                                                                            >  1:56267/1‑150 (MQ=255)

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

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

                                              aaaaCCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTGCAGTTGAGTGGTCATGCGCaa                                                                                                         >  2:408588/1‑150 (MQ=255)

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

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

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

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

                                                                gCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCCCCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCgt                                                                                       >  1:1613188/1‑150 (MQ=255)

                                                                    tGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTag                                                                                   <  1:506118/150‑1 (MQ=255)

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

                                                                       cagcagTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGt                                                                                >  2:1988906/1‑150 (MQ=255)

                                                                       cagcagTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGt                                                                                >  2:1395644/1‑150 (MQ=255)

                                                                                aTTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTgcg                                                                       >  1:513734/1‑150 (MQ=255)

                                                                                aTTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTgcg                                                                       <  1:1510892/150‑1 (MQ=255)

                                                                                   gAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGccc                                                                    >  2:446080/1‑149 (MQ=255)

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

                                                                                         ggTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGt                                                              <  1:1917977/150‑1 (MQ=255)

                                                                                         cgTAAACGATCAAGTTGGAGCAGCGCCTCCTCGACATAAATCCCCACCGCCGCCAAACGCAGATTCGCTAGCTCCCAGCCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGt                                                              >  2:1122678/2‑150 (MQ=255)

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

                                                                                          gTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTg                                                             >  1:2075089/1‑150 (MQ=255)

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

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

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

                                                                                               cGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTATGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCaa                                                        >  2:593246/1‑150 (MQ=255)

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

                                                                                                    aaGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCACGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCa                                                   <  2:56267/150‑1 (MQ=255)

                                                                                                     aGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTTAGTGGTCATGCGCAAAGCTACTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGTCAATTTCAc                                                  <  2:1288934/150‑1 (MQ=255)

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

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

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

                                                                                                              aGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGAt                                         <  1:2041010/150‑1 (MQ=255)

                                                                                                               gggcCTCCTCGTCATCTTTGCCCAGCTCCGCGAAACGCTGTTTCGCTAGGTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGAtt                                        <  1:2075840/148‑1 (MQ=255)

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

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

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

                                                                                                                         gACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTAAAg                              >  2:1830939/1‑150 (MQ=255)

                                                                                                                          aCATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTAAAgg                             <  2:841596/150‑1 (MQ=255)

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

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

                                                                                                                                aTCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTAAAGGTGGTCa                       <  1:41142/150‑1 (MQ=255)

                                                                                                                                aTCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTAAAGGTGGTCa                       <  1:1958650/150‑1 (MQ=255)

                                                                                                                                  ccccACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTAAAGGTGGTCAGa                     >  2:1456486/1‑150 (MQ=255)

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

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

                                                                                                                                       agctagcgAAACAGCGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTAAAGGTGGTCAGATTcgc                >  1:1764972/6‑150 (MQ=255)

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

                                                                                                                                        cgccgcGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTAAAGGTGGTCAGATTcgcg               <  2:665136/150‑1 (MQ=255)

                                                                                                                                         ctagcgAAACAGCGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTAAAGGTGGTCAGATTCGCgg              >  1:932509/4‑150 (MQ=255)

                                                                                                                                         cccgcGAAACGCTGTTTCGCTAGCTCCCATGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTAAAGGTGGTCAGATTCGCgg              <  2:1597355/149‑1 (MQ=255)

                                                                                                                                          tagcgAAACAGCGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTAAAGGTGGTCAGATTCGCggt             >  2:15560/3‑150 (MQ=255)

                                                                                                                                           agcgAAACAGCGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTAAAGGTGGTCAGATTTGCggtg            >  2:991122/2‑150 (MQ=255)

                                                                                                                                            gcgAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTAAAGGTGGTCAGATTCGCggtgg           <  1:756457/150‑1 (MQ=255)

                                                                                                                                              gAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTAAAGGTGGTCAGATTCGCGgtggtg         >  1:1648132/1‑150 (MQ=255)

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

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

                                                                                                                                                   gCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTAAAGGTGGTCAGATTCGCGGTGGTGCTgac    <  1:2018309/150‑1 (MQ=255)

                                                                                                                                                   gCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTAAAGGTGGTCAGATTCGCGGTGGTGCTgac    <  1:26966/150‑1 (MQ=255)

                                                                                                                                                    cTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTAAAGGTGGTCAGATTCGCGGTGGTGCTgacg   <  1:1988906/150‑1 (MQ=255)

                                                                                                                                                     tGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCGCATAGTGGGCAATTTCACTGTCAGGATTAGGGGTAAAGGTGGTCAGATTCGCGGTGGTGCTgacga  >  2:1313638/1‑150 (MQ=255)

TGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGCAAAGCTCCTTTGTCTGTCGTGTAGAGTGAATCTGCGCCACATAGTGGGCAATTTCACTGTCAGGATTAGGGGTAAAGGTGGTCAGATTCGCGGTGGTGCTGACGA  >  NC_000913/4095793‑4096091

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