BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,476,619	A→G	V91A (GTG→GCG)	b4256 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,476,619	0	A	G	99.6%	784.2 / ‑5.2	245	V91A (GTG→GCG)	b4256	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (113/130); total (113/132)
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.37e-01

Only 100 of 252 total aligned reads displayed.

ATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAA  >  NC_000913/4476529‑4476707

 tCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCtt                                                                                       >  1:3223805/1‑93 (MQ=255)

  ccGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTc                                                                                      <  2:5822467/93‑1 (MQ=255)

  ccGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTc                                                                                      >  2:6228779/1‑93 (MQ=255)

   cGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTGAGGTTCGCGCAACATCCGCGCTTca                                                                                     >  1:163246/1‑93 (MQ=255)

   cGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTca                                                                                     >  2:4265018/1‑93 (MQ=255)

    ggCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTcac                                                                                    >  2:787220/1‑93 (MQ=255)

    ggCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTcac                                                                                    >  1:5148085/1‑93 (MQ=255)

      cATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACaa                                                                                  >  1:4159038/1‑93 (MQ=255)

        tAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATc                                                                                <  2:1472144/93‑1 (MQ=255)

        tAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATc                                                                                >  2:3916037/1‑93 (MQ=255)

        tAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATc                                                                                >  1:5923815/1‑93 (MQ=255)

        tAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATc                                                                                <  2:3475191/93‑1 (MQ=255)

        tAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATc                                                                                <  2:2405711/93‑1 (MQ=255)

        tAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATc                                                                                >  2:4934864/1‑93 (MQ=255)

        tAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATc                                                                                >  2:2610173/1‑93 (MQ=255)

          aaCAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCgc                                                                              >  2:4395753/1‑92 (MQ=255)

          aaCAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGa                                                                              >  1:2724639/1‑93 (MQ=255)

            cAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACa                                                                            <  2:3693730/93‑1 (MQ=255)

            cAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACa                                                                            >  2:285782/1‑93 (MQ=255)

             aaaaCGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACAt                                                                           >  1:3052519/1‑93 (MQ=255)

             aaaaCGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACAt                                                                           <  2:5759922/93‑1 (MQ=255)

             aaaaCGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACAt                                                                           >  1:2466010/1‑93 (MQ=255)

              aaaCGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATg                                                                          <  1:5756691/93‑1 (MQ=255)

              aaaCGTACGTGGTTAATGTGCGAAGGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATg                                                                          <  1:4424426/93‑1 (MQ=255)

                acgtacgtTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGc                                                                        <  1:3124961/93‑1 (MQ=255)

                acgtacgtTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGc                                                                        <  2:3182307/93‑1 (MQ=255)

                 cgtacgtTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGcc                                                                       >  1:1015631/1‑93 (MQ=255)

                 cgtacgtTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGcc                                                                       <  2:2421127/93‑1 (MQ=255)

                    acgtTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGt                                                                    <  2:5505967/93‑1 (MQ=255)

                     cgtTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTg                                                                   <  1:2665215/93‑1 (MQ=255)

                     cgtTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTg                                                                   <  2:229428/93‑1 (MQ=255)

                      gttgttAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGc                                                                  >  1:516703/1‑93 (MQ=255)

                       ttgttAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCa                                                                 <  1:5267552/93‑1 (MQ=255)

                            aaTGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCaa                                                            <  1:4805234/93‑1 (MQ=255)

                            aaTGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGACCGTGCAGCCaa                                                            >  2:5250908/1‑93 (MQ=255)

                             aTGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAg                                                           <  1:1958089/93‑1 (MQ=255)

                             aTGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAg                                                           <  2:3950802/93‑1 (MQ=255)

                             aTGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAg                                                           >  2:296870/1‑93 (MQ=255)

                              tGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCCAgc                                                          >  1:5459690/1‑93 (MQ=255)

                              tGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAgc                                                          >  1:809252/1‑93 (MQ=255)

                               gTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAgcg                                                         <  1:5410793/93‑1 (MQ=255)

                                   gAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATa                                                     <  1:4296257/93‑1 (MQ=255)

                                   gAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATa                                                     <  2:2100455/93‑1 (MQ=255)

                                   gAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATa                                                     >  1:6261944/1‑93 (MQ=255)

                                    aaTGCCCGCGAGTTGAACCCGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAg                                                    <  1:1526087/93‑1 (MQ=255)

                                     aTGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGt                                                   <  1:3606552/93‑1 (MQ=255)

                                     aTGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGt                                                   <  2:1050160/93‑1 (MQ=255)

                                     aTGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGt                                                   <  2:4648037/93‑1 (MQ=255)

                                      tGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTc                                                  >  1:86900/1‑93 (MQ=255)

                                        cccGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGa                                                <  2:6274184/93‑1 (MQ=255)

                                         ccGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGat                                               >  1:1910413/1‑93 (MQ=255)

                                         ccGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGat                                               <  1:2331804/93‑1 (MQ=255)

                                          cgcgAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGata                                              <  1:5514078/93‑1 (MQ=255)

                                           gcgAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGatat                                             <  1:6050675/93‑1 (MQ=255)

                                           gcgAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGatat                                             <  2:5888400/93‑1 (MQ=255)

                                           gcgAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGatat                                             <  2:6141454/93‑1 (MQ=255)

                                             gAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGc                                           >  2:3771394/1‑93 (MQ=255)

                                                 tGACCACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCaa                                       <  1:4312414/93‑1 (MQ=255)

                                                 tGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCaa                                       >  1:2421118/1‑93 (MQ=255)

                                                 tGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCaa                                       <  1:4192553/93‑1 (MQ=255)

                                                   aaCACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAg                                     >  1:2222553/1‑93 (MQ=255)

                                                    acacGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGc                                    >  1:6114172/1‑93 (MQ=255)

                                                    acacGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGc                                    <  2:5451510/93‑1 (MQ=255)

                                                    acacGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGc                                    <  1:2786190/93‑1 (MQ=255)

                                                       cgcgGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCa                                 <  2:3852382/93‑1 (MQ=255)

                                                        gcgGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCaa                                >  2:60682/1‑93 (MQ=255)

                                                        gcgGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCaa                                >  2:4813156/1‑93 (MQ=255)

                                                         cgGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCaaa                               <  1:690084/93‑1 (MQ=255)

                                                           ggCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATg                             >  1:1402677/1‑93 (MQ=255)

                                                            gCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGc                            <  2:1803856/93‑1 (MQ=255)

                                                            gCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCAATATGCTCAAAGCCCAAATGc                            <  1:564620/93‑1 (MQ=255)

                                                             cATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCt                           <  1:5988897/93‑1 (MQ=255)

                                                               tttttttAGAGTTCGCGCAACATCCGCTCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCa                         <  2:2735032/93‑1 (MQ=255)

                                                               tttttttAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCa                         <  2:454722/93‑1 (MQ=255)

                                                               tttttttAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCa                         <  2:4208209/93‑1 (MQ=255)

                                                                ttttttAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCAt                        >  1:3904895/1‑93 (MQ=255)

                                                                   tttAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATaaa                     <  2:1412423/93‑1 (MQ=255)

                                                                       gagTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAgcgc                 <  1:1181545/93‑1 (MQ=255)

                                                                        agTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCa                <  2:4429786/93‑1 (MQ=255)

                                                                         gTTCGCGCACCATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCaa               <  1:739620/93‑1 (MQ=255)

                                                                         gTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCaa               <  1:338815/93‑1 (MQ=255)

                                                                         gTTCGCGCAACATCAGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCaa               <  1:2959906/93‑1 (MQ=255)

                                                                          ttCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAAt              <  1:5160334/93‑1 (MQ=255)

                                                                          ttCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAAt              >  1:4645901/1‑93 (MQ=255)

                                                                          ttCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAAt              <  2:2059867/93‑1 (MQ=255)

                                                                           tCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATg             <  2:235780/93‑1 (MQ=255)

                                                                           tCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATg             >  1:4001413/1‑93 (MQ=255)

                                                                            cgcgCAACAGCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATgg            <  1:1288449/93‑1 (MQ=255)

                                                                             gcgcAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGc           <  2:1177796/93‑1 (MQ=255)

                                                                             gcgcAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGc           <  2:5640424/93‑1 (MQ=255)

                                                                               gcAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCtt         >  2:3085619/1‑93 (MQ=255)

                                                                                cAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTc        >  1:1663867/1‑93 (MQ=255)

                                                                                 aaCATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTcc       <  2:3558939/93‑1 (MQ=255)

                                                                                  aCATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCt      <  1:2926159/93‑1 (MQ=255)

                                                                                  aCATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCt      >  1:1101637/1‑93 (MQ=255)

                                                                                   cATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCtt     >  1:4505757/1‑93 (MQ=255)

                                                                                   cATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCtt     >  2:1826220/1‑93 (MQ=255)

                                                                                   aaTCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCtt     <  1:4566788/92‑1 (MQ=255)

                                                                                    aTCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCttt    >  2:5038578/1‑93 (MQ=255)

                                                                                      ccGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCACAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTaa  >  1:2597147/1‑93 (MQ=255)

ATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAA  >  NC_000913/4476529‑4476707

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