BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,638,193	(C)_6→7	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.4%	368.5 / 9.9	109	coding (10/1353 nt)	creD	inner membrane protein
Reads supporting (aligned to +/- strand): ref base . (2/3); new base C (62/42); total (64/45)
Fisher's exact test for biased strand distribution p-value = 6.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 120 total aligned reads displayed.

GCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑CCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGC  >  NC_000913/4638088‑4638281

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

 ncTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCCTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑c                                                                                                >  1:991475/2‑100 (MQ=255)

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

   tCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑ccc                                                                                              >  1:1461354/1‑100 (MQ=255)

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

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

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

         cGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑ccc‑cccctg                                                                                       >  2:1309542/1‑97 (MQ=255)

          gACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑ccc‑cccctgt                                                                                      >  1:2270086/1‑96 (MQ=255)

          gACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑ccc‑cccctgt                                                                                      >  1:219700/1‑96 (MQ=255)

             ttCACCGTCACTGCACATAGCTTCAAATTCTTCCCACCTAGTCTTCGTCTCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTc                                                                                     >  2:2702869/1‑99 (MQ=255)

              tCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTg                                                                                   <  1:2046994/100‑1 (MQ=255)

               cACCTTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTgg                                                                                  <  2:2667884/100‑1 (MQ=255)

               cACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTgg                                                                                  >  1:892382/1‑100 (MQ=255)

               cACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑CCC‑CCCTGTTCTGGa                                                                                 >  1:1749723/1‑100 (MQ=255)

                aCCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGa                                                                                 <  2:1973093/100‑1 (MQ=255)

                aCCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTCCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGa                                                                                 <  2:1493728/100‑1 (MQ=255)

                  cGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGaaa                                                                               >  1:2336489/1‑100 (MQ=255)

                   gTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGaaaa                                                                              >  1:1728732/1‑100 (MQ=255)

                   gTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGaaaa                                                                              <  2:3012974/100‑1 (MQ=255)

                     cACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATg                                                                            >  2:264115/1‑100 (MQ=255)

                     cACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATg                                                                            >  2:484820/1‑100 (MQ=255)

                      aCTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGa                                                                           >  2:954998/1‑100 (MQ=255)

                      aCTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGa                                                                           >  1:2283082/1‑100 (MQ=255)

                       cTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGAc                                                                          <  2:1870222/100‑1 (MQ=255)

                       cTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGAc                                                                          >  1:1356801/1‑100 (MQ=255)

                       cTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGAc                                                                          >  2:1196109/1‑100 (MQ=255)

                         tCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTa                                                                        >  2:2330406/1‑100 (MQ=255)

                          cacaTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAg                                                                       >  2:1366247/1‑100 (MQ=255)

                            caTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGcc                                                                     <  2:1460644/100‑1 (MQ=255)

                            caTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGcc                                                                     >  1:462410/1‑100 (MQ=255)

                               aGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGt                                                                  <  1:1972091/100‑1 (MQ=255)

                                 cTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGttt                                                                >  2:1425303/1‑100 (MQ=255)

                                 cTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGttt                                                                >  2:1819048/1‑100 (MQ=255)

                                  ttccgatctttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTg                                                               <  2:2731756/91‑1 (MQ=255)

                                  ttCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTg                                                               >  1:66951/1‑100 (MQ=255)

                                   tCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTgg                                                              >  1:3018376/1‑100 (MQ=255)

                                   tCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTgg                                                              <  1:2812564/100‑1 (MQ=255)

                                    cAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGt                                                             >  1:2935503/1‑100 (MQ=255)

                                        ttcttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAg                                                         <  1:278687/100‑1 (MQ=255)

                                         tcttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGt                                                        <  1:569878/100‑1 (MQ=255)

                                           ttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTgagatcggaa                                                      >  1:2731793/1‑91 (MQ=255)

                                             cccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTg                                                    >  2:1107693/1‑100 (MQ=255)

                                             cccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTg                                                    >  1:265921/1‑100 (MQ=255)

                                              ccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGc                                                   >  2:906241/1‑100 (MQ=255)

                                              ccACATAGTCTTCGTATCCTGCTGCCATTACAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGc                                                   <  1:1024553/100‑1 (MQ=255)

                                               cacaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCt                                                  >  1:1727969/1‑100 (MQ=255)

                                               cacaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCt                                                  >  2:2843816/1‑100 (MQ=255)

                                                acaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtg                                                 <  2:1715242/100‑1 (MQ=255)

                                                     gTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTgttgtt                                            <  1:1712895/100‑1 (MQ=255)

                                                     gTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTgttgtt                                            <  2:437220/100‑1 (MQ=255)

                                                     gTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTgttgtt                                            >  2:335533/1‑100 (MQ=255)

                                                     gTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTgttgtt                                            >  2:1224417/1‑100 (MQ=255)

                                                      tCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGttgttg                                           <  1:2036482/100‑1 (MQ=255)

                                                       cTTCGTATCCTGCTGCCATTGAAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGttgttg                                           >  2:1273686/1‑99 (MQ=255)

                                                          cGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTc                                       <  1:598363/100‑1 (MQ=255)

                                                          cGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTc                                       >  2:3053939/1‑100 (MQ=255)

                                                           gTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTcc                                      <  1:2945328/100‑1 (MQ=255)

                                                           gTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTcc                                      <  1:552316/100‑1 (MQ=255)

                                                           gTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGGATTGCTGTTGTTGATTcc                                      >  2:291306/1‑100 (MQ=255)

                                                             aTCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGa                                    <  2:1401655/100‑1 (MQ=255)

                                                               ccTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTccaa                                    >  2:2508146/1‑96 (MQ=255)

                                                               ccTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATa                                  >  2:1842850/1‑100 (MQ=255)

                                                                 tgctgcCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAAt                                <  1:1123555/100‑1 (MQ=255)

                                                                  gctgcCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATg                               <  1:113535/100‑1 (MQ=255)

                                                                  gctgcCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATg                               >  2:1296054/1‑100 (MQ=255)

                                                                  gctgcCATTGCAAAGGAGAAGACTATGTTGAAAT‑CCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGc                              <  2:1496790/100‑1 (MQ=255)

                                                                   cggcCATTGCAAAGGAGAAGACTATGTTGAAAT‑CCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCt                             <  1:393642/98‑1 (MQ=255)

                                                                    tgcCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCt                             <  1:585233/100‑1 (MQ=255)

                                                                    tgcCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCt                             <  1:1096508/100‑1 (MQ=255)

                                                                    tgcCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCt                             <  1:488330/100‑1 (MQ=255)

                                                                     gcCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTg                            >  2:2646840/1‑100 (MQ=255)

                                                                     gcCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTg                            >  2:1106608/1‑100 (MQ=255)

                                                                     gcCATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTg                            >  2:1212171/1‑100 (MQ=255)

                                                                      ccATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGa                           >  1:831583/1‑100 (MQ=255)

                                                                      ccATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGa                           >  2:455858/1‑100 (MQ=255)

                                                                      ccATTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGa                           >  2:1784067/1‑100 (MQ=255)

                                                                        aTTGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGAtt                         >  2:2287624/1‑100 (MQ=255)

                                                                         ttGCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTc                        <  1:409316/100‑1 (MQ=255)

                                                                           gCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCgg                      >  2:1539215/1‑100 (MQ=255)

                                                                           gCAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCgg                      >  2:1436532/1‑100 (MQ=255)

                                                                            cAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGc                     <  2:345958/100‑1 (MQ=255)

                                                                            cAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGc                     >  1:866063/1‑100 (MQ=255)

                                                                            cAAAGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGc                     <  1:2686136/100‑1 (MQ=255)

                                                                             aaaGGAGAAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCa                    <  2:254671/100‑1 (MQ=255)

                                                                                 gagaAGACTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTg                >  1:2024653/1‑100 (MQ=255)

                                                                                     agaCTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTg            >  1:2326707/1‑100 (MQ=255)

                                                                                     agaCTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTg            >  1:2934446/1‑100 (MQ=255)

                                                                                     agaCTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTg            >  1:2935585/1‑100 (MQ=255)

                                                                                      gaCTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGATTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGt           >  2:109319/1‑100 (MQ=255)

                                                                                       aCTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTc          <  2:2695023/100‑1 (MQ=255)

                                                                                       aCTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTc          <  1:3053982/100‑1 (MQ=255)

                                                                                        cTATGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCg         <  1:249678/100‑1 (MQ=255)

                                                                                          aTGTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGaa       >  1:1568295/1‑100 (MQ=255)

                                                                                            gTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACg     >  2:2988499/1‑100 (MQ=255)

                                                                                            gTTGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACg     <  1:1061059/100‑1 (MQ=255)

                                                                                            gTTGAAAG‑CCCGCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACg     >  2:1601854/1‑100 (MQ=255)

                                                                                             ttGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGt    >  2:2122414/1‑100 (MQ=255)

                                                                                              tGAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTg   >  1:1545448/1‑100 (MQ=255)

                                                                                               gAAATCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATACCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGc  <  2:617371/100‑1 (MQ=255)

GCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑CCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGC  >  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: `‑`