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	97.4%	404.8 / 4.2	114	coding (10/1353 nt)	creD	inner membrane protein
Reads supporting (aligned to +/- strand): ref base . (3/0); new base C (55/56); total (58/56)
Fisher's exact test for biased strand distribution p-value = 2.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

Only 100 of 129 total aligned reads displayed.

GCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGC  >  NC_000913/4638088‑4638281

gCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAt                                                                                                  <  2:809716/100‑1 (MQ=255)

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

 ccTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑c                                                                                                <  1:176153/100‑1 (MQ=255)

    cGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑cccc                                                                                             >  2:2699459/1‑100 (MQ=255)

     gCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑ccccc                                                                                            >  1:515609/1‑100 (MQ=255)

      cTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑ccccc‑c                                                                                          >  1:2915869/1‑100 (MQ=255)

       ttCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑ccccc‑cc                                                                                         <  2:1446018/100‑2 (MQ=255)

       ttCGAATTCACCGTCACTTCACCTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑ccccc‑cc                                                                                         <  1:2515741/100‑2 (MQ=255)

        tCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑ccccc‑cct                                                                                        <  2:217740/100‑3 (MQ=255)

          gACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑ccccc‑cctgt                                                                                      >  2:2698352/1‑96 (MQ=255)

           aCTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑ccccc‑cctgtt                                                                                     <  1:638340/100‑6 (MQ=255)

            cTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTc                                                                                     >  1:498435/1‑100 (MQ=255)

            cTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTc                                                                                     >  1:361720/1‑100 (MQ=255)

            cTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTc                                                                                     <  1:620785/100‑1 (MQ=255)

            cTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTc                                                                                     <  2:1696666/100‑1 (MQ=255)

            cTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑CCCCCACTGTTc                                                                                     >  2:2754825/1‑100 (MQ=255)

             ttCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCt                                                                                    <  1:3069440/100‑1 (MQ=255)

             ttCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCt                                                                                    <  1:205689/100‑1 (MQ=255)

              tCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTg                                                                                   <  1:2493381/100‑1 (MQ=255)

              tCACCGGCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTg                                                                                   <  2:902967/100‑1 (MQ=255)

               cACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTgg                                                                                  >  1:1516216/1‑100 (MQ=255)

                aCCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTGGa                                                                                 <  2:3192951/100‑1 (MQ=255)

                 ccGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTGGaa                                                                                <  1:534381/100‑1 (MQ=255)

                   gTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTTTGGaaaa                                                                              <  1:1716502/100‑1 (MQ=255)

                   gTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTGGaaaa                                                                              >  1:3148927/1‑100 (MQ=255)

                   gTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTGGaaaa                                                                              <  1:1787491/100‑1 (MQ=255)

                    tCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTGGAAAAt                                                                             <  2:182044/100‑1 (MQ=255)

                     cACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTGGAAAATg                                                                            <  2:2563991/100‑1 (MQ=255)

                     cACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTGGAAAATg                                                                            <  2:2763774/100‑1 (MQ=255)

                     cACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTGGAAAATg                                                                            <  2:2850205/100‑1 (MQ=255)

                      aCTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTGGAAAATGa                                                                           >  1:876130/1‑100 (MQ=255)

                       cTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTGGAAAATGAc                                                                          <  2:125002/100‑1 (MQ=255)

                        ttCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTGGAAAATGACt                                                                         <  2:1637440/100‑1 (MQ=255)

                            caTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTGGAAAATGACTAGcc                                                                     <  1:1928173/100‑1 (MQ=255)

                            caTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTGGAAAATGACTAGcc                                                                     >  1:31861/1‑100 (MQ=255)

                            caTAGCTTCAAAGTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTGGAAAATGACTAGcc                                                                     >  1:2153869/1‑100 (MQ=255)

                                gCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGtt                                                                 >  2:1101062/1‑100 (MQ=255)

                                 cTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGttt                                                                >  1:241439/1‑100 (MQ=255)

                                   tCAAATTCTTCCCCCATAGTCTTCGTATCCTGCTGCCACTGCAAAGGAGAAGACTATGTTGCAATCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTgg                                                              >  2:2021725/1‑100 (MQ=255)

                                      aaTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGc                                                           <  2:280585/100‑1 (MQ=255)

                                      aaTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGc                                                           >  1:1648265/1‑100 (MQ=255)

                                       aTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCa                                                          <  1:461840/100‑1 (MQ=255)

                                       aTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCa                                                          >  2:2005871/1‑100 (MQ=255)

                                        ttcttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAg                                                         <  1:2145732/100‑1 (MQ=255)

                                        ttcttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAg                                                         <  2:1207078/100‑1 (MQ=255)

                                           ttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTAt                                                      >  2:544231/1‑100 (MQ=255)

                                             cccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTg                                                    <  1:2851125/100‑1 (MQ=255)

                                              ccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGc                                                   <  1:1047544/100‑1 (MQ=255)

                                              ccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGc                                                   <  2:1190062/100‑1 (MQ=255)

                                              ccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGc                                                   <  2:1565699/100‑1 (MQ=255)

                                               cacaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCt                                                  >  1:2181646/1‑100 (MQ=255)

                                                acaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtg                                                 >  2:2894152/1‑100 (MQ=255)

                                                acaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtg                                                 <  1:394380/100‑1 (MQ=255)

                                                 caTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgt                                                <  2:692899/100‑1 (MQ=255)

                                                  aTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgtt                                               <  2:673775/100‑1 (MQ=255)

                                                  aTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgtt                                               >  1:2575083/1‑100 (MQ=255)

                                                   tAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgttg                                              <  2:1738447/100‑1 (MQ=255)

                                                     gTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTgttgtt                                            >  2:3089377/1‑100 (MQ=255)

                                                     gTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTgttgtt                                            <  1:216254/100‑1 (MQ=255)

                                                     gTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTgttgtt                                            >  2:1826322/1‑100 (MQ=255)

                                                     gTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTgttgtt                                            <  1:2877840/100‑1 (MQ=255)

                                                     gTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTgttgtt                                            <  1:919985/100‑1 (MQ=255)

                                                     gTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGATGCAGTATTGCtgttgt                                             >  2:630467/1‑99 (MQ=255)

                                                      tCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGttgttg                                           <  1:2655857/100‑1 (MQ=255)

                                                          cGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTc                                       >  2:2448621/1‑100 (MQ=255)

                                                           gTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTcc                                      <  2:384520/100‑1 (MQ=255)

                                                           gTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTcc                                      <  1:2604129/100‑1 (MQ=255)

                                                             aTCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGa                                    <  1:978093/100‑1 (MQ=255)

                                                                ctgctgCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGAt                                   <  2:850457/98‑1 (MQ=255)

                                                                 tgctgcCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAAt                                >  2:2575390/1‑100 (MQ=255)

                                                                  gctgcCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATg                               <  1:2503459/100‑1 (MQ=255)

                                                                  gctgcCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATg                               >  2:426086/1‑100 (MQ=255)

                                                                  gctgcCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATg                               <  2:129296/100‑1 (MQ=255)

                                                                   ctgcCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGc                              >  1:1979257/1‑100 (MQ=255)

                                                                    tgcCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCt                             <  1:408288/100‑1 (MQ=255)

                                                                    tgcCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCt                             >  2:2679213/1‑100 (MQ=255)

                                                                     gcCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTg                            >  2:1112841/1‑100 (MQ=255)

                                                                      ccATTGCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGa                           >  2:2469657/1‑100 (MQ=255)

                                                                      ccATTGCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTGGAAAATGACTAGCCCGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGa                           >  2:630634/1‑100 (MQ=255)

                                                                      ccATTGCAAAGGAGAAGACTATGTTGAAAT‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGAt                          >  2:2440540/1‑100 (MQ=255)

                                                                        aTTGCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGAtt                         <  2:1945123/100‑1 (MQ=255)

                                                                         ttGCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTc                        >  1:2046435/1‑100 (MQ=255)

                                                                          tGCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCg                       <  2:1980736/100‑1 (MQ=255)

                                                                           gCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCgg                      >  1:2267572/1‑100 (MQ=255)

                                                                           gCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCgg                      >  2:504830/1‑100 (MQ=255)

                                                                           gCAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCgg                      >  2:2303672/1‑100 (MQ=255)

                                                                            cAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGc                     >  2:2807226/1‑100 (MQ=255)

                                                                            cAAAGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGc                     >  1:708685/1‑100 (MQ=255)

                                                                             aaaGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCa                    <  1:2807842/100‑1 (MQ=255)

                                                                             aaaGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCa                    <  1:1827687/100‑1 (MQ=255)

                                                                             aaaGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCa                    <  2:847238/100‑1 (MQ=255)

                                                                              aaGGAGAAGACTATGTTGAAATCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAg                   <  1:1102951/100‑1 (MQ=255)

                                                                                 gagaAGACTATGTTGAAATCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTg                >  2:2167977/1‑100 (MQ=255)

                                                                                     agaCTATGTTGAAATCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTg            <  2:627963/100‑1 (MQ=255)

                                                                                     agaCTATGTTGAAATCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTg            >  1:1243354/1‑100 (MQ=255)

                                                                                      gaCTATGTTGAAATCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGt           >  1:944679/1‑100 (MQ=255)

                                                                                        cTATGTTGAAATCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCg         >  2:109049/1‑100 (MQ=255)

                                                                                            gTTGAAATCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACg     >  2:3061370/1‑100 (MQ=255)

                                                                                            gTTGAAATCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACg     >  1:105600/1‑100 (MQ=255)

                                                                                               gAAATCCCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGc  >  1:1941758/1‑100 (MQ=255)

GCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑CCCCC‑CTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGC  >  NC_000913/4638088‑4638281

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