BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,314,600	C→T	R66H (CGC→CAC)	phnM ←	PhnM

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,314,600	0	C	T	87.3%	678.6 / 61.1	244	R66H (CGC→CAC)	phnM	PhnM
Reads supporting (aligned to +/- strand): ref base C (9/22); new base T (82/131); total (91/153)
Fisher's exact test for biased strand distribution p-value = 4.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.16e-01

Only 100 of 253 total aligned reads displayed.

ACGGCATCCAGTACGGTGGTGATGCCGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGG  >  NC_000913/4314510‑4314692

aCGGCATCCAGTACGGTGGTGATGCCGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGgtgc                                                                                            <  2:6143560/93‑1 (MQ=255)

 cGGCATCCAGTACGGTGGTGATGCCGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGtgcg                                                                                           >  1:4202698/1‑93 (MQ=255)

 cGGCATCCAGTACGGTGGTGATGCCGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGtgcg                                                                                           >  1:3896258/1‑93 (MQ=255)

 cGGCATCCAGTACGGTGGTGATGCCGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGcgcg                                                                                           >  2:3884419/1‑93 (MQ=255)

  ggCATCCAGTACGGTGGTGATGCCGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGTGCgg                                                                                          <  1:2800502/93‑1 (MQ=255)

  ggCATCCAGTACGGTGGTGATGCCGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGCGCgg                                                                                          >  1:5423212/1‑93 (MQ=255)

   gCATCCAGTACGGTGGTGATGCCGCTCGCCACCATCAGCGCGTCGTTGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGCGCggg                                                                                         <  1:1289663/93‑1 (MQ=255)

   gCATCCAGTACGGTGGTGATGCCGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGTGCggg                                                                                         <  1:2519288/93‑1 (MQ=255)

    cATCCAGTACGGTGGTGATGCCGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGTGCgggg                                                                                        <  2:2639312/93‑1 (MQ=255)

    cATCCAGTACGGTGGTGATGCCGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGTGCgggg                                                                                        >  2:3363295/1‑93 (MQ=255)

     aTCCAGTACGGTGGTGATGCCGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGTGCGGGGt                                                                                       <  1:4937164/93‑1 (MQ=255)

      tCCAGTACGGTGGTGATGCCGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGTGCGGGGTg                                                                                      >  2:6068842/1‑93 (MQ=255)

      tCCAGTACGGTGGTGATGCCGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGTGCGGGGTg                                                                                      >  1:1131432/1‑93 (MQ=255)

        cAGTACGGTGGTGATGCCGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGTGCGGGGTgaa                                                                                    <  1:4843882/93‑1 (MQ=255)

        cAGTACGGTGGTGATGCCGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGTGCGGGGTgaa                                                                                    <  2:6935520/93‑1 (MQ=255)

        cAGTACGGTGGTGATGCCGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGTGCGGGGTgaa                                                                                    <  1:2650425/93‑1 (MQ=255)

         aGTACGGTGGTGATGCCGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGTGCGGGGCgaag                                                                                   >  1:6396761/1‑93 (MQ=255)

         aGTACGGTGGTGATGCCGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGTGCGGGGAgaag                                                                                   <  2:1815848/93‑1 (MQ=255)

          gTACGGTGGTGATGCCGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGTGCGGGGTgaaga                                                                                  <  1:6792598/93‑1 (MQ=255)

           tACGGTGGTGATGCCGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGTGCGGGGTgaagaa                                                                                 >  1:1902246/1‑93 (MQ=255)

           tACGGTGGTGATGCCGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTgaagaa                                                                                 >  2:4167531/1‑93 (MQ=255)

            aCGGTGGTGATGCCGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGTGCGGGGTGAAGAAt                                                                                >  1:1805168/1‑93 (MQ=255)

             cGGTGGTGATGCCGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGTGCGGGGTGAAGAAtt                                                                               <  2:5647116/93‑1 (MQ=255)

             cGGTGGTGATGCCGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGTGCGGGGTGAAGAAtt                                                                               <  2:2803930/93‑1 (MQ=255)

               gtggtgATGCCGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGTGCGGGGTGAAGAATTTa                                                                             >  2:1842814/1‑93 (MQ=255)

               gtggtgATGCCGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGTGCGGGGTGAAGAATTTa                                                                             >  1:5512445/1‑93 (MQ=255)

               gtggtgATGCCGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTa                                                                             <  2:5474847/93‑1 (MQ=255)

                tggtgATGCCGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGTGCGGGGTGAAGAATTTAt                                                                            >  1:6227668/1‑93 (MQ=255)

                tggtgATGCCGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGTGCGGGGTGAAGAATTTAt                                                                            >  2:2970012/1‑93 (MQ=255)

                 ggtgATGCCGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGTGCGGGGTGAAGAATTTATc                                                                           <  2:5733876/93‑1 (MQ=255)

                  gtgATGCCGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGTGCGGGGTGAAGAATTTATcc                                                                          <  2:3468079/93‑1 (MQ=255)

                  gtgATGCCGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGTGCGGGGTGAAGAATTTATcc                                                                          <  1:283087/93‑1 (MQ=255)

                  gtgATGCCGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGTGCGGGGTGAAGAATTTATcc                                                                          >  2:5040977/1‑93 (MQ=255)

                   tgATGCCGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGTGCGGGGTGAAGAATTTATCCa                                                                         <  2:5861925/93‑1 (MQ=255)

                   tgATGCCGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGTGCGGGGTGAAGAATTTATCCa                                                                         <  1:5066553/93‑1 (MQ=255)

                    gATGCCGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGTGCGGGGTGAAGAATTTATCCAg                                                                        >  2:5255193/1‑93 (MQ=255)

                     aTGCCGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGa                                                                       <  2:1435529/93‑1 (MQ=255)

                      tGCCGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGTGCGGGGTGAAGAATTTATCCAGAt                                                                      <  2:6315191/93‑1 (MQ=255)

                        ccGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGTGCGGGGTGAAGAATTTATCCAGATTa                                                                    <  1:5763657/93‑1 (MQ=255)

                         cGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTAt                                                                   <  2:7110777/93‑1 (MQ=255)

                          nntCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGTGCGGGGTGAAGAATTTATCCAGATTATc                                                                  <  2:2084866/91‑1 (MQ=255)

                          gCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGTGCGGGGTGAAGAATTTATCCAGATTATc                                                                  <  1:494675/93‑1 (MQ=255)

                          gCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATc                                                                  <  1:5310808/93‑1 (MQ=255)

                           cTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGTGCGGGGTGAAGAATTTATCCAGATTATCg                                                                 <  2:7250338/93‑1 (MQ=255)

                           cTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGTGCGGGGTGAAGAATTTATCCAGATTATCg                                                                 <  2:4361218/93‑1 (MQ=255)

                           cTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGTGCGGGGTGAAGAATTTATCCAGATTATCg                                                                 <  2:6820522/93‑1 (MQ=255)

                           cTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCg                                                                 <  2:6305507/93‑1 (MQ=255)

                           cTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCg                                                                 <  1:4281976/93‑1 (MQ=255)

                            tCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGTGCGGGGTGAAGAATTTATCCAGATTATCgg                                                                <  2:6722034/93‑1 (MQ=255)

                            tCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGTGCGGGGTGAAGAATTTATCCAGATTATCgg                                                                >  2:136242/1‑93 (MQ=255)

                             cGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGTGCGGGGTGAAGAATTTATCCAGATTATCGGt                                                               <  1:7217767/93‑1 (MQ=255)

                               ccaccaTCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGTGCGGGGTGAAGAATTTATCCAGATTATCGGTAt                                                             <  1:1477596/93‑1 (MQ=255)

                                caccaTCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGTGCGGGGTGAAGAATTTATCCAGATTATCGGTATg                                                            >  1:7050228/1‑93 (MQ=255)

                                 accaTCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGc                                                           >  1:7243945/1‑93 (MQ=255)

                                 accaTCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTAAGc                                                           >  2:5265640/1‑93 (MQ=255)

                                  ccaTCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGTGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCa                                                          >  2:4757294/1‑93 (MQ=255)

                                  ccaTCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGTGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCa                                                          <  1:4409992/93‑1 (MQ=255)

                                  ccaTCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGTGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCa                                                          <  2:5398861/93‑1 (MQ=255)

                                    aTCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGTGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGc                                                        >  2:3608036/1‑93 (MQ=255)

                                    aTCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGTGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGc                                                        <  2:2713152/93‑1 (MQ=255)

                                     tCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGTGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCt                                                       >  1:194709/1‑93 (MQ=255)

                                       aGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGTGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCa                                                     <  2:4849855/93‑1 (MQ=255)

                                       aGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGTGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCa                                                     >  2:1509324/1‑93 (MQ=255)

                                        gcgcGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGTGCGGGGGGAAGAATTTATCCAGATTATCGGTATGCAGCTCaa                                                    >  1:74932/1‑93 (MQ=255)

                                         cgcgTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGTGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAAt                                                   <  1:1504387/93‑1 (MQ=255)

                                          gcgTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGTGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATc                                                  <  2:7215587/93‑1 (MQ=255)

                                          gcgTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATc                                                  <  1:2987165/93‑1 (MQ=255)

                                           cgtcgtGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGTGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCa                                                 <  2:1691720/93‑1 (MQ=255)

                                           cgtcgtGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGTGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCa                                                 <  2:5978502/93‑1 (MQ=255)

                                           cgtcgtGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCa                                                 >  2:1078461/1‑93 (MQ=255)

                                            gtcgtGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGTGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAg                                                <  2:6432303/93‑1 (MQ=255)

                                            gtcgtGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGTGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAg                                                >  1:5136807/1‑93 (MQ=255)

                                                    tgctCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGTGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGgca                                        <  1:6043395/93‑1 (MQ=255)

                                                     gctCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGTGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGgcag                                       <  1:5840572/93‑1 (MQ=255)

                                                          tCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGTGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGcagcc                                  >  2:5146839/1‑93 (MQ=255)

                                                           cGCCGAGTGGGCAGGCCAGTCAACTTTCGGGTGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGcagcca                                 <  1:2673701/93‑1 (MQ=255)

                                                             ccGAGTGGGCAGGCCAGTCAACTTTCGGGTGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGcagccagc                               >  1:2229767/1‑93 (MQ=255)

                                                             ccGAGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGcagccagc                               <  2:7440459/93‑1 (MQ=255)

                                                              cGAGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCagccagcc                              <  1:206663/93‑1 (MQ=255)

                                                                aGTGGGCAGGCCAGTCAACTTTCGGGTGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAgccgc                            <  1:1797106/93‑1 (MQ=255)

                                                                 gTGGGCAGGCCAGTCAACTTTCGGGTGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAgccgcc                           >  2:6335119/1‑93 (MQ=255)

                                                                  tGGGCAGGCCAGTCAACTTTCGGGTGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCt                          >  1:3838738/1‑93 (MQ=255)

                                                                     gCAGGCCAGTCAACTTTCGGGTGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCg                       <  2:3305999/93‑1 (MQ=255)

                                                                       aGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGcc                     <  1:4454860/93‑1 (MQ=255)

                                                                        ggCCAGTCAACTTTCGGGTGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCg                    >  2:5723811/1‑93 (MQ=255)

                                                                         gCCAGTCAACTTTCGGGTGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGt                   <  1:3636194/93‑1 (MQ=255)

                                                                            aGTCAACTTTCGGGTGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCa                <  1:2424671/93‑1 (MQ=255)

                                                                             gTCAACTTTCGGGTGCGGGGTGAAGAATTTATCCAGATTGTCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCAt               <  1:81664/93‑1 (MQ=255)

                                                                             gTCAACTTTCGGGTGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCAt               <  1:4621252/93‑1 (MQ=255)

                                                                             gTCAACTTTCGGGTGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCAt               >  1:3188737/1‑93 (MQ=255)

                                                                                 aCTTTCGGGTGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGcc           <  1:659189/93‑1 (MQ=255)

                                                                                  cTTTCGGGTGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCt          <  2:42831/93‑1 (MQ=255)

                                                                                   tttCGGGTGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTc         <  1:5045680/93‑1 (MQ=255)

                                                                                   tttCGGGTGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTc         <  2:2048875/93‑1 (MQ=255)

                                                                                    ttCGGGTGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTcc        <  2:2260364/93‑1 (MQ=255)

                                                                                    ttCGGGTGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTcc        >  1:5601027/1‑93 (MQ=255)

                                                                                    ttCGGGTGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTcc        >  2:2929532/1‑93 (MQ=255)

                                                                                        ggCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCa    <  2:4012987/93‑1 (MQ=255)

                                                                                         gcgcGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAg   <  1:5192952/93‑1 (MQ=255)

                                                                                          cgcgGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAgg  <  1:6291397/93‑1 (MQ=255)

ACGGCATCCAGTACGGTGGTGATGCCGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGG  >  NC_000913/4314510‑4314692

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