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.9%	445.9 / 7.6	127	coding (10/1353 nt)	creD	inner membrane protein
Reads supporting (aligned to +/- strand): ref base . (2/2); new base C (60/63); total (62/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.68e-01

Only 100 of 144 total aligned reads displayed.

CCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑CCC‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGC  >  NC_000913/4638089‑4638281

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

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

 cTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAGAT‑cg                                                                                                >  2:3021355/1‑99 (MQ=255)

 cTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑cc                                                                                                <  1:749345/100‑1 (MQ=255)

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

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

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

     cTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑ccc‑ccc                                                                                           >  2:1959655/1‑100 (MQ=255)

      ttCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑ccc‑cccc                                                                                          <  1:2149075/100‑2 (MQ=255)

       tCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑ccc‑cccct                                                                                         >  2:1949559/1‑98 (MQ=255)

         gACTTCACCGTCACTTCACACAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑ccc‑cccctgt                                                                                       >  2:3250660/1‑96 (MQ=255)

          aCTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑ccc‑cccctgtt                                                                                      >  1:3424210/1‑95 (MQ=255)

           cTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCC‑TGTTc                                                                                     >  1:369569/1‑100 (MQ=255)

           cTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCC‑TGTTc                                                                                     >  1:703926/1‑100 (MQ=255)

             tCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCC‑TGTTCTg                                                                                   <  2:1881549/100‑1 (MQ=255)

             tCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCC‑TGTTCTg                                                                                   <  2:1120136/100‑1 (MQ=255)

              cACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCC‑TGTTCTgg                                                                                  <  1:1581533/100‑1 (MQ=255)

              cACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCC‑TGTTCTgg                                                                                  <  1:954671/100‑1 (MQ=255)

               aCCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCC‑TGTTCTGGa                                                                                 >  1:2848957/1‑100 (MQ=255)

               aCCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCC‑TGTTCTGGa                                                                                 <  1:1624533/100‑1 (MQ=255)

                  gTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCC‑TGTTCTGGaaaa                                                                              <  1:1959639/100‑1 (MQ=255)

                  gTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCC‑TGTTCTGGaaaa                                                                              >  2:1601559/1‑100 (MQ=255)

                  gTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCC‑TGTTCTGGaa                                                                                <  2:1379584/98‑1 (MQ=255)

                   tCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCC‑TGTTCTGGAAAAt                                                                             <  2:1252182/100‑1 (MQ=255)

                    cACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCC‑TGTTCTGGAAAATg                                                                            <  1:1125883/100‑1 (MQ=255)

                    cACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑CCC‑CCC‑TGTTCTGGAAAATGa                                                                           >  2:2033025/1‑100 (MQ=255)

                     aCTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCC‑TGTTCTGGAAAATGa                                                                           >  1:1746371/1‑100 (MQ=255)

                     aCTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCC‑TGTTCTGGAAAATGa                                                                           >  2:476817/1‑100 (MQ=255)

                        tCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCC‑TGTTCTGGAAAATGACTa                                                                        >  1:2201199/1‑100 (MQ=255)

                         cacaTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCC‑TGTTCTGGAAAATGACTAg                                                                       >  1:2062466/1‑100 (MQ=255)

                          acaTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCC‑TGTTCTGGAAAATGACTAGc                                                                      <  1:1240054/100‑1 (MQ=255)

                             tAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCC‑TGTTCTGGAAAATGACTAGCCTg                                                                   <  2:3310732/100‑1 (MQ=255)

                               gCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCC‑TGTTCTGGAAAATGACTAGCCTGtt                                                                 >  1:2294638/1‑100 (MQ=255)

                                cTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCC‑TGTTCTGGAAAATGACTAGCCTGttt                                                                >  2:2117377/1‑100 (MQ=255)

                                 ttCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTg                                                               <  2:1991302/100‑1 (MQ=255)

                                  tCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCC‑TGTTATGGAAAATGACTAGACTgcctgg                                                              >  2:635137/1‑95 (MQ=255)

                                    aaaTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTg                                                            >  2:914808/1‑100 (MQ=255)

                                    aaaTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTg                                                            <  2:8824/100‑1 (MQ=255)

                                        tcttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGt                                                        >  2:2262751/1‑100 (MQ=255)

                                         cttcCCACATAGTATTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑CCCGCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAg                                                         >  1:1519001/1‑98 (MQ=255)

                                          ttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCC‑TGTTCTGGCAAATGACTAGCCTGTTTGGTGCAGTAt                                                      >  1:2581901/1‑100 (MQ=255)

                                            cccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTg                                                    <  1:1443868/100‑1 (MQ=255)

                                             ccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGc                                                   <  2:1708707/100‑1 (MQ=255)

                                             ccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGc                                                   <  2:3208264/100‑1 (MQ=255)

                                              cacaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTTCt                                                  >  1:1481089/1‑100 (MQ=255)

                                              cacaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCt                                                  >  1:1834450/1‑100 (MQ=255)

                                              cacaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCt                                                  >  2:299028/1‑100 (MQ=255)

                                               acaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑CCC‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgt                                                <  2:2252284/100‑1 (MQ=255)

                                                caTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgt                                                >  2:2533177/1‑100 (MQ=255)

                                                    gTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTgttgtt                                            >  2:1884734/1‑100 (MQ=255)

                                                    gTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTgttgtt                                            >  2:235463/1‑100 (MQ=255)

                                                     tCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGttgttg                                           <  2:1200639/100‑1 (MQ=255)

                                                     tCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGttgttg                                           >  2:2735447/1‑100 (MQ=255)

                                                     tCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGttgttg                                           <  2:721053/100‑1 (MQ=255)

                                                       ttCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGAt                                         <  1:561447/100‑1 (MQ=255)

                                                        tCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGAtt                                        <  1:3207496/100‑1 (MQ=255)

                                                        tCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGAtt                                        <  1:2487343/100‑1 (MQ=255)

                                                          gTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTcc                                      <  2:2839345/100‑1 (MQ=255)

                                                          gTATCCTGCGGCCATTGCAAAGGATAAGTGGATGTTGAAATCCCC‑CCC‑TGTTCTGGAAAATGACTAGCCGTTTTGGTGCAGTATTGCTGTTGTTGATTcc                                      <  1:3315900/100‑1 (MQ=255)

                                                           tATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCC‑GGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATGGCTGTTGTTGATTCCg                                     >  2:555179/1‑100 (MQ=255)

                                                              ccTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATa                                  <  1:3370300/100‑1 (MQ=255)

                                                               ctgctgCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATaa                                 >  1:905802/1‑100 (MQ=255)

                                                                  ctgcCATTGCAAAGCAGAAGACTATGTTGAGAT‑CCC‑CCCTTGTTCTGGAAAATGACTAGCCTGTTTTGTGCAGTATTGCTGTTGTTGATTCCTATAATGc                              <  1:1703408/100‑1 (MQ=255)

                                                                   tgcCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCt                             >  1:2561434/1‑100 (MQ=255)

                                                                    gcCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTg                            <  1:2117360/100‑1 (MQ=255)

                                                                     ccATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGa                           >  2:397106/1‑100 (MQ=255)

                                                                      cATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGAt                          >  2:3044046/1‑100 (MQ=255)

                                                                      cATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGAt                          <  1:1134155/100‑1 (MQ=255)

                                                                      cATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCC‑TGTTCTGGAAAATGACTAGCATGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGAt                          >  1:2214931/1‑100 (MQ=255)

                                                                         tGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCg                       <  1:2735426/100‑1 (MQ=255)

                                                                          gCAAAGGAGAAGACTATGTTGAAATCCCC‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCgg                      >  2:2790118/1‑100 (MQ=255)

                                                                          gCAAAGGAGAAGACTATGTTGAAATCCCC‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCgg                      >  2:325165/1‑100 (MQ=255)

                                                                           cAAAGGAGAAGACTATGTTGAAATCCCC‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGc                     >  1:250983/1‑100 (MQ=255)

                                                                           cAAAGGAGAAGACTATGTTGAAATCCCC‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGATGATTTgg                      >  1:3381250/1‑99 (MQ=255)

                                                                            aaaGGAGAAGACTATGTTGAAATCCCC‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCa                    <  2:2294656/100‑1 (MQ=255)

                                                                            aaaGGAGAAGACTATGTTGAAATCCCC‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCa                    <  1:914801/100‑1 (MQ=255)

                                                                             aaGGAGAAGACTATGTTGAAATCCCC‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAg                   >  2:3146699/1‑100 (MQ=255)

                                                                             aaGGAGAAGACTATGTTGAAATCCCC‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAg                   <  1:1854534/100‑1 (MQ=255)

                                                                              aGGAGAAGACTATGTTGAAATCCCC‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAgg                  >  1:1002913/1‑100 (MQ=255)

                                                                              aGGAGAAGACTATGTTGAAATCCCC‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAgg                  >  1:3229917/1‑100 (MQ=255)

                                                                              aGGAGAAGACTATGTTGAAATCCCC‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAgg                  <  1:3250636/100‑1 (MQ=255)

                                                                               ggAGAAGACTATGTTGAAATCCCC‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGt                 >  2:1608951/1‑100 (MQ=255)

                                                                               ggAGAAGACTATGTTGAAATCCCC‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGt                 >  2:2554595/1‑100 (MQ=255)

                                                                                gagaAGACTATGTTGAAATCCCC‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTg                <  1:3309708/100‑1 (MQ=255)

                                                                                  gaagaCTATGTTGAAATCCCC‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGAt              <  2:294099/100‑1 (MQ=255)

                                                                                  gaagaCTATGTTGAAATCCCC‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGAt              >  2:1733853/1‑100 (MQ=255)

                                                                                  gaagaCTATGTTGAAATCCCC‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGAt              >  2:2374663/1‑100 (MQ=255)

                                                                                    agaCTATGTTGAAATCCCC‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGAtt             >  1:2916081/1‑99 (MQ=255)

                                                                                    agaCTATGTTGAAATCCCC‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTg            <  2:369572/100‑1 (MQ=255)

                                                                                    agaCTATGTTGAAATCCCC‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTg            >  2:1562726/1‑100 (MQ=255)

                                                                                     gaCTATGTTGAAATCCCC‑CCC‑TGTTCTGGAAAAGGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGACGGTGATTGt           >  1:2265463/1‑100 (MQ=255)

                                                                                      aCTATGTTGAAATCCCC‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTc          <  1:973347/100‑1 (MQ=255)

                                                                                      aCTATGTTGAAATCCCC‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTc          >  1:2624584/1‑100 (MQ=255)

                                                                                       cTATGTTGAAATCCCC‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCg         >  1:1945546/1‑100 (MQ=255)

                                                                                           gTTGAAATCCCC‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACg     >  1:1264019/1‑100 (MQ=255)

                                                                                           gTTGAAATCCCC‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACg     >  1:3220394/1‑100 (MQ=255)

                                                                                           gTTGAAATCCCC‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACg     >  2:1184605/1‑100 (MQ=255)

                                                                                             tGAAATCCCC‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTg   <  1:2468319/100‑1 (MQ=255)

                                                                                              gAAATCCCC‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGc  >  2:799874/1‑100 (MQ=255)

                                                                                              gAAATCCCC‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGc  <  1:235461/100‑1 (MQ=255)

CCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑CCC‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGC  >  NC_000913/4638089‑4638281

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