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	96.8%	315.3 / 2.5	95	coding (10/1353 nt)	creD	inner membrane protein
Reads supporting (aligned to +/- strand): ref base . (1/2); new base C (42/50); total (43/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 106 total aligned reads displayed.

GCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACG  >  NC_000913/4638088‑4638278

gCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑t                                                                                              >  1:2188459/1‑100 (MQ=255)

  cTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑cc                                                                                           <  2:1800074/100‑1 (MQ=255)

   tCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑ccc                                                                                          <  2:778318/100‑1 (MQ=255)

    ttccgatctttCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑cccc                                                                                         <  1:2258831/91‑1 (MQ=255)

     gCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑ccccc                                                                                        >  1:2347788/1‑100 (MQ=255)

     gCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑ccccc                                                                                        >  2:3635/1‑100 (MQ=255)

       ttCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑ccccccc                                                                                      >  2:1472366/1‑99 (MQ=255)

       ttCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑ccccccc                                                                                      <  1:3635/100‑2 (MQ=255)

          gACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGATAAGACTATGTTGAAA‑T‑ccccccctgt                                                                                   <  2:1571831/100‑5 (MQ=255)

          gACTTCACCGTCACTTCACATAGCTCCAAATTCTTCCCACCTACCCTTCGTATCCTGCTGCCATTGCTAAGGAGAAGACTATGTTGAAA‑T‑ccccccctgt                                                                                   >  1:1163545/1‑96 (MQ=255)

            cTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTc                                                                                  >  1:2487237/1‑100 (MQ=255)

             ttCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCt                                                                                 <  1:962895/100‑1 (MQ=255)

             ttCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑ccccagatcggaa                                                                                >  2:2258846/1‑91 (MQ=255)

               cACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTgg                                                                               <  2:1694026/100‑1 (MQ=255)

                    tCACTTCACATAGCTTCAAATTCTTCCCACCTAGTCTCCGCATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGaaa                                                                            <  1:1817088/98‑1 (MQ=255)

                    tCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAAt                                                                          <  2:2544134/100‑1 (MQ=255)

                      aCTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGa                                                                        >  1:1771396/1‑100 (MQ=255)

                      aCTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGa                                                                        >  1:699972/1‑100 (MQ=255)

                       cTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGAc                                                                       <  1:739586/100‑1 (MQ=255)

                         tCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTa                                                                     <  1:1646444/100‑1 (MQ=255)

                         tCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTa                                                                     <  2:8648/100‑1 (MQ=255)

                          cacaTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAg                                                                    <  1:586539/100‑1 (MQ=255)

                            caTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGcc                                                                  <  1:263277/100‑1 (MQ=255)

                             aTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCt                                                                 >  2:1626630/1‑100 (MQ=255)

                              tAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCCTTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTg                                                                >  1:1271226/1‑100 (MQ=255)

                                gCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGc                                                                   >  1:1253610/1‑95 (MQ=255)

                                gCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGtt                                                              >  2:1372130/1‑100 (MQ=255)

                                gCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGtt                                                              >  2:1276962/1‑100 (MQ=255)

                                  ttCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTg                                                            >  2:1315330/1‑100 (MQ=255)

                                   tCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTgg                                                           <  2:1379387/100‑1 (MQ=255)

                                   tCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTgg                                                           >  2:1371342/1‑100 (MQ=255)

                                   acAAATTCTTGCCACATAGACCTCGTATCTTGCTGCGATTGCAAAGAAGAAGACTATCTTAAAATT‑CCCCCGTATTCTGGAAAATGACTAGCCTGt                                                               >  1:240600/2‑96 (MQ=255)

                                     aaaTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTg                                                         <  2:1638917/100‑1 (MQ=255)

                                     aaaTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTg                                                         >  1:465879/1‑100 (MQ=255)

                                      aaTTCTTCCCTGATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACGATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGGCTAGCCTGTTTGGTg                                                         <  1:2338333/99‑1 (MQ=255)

                                      aaTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGc                                                        <  1:2321435/100‑1 (MQ=255)

                                        ccgatctcacaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAg                                                      <  2:457352/93‑1 (MQ=255)

                                         tcttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGt                                                     <  2:702929/100‑1 (MQ=255)

                                         tcttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGt                                                     >  2:1087262/1‑100 (MQ=255)

                                          cttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTAt                                                   <  2:96136/100‑1 (MQ=255)

                                           ttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCGTGTTTGGTGCAGTAt                                                   >  1:2350229/1‑100 (MQ=255)

                                           ttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTAt                                                   <  2:2487256/100‑1 (MQ=255)

                                           ttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTAt                                                   <  1:2080234/100‑1 (MQ=255)

                                              ccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGc                                                <  1:1565516/100‑1 (MQ=255)

                                              ccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGc                                                <  1:1147759/100‑1 (MQ=255)

                                              ccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGCTTGGTGCAGTATTGc                                                >  1:20884/1‑100 (MQ=255)

                                               cacaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAgagatcgg                                               >  1:457350/1‑93 (MQ=255)

                                                acaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtg                                              <  2:544243/100‑1 (MQ=255)

                                                   tAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgttg                                           >  1:158558/1‑100 (MQ=255)

                                                    aGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgttgt                                          <  1:1985615/100‑1 (MQ=255)

                                                     gTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTgttgtt                                         >  1:538279/1‑100 (MQ=255)

                                                      tCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGttgttg                                        <  1:99316/100‑1 (MQ=255)

                                                      tCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGttgttg                                        >  2:1348776/1‑100 (MQ=255)

                                                      tCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGttgttg                                        <  2:1262924/100‑1 (MQ=255)

                                                       cTTCGTATCCTGCTGCCATTGCAAAGGAGAAGCCTACGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTAATTGGTGCAGTATTGCGgttgtt                                         >  1:2311474/1‑98 (MQ=255)

                                                       cTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGa                                       <  2:143629/100‑1 (MQ=255)

                                                       cTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGa                                       <  2:1174668/100‑1 (MQ=255)

                                                        ttCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGAt                                      <  1:28235/100‑1 (MQ=255)

                                                        ttCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCCGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGAt                                      <  1:146427/100‑1 (MQ=255)

                                                         tCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGAtt                                     >  1:619507/1‑100 (MQ=255)

                                                           gTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTcc                                   >  2:2251636/1‑100 (MQ=255)

                                                           gTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTcc                                   <  2:942427/100‑1 (MQ=255)

                                                            tATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTcc                                   <  1:1573660/99‑1 (MQ=255)

                                                             aTCCTGCTGCCATTGAAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTTGAAGATGACTAGCCTGTTTGGTGCAGTATTTTTGTTGGTGATTCCGa                                 <  1:414435/100‑1 (MQ=255)

                                                              tCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGAt                                <  2:1337296/100‑1 (MQ=255)

                                                              tCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGAt                                <  1:1001059/100‑1 (MQ=255)

                                                              tCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGAt                                <  2:487154/100‑1 (MQ=255)

                                                              tCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGAt                                <  2:17089/100‑1 (MQ=255)

                                                                ctgctgCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATaa                              >  1:2026013/1‑100 (MQ=255)

                                                                ctgctgCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATaa                              <  2:1603960/100‑1 (MQ=255)

                                                                  gctgcCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATg                            <  2:1230131/100‑1 (MQ=255)

                                                                  gctgcCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATg                            >  1:339396/1‑100 (MQ=255)

                                                                   ctgcCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGc                           <  2:2585907/100‑1 (MQ=255)

                                                                   ctgcCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGCGCAGTATTGCTGTTGTTGATTCCGATAATGc                           <  1:765903/100‑1 (MQ=255)

                                                                     gcCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTg                         <  1:2251621/100‑1 (MQ=255)

                                                                       cATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGAt                       <  1:381284/100‑1 (MQ=255)

                                                                       cATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGAt                       >  2:2285573/1‑100 (MQ=255)

                                                                          tGCAAAGGAGAAGATTATGTTGAAATT‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCg                    <  1:186823/100‑1 (MQ=255)

                                                                          tGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCg                    >  2:2255034/1‑100 (MQ=255)

                                                                           gCAAAGGAGAAGACTCTGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTAGGTTCAGTATTGCTGTTGTTGATACCGATAATGCTGttt                      >  1:533605/1‑97 (MQ=255)

                                                                           gCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCgg                   >  2:870746/1‑100 (MQ=255)

                                                                           gCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCgg                   >  1:2083747/1‑100 (MQ=255)

                                                                            cAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGc                  >  1:2578187/1‑100 (MQ=255)

                                                                            cAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGc                  <  2:2395264/100‑1 (MQ=255)

                                                                            cAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGc                  >  1:919520/1‑100 (MQ=255)

                                                                             aaaGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCa                 <  2:1255848/100‑1 (MQ=255)

                                                                              aaGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAg                <  1:1372122/100‑1 (MQ=255)

                                                                                ggAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTAGTGCTGTTGTTGATTCCGATAATGCTGAGTCGGCAGGt              <  1:1768759/100‑1 (MQ=255)

                                                                                       aCTGTGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTc       <  1:2149799/100‑1 (MQ=255)

                                                                                       aCTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTc       >  1:2122631/1‑100 (MQ=255)

                                                                                        cTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCg      >  2:776761/1‑100 (MQ=255)

                                                                                        cTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCg      >  2:973448/1‑100 (MQ=255)

                                                                                        cTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCg      >  2:1432350/1‑100 (MQ=255)

                                                                                        cTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCg      >  1:1013775/1‑100 (MQ=255)

                                                                                          aTGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGaa    >  2:1538243/1‑100 (MQ=255)

                                                                                           tGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAAc   >  2:609763/1‑100 (MQ=255)

                                                                                           tGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAAc   >  2:280255/1‑100 (MQ=255)

                                                                                           tGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAAc   <  2:206695/100‑1 (MQ=255)

                                                                                           tGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAAc   >  1:1913210/1‑100 (MQ=255)

                                                                                            gTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACg  <  2:1825598/100‑1 (MQ=255)

GCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACG  >  NC_000913/4638088‑4638278

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