BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,183,028	T→C	V70A (GTC→GCC)	b3042 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,183,028	0	T	C	99.5%	601.6 / ‑5.0	188	V70A (GTC→GCC)	b3042	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (81/106); total (82/106)
Fisher's exact test for biased strand distribution p-value = 4.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.47e-01

Only 100 of 197 total aligned reads displayed.

ATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTGAGATCAAAAAGCAGCCA  >  NC_000913/3182938‑3183116

aTGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGccc                                                                                        <  1:2712354/93‑1 (MQ=255)

aTGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGccc                                                                                        <  2:1392922/93‑1 (MQ=255)

aTGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGccc                                                                                        <  1:382047/93‑1 (MQ=255)

 tgtgGAGAAAAAAATCCGCCAACCCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCt                                                                                       <  2:2665460/93‑1 (MQ=255)

 tgtgGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCt                                                                                       >  2:3269718/1‑93 (MQ=255)

  gtgGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTg                                                                                      <  2:4494023/93‑1 (MQ=255)

  gtgAAGAAAAAAATCCGCCAACCCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTg                                                                                      <  2:4436157/93‑1 (MQ=255)

   ttgagaAAAAACTCCGCCAACCCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGt                                                                                     <  1:1812613/91‑1 (MQ=255)

    ggAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGtt                                                                                    <  2:514440/93‑1 (MQ=255)

    ggAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGtt                                                                                    >  1:2297941/1‑93 (MQ=255)

     gagaAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTa                                                                                   <  1:561293/93‑1 (MQ=255)

     gagaAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTa                                                                                   <  1:2587876/93‑1 (MQ=255)

       gaAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACg                                                                                 >  1:3610768/1‑93 (MQ=255)

        aaaaaaaTCCGCCAAACCCTACAAGCGCAGCTGGCGCGCCTCGACCCTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACgg                                                                                >  2:276276/1‑93 (MQ=255)

         aaaaaaTCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGa                                                                               <  1:730902/93‑1 (MQ=255)

         aaaaaaTCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGa                                                                               <  2:2860280/93‑1 (MQ=255)

         aaaaaaCCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGa                                                                               <  2:2960413/93‑1 (MQ=255)

           aaaaTCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACt                                                                             <  2:3777145/93‑1 (MQ=255)

           aaaaTCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACt                                                                             <  2:3922755/93‑1 (MQ=255)

            aaaCCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTc                                                                            <  2:2044340/93‑1 (MQ=255)

             aaTCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCg                                                                           >  1:2225817/1‑93 (MQ=255)

             aaTCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCg                                                                           <  1:820407/93‑1 (MQ=255)

              aTCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAATCCCTGTTACGGACTCGt                                                                          >  2:4204828/1‑93 (MQ=255)

              aTCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGt                                                                          >  1:3560745/1‑93 (MQ=255)

                  cccAAACCCTACAGGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGaaa                                                                      <  1:3180718/92‑1 (MQ=255)

                    cAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGaaaaa                                                                    >  2:631328/1‑93 (MQ=255)

                      aaCCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACt                                                                  >  2:3684449/1‑93 (MQ=255)

                       aCCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTg                                                                 <  2:136718/93‑1 (MQ=255)

                        cccTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTgg                                                                <  1:1106455/93‑1 (MQ=255)

                        cccTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTgg                                                                <  2:1171826/93‑1 (MQ=255)

                        cccTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTgg                                                                >  2:3426827/1‑93 (MQ=255)

                        cccTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTgg                                                                >  1:1064249/1‑93 (MQ=255)

                         ccTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGgc                                                               <  2:1948404/93‑1 (MQ=255)

                         ccTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGgc                                                               <  1:2416048/93‑1 (MQ=255)

                         ccTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGgc                                                               <  1:2163821/93‑1 (MQ=255)

                         ccTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGgc                                                               >  2:4710284/1‑93 (MQ=255)

                            aCAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCgct                                                            <  1:2267216/93‑1 (MQ=255)

                             cAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCgctg                                                           <  2:3230128/93‑1 (MQ=255)

                              aaGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCACTGTTACGGACTCGTGAAAAACTCGCgctgc                                                          >  1:4480949/1‑93 (MQ=255)

                               aGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCgctggt                                                         >  2:947531/1‑91 (MQ=255)

                               aGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCgctgct                                                         >  2:508200/1‑93 (MQ=255)

                               aGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCgctgct                                                         >  1:4218820/1‑93 (MQ=255)

                               aGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCgctgct                                                         >  1:3622658/1‑93 (MQ=255)

                                gcgcAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCtt                                                        >  1:891022/1‑93 (MQ=255)

                                gcgcAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCtt                                                        <  1:2905184/93‑1 (MQ=255)

                                  gcAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGa                                                      >  2:1279743/1‑93 (MQ=255)

                                  gcAACTGCCGCGCCTCGACCTTGTAAGCCGTGAAGCGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCATGa                                                      >  1:2572239/1‑93 (MQ=255)

                                    aGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAAc                                                    >  1:469184/1‑93 (MQ=255)

                                    aGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAAc                                                    >  2:1541101/1‑93 (MQ=255)

                                    aGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAAc                                                    >  2:2479638/1‑93 (MQ=255)

                                    aGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAAAGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCCCTGCTTGAAc                                                    >  1:3140712/1‑93 (MQ=255)

                                     gCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACa                                                   >  2:1124056/1‑93 (MQ=255)

                                       tGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACaaa                                                 >  2:2711266/1‑93 (MQ=255)

                                       tGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACaaa                                                 <  2:1690340/93‑1 (MQ=255)

                                           gcgcCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAAt                                             >  2:4106667/1‑93 (MQ=255)

                                           gcgcCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAAt                                             <  2:1750524/93‑1 (MQ=255)

                                            cgcCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATc                                            <  2:2628593/93‑1 (MQ=255)

                                            cgcCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATc                                            <  1:3335300/93‑1 (MQ=255)

                                             gcCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAACAACTGGCGCTGCTTGAACAAAGAATCa                                           >  1:705567/1‑93 (MQ=255)

                                             gcCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCa                                           <  2:2547186/93‑1 (MQ=255)

                                              ccTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAg                                          >  1:508407/1‑93 (MQ=255)

                                              ccTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAg                                          >  2:1458772/1‑93 (MQ=255)

                                              ccTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAg                                          <  2:1238780/93‑1 (MQ=255)

                                                tCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCg                                        >  1:1899428/1‑93 (MQ=255)

                                                 cGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGa                                       <  2:4295949/93‑1 (MQ=255)

                                                  gACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGaa                                      >  2:3049453/1‑93 (MQ=255)

                                                   aCCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAAc                                     >  2:1722364/1‑93 (MQ=255)

                                                   aCCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAAc                                     >  1:1340834/1‑93 (MQ=255)

                                                    ccTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACt                                    <  1:2755658/93‑1 (MQ=255)

                                                     cTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACACAGAATCAGCGAACTg                                   >  2:1692316/1‑93 (MQ=255)

                                                     cTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTg                                   >  2:934329/1‑93 (MQ=255)

                                                     cTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTg                                   >  1:1198511/1‑93 (MQ=255)

                                                      ttGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTgg                                  >  1:678868/1‑93 (MQ=255)

                                                      ttGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTgg                                  >  2:3256990/1‑93 (MQ=255)

                                                         tAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGaaa                               <  2:2387858/93‑1 (MQ=255)

                                                          aaGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGaaaa                              <  2:363094/93‑1 (MQ=255)

                                                            ggcGTGAAGAGTTCGACGTCCAAACGCAACCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATc                            <  2:3569325/91‑1 (MQ=255)

                                                            gCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATc                            <  1:1398020/93‑1 (MQ=255)

                                                            accGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATc                            <  2:3621772/92‑1 (MQ=255)

                                                               gTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCa                         <  1:1958039/93‑1 (MQ=255)

                                                               gTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCa                         <  2:2450358/93‑1 (MQ=255)

                                                               gTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCa                         <  2:228551/93‑1 (MQ=255)

                                                                tGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAg                        >  1:1760884/1‑93 (MQ=255)

                                                                 gAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGt                       >  2:1859719/1‑93 (MQ=255)

                                                                   agagTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTAc                     <  1:248873/93‑1 (MQ=255)

                                                                    gagTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACt                    <  2:1442195/93‑1 (MQ=255)

                                                                     agTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTg                   <  1:4150196/93‑1 (MQ=255)

                                                                      gTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTga                  <  1:433906/93‑1 (MQ=255)

                                                                      gTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTga                  <  1:2700841/93‑1 (MQ=255)

                                                                      gTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTga                  >  1:3715525/1‑93 (MQ=255)

                                                                        tCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTgaga                <  1:1596071/93‑1 (MQ=255)

                                                                        tCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTgaga                <  1:1411293/93‑1 (MQ=255)

                                                                          gACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTGAGATc              >  1:4723849/1‑93 (MQ=255)

                                                                           aCGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTGAGATCa             >  2:1798723/1‑93 (MQ=255)

                                                                            cGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTGAGATCaa            <  1:3049260/93‑1 (MQ=255)

                                                                                cAAACGCAAGCGCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTGAGATCAAAAag        <  1:2526902/93‑1 (MQ=255)

                                                                                    cGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTGAGATCAAAAagcagc    <  2:695131/93‑1 (MQ=255)

                                                                                     gCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTGAGATCAAAAAGCAGcc   <  1:1825712/93‑1 (MQ=255)

                                                                                     gCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTGAGATCAAAAAGCAGcc   <  2:4039879/93‑1 (MQ=255)

                                                                                     gCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTGAGATCAAAAAGCAGcc   <  1:1142376/93‑1 (MQ=255)

ATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTGAGATCAAAAAGCAGCCA  >  NC_000913/3182938‑3183116

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