BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,638,193:1	(C)_6→7	100%	coding (16/1353 nt)	creD →	inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,638,187	1	.	C	95.9%	499.7 / 13.6	146	coding (10/1353 nt)	creD	inner membrane protein
Reads supporting (aligned to +/- strand): ref base . (3/3); new base C (78/61); total (81/65)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 162 total aligned reads displayed.

GCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑CC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGC  >  NC_000913/4638088‑4638281

 ccTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑c                                                                                                >  2:3786808/1‑100 (MQ=255)

  cTCGCTTCGACTGCACCGGCACTTCACATAGCTGCAAATTCTTCCCACATAGTCTTCGTATCGTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑cc                                                                                               <  1:221220/100‑1 (MQ=255)

    cGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑cc‑cc                                                                                            <  2:2871721/100‑1 (MQ=255)

     gCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑cc‑ccc                                                                                           <  1:601773/100‑1 (MQ=255)

     gCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑cc‑ccc                                                                                           <  1:423126/100‑1 (MQ=255)

     gCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑cc‑ccc                                                                                           >  2:1101118/1‑100 (MQ=255)

     gCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑cc‑cc                                                                                            <  1:661883/99‑1 (MQ=255)

      cTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑cc‑cccc                                                                                          >  1:1548891/1‑100 (MQ=255)

          gACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCACAGGAGAAGACTATGTTGAAAT‑cc‑ccccctgt                                                                                      >  1:3320407/1‑96 (MQ=255)

          gACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑cc‑ccccctgt                                                                                      <  2:514484/100‑5 (MQ=255)

            cTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTc                                                                                     >  2:1689350/1‑100 (MQ=255)

              tCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTg                                                                                   <  1:2781265/100‑1 (MQ=255)

              tCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTg                                                                                   >  1:64733/1‑100 (MQ=255)

               cACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTgg                                                                                  >  2:3340912/1‑100 (MQ=255)

                 ccGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGaa                                                                                >  2:599019/1‑100 (MQ=255)

                 ccGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGaa                                                                                >  2:3508670/1‑100 (MQ=255)

                    ttaCTTCACATAGATTCAAATTCTTCCCACATAGTCTTCGTACCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAAt                                                                             >  2:2571138/3‑100 (MQ=255)

                    tCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAAt                                                                             <  1:3497067/100‑1 (MQ=255)

                    tCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAAt                                                                             <  2:3419961/100‑1 (MQ=255)

                       cTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGAc                                                                          <  2:3373384/100‑1 (MQ=255)

                       cTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGAc                                                                          >  1:1406212/1‑100 (MQ=255)

                          cacaTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAg                                                                       >  1:3607455/1‑100 (MQ=255)

                            caTAGCTTCAAATTCTTCCCACATAGTCCTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGcc                                                                     <  2:584559/100‑1 (MQ=255)

                             aTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTTTGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCt                                                                    <  1:2334830/100‑1 (MQ=255)

                                gCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGtt                                                                 >  1:1337597/1‑100 (MQ=255)

                                gCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGtt                                                                 >  2:1967840/1‑100 (MQ=255)

                                 cTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTCTGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGt                                                                  >  2:3117458/1‑98 (MQ=255)

                                 cTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGttt                                                                >  1:1211686/1‑100 (MQ=255)

                                 cTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGttt                                                                >  2:2764689/1‑100 (MQ=255)

                                 cTTAAAATTCTTCCCACTTACTCCTCATATCCTGCTGCCATTGAAGAGGAGAAGACTATGTCGAAATCCC‑CCCCTGTTCTGGAAAACGACTAGCCt                                                                    >  1:487552/1‑96 (MQ=255)

                                  ttCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTg                                                               <  2:3796839/100‑1 (MQ=255)

                                   tCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTgg                                                              >  2:2900568/1‑100 (MQ=255)

                                      aaTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGc                                                           <  1:152499/100‑1 (MQ=255)

                                        ttcttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAg                                                         <  2:1029284/100‑1 (MQ=255)

                                          cttcCCACATAGTCTTCGTATCCTGCTGCCAGTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTa                                                       <  1:1819459/100‑1 (MQ=255)

                                           ttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTAt                                                      >  2:3295019/1‑100 (MQ=255)

                                           ttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTAt                                                      >  2:3164769/1‑100 (MQ=255)

                                           ttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTAt                                                      >  1:2516143/1‑100 (MQ=255)

                                           ttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTAt                                                      >  1:2378375/1‑100 (MQ=255)

                                             cccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTg                                                    <  2:3289188/100‑1 (MQ=255)

                                             cccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTg                                                    >  2:544954/1‑100 (MQ=255)

                                             cccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTg                                                    <  1:2203298/100‑1 (MQ=255)

                                              ccACATAGTCTTCGTATCCTGCTGCCATTGCACAGGAGATGAGTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGc                                                   <  1:2669231/100‑1 (MQ=255)

                                              ccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGc                                                   >  1:3697596/1‑100 (MQ=255)

                                              ccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGc                                                   >  2:2361483/1‑100 (MQ=255)

                                              ccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGc                                                   <  1:2136225/100‑1 (MQ=255)

                                               cacaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCt                                                  >  2:978515/1‑100 (MQ=255)

                                               cacaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCt                                                  <  1:3313427/100‑1 (MQ=255)

                                               cacaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCt                                                  >  2:547229/1‑100 (MQ=255)

                                                acaTAGTCTTCGTATGCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtg                                                 <  2:3743437/100‑1 (MQ=255)

                                                acaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtg                                                 <  1:290439/100‑1 (MQ=255)

                                                acaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtg                                                 <  2:1236775/100‑1 (MQ=255)

                                                  aTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgtt                                               >  2:1535432/1‑100 (MQ=255)

                                                   tAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgttg                                              <  1:993650/100‑1 (MQ=255)

                                                    aGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgttgt                                             <  2:2721331/100‑1 (MQ=255)

                                                     gTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTgttgtt                                            <  1:3505998/100‑1 (MQ=255)

                                                      tCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGttgttg                                           <  1:1409660/100‑1 (MQ=255)

                                                       cTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGa                                          >  2:3723552/1‑100 (MQ=255)

                                                       cTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGa                                          >  2:1732587/1‑100 (MQ=255)

                                                         tCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGAtt                                        >  1:2568563/1‑100 (MQ=255)

                                                         acGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑CCACCCCTGATCTGGACAATGACTAGACGGTGTGGTGCAGTAGTGCTGttgttggt                                         >  1:2525804/2‑97 (MQ=255)

                                                           gtttCCAGCTGGCATTGCAAAGGAGAAGACTCTGTTGAAATCCC‑CCCCTGTTCTGGAAAATGCCTAGCCTGTTTGGTGCAGTATTGCCGTTGTTGATTcc                                      <  1:943233/97‑1 (MQ=255)

                                                            tATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCg                                     >  2:1353528/1‑100 (MQ=255)

                                                             aTCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGa                                    <  1:544951/100‑1 (MQ=255)

                                                               ccTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATa                                  >  2:1512547/1‑100 (MQ=255)

                                                               ccTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATa                                  <  2:401847/100‑1 (MQ=255)

                                                               ccTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGGTCTGGAAAATGACTATCCTTTTTGGTGCAGTATTGCTGTTGTTGATTCCGATa                                  >  1:3222307/1‑100 (MQ=255)

                                                                ctgctgCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTGTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATa                                  >  2:1796181/1‑99 (MQ=255)

                                                                ctgctgCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATTACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATaa                                 >  1:2424695/1‑100 (MQ=255)

                                                                ctgctgCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATaa                                 >  1:802384/1‑100 (MQ=255)

                                                                 tgctgcCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAAt                                >  2:1569554/1‑100 (MQ=255)

                                                                 tgctgcCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAAt                                <  1:1103181/100‑1 (MQ=255)

                                                                 tgctgcCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAAt                                >  2:792463/1‑100 (MQ=255)

                                                                  gctgcCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATg                               >  1:801751/1‑100 (MQ=255)

                                                                  gctgcCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATg                               <  2:3536585/100‑1 (MQ=255)

                                                                    tgcCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCt                             >  2:686811/1‑100 (MQ=255)

                                                                    tgcCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCt                             <  1:1252921/100‑1 (MQ=255)

                                                                           gCAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCgg                      >  1:3079303/1‑100 (MQ=255)

                                                                            cAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGc                     >  2:2788155/1‑100 (MQ=255)

                                                                            cAAAGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGc                     >  1:3008411/1‑100 (MQ=255)

                                                                             aaaGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCa                    <  2:2899864/100‑1 (MQ=255)

                                                                              aaGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTATCCTGTTTGTTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAg                   <  2:929343/100‑1 (MQ=255)

                                                                              aaGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAg                   >  2:2285417/1‑100 (MQ=255)

                                                                              aaGGAGAAGACTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAg                   >  1:3470596/1‑100 (MQ=255)

                                                                                   gaagaCTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGAt              >  1:828861/1‑100 (MQ=255)

                                                                                    aagaCTATGTTGAAATCCC‑CCCCTGTTTTGCAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCTGATAATGCTGATTCGGCAGGTGa               <  1:2798262/98‑1 (MQ=255)

                                                                                    aagaCTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGAtt             <  2:1139794/100‑1 (MQ=255)

                                                                                     agaCTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCCGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTg            >  1:3155487/1‑100 (MQ=255)

                                                                                     agaCTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTg            <  1:1072046/100‑1 (MQ=255)

                                                                                        cTATGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCg         >  1:2627386/1‑100 (MQ=255)

                                                                                          aTGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGaa       >  1:1698333/1‑100 (MQ=255)

                                                                                           tGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAAc      >  1:532528/1‑100 (MQ=255)

                                                                                           tGTTGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAAc      >  2:2256199/1‑100 (MQ=255)

                                                                                             ttGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGt    >  2:2998443/1‑100 (MQ=255)

                                                                                             ttGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGt    >  2:199464/1‑100 (MQ=255)

                                                                                              tGAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGCGATTGTCGAACGTg   >  1:492769/1‑100 (MQ=255)

                                                                                               gAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGc  <  1:2983640/100‑1 (MQ=255)

                                                                                               gAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGc  >  2:421358/1‑100 (MQ=255)

                                                                                               gAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGc  <  1:1731498/100‑1 (MQ=255)

                                                                                               gAAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGc  <  1:686807/100‑1 (MQ=255)

GCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑CC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGC  >  NC_000913/4638088‑4638281

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