BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,312,334	A→G	L116P (CTG→CCG)	phnP ←	PhnP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,312,334	0	A	G	99.6%	808.3 / ‑5.0	236	L116P (CTG→CCG)	phnP	PhnP
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (126/109); total (126/110)
Fisher's exact test for biased strand distribution p-value = 4.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.17e-01

Only 100 of 245 total aligned reads displayed.

AGTTTTGAGTGGTTGAGCGGCAGGGGCGTGACCTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCAGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCGGGATCGGAT‑CGCCAACGCCCCAGCGCAGCGGA  >  NC_000913/4312244‑4312423

aGTTTTGAGTGGTTGAGCGGCAGGGGCGTGACCTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCggg                                                                                          <  1:6328211/93‑1 (MQ=255)

   tttGTGTGGTTGAGCGACAGGGACGTGACCTGTAACCCCTAAAAATCAAACACCACAAACGTTTCCACCGTGTGGCTGAAATCAAGCGGGccc                                                                                       <  2:2028578/93‑1 (MQ=255)

    ttGAGTGGTTGAGCGGCAGGGGCGTGACCTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCg                                                                                      <  2:4205491/93‑1 (MQ=255)

      gAGTGGTTGAGCGGCAGGGGCGTGACCTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGa                                                                                    <  1:1529970/93‑1 (MQ=255)

         tGGTTGAGCGGCAGGGGCGTGACCTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGt                                                                                 <  1:4190722/93‑1 (MQ=255)

           gTTGAGCGGCAGGGGCGTGACCTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGttt                                                                               <  1:5083234/93‑1 (MQ=255)

            ttGAGCGGCAGGGGCGTGACCTGTAACCCCTGCAAATCAAACACCACAAACGGATCCACCGTGTGGCTGAAATCCAGAGGGTCCGGCAGTTTa                                                                              >  2:4670482/1‑93 (MQ=255)

             tGAGCGGCAGGGGCGTGACCTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTaa                                                                             <  2:2719941/93‑1 (MQ=255)

               aGCGGCAGGGGCGTGACCTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAAc                                                                           >  1:5938310/1‑93 (MQ=255)

               aGCGGCAGGGGCGTGACCTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAAc                                                                           <  1:2045594/93‑1 (MQ=255)

                gCGGCAGGGGCGTGACCTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACa                                                                          <  2:2790692/93‑1 (MQ=255)

                gCGGCAGGGGCGTGACCTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACa                                                                          <  2:3872628/93‑1 (MQ=255)

                   gCAGGGGCGTGACCTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGAt                                                                       >  2:5611258/1‑93 (MQ=255)

                   gCAGGGGCGTGACCTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCAGGCCCGGATGTTTAAACAGAt                                                                       <  2:23270/93‑1 (MQ=255)

                    gaGGGGCGTGACCTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGAtc                                                                      <  1:2579806/92‑1 (MQ=255)

                    cAGGGGCGTGACCTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGAtc                                                                      <  1:4093511/93‑1 (MQ=255)

                    cAGGGGCGTGACCTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGAtc                                                                      <  1:5395974/93‑1 (MQ=255)

                     aGGGGCGTGACCTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGAtcg                                                                     >  2:2175048/1‑93 (MQ=255)

                      ggggCGTGACCTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGAtcgt                                                                    >  2:4929157/1‑93 (MQ=255)

                      ggggCGTGACCTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGAtcgt                                                                    >  2:5716602/1‑93 (MQ=255)

                      ggggCGTGACCTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGAtcgt                                                                    <  2:4295586/93‑1 (MQ=255)

                        ggCGTGACCTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGAtcgtcg                                                                  <  2:3290346/93‑1 (MQ=255)

                         gCGTGACCTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGc                                                                 <  1:4884340/93‑1 (MQ=255)

                         gCGTGACCTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGc                                                                 <  2:1060/93‑1 (MQ=255)

                            tGACCTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGc                                                              <  1:3837444/93‑1 (MQ=255)

                             gACCTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGcc                                                             >  1:4815510/1‑93 (MQ=255)

                              aCCTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGccc                                                            <  1:6051192/93‑1 (MQ=255)

                              aCCTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGccc                                                            >  2:5312834/1‑93 (MQ=255)

                              aCCTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGccc                                                            <  1:2660544/93‑1 (MQ=255)

                               ccTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCt                                                           >  1:1832033/1‑93 (MQ=255)

                               ccTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCt                                                           <  1:1702984/93‑1 (MQ=255)

                                cTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTg                                                          >  2:3704684/1‑93 (MQ=255)

                                cTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTg                                                          <  2:450242/93‑1 (MQ=255)

                                  gTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGtt                                                        >  1:2804068/1‑93 (MQ=255)

                                  gTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGtt                                                        <  1:4587535/93‑1 (MQ=255)

                                   tAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTc                                                       <  2:439229/93‑1 (MQ=255)

                                    aaCCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCa                                                      >  1:4660044/1‑93 (MQ=255)

                                      ccccTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATc                                                    >  2:759154/1‑93 (MQ=255)

                                      ccccTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATc                                                    >  2:6204349/1‑93 (MQ=255)

                                      ccccTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATc                                                    <  1:5694188/93‑1 (MQ=255)

                                      ccccTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATc                                                    <  1:3486500/93‑1 (MQ=255)

                                       cccTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATcc                                                   <  1:2607511/93‑1 (MQ=255)

                                       cccTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGAGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATcc                                                   >  2:4106825/1‑93 (MQ=255)

                                          tGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCcggc                                                <  2:3033896/93‑1 (MQ=255)

                                          tGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCcggc                                                <  2:2528339/93‑1 (MQ=255)

                                          tGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCcggc                                                >  1:5149161/1‑93 (MQ=255)

                                            cAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCcggcgg                                              <  1:5915742/93‑1 (MQ=255)

                                            cAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCcggcgg                                              <  1:1208303/93‑1 (MQ=255)

                                            cAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCcggcgg                                              >  1:246744/1‑93 (MQ=255)

                                             aaaTCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATccggccgg                                             >  1:2181344/1‑93 (MQ=255)

                                              aaTCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGc                                            >  2:585754/1‑93 (MQ=255)

                                               aTCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGcc                                           >  2:3558019/1‑93 (MQ=255)

                                               aTCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGcc                                           >  1:4641010/1‑93 (MQ=255)

                                               aTCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAAAAGATCGTCGCAGCCCTGTTCATCCGGCGGGcc                                           >  2:6175453/1‑93 (MQ=255)

                                                 cAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGt                                         >  2:4446424/1‑93 (MQ=255)

                                                 cAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGGTTAAACAGATCGTCGCATCCCTGTTCATCCGGCGGGCCGt                                         >  1:284377/1‑93 (MQ=255)

                                                  aaaCACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGCTCATCCGGCGGGCCGTa                                        >  2:578737/1‑93 (MQ=255)

                                                   aaCACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTaa                                       >  2:390017/1‑93 (MQ=255)

                                                    acacCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTaaa                                      >  2:4331622/1‑93 (MQ=255)

                                                     caccacAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAAc                                     <  2:1799916/93‑1 (MQ=255)

                                                       ccacAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCg                                   >  2:2591252/1‑93 (MQ=255)

                                                       ccacAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCg                                   >  2:13738/1‑93 (MQ=255)

                                                       ccacAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCg                                   >  1:3246760/1‑93 (MQ=255)

                                                         acaAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCggg                                 >  2:651932/1‑93 (MQ=255)

                                                          caAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCGgga                                <  1:6429546/93‑1 (MQ=255)

                                                           aaaCGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCGggat                               <  1:1191391/93‑1 (MQ=255)

                                                           aaaCGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCGggat                               >  1:1043180/1‑93 (MQ=255)

                                                            aaCGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGCTGTTTAAACAGATCGTCGCAGCCCTGTTCATTCCGCGGGCCGTAAACCGggatc                              >  1:3089692/1‑93 (MQ=255)

                                                            aaCGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCGggatc                              <  2:2969911/93‑1 (MQ=255)

                                                            aaCGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCGggatc                              >  2:5495016/1‑93 (MQ=255)

                                                            aaCGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCGggatc                              >  1:2928867/1‑93 (MQ=255)

                                                            aaCGGTTCCACCGTGTGGCGGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCGggatc                              <  1:3364282/93‑1 (MQ=255)

                                                             aCGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCGggatcg                             <  2:2877865/93‑1 (MQ=255)

                                                             aCGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCGggatcg                             >  2:5145594/1‑93 (MQ=255)

                                                                gTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCGggatcggat                          >  2:395257/1‑93 (MQ=255)

                                                                gTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCGggatcggat                          <  1:1387967/93‑1 (MQ=255)

                                                                 ttCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCGggatcggat‑c                        <  1:5399845/93‑1 (MQ=255)

                                                                   ccACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCGGGATCGGAT‑CGc                      >  1:854170/1‑93 (MQ=255)

                                                                   ccACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCGGGATCGGAT‑CGc                      >  1:5192640/1‑93 (MQ=255)

                                                                    cACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCGGGATCGGAT‑CGcc                     <  2:4604565/93‑1 (MQ=255)

                                                                     aCCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCGGGATCGGAT‑CGCCa                    >  2:879139/1‑93 (MQ=255)

                                                                     aCCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCGGGATCGGAT‑CGCCa                    >  2:3056912/1‑93 (MQ=255)

                                                                      ccGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCGGGATCGGAT‑CGCCaa                   >  1:219641/1‑93 (MQ=255)

                                                                       cGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCGGGATCGGAT‑CGCCAAc                  >  1:6256380/1‑93 (MQ=255)

                                                                        gtgtGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGGCGTAAACCGGGCATGGATCCGCCAAc                  >  1:157611/1‑93 (MQ=255)

                                                                           tgGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCGGGATCGGAT‑CGCCAACcccc              >  1:5276471/1‑93 (MQ=255)

                                                                            ggCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCGGGATCGGAT‑CGCCAACGcccc             <  2:2601863/93‑1 (MQ=255)

                                                                            ggCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCGGGATCGGAT‑CGCCAACGcccc             >  2:1671567/1‑93 (MQ=255)

                                                                             gCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCGGGATCGGAT‑CGCCAACGCCCca            <  1:2525137/93‑1 (MQ=255)

                                                                              cTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCGGGATCGGAT‑CGCCCACGCCCcag           >  1:686523/1‑93 (MQ=255)

                                                                              cTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCGGGATCGGAT‑CGCCAACGCCCcag           <  1:1250609/93‑1 (MQ=255)

                                                                              cTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCGGGATCGGAT‑CGCCAACGCCCcag           <  1:4040199/93‑1 (MQ=255)

                                                                              cTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCGGGATCGGAT‑CGCCAACGCCCcag           >  2:477956/1‑93 (MQ=255)

                                                                                 aaaTCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCGGGATCGGAT‑CGCCAACGCCCcagcgc        <  1:4109459/93‑1 (MQ=255)

                                                                                  aaTCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCGGGATCGGAT‑CGCCAACGCCCcagcgca       >  2:6365251/1‑93 (MQ=255)

                                                                                   aTCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCGGGATCGGAT‑CGCCAACGCCCcagcgcag      <  1:3646862/93‑1 (MQ=255)

                                                                                     cAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCGGGATCGGAT‑CGCCAACGCCCcagcgcagcg    <  1:4281970/93‑1 (MQ=255)

                                                                                      aaGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCGGGATCGGAT‑CGCCAACGCCCCAGCGCAGCgg   <  1:5062919/93‑1 (MQ=255)

                                                                                      aaGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCGGGATCGGAT‑CGCCAACGCCCCAGCGCAGCgg   >  2:248331/1‑93 (MQ=255)

                                                                                      aaGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCGGGATCGGAT‑CGCCAACGCCCCAGCGCAGCgg   >  2:3362721/1‑93 (MQ=255)

AGTTTTGAGTGGTTGAGCGGCAGGGGCGTGACCTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCAGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCGGGATCGGAT‑CGCCAACGCCCCAGCGCAGCGGA  >  NC_000913/4312244‑4312423

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