BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,965,001	T→C	M486V (ATG→GTG)	ptsP ←	PtsP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,965,001	0	T	C	99.7%	1156.9 / ‑5.9	368	M486V (ATG→GTG)	ptsP	PtsP
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (191/176); total (191/177)
Fisher's exact test for biased strand distribution p-value = 4.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 385 total aligned reads displayed.

CTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCAT  >  NC_000913/2964911‑2965090

ctTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCAccc                                                                                         <  2:2634111/93‑1 (MQ=255)

 ttCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCAcccc                                                                                        <  2:287064/93‑1 (MQ=255)

  tCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACcccc                                                                                       >  2:4442044/1‑93 (MQ=255)

  tCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACcccc                                                                                       <  1:7149132/93‑1 (MQ=255)

   cGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCt                                                                                      >  2:3587594/1‑93 (MQ=255)

   cGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCt                                                                                      <  1:5311470/93‑1 (MQ=255)

    gCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTg                                                                                     >  1:2175665/1‑93 (MQ=255)

     cTGATCGGCATGTAAGGCATCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGATCATTTGCAGCACCCCCTGa                                                                                    >  1:7033933/1‑93 (MQ=255)

      tGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGAt                                                                                   >  1:4753904/1‑93 (MQ=255)

       gATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATa                                                                                  >  2:4741834/1‑93 (MQ=255)

       gATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATa                                                                                  >  2:7638518/1‑93 (MQ=255)

        aTCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATAc                                                                                 <  2:5357434/93‑1 (MQ=255)

         tCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACt                                                                                >  1:45347/1‑93 (MQ=255)

         tCGGCATGGAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACt                                                                                <  1:374261/93‑1 (MQ=255)

          cGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTg                                                                               <  2:2272263/93‑1 (MQ=255)

           ggCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTgc                                                                              <  2:4535014/93‑1 (MQ=255)

            gCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTgcg                                                                             >  2:7883533/1‑93 (MQ=255)

            gCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTgcg                                                                             >  1:177592/1‑93 (MQ=255)

             cATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGCTACTgcgc                                                                            >  2:1585391/1‑93 (MQ=255)

               tGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCa                                                                          >  2:4192319/1‑93 (MQ=255)

               tGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCa                                                                          >  2:3344548/1‑93 (MQ=255)

               tGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCa                                                                          >  2:6494700/1‑93 (MQ=255)

                 tAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCAcc                                                                        >  1:6961592/1‑93 (MQ=255)

                   aGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTg                                                                      <  2:7336277/93‑1 (MQ=255)

                   aGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTg                                                                      >  2:5453184/1‑93 (MQ=255)

                   aGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTg                                                                      <  1:7749936/93‑1 (MQ=255)

                    ggCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGt                                                                     <  2:1968335/93‑1 (MQ=255)

                    ggCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGt                                                                     >  1:6664815/1‑93 (MQ=255)

                      cAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGttc                                                                   <  1:1261774/93‑1 (MQ=255)

                      cAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGttc                                                                   <  2:1663105/93‑1 (MQ=255)

                        gctgctTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGttctt                                                                 >  2:6761736/1‑93 (MQ=255)

                          tgctTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGttcttct                                                               >  1:5882273/1‑93 (MQ=255)

                           gctTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGttcttctt                                                              >  1:5581939/1‑93 (MQ=255)

                           gctTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGttcttctt                                                              >  1:4639632/1‑93 (MQ=255)

                           gctTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGttcttctt                                                              >  1:6535690/1‑93 (MQ=255)

                            ctTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGttcttcttc                                                             <  1:1672051/93‑1 (MQ=255)

                             ttATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTcc                                                            <  2:5580746/93‑1 (MQ=255)

                             ttATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTcc                                                            <  1:5185161/93‑1 (MQ=255)

                             ttATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTcc                                                            >  1:855178/1‑93 (MQ=255)

                                tCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGAc                                                         >  1:2961056/1‑93 (MQ=255)

                                tCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGAc                                                         >  2:5215149/1‑93 (MQ=255)

                                tCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCACTGCTACTGCGCCACCTGTTCTTCTTCCgcc                                                         >  2:2117412/1‑91 (MQ=255)

                                 cTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACg                                                        <  2:712843/93‑1 (MQ=255)

                                 cTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACg                                                        >  2:6396265/1‑93 (MQ=255)

                                   gCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACggg                                                      >  2:5682963/1‑93 (MQ=255)

                                     tCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGaa                                                    >  1:2286412/1‑93 (MQ=255)

                                     tCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCCACGGGaa                                                    >  1:2680759/1‑93 (MQ=255)

                                       cGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAAc                                                  <  1:2058664/93‑1 (MQ=255)

                                        gACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAAcc                                                 <  2:3493025/93‑1 (MQ=255)

                                        gACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAAcc                                                 >  1:360898/1‑93 (MQ=255)

                                        gACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAAcc                                                 >  2:2443118/1‑93 (MQ=255)

                                         gcATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCa                                                <  2:7761239/92‑1 (MQ=255)

                                           aTCCAGCGTACGCAGGGTGACGGGTTTATCAGTGAACATTTGCAGCACCCCCTGATTCTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACt                                              <  2:692902/93‑1 (MQ=255)

                                           aTCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACt                                              >  2:4111213/1‑93 (MQ=255)

                                            tCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACtt                                             >  1:2450719/1‑93 (MQ=255)

                                            tCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACtt                                             >  2:178313/1‑93 (MQ=255)

                                            tCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACtt                                             >  1:42174/1‑93 (MQ=255)

                                              cAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTg                                           <  1:6246029/93‑1 (MQ=255)

                                              cAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTg                                           >  2:6315417/1‑93 (MQ=255)

                                                  gTACGCGAGGTGACGGGTTTATCATTGAACATTTGCAGCCCCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTgcagc                                       <  1:5815304/93‑1 (MQ=255)

                                                  gTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTgcagc                                       >  2:4582815/1‑93 (MQ=255)

                                                  gTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTgcagc                                       >  2:2151946/1‑93 (MQ=255)

                                                    aCGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCAt                                     <  2:324945/93‑1 (MQ=255)

                                                    aCGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCAt                                     >  2:4518250/1‑93 (MQ=255)

                                                     cGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATg                                    <  1:2914843/93‑1 (MQ=255)

                                                      gCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGa                                   >  1:1384243/1‑93 (MQ=255)

                                                       cAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGaa                                  <  2:3847279/93‑1 (MQ=255)

                                                        aaGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAAt                                 <  1:2388087/93‑1 (MQ=255)

                                                        aaGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAAt                                 <  2:5953595/93‑1 (MQ=255)

                                                         aGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATg                                <  1:1562909/93‑1 (MQ=255)

                                                          ggTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATgg                               <  2:6140491/93‑1 (MQ=255)

                                                            tGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGa                             >  1:7586297/1‑93 (MQ=255)

                                                             gACGGTTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGAt                            <  2:52582/93‑1 (MQ=255)

                                                             gACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGACTGGGAt                            >  2:426555/1‑93 (MQ=255)

                                                             gACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGAt                            >  1:1821675/1‑93 (MQ=255)

                                                               cGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGAttt                          >  1:7439846/1‑93 (MQ=255)

                                                                gggTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTc                         <  1:1035434/93‑1 (MQ=255)

                                                                 ggTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCa                        >  2:1995869/1‑93 (MQ=255)

                                                                   tttATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGt                      >  1:3944704/1‑93 (MQ=255)

                                                                      aTCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCg                   >  2:1182519/1‑93 (MQ=255)

                                                                       tCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGa                  <  1:3485806/93‑1 (MQ=255)

                                                                       tCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGa                  >  2:3046633/1‑93 (MQ=255)

                                                                       tCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGa                  <  2:5318226/93‑1 (MQ=255)

                                                                        cATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGAt                 >  2:1396173/1‑93 (MQ=255)

                                                                        cATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGAt                 >  1:6497211/1‑93 (MQ=255)

                                                                        cATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGAt                 <  1:326713/93‑1 (MQ=255)

                                                                         aTTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGTGATTTCAGTGCGATa                >  2:6502321/1‑93 (MQ=255)

                                                                          ttGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATaa               >  2:1574517/1‑93 (MQ=255)

                                                                          ttGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATaa               >  1:7814067/1‑93 (MQ=255)

                                                                            gAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAg             <  1:5444123/93‑1 (MQ=255)

                                                                              aCATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGAc           <  2:107614/93‑1 (MQ=255)

                                                                                 tttGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGa        >  2:5088446/1‑93 (MQ=255)

                                                                                 tttGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGa        >  2:6636937/1‑93 (MQ=255)

                                                                                 tttGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGa        <  1:5683449/93‑1 (MQ=255)

                                                                                 tttGCAGCACCCCCTGATACTGCGCCACCGGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGa        <  1:5836864/93‑1 (MQ=255)

                                                                                  ttGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGAt       <  2:5390895/93‑1 (MQ=255)

                                                                                    gcagcaCCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGc     <  1:213287/93‑1 (MQ=255)

                                                                                    gcagcaCCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGc     >  1:1907739/1‑93 (MQ=255)

                                                                                      agcaCCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCa   >  1:4968173/1‑93 (MQ=255)

                                                                                      agcaCCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCa   <  1:1058472/93‑1 (MQ=255)

CTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCAT  >  NC_000913/2964911‑2965090

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