BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,736,520	A→G	Q218R (CAG→CGG)	purR →	PurR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,736,520	0	A	G	99.1%	372.7 / ‑4.8	117	Q218R (CAG→CGG)	purR	PurR
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (66/49); total (68/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-01

Only 100 of 129 total aligned reads displayed.

TGGAACGTAACACCGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCAGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCG  >  NC_000913/1736430‑1736609

tGGAACGTAACACCGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCggg                                                                                         >  1:2839123/1‑93 (MQ=255)

tGGAACGTAACACCGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCggg                                                                                         <  2:1653544/93‑1 (MQ=255)

 ggAACGTAACACCGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCgggg                                                                                        >  2:2715776/1‑93 (MQ=255)

 ggAACGTAACACCGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCgggg                                                                                        >  1:2902623/1‑93 (MQ=255)

  gAACGTAACACCGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCggggg                                                                                       >  1:2282051/1‑92 (MQ=255)

   aaCGTAACACCGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTg                                                                                      >  2:533741/1‑93 (MQ=255)

   aaCGTAACACCGGCGCAGGCCGCCTTGCCGGTCTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTg                                                                                      >  2:4484384/1‑93 (MQ=255)

   aaCATAACACCGGCGAAGGCCGCATTGCCGGTATTATAAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTg                                                                                      <  2:4339689/93‑1 (MQ=255)

    aCGTAACACCGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGa                                                                                     >  2:2242875/1‑93 (MQ=255)

    aCGTAACACCGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGa                                                                                     >  1:750549/1‑93 (MQ=255)

     cGTAACACCGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGAc                                                                                    <  1:4253011/93‑1 (MQ=255)

      gTAACACCGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACt                                                                                   >  1:1111630/1‑93 (MQ=255)

           acCGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGaa                                                                              >  2:3560746/1‑93 (MQ=255)

                 gcAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGaa                                                                        >  1:3456548/1‑93 (MQ=255)

                 gcAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGaa                                                                        >  2:3851357/1‑93 (MQ=255)

                 gcAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGGGACTTTGAACCTGaa                                                                        >  2:3618040/1‑93 (MQ=255)

                  cAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAAt                                                                       <  2:229374/93‑1 (MQ=255)

                   aGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATc                                                                      <  1:3915979/93‑1 (MQ=255)

                    ggCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATcc                                                                     >  2:1406825/1‑93 (MQ=255)

                    ggCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATcc                                                                     >  2:1919565/1‑93 (MQ=255)

                    ggCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATcc                                                                     >  1:3200889/1‑93 (MQ=255)

                     gccgccTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCg                                                                    <  2:701092/93‑1 (MQ=255)

                        gccTTGGCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGtt                                                                 <  1:2957355/93‑1 (MQ=255)

                        gccTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCTGGGTGACTTTGAACCTGAATCCGGtt                                                                 >  1:864064/1‑93 (MQ=255)

                        gccTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGtt                                                                 >  1:2316439/1‑93 (MQ=255)

                          cTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTAt                                                               >  2:1148053/1‑93 (MQ=255)

                          cTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTAt                                                               <  1:931325/93‑1 (MQ=255)

                           ttGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATc                                                              <  2:1035228/93‑1 (MQ=255)

                           ttGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATc                                                              <  2:2545457/93‑1 (MQ=255)

                           ttGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATc                                                              >  1:168733/1‑93 (MQ=255)

                           ttGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATc                                                              >  2:4547847/1‑93 (MQ=255)

                           ttGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATc                                                              >  2:1953102/1‑93 (MQ=255)

                           ttGCCGGTTTTATGAAAGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTAtg                                                              >  2:1981698/1‑92 (MQ=255)

                            tGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGGGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATcg                                                             >  1:1421229/1‑93 (MQ=255)

                             gCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATcgc                                                            <  1:4421515/93‑1 (MQ=255)

                               cGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCgcgc                                                          >  1:2472170/1‑93 (MQ=255)

                               cGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCgcgc                                                          >  2:3734389/1‑93 (MQ=255)

                                ggTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGcc                                                         <  1:3733492/93‑1 (MQ=255)

                                ggTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGcc                                                         >  2:1118540/1‑93 (MQ=255)

                                 gTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCa                                                        >  2:4293740/1‑93 (MQ=255)

                                 gTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCa                                                        <  2:868511/93‑1 (MQ=255)

                                  ttttATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCAt                                                       <  2:1501062/93‑1 (MQ=255)

                                    ttattaAGGCGATGGAAGAAGCGATGAGCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATgc                                                     <  2:1528776/93‑1 (MQ=255)

                                     tATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATgca                                                    <  1:4062253/93‑1 (MQ=255)

                                      aTGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATgcag                                                   >  2:323898/1‑93 (MQ=255)

                                      aTGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATgcag                                                   >  1:594746/1‑93 (MQ=255)

                                       tGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATgcagc                                                  <  2:774570/93‑1 (MQ=255)

                                       tGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATgcagc                                                  >  1:1831501/1‑93 (MQ=255)

                                        gAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATgcagca                                                 >  1:1386109/1‑93 (MQ=255)

                                        gAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATgcagca                                                 >  1:1007774/1‑93 (MQ=255)

                                         aaGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCaa                                                >  2:1561493/1‑93 (MQ=255)

                                          aGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCaaa                                               <  1:1084122/93‑1 (MQ=255)

                                           ggCGATGGAAGATGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTAAAGCGCCATGCAGCaaaa                                              <  2:4212873/93‑2 (MQ=255)

                                           ggCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAAt                                              >  2:964484/1‑93 (MQ=255)

                                           ggCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAAt                                              >  1:2731727/1‑93 (MQ=255)

                                            gCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATc                                             >  2:3369764/1‑93 (MQ=255)

                                             cGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATcc                                            <  1:599436/93‑1 (MQ=255)

                                              gATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCt                                           >  2:3840931/1‑93 (MQ=255)

                                              gATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCt                                           <  1:4063560/93‑1 (MQ=255)

                                              gATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCt                                           <  1:4063559/93‑1 (MQ=255)

                                                tGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGt                                         <  2:4019780/93‑1 (MQ=255)

                                                 ggAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTc                                        <  2:1491752/93‑1 (MQ=255)

                                                 ggAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTc                                        >  1:2914410/1‑93 (MQ=255)

                                                  gaagaaGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGACTCCGGTTATCGCGCCATGCAGCAAATCCTGTCg                                       >  1:228221/1‑93 (MQ=255)

                                                  gaagaaGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCg                                       <  1:1112922/93‑1 (MQ=255)

                                                  gaagaaGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCg                                       >  2:2201397/1‑93 (MQ=255)

                                                    agaagCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCa                                     <  2:3264021/93‑1 (MQ=255)

                                                    agaagCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCa                                     <  1:1539198/93‑1 (MQ=255)

                                                    agaagCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCa                                     >  1:3084610/1‑93 (MQ=255)

                                                     gaagCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAg                                    >  2:1202204/1‑93 (MQ=255)

                                                         cGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGc                                <  1:369478/93‑1 (MQ=255)

                                                          gatgatCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCa                               <  1:179646/93‑1 (MQ=255)

                                                           atgatCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCAt                              >  1:2262199/1‑93 (MQ=255)

                                                             gatCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCg                            >  2:118966/1‑93 (MQ=255)

                                                              atCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGc                           >  2:4489416/1‑93 (MQ=255)

                                                              atCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGc                           >  1:4423024/1‑93 (MQ=255)

                                                               tCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCAATGCAGCAAATCCTGTCGCACCCGCATCGcc                          >  1:3759834/1‑93 (MQ=255)

                                                                   ggggCCGGAAAGCTGGATTGTGCAGGGTGACTTTGAACCTGAATCCGTTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTAc                      <  2:2729567/90‑1 (MQ=255)

                                                                    gTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACt                     <  1:3302657/93‑1 (MQ=255)

                                                                     tGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTg                    >  1:1684753/1‑93 (MQ=255)

                                                                      gCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGc                   <  2:1603110/93‑1 (MQ=255)

                                                                       ccGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGca                  >  2:2995152/1‑92 (MQ=255)

                                                                         ggAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGt                >  2:3507875/1‑93 (MQ=255)

                                                                         ggAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGt                >  1:3811971/1‑93 (MQ=255)

                                                                              gCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGtcttct           >  1:1434855/1‑93 (MQ=255)

                                                                              gCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGtcttct           >  1:3698031/1‑93 (MQ=255)

                                                                               cTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCtg          <  1:2617646/93‑1 (MQ=255)

                                                                                 ggATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCtgtg        >  1:2490749/1‑93 (MQ=255)

                                                                                 ggATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCtgtg        >  1:4487202/1‑93 (MQ=255)

                                                                                 ggATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCtgtg        <  2:2722016/93‑1 (MQ=255)

                                                                                 ggATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCtgtg        >  1:916418/1‑93 (MQ=255)

                                                                                  gATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTgtgg       <  1:994946/93‑1 (MQ=255)

                                                                                   aTTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTgtggt      <  1:61203/93‑1 (MQ=255)

                                                                                    ttGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTgtggtg     >  2:671777/1‑93 (MQ=255)

                                                                                    ttGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTgtggtg     <  1:78720/93‑1 (MQ=255)

                                                                                    ttGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTgtggtg     <  2:889842/93‑1 (MQ=255)

                                                                                     tgtgCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGtggtgg    >  1:1700379/1‑93 (MQ=255)

                                                                                      gtgCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGc   <  1:2084938/93‑1 (MQ=255)

                                                                                      gtgCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGc   >  1:3066675/1‑93 (MQ=255)

                                                                                       tgCAGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCg  >  1:3258577/1‑93 (MQ=255)

TGGAACGTAACACCGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCAGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCG  >  NC_000913/1736430‑1736609

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