BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,619,453	T→G	V39G (GTG→GGG)	yjjJ →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,619,453	0	T	G	99.4%	566.5 / ‑5.1	174	V39G (GTG→GGG)	yjjJ	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/1); new base G (80/93); total (80/94)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.37e-01

Only 100 of 175 total aligned reads displayed.

TGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGAC  >  NC_000913/4619363‑4619541

tGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGTCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGGGa                                                                                        <  1:3629741/93‑1 (MQ=255)

tGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGGGa                                                                                        >  2:3487928/1‑93 (MQ=255)

tGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGGGa                                                                                        >  2:2365653/1‑93 (MQ=255)

 gCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGGGAt                                                                                       <  2:622566/93‑1 (MQ=255)

 gCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGGGAt                                                                                       >  1:2850594/1‑93 (MQ=255)

   aaGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGGGATTc                                                                                     <  2:3369660/93‑1 (MQ=255)

    aGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGGGATTCg                                                                                    >  2:3983431/1‑93 (MQ=255)

    aGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGGGATTCg                                                                                    >  1:2948837/1‑93 (MQ=255)

        ccGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGGGATTCGCttt                                                                                >  1:392641/1‑93 (MQ=255)

         cgcgTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGGGATTCGCTTTg                                                                               <  1:3192713/93‑1 (MQ=255)

          gcgTTCTGCCCCGGAATGGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGGGATTCGCTTTgg                                                                              <  1:1973153/93‑1 (MQ=255)

           cgTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGGGATTCGCTTTGGt                                                                             <  2:3172697/93‑1 (MQ=255)

           cgTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGGGATTCGCTTTGGt                                                                             >  2:3041447/1‑93 (MQ=255)

           cgTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGGGATTCGCTTTGGt                                                                             <  1:732426/93‑1 (MQ=255)

           cgTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGGGATTCGCTTTGGt                                                                             >  2:4638851/1‑93 (MQ=255)

            gTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGGGATTCGCTTTGGTa                                                                            >  1:4709862/1‑93 (MQ=255)

            gTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGGGATTCGCTTTGGTa                                                                            >  1:255884/1‑93 (MQ=255)

            gTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGGGATTCGCTTTGGTa                                                                            >  1:2248618/1‑93 (MQ=255)

             ttCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGGGATTCGCTTTGGTaa                                                                           >  1:2211470/1‑93 (MQ=255)

              tCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGGGATTCGCTTTGGTaaa                                                                          <  1:2836327/93‑1 (MQ=255)

                tGCCCCGGAATTGCGCCATCGTCGGGCAATCAGTCATGCGACGTTCTCACGCCTTGTTGCCAGCGAACATCGGGGGATTCGCTTTGGTAAAGc                                                                        <  2:338572/93‑1 (MQ=255)

                 gCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGGGATTCGCTTTGGTAAAGCa                                                                       <  2:2933947/93‑1 (MQ=255)

                  ccccGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGGGATTCGCTTTGGTAAAGCAc                                                                      <  2:449985/93‑1 (MQ=255)

                  ccccGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGGGATTCGCTTTGGTAAAGCAc                                                                      <  1:602843/93‑1 (MQ=255)

                   cccGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCACg                                                                     <  2:2079997/93‑1 (MQ=255)

                    ccGGTAATGCGCCAGAATCAGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGGGATTAGCTTTAGTCAGGCCCgg                                                                    <  2:4226243/93‑1 (MQ=255)

                     cGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGGGATTCGCTTTGGTAAAGCACggg                                                                   <  1:924087/93‑1 (MQ=255)

                      ggAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGGGATTCGCTTTGGTAAAGCACGGGc                                                                  >  1:2564539/1‑93 (MQ=255)

                       gAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGGGATTCGCTTTGGTAAAGCACGGGCa                                                                 >  1:120244/1‑93 (MQ=255)

                       gAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGGGATTCGCTTTGGTAAAGCACGGGCa                                                                 >  1:3530555/1‑93 (MQ=255)

                        aaTTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGGGATTCGCTTTGGTAAAGCACGGGCaa                                                                <  2:185226/93‑1 (MQ=255)

                        aaTTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGGGATTCGCTTTGGTAAAGCACGGGCaa                                                                <  2:4141172/93‑1 (MQ=255)

                         aTTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGGGATTCGCTTTGGTAAAGCACGGGCAAc                                                               <  1:1695111/93‑1 (MQ=255)

                          ttGCGCCAGCGTCTGGCGATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGGGATTCGCTTTGGTAAAGCACGGGCAAcg                                                              <  2:2804287/93‑1 (MQ=255)

                          ttGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGGGATTCGCTTTGGTAAAGCACGGGCAAcg                                                              >  1:1811086/1‑93 (MQ=255)

                            gcgcCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGGGATTCGCTTTGGTAAAGCACGGGCAAcgcg                                                            <  2:1559284/93‑1 (MQ=255)

                             cgcCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGGGATTCGCTTTGGTAAAGCACGGGCAACGCGa                                                           >  1:4163974/1‑93 (MQ=255)

                              gcCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGGGATTCGCTTTGGTAAAGCACGGGCAACGCGat                                                          <  1:2800234/93‑1 (MQ=255)

                               ccAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGGGATTCGCTTTGGTAAAGCACGGGCAACGCGata                                                         <  1:1439991/93‑1 (MQ=255)

                                cAGCGTCAGGCAATCAGTCAGGCGACGTGCTCACGCCTTGTTGCCAGAGAAGATCGGGGGATTCGCTTTGGTAAAGCACGGGCAACGCGatat                                                        <  1:2021032/93‑1 (MQ=255)

                                 aGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGGGATTCGCTTTGGTAAAGCACGGGCAACGCGATATg                                                       >  2:2426415/1‑93 (MQ=255)

                                  gCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGGGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGc                                                      <  2:2544631/93‑1 (MQ=255)

                                   cGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGGGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCa                                                     <  2:866741/93‑1 (MQ=255)

                                   cGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGGGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCa                                                     >  2:2901962/1‑93 (MQ=255)

                                    gTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGGGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCAc                                                    >  2:1707337/1‑93 (MQ=255)

                                    gTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGGGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCAc                                                    <  1:2127233/93‑1 (MQ=255)

                                    gTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGGGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCAc                                                    >  2:438104/1‑93 (MQ=255)

                                       tGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGGGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCActgc                                                 <  1:3972612/93‑1 (MQ=255)

                                        ggCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGGGATTCGCTTTGGTAAAGCACGGGCCACGCGATATGCActgct                                                >  2:610371/1‑93 (MQ=255)

                                        ggCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGGGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCActgct                                                <  2:4229394/93‑1 (MQ=255)

                                            aTCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGGGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGt                                            >  1:1234623/1‑93 (MQ=255)

                                              cagtcagGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGGGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTcc                                          >  1:1538261/1‑93 (MQ=255)

                                               agtcagGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGGGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCt                                         <  2:3789850/93‑1 (MQ=255)

                                                 tcagGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGGGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTa                                       <  2:1137224/93‑1 (MQ=255)

                                                  cagGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGGGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTAt                                      <  2:1388998/93‑1 (MQ=255)

                                                   agGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGGGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATc                                     >  1:3773498/1‑93 (MQ=255)

                                                   agGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGGGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATc                                     >  2:827072/1‑93 (MQ=255)

                                                   agGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGGGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATc                                     <  2:2743998/93‑1 (MQ=255)

                                                    ggCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGGGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATcg                                    >  1:3262214/1‑93 (MQ=255)

                                                    ggCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGGGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATcg                                    <  1:2391367/93‑1 (MQ=255)

                                                     gCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGGGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATcgc                                   <  1:2593899/93‑1 (MQ=255)

                                                       gACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGGGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCgg                                 >  2:1716599/1‑93 (MQ=255)

                                                       gACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGGGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCgg                                 <  2:2770230/93‑1 (MQ=255)

                                                       gACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGGGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCgg                                 >  2:1014809/1‑93 (MQ=255)

                                                        aCGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGGGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGa                                >  1:2364936/1‑93 (MQ=255)

                                                          gTTCTCACGCCTTGTTGCCAGAGAAGATCGGGGGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAAt                              <  1:4546551/93‑1 (MQ=255)

                                                          gTTCTCACGCCTTGTTGCCAGAGAAGATCGGGGGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAAt                              >  1:3128242/1‑93 (MQ=255)

                                                           ttCTCACGCCTTGTTGCCAGAGAAGATCGGGGGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAAtt                             >  2:793163/1‑93 (MQ=255)

                                                            tctcACGCCTTGTTGCCAGAGAAGATCGGGGGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTg                            >  1:4274489/1‑93 (MQ=255)

                                                             ctcACGCCTTGTTGCCAGAGAAGATCGGGGGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGa                           <  1:3851735/93‑1 (MQ=255)

                                                             ctcACGCCTTGTTGCCAGAGAAGATCGGGGGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGa                           <  2:1811329/93‑1 (MQ=255)

                                                              tcACGCCTTGTTGCCAGAGAAGATCGGGGGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAg                          <  2:2694489/93‑1 (MQ=255)

                                                              tcACGCCTTGTTGCCAGAGAAGATCGGGGGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAg                          >  1:1894453/1‑93 (MQ=255)

                                                               cACGCCTTGTTGCCAGAGAAGATCGGGGGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGc                         <  2:4471394/93‑1 (MQ=255)

                                                               cACGCCTTGTTGCCAGAGAAGATCGGGGGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGc                         <  1:387708/93‑1 (MQ=255)

                                                               cACGCCTTGTTGCCAGAGAAGATCGGGGGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGc                         >  2:1235594/1‑93 (MQ=255)

                                                               cACGCCTTGTTGCCAGAGAAGATCGGGGGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGc                         <  2:376188/93‑1 (MQ=255)

                                                                  gCCTTGTTGCCAGAGAAGATCGGGGGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTa                      <  2:196203/93‑1 (MQ=255)

                                                                     ttgttgCCAGAGAAGATCGGGGGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTc                   >  2:3627540/1‑93 (MQ=255)

                                                                     ttgttgCCAGAGAAGATCGGGGGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTc                   <  1:404566/93‑1 (MQ=255)

                                                                        ttgCCAGAGAAGATCGGGGGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCg                >  2:3624979/1‑93 (MQ=255)

                                                                         tgCCAGAGAAGATCGGGGGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGt               <  2:898910/93‑1 (MQ=255)

                                                                         tgCCAGAGAAGATCGGGGGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCC‑TATCGCGGAATTGAGCGTATTCCCGTc              <  1:4602141/93‑1 (MQ=255)

                                                                          gCCAGAGAAGATCGGGGGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTc              <  1:141437/93‑1 (MQ=255)

                                                                          gCCAGAGAAGATCGGGGGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTc              <  2:573613/93‑1 (MQ=255)

                                                                          gCCAGAGAAGATCGGGGGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTc              <  2:3671894/93‑1 (MQ=255)

                                                                           ccAGAGAAGATCGGGGGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCt             <  1:3432773/93‑1 (MQ=255)

                                                                           ccAGAGAAGATCGGGGGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCt             <  2:4311605/93‑1 (MQ=255)

                                                                            cAGAGAAGATCGGGGTATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCTCGGAATGGAGCGTATTCCCGTCTg            >  2:2522355/1‑93 (MQ=255)

                                                                             agagAAGATCGGGGGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTgg           <  2:2695461/93‑1 (MQ=255)

                                                                             agagAAGATCGGGGGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTgg           >  1:4077384/1‑93 (MQ=255)

                                                                              gagaAGATCGGGGGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGc          <  2:753274/93‑1 (MQ=255)

                                                                              gagaAGATCGGGGGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGc          <  1:706203/93‑1 (MQ=255)

                                                                               agaagaTCGGGGGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCg         >  1:4499432/1‑93 (MQ=255)

                                                                                  agaTCGGGGGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGt      <  2:844/93‑1 (MQ=255)

                                                                                   gaTCGGGGGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCTTCTGGCGGGTg     >  2:3971719/1‑93 (MQ=255)

                                                                                   gaTCGGGGGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTg     >  1:3194989/1‑93 (MQ=255)

                                                                                    aTCGGGGGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTgg    >  1:596943/1‑93 (MQ=255)

                                                                                     tCGGGGGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGa   <  1:1545265/93‑1 (MQ=255)

                                                                                      cGGGGGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGag  >  1:3530882/1‑92 (MQ=255)

TGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGAC  >  NC_000913/4619363‑4619541

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