BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,484,261	T→C	L159P (CTG→CCG)	yjgP →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,484,261	0	T	C	99.4%	972.1 / ‑5.0	308	L159P (CTG→CCG)	yjgP	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (2/0); new base C (160/146); total (162/146)
Fisher's exact test for biased strand distribution p-value = 5.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 324 total aligned reads displayed.

ATGAAGTGTTAGCAGAAGCGAAAGCGAACCCTGGCATGGCGGCGCTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCTGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCAAAAGGTAATGCACGTCCTTCTGTGGTGGTG  >  NC_000913/4484171‑4484352

   aaGTGTTAGCAGAAGCGAAAGCGAACCCTGGCATGGCGGCGCTGGCGCAAGGGTAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCCGTTCa                                                                                        <  2:4368573/93‑1 (MQ=255)

    aGTGTTAGCAGAAGCGAAAGCGAACCCTGGCATGGCGGCGCTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCCGTTCAt                                                                                       >  2:7708045/1‑93 (MQ=255)

        ttAGCAGAAGCGAAAGCGAACCCTGGCATGGCGGCGCTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCCGTTCATCGaa                                                                                   >  2:415730/1‑93 (MQ=255)

          aGCAGAAGCGAAAGCGAACCCTGGCATGGCGGCGCTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCCGTTCATCGAAAg                                                                                 <  1:940753/93‑1 (MQ=255)

          aGCAGAAGCGAAAGCGAACCCTGGCATGGCGGCGCTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCCGGTCATCGAAAg                                                                                 >  1:2013118/1‑93 (MQ=255)

            cAGAAGCGAAAGCGAACCCTGGCATGGCGGCGCTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCCGTTCATCGAAAGCg                                                                               <  1:2736787/93‑1 (MQ=255)

             aGAAGCGAAAGCGAACCCTGGCATGGCGGCGCTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCCGTTCATCGCAAGCGt                                                                              >  2:2100676/1‑93 (MQ=255)

             aGAAGCGAAAGCGAACCCTGGCATGGCGGCGCTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCCGTTCATCGAAAGCGt                                                                              >  1:4537789/1‑93 (MQ=255)

              gAAGCGAAAGCGAACCCTGGCATGGCGGCGCTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCCGTTCATCGAAAGCGtt                                                                             >  2:1529179/1‑93 (MQ=255)

              gAAGCGAAAGCGAACCCTGGCATGGCGGCGCTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCCGTTCATCGAAAGCGtt                                                                             >  1:5085488/1‑93 (MQ=255)

              gAAGCGAAAGCGAACCCTGGCATGGCGGCGCTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCCGTTCATCGAAAGCGtt                                                                             >  2:3772065/1‑93 (MQ=255)

                agcgaaagcgaaCCCTGGCATGGCGGCGCTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCCGTTCATTGAAAGCGTTGa                                                                           >  2:7111904/1‑93 (MQ=255)

                agcgaaagcgaaCCCTGGCATGGCGGCGCTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCCGTTCATCGAAAGCGTTGa                                                                           <  2:6674558/93‑1 (MQ=255)

                agcgaaagcgaaCCCTGGCATGGCGGCGCTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCCGTTCATCGAAAGCGTTGa                                                                           >  2:755864/1‑93 (MQ=255)

                 gcgaaagcgaaCCCTGGCATGGCGGCGCTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCCGTTCATCGAAAGCGTTGAc                                                                          <  1:6516362/93‑1 (MQ=255)

                   gaaagcgaaCCCTGGCATGGCGGCGCTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCCGTTCATCGAAAGCGTTGACgg                                                                        <  1:3580737/93‑1 (MQ=255)

                    aaagcgaaCCCTGGCATGGCGGCGCTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCCGTTCATCGAAAGCGTTGACGGc                                                                       >  2:6347636/1‑93 (MQ=255)

                     aagcgaaCCCTGGCATGGCGGCGCTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCCGTTCATCGAAAGCGTTGACGGCa                                                                      >  2:6856385/1‑93 (MQ=255)

                     aagcgaaCCCTGGCATGGCGGCGCTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCCGTTCATCGAAAGCGTTGACGGCa                                                                      >  2:3405145/1‑93 (MQ=255)

                     aagcgaaCCCTGGCATGGCGGCGCTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCCGTTCATCGAAAGCGTTGACGGCa                                                                      >  1:7041603/1‑93 (MQ=255)

                       gcgaaCCCTGGCATGGCGGCGCTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCCGTTCATCGAAAGCGTTGACGGCAGc                                                                    <  1:1819348/93‑1 (MQ=255)

                       gcgaaCCCTGGCATGGCGGCGCTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCCGTTCATCGAAAGCGTTGACGGCAGc                                                                    >  2:3235999/1‑93 (MQ=255)

                       gcgaaCCCTGGCATGGCGGCGCTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCCGTGACGTTCATCGAAAGCGTTGACGGAAGc                                                                    >  1:48981/1‑93 (MQ=255)

                          aaCCCTGGCATGGCGGCGCTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCCGTTCATCGAAAGCGTTGACGGCAGCGAt                                                                 >  1:5838313/1‑93 (MQ=255)

                           aCCCTGGCATGGCGGCGCTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCCGTTCATCGAAAGCGTTGACGGCAGCGAtt                                                                >  1:4411920/1‑93 (MQ=255)

                           aCCCTGGCATGGCGGCGCTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCCGTTCATCGAAAGCGTTGACGGCAGCGAtt                                                                >  1:7373721/1‑93 (MQ=255)

                           aCCCTGGCATGGCGGCGCTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCCGTTCATCGAAAGCGTTGACGGCAGCGAtt                                                                <  2:2239983/93‑1 (MQ=255)

                             ccTGGCATGGCGGCGCTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCCGTTCATCGAAAGCGTTGACGGCAGCGATTTc                                                              <  2:1948067/93‑1 (MQ=255)

                              cTGGCATGGCGGCGCTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCCGTTCATCGAAAGCGTTGACGGCAGCGATTTCa                                                             >  2:5208560/1‑93 (MQ=255)

                              cTGGCATGGCGGCGCTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCCGTTCATCGAAAGCGTTGACGGCAGCGATTTCa                                                             >  1:7930465/1‑93 (MQ=255)

                               tGGCATGGCGGCGCTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCCGTTCATCGAAAGCGTTGACGGCAGCGATTTCaa                                                            >  2:3634718/1‑93 (MQ=255)

                                ggCATGGCGGCGCTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCCGTTCATCGAAAGCGTTGACGGCAGCGATTTCaaa                                                           <  1:3266845/93‑1 (MQ=255)

                                 gCATGGCGGCGCTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTTGGTGCCGTTCATCGAAATCGTTGACGGCAGCGATTTCAAAg                                                          <  2:6955323/93‑1 (MQ=255)

                                 gCATGGCGGCGCTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCCGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAg                                                          <  1:1559927/93‑1 (MQ=255)

                                 gCATGGCGGCGCTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCCGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAg                                                          <  1:6779779/93‑1 (MQ=255)

                                  cATGGCGGCGCTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCCGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGa                                                         <  2:533436/93‑1 (MQ=255)

                                  cATGGCGGCGCTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCCGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGa                                                         >  1:639986/1‑93 (MQ=255)

                                    tGGCGGCGCTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCCGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGAtg                                                       <  2:2651153/93‑1 (MQ=255)

                                    tGGCGGCGCTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCCGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGAtg                                                       <  2:6143872/93‑1 (MQ=255)

                                    tGGCGGCGCTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCCGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGAtg                                                       <  2:3005656/93‑1 (MQ=255)

                                    tGGCGGCGCTGGCGCAAGGGCAATTACAGCAAGCGACTAATGGCAGCTCGGTGAAGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGAtg                                                       <  1:4852948/93‑1 (MQ=255)

                                     ggcggcGCTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCCGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGAtgt                                                      <  1:7229978/93‑1 (MQ=255)

                                      gcggcgCTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCCGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGAtgtg                                                     >  1:5086057/1‑93 (MQ=255)

                                       gggcgCTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCCGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATgtgt                                                    <  2:1488612/92‑1 (MQ=255)

                                        ggcgCTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCCGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGtt                                                   >  1:6719717/1‑93 (MQ=255)

                                        ggcgCTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCCGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGtt                                                   >  2:6422031/1‑93 (MQ=255)

                                         gcgcTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCCGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTc                                                  <  1:2494127/93‑1 (MQ=255)

                                         gcgcTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCCGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTc                                                  >  2:1225642/1‑93 (MQ=255)

                                           gcTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCCGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCt                                                >  1:5659804/1‑93 (MQ=255)

                                            cTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCCGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTc                                               <  2:2958002/93‑1 (MQ=255)

                                              ggCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCCGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTcgc                                             >  2:4823921/1‑93 (MQ=255)

                                                   aaGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCCGTTCATCGCAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCaaa                                        >  2:2386131/1‑93 (MQ=255)

                                                   aaGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCCGTTCATCGAAAGCGTTGACGGCAGCGCTTTCAAAGATGTGTTCCTCGCGCaaa                                        >  1:7600863/1‑93 (MQ=255)

                                                    aGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCCGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAAt                                       <  1:801285/93‑1 (MQ=255)

                                                    aGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCCGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAAt                                       >  2:7854364/1‑93 (MQ=255)

                                                     gggCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCCGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAAtt                                      <  1:5774410/93‑1 (MQ=255)

                                                     gggCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCCGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAAtt                                      <  2:3050598/93‑1 (MQ=255)

                                                     gggCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCCGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAAtt                                      >  1:5988517/1‑93 (MQ=255)

                                                     gggCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCCGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAAtt                                      <  1:5548224/93‑1 (MQ=255)

                                                     gggCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCCGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAAtt                                      <  2:5918202/93‑1 (MQ=255)

                                                     gggCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCCGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAAtt                                      <  1:1062956/93‑1 (MQ=255)

                                                      ggCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCCGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTc                                     >  2:1021170/1‑93 (MQ=255)

                                                             ccAGCAAGCGACTAATGGCAGCTCGGTGCCGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCaaa                              >  2:3991329/1‑93 (MQ=255)

                                                             ccAGCAAGCGACTAATGGCAGCTCGGTGCCGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCaaa                              >  1:1209820/1‑93 (MQ=255)

                                                             ccAGCAAGCGACTAATGGCAGCTCGGTGCCGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCaaa                              <  1:7727092/93‑1 (MQ=255)

                                                             ccAGCAAGCGACTAATGGCAGCTCGGTGCCGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCaaa                              >  1:7397813/1‑93 (MQ=255)

                                                              cAGCAAGCGACTAATGGCAGCTCGGTGCCGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCaaaa                             <  1:4097916/93‑1 (MQ=255)

                                                               aGCAAGCGACTAATGGCAGCTCGGTGCCGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCAAAAg                            <  1:6652184/93‑1 (MQ=255)

                                                               aGCAAGCGACTAATGGCAGCTCGGTGCCGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCAAAAg                            <  1:5241465/93‑1 (MQ=255)

                                                                gCAAGCGACTAATGGCAGCTCGGTGCCGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCAAAAgg                           <  2:5042237/93‑1 (MQ=255)

                                                                 cAAGCGACTAATGGCAGCTCGGTGCCGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCAAAAGGt                          >  1:34929/1‑93 (MQ=255)

                                                                 cAAGCGACTAATGGCAGCTCGGTGCCGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCAAAAGGt                          <  1:2268746/93‑1 (MQ=255)

                                                                  aaGCGACTAATGGCAGCTCGGTGCCGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCAAAAGGTa                         <  2:7317417/93‑1 (MQ=255)

                                                                   aGCGACTAATGGCAGCTCGTTGCCGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCAAAAGGTaa                        <  1:6578411/93‑1 (MQ=255)

                                                                   aGCGACTAATGGCAGCTCGGTGCCGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCAAAAGGTaa                        <  1:4719251/93‑1 (MQ=255)

                                                                     cGACTAATGGCAGCTCGGTGCCGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCAAAAGGTAATg                      <  1:1352166/93‑1 (MQ=255)

                                                                      gACTAATGGCAGCTCGGTGCCGTTCATCGAAAGCGTTGACGGCAGCGCGTTCAAAGATGTGGTCCTCGCGCAAATTCGACCAAAAGGGAATGc                     >  2:4576551/1‑93 (MQ=255)

                                                                      gACTAATGGCAGCTCGGTGCCGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCAAAAGGTAATGc                     >  2:7895935/1‑93 (MQ=255)

                                                                      gACTAATGGCAGCTCGGTGCCGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCAAAAGGTAATGc                     <  2:5735438/93‑1 (MQ=255)

                                                                      gACTAATGGCAGCTCGGTGCCGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCAAAAGGTAATGc                     <  2:6560919/93‑1 (MQ=255)

                                                                        cTAATGGCAGCTCGGTGCCGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCAAAAGGTAATGCAc                   <  1:40146/93‑1 (MQ=255)

                                                                        cTAATGGCAGCTCGGTGCCGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCAAAAGGTAATGCAc                   >  2:1682790/1‑93 (MQ=255)

                                                                         tAATGGCAGCTCGGTGCTGCTCATCGAAAGGGTTGACGGCAGCGATTTCAAGGATGTGTTCCGCGCGCACATCGCTCCAAAGGGAAATGcgct                  >  1:2691108/1‑90 (MQ=255)

                                                                         tAATGGCAGCTAGGTGCCGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCCAATTCGACCAAAAGGTAATGCACg                  >  2:839765/1‑93 (MQ=255)

                                                                           aTGGCAGCTCGGTGCCGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCAAAAGGTAATGCACGTc                <  1:2532698/93‑1 (MQ=255)

                                                                           aTGGCAGCTCGGTGCCGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCAAAAGGTAATGCACGTc                >  2:5404638/1‑93 (MQ=255)

                                                                           aTGGCAGCTCGGTGCCGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCAAAAGGTAATGCACGTc                <  2:3872461/93‑1 (MQ=255)

                                                                            tGGCAGCTCGGTGCCGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCAAAAGGTAATGCACGTcc               >  1:3367798/1‑93 (MQ=255)

                                                                             ggCAGCTCGGTGCCGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCAAAAGGTAATGCACGTCCt              <  1:78346/93‑1 (MQ=255)

                                                                               cAGCTCGGTGCCGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCAAAAGGTAATGCACGTCCTTc            <  1:292126/93‑1 (MQ=255)

                                                                                 gCTCGGTGCCGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCAAAAGGTAATGCACGTCCTTCtg          >  1:2364994/1‑93 (MQ=255)

                                                                                  cTCGGTGCCGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCAAAAGGTAATGCACGTCCTTCtgt         <  2:4187809/93‑1 (MQ=255)

                                                                                   tCGGTGCCGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGCCCAAACGGTAATGCACGTCCTTCtgtg        >  1:2126346/1‑93 (MQ=255)

                                                                                   tCGGTGCCGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCAAAAGGTAATGCACGTCCTTCtgtg        >  1:355318/1‑93 (MQ=255)

                                                                                    cGGTGCCGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCAAAAGGTAATGCACGTCCTTCTgtgg       >  2:4755189/1‑93 (MQ=255)

                                                                                     ggTGCCGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCAAAAGGTAATGCACGTCCTTCTgtggt      >  2:6735678/1‑93 (MQ=255)

                                                                                      gTGCCGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCAAAAGGTAATGCACGTCCTTCTgtggtg     >  1:5777353/1‑93 (MQ=255)

                                                                                      gTGCCGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCAAAAGGTAATGCACGTCCTTCTgtggtg     >  1:7479010/1‑93 (MQ=255)

                                                                                       tGCCGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCAAAAGGTAATGCACGTCCTTCTgtggtgg    >  1:3150962/1‑93 (MQ=255)

                                                                                         cTGTTCATCGAAAGCGTGGTCGGCGGCGAATCCAAAGATGTGTTCCGGGCGCAAATTCGACCAAAATGTAAGGCGCGTCGTTCTgtggtggtg  <  2:5017812/93‑1 (MQ=255)

ATGAAGTGTTAGCAGAAGCGAAAGCGAACCCTGGCATGGCGGCGCTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCTGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGATGTGTTCCTCGCGCAAATTCGACCAAAAGGTAATGCACGTCCTTCTGTGGTGGTG  >  NC_000913/4484171‑4484352

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