BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,654,651	T→C	I36V (ATT→GTT)	yfhJ ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,654,651	0	T	C	99.4%	540.1 / NA	173	I36V (ATT→GTT)	yfhJ	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (90/82); total (90/83)
Fisher's exact test for biased strand distribution p-value = 4.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 179 total aligned reads displayed.

CGGCCTCGTCCAGCCAGACTAACAAAATCGCTTCGAGGATTTTCTCGTTGGATGCCTGCGGGTCGTCGTCGAAATCTTCCAGATCGCAAATCCACTGATGCATATCGGTGAATCGAACCGTTTTCGGATCAAGATCGGGATACGCATCGTACAGTGCTTCGCCAATTTCGCGGCTATCGG  >  NC_000913/2654561‑2654740

cGGCCTCGTCCAGCCAGACTAACAAAATCGCTTCGAGGATTTTCTCGTTGGATGCCTGCGGGTCGTCGTCGAAATCTTCCAGATCGCAAAccc                                                                                         <  1:107416/93‑1 (MQ=255)

    cTCGTCCAGCCAGACTAACAAAATCGCTTCGAGGATTTTCTCGTTGGATGCCTGCGGGTCGTCGTCGAAATCTTCCAGATCGCAAACCCACTg                                                                                     >  2:3000903/1‑93 (MQ=255)

    cTCGTCCAGCCAGACTAACAAAATCGCTTCGAGGATTTTCTCGTTGGATGCCTGCGGGTCGTCGTCGAAATCTTCCAGATCGCAAACCCACTg                                                                                     >  2:4508279/1‑93 (MQ=255)

     tCGTCCAGCCAGACTAACAAAATCGCTTCGAGGATTTTCTCGTTGGATGCCTGCGGGTCGTCGTCGAAATCTTCCAGATCGCAAACCCACTGa                                                                                    <  1:16452/93‑1 (MQ=255)

     ccGTCCAGCCAGACTAACAAAATCGCTTCGAGGATTTTCTCGTTGGATGCCTGCGGGTCGTCGTCGAAATCTTCCAGATCGCAAACCCACTGa                                                                                    <  1:3523469/92‑1 (MQ=255)

         ccagccagACTAACAAAATCGCTTCGAGGATTTTCTCGTTGGATGCCTGCGGGTCGTCGTCGAAATCTTCCAGATCGCAAACCCACTGATGCa                                                                                <  1:4264133/93‑1 (MQ=255)

          cagccagACTAACAAAATCGCTTCGAGGATTTTCTCGTTGGATGCCTGCGGGTCGTCGTCGAAATCTTCCAGATCGCAAACCCGCTGATGCat                                                                               <  2:735949/93‑1 (MQ=255)

           agccagACTAACAAAATCGCTTCGAGGATTTTCTCGTTGGATGCCTGCGGGTCGTCGTCGAAATCTTCCAGATCGCAAACCCACTGATGCata                                                                              >  2:603610/1‑93 (MQ=255)

            gccagACTAACAAAATCGCTTCGAGGATTTTCTCGTTGGATGCCTGCGGGTCGTCGTCGAAATCTTCCAGATCGCAAACCCACTGATGCatat                                                                             >  2:2835243/1‑93 (MQ=255)

              cagACTAACAAAATCGCTTCGAGGATTTTCTCGTTGGATGCCTGCGGGTCGTCGTCGAAATCTTCCAGATCGCAAACCCACTGATGCATATCg                                                                           >  1:702196/1‑93 (MQ=255)

              cagACTAACAAAATCGCTTCGAGGATTTTCTCGTTGGATGCCTGCGGGTCGTCGTCGAAATCTTCCAGATCGCAAACCCACTGATGCATATCg                                                                           >  2:477218/1‑93 (MQ=255)

               agACTAACAAAATCGCTTCGAGGATTTTCTCGTTGGATGCCTGCGGGTCGTCGTCGAAATCTTCCAGATCGCAAACCCACTGATGCATATCgg                                                                          >  2:3968992/1‑93 (MQ=255)

                gACTAACAAAATCGCTTCGAGGATTTTCTCGTTGGATGCCTGCGGGTCGTCGTCGAAATCTTCCAGATCGCAAACCCACTGATGCATATCGGt                                                                         >  1:3775513/1‑93 (MQ=255)

                gACTAACAAAATCGCTTCGAGGATTTTCTCGTTGGATGCCTGCGGGTCGTCGTCGAAATCTTCCAGATCGCAAACCCACTGATGCATATCGGt                                                                         >  2:2012400/1‑93 (MQ=255)

                gACTAACAAAATCGCTTCGAGGATTTTCTCGTTGGATGCCTGCGGGTCGTCGTCGAAATCTTCCAGATCGCAAACCCACTGATCCATATCggg                                                                         >  1:1486900/1‑92 (MQ=255)

                  cTAACAAAATCGCTTCGAGGATTTTCTCGTTGGATGCCTGCGGGTCGTTGTCGAAATCTTCCAGATCGCAAACCCACTGATGCATATCGGTGa                                                                       >  2:2506228/1‑93 (MQ=255)

                   tAACAAAATCGCTTCGAGGATTTTCTCGTTGGATGCCTGCGGGTCGTCGTCGAAATCTTCCAGATCGCAAACCCACTGATGCATATCGGTGaa                                                                      <  1:920021/93‑1 (MQ=255)

                   tAACAAAATCGCTTCGAGGATTTTCTCGTTGGATGCCTGCGGGTCGTCGTCGAAATCTTCCAGATCGCAAACCCACTGATGCATATCGGTGaa                                                                      <  2:4290277/93‑1 (MQ=255)

                   tAACAAAATCGCTTCGAGGATTTTCTCGTTGGATGCCTGCGGGTCGTCGTCGAAATCTTCCAGATCGCAAACCCACTGATGCATATCGGTGaa                                                                      <  2:3841137/93‑1 (MQ=255)

                       aaaaTCGCTTCGAGGATTTTCTCGTTGGATGCCTGCGGGTCGTCGTCGAAATCTTCCAGATCGCAAACCCACTGATGCATATCGGTGAATCGa                                                                  <  2:17372/93‑1 (MQ=255)

                        aaaTCGCTTCGAGGATTTTCTCGTTGGATGCCTGCGGGTCGTCGTCGAAATCTTCCAGATCGCAAACCCACTGATGCATATCGGTGAATCGaa                                                                 <  1:4155540/93‑1 (MQ=255)

                        aaaTCGCTTCGAGGATTTTCTCGTTGGATGCCTGCGGGTCGTCGTCGAAATCTTCCAGATCGCAAACCCACTGATGAATATCGGTGAATCGaa                                                                 >  2:507536/1‑93 (MQ=255)

                         aaTCGCTTCGAGGATTTTCTCGTTGGATGCCTGCGGGTCGTCGTCGAAATCTTCCAGATCGCAAACCCACTGATGCATATCGGTGAATCGAAc                                                                <  1:1954169/93‑1 (MQ=255)

                          aTCGCTTCGAGGATTTTCTAGTTGGATGCCTGCGGGTCGTCGTCGAAATCTTCCAGATCGCAAACCCACTGATGCATATCGGTGAATCGAAcc                                                               >  2:4584352/1‑93 (MQ=255)

                            cGCTTCGAGGATTTTCTCGTTGGATGCCTGCGGGTCGTCGTCGAAATCTTCCAGATCGCAAACCCACTGATGCATATCGGTGAATCGAACCGt                                                             >  1:1206689/1‑93 (MQ=255)

                            cGCTTCGAGGATTTTCTCGTTGGATGCCTGCGGGTCGTCGTCGAAATCTTCCAGATCGCAAACCCACTGATGCATATCGGTGAATCGAACCGt                                                             >  2:2470331/1‑93 (MQ=255)

                            cGCTTCGAGGATTTTCTCGTTGGATGCCTGCGGGTCGTCGTCGAAATCTTCCAGATCGCAAACCCACTGATGCATATCGGTGAATCGAACCGt                                                             <  1:3701838/93‑1 (MQ=255)

                              cTTCGAGGATTTTCTCGTTGGATGCCTGCGGGTCGTCGTCGAAATCTTCCAGATCGCAAACCCACTGATGCATATCGGTGAATCGAACCGttt                                                           <  2:1972158/93‑1 (MQ=255)

                               ttCGAGGATTTTCTCGTTGGATGCCTGCGGGTCGTCGTCGAAATCTTCCAGATCGCAAACCCACTGATGCATATCGGTGAATCGAACCGtttt                                                          >  1:2570626/1‑93 (MQ=255)

                                tCGAGGATTTTCTCGTTGGATGCCTGCGGGTCGTCGTCGAAATCTTCCAGATCGCAAACCCACTGATGCATATCGGTGAATCGAACCGTTTTc                                                         <  1:1057120/93‑1 (MQ=255)

                                    ggATTTTCTCGTTGGATGCCTGCGGGTCGTCGTCGAAATCTTCCAGATCGCAAACCCACTGATGCATATCGGTGAATCGAACCGTTTTCGGAt                                                     >  1:353963/1‑93 (MQ=255)

                                    ggATTTTCTCGTTGGATGCCTGCGGCTCGTCGTCGAAATCTTCCAGATCGCAAACCCACTGATGCATATCGGTGAATCGAACCGTTTTCGGAt                                                     >  2:2190368/1‑93 (MQ=255)

                                    ggATTTTCTCGTTGGATGCCTGAGGGTCGTCGTCGAAATCTTCCAGATCGCAAACCCACTGATGCATATCGGTGCATCGAACAGTTGTCGGAt                                                     >  1:3353459/1‑93 (MQ=255)

                                     gATTTTCTCGTTGGATGCCTGCGGGTCGTCGTCGAAATCTTCCAGATCGCAAACCCACTGATGCATATCGGTGAATCGAACCGTTTTCGGAta                                                    >  2:2321571/1‑92 (MQ=255)

                                     gATTTTCTCGTTGGATGCCTGCGGGTCGTCGTCGAAATCTTCCAGATCGCAAACCCACTGATGCATATCGGTGAATCGAACCGTTTTCGGATc                                                    <  1:2290906/93‑1 (MQ=255)

                                     gATTTTCTCGTTGGATGCCTGCGGGTCGTCGTCGAAATCTTCCAGATCGCAAACCCACTGATGCATATCGGTGAATCGAACCGTTTTCGGATc                                                    <  2:1885042/93‑1 (MQ=255)

                                      aTTTTCTCGTTGGATGCCTGCGGGTCGTCGTCGAAATCTTCCAGATCGCAAACCCACTGATGCATATCGGTGAATCGAACCGTTTTCGGATCa                                                   <  1:2779247/93‑1 (MQ=255)

                                       ttttCTCGTTGGATTCCTGCGGGTCGTCGTCGAAATCTTCCAGATCGCAAACCCACTGATGCATATCGGTGAATCGAACCGTTTTCGGATCaa                                                  <  1:1208626/93‑1 (MQ=255)

                                       ttttCTCGTTGGATGCCTGCGGGTCGTCGTCGAAATCTTCCAGATCGCAAACCCACTGATGCATATCGGTGAATCGAACCGTTTTCGGATCaa                                                  <  1:57850/93‑1 (MQ=255)

                                        tttCTCGTTGGATGCCTGCGGGTCGTCGTCGAAATCTTCCAGATCGCAAACCCACTGATGCATATCGGTGAATCGAACCGTTTTCGGATCAAg                                                 >  1:1262471/1‑93 (MQ=255)

                                         ttCTCGTTGGATGCCTGCGGGTCGTCGTCGAAATCTTCCAGATCGCAAACCCACTGATGCATATCGGTGAATCGAACCGTTTTCGGATCAAGa                                                <  2:2089048/93‑1 (MQ=255)

                                         ttCTCGTTGGATGCCTGCGGGTCGTCGTCGAAATCTTCCAGATCGCAAACCCACTGATGCATATCGGTGAATCGAACCGTTTTCGGATCAAGa                                                >  2:1322994/1‑93 (MQ=255)

                                          tctcGTTGGATGCCTGCGGGTCGTCGTCGAAATCTTCCAGATCGCAAACCCACTGATGCATATCGGTGAATCGAACCGTTTTCGGATCAAGAt                                               <  2:4407294/93‑1 (MQ=255)

                                          tctcGTTGGATGCCTGCGGGTCGTCGTCGAAATCTTCCAGATCGCAAACCCACTGATGCATATCGGTGAATCGAACCGTTTTCGGATCAAGAt                                               <  1:80914/93‑1 (MQ=255)

                                           ctcGTTGGATGCCTGCGGGTCGTCGTCGAAATCTTCCAGATCGCAAACCCACTGATGCATATCGGTGAATCGAACCGTTTTCGTATCAAGATc                                              <  1:2034850/93‑1 (MQ=255)

                                           ctcGTTGGATGCCTGCGGGTCGTCGTCGAAATCTTCCAGATCGCAAACCCACTGATGCATATCGGTGAATCGAACCGTTTTCGGATCAAGATc                                              <  2:1026993/93‑1 (MQ=255)

                                            tcGTTGGATGCCTGCGGGTCGTCGTCGAAATCTTCCAGATCGCAAACCCACTGATGCATATCGGTGAATCGAACCGTTTTCGGATCAAGATCg                                             <  2:2460376/93‑1 (MQ=255)

                                             cGTTGGATGCCTGCGGGTCGTCGTCGAAATCTTCCAGATCGCAAACCCACTGATGCATATCGGTGAATCGAACCGTTTTCGGATCAAGATCgg                                            <  1:2698486/93‑1 (MQ=255)

                                             cGTTGGATGCCTGCGGGTCGTCGTCGAAATCTTCCAGATCGCAAACCCACTGATGCATATCGGTGAATCGAACCGTTTTCGGATCAAGATCgg                                            <  1:3085647/93‑1 (MQ=255)

                                              gTTGGATGCCTGCGGGTCGTCGTCGAAATCTTCCAGATCGCAAACCCACTGATGCATATCGGTGAATCGAACCGTTTTCGGATCAAGATCggg                                           >  2:1871347/1‑93 (MQ=255)

                                              gTTGGATGCCTGCGGGTCGTCGTCGAAATCTTCCAGATCGCAAACCCACTGATGCATATCGGTGAATCGAACCGTTTTCGGATCAAGATCggg                                           >  2:1217046/1‑93 (MQ=255)

                                              gTTGGATGCCTGCGGGTCGTCGTCGAAATCTTCCAGATCGCAAACCCACTGATGCATATCGGTGAATCGAACCGTTTTCGGATCAAGATCggg                                           >  1:646697/1‑93 (MQ=255)

                                              gTTGGATGCCTGCGGGTCGTCGTCGAAATCTTCCAGATCGCAAACCCACTGATGCATATCGGTGAATCGAACCGCGTTCAGATTAAGATCggg                                           >  2:2486106/1‑93 (MQ=255)

                                               ttGGATGCCTGCGGGTCGTCGTCGAAATCTTCCAGATCGCAAACCCACTGATGCATATCGGTGAATCGAACCGTTTTCGGATCAAGATCGGGa                                          >  2:3899321/1‑93 (MQ=255)

                                                tGGATGCCTGCGGGTCGTCGTCGAAATCTTCCAGATCGCAAACCCACTGATGCATATCGGTGAATCGAACCGTTTTCGGATCAAGATCGGGAt                                         <  2:3258957/93‑1 (MQ=255)

                                                 ggATGCCTGCGGGTCGTCGTCGAAATCTTCCAGATCGCAAACCCACTGATGCATATCGGTGAATCGAACCGTTTTCGGATCAAGATCGGGATa                                        >  1:3847022/1‑93 (MQ=255)

                                                  gATGCCTGCGGGTCGTCGTCGAACTCTTCCAGATCGCAAACCCACTGATGCATATCGGTGATTCGAACCGTTTTCGGATCAAGACCGGGATAc                                       >  2:3582625/1‑93 (MQ=255)

                                                  gATGCCTGCGGGTCGTCGTCGAAATCTTCCAGATCGCAAACCCACTGATGCATATCGGTGAATCGAACCGTTTTCGGATCAAGATCGGGATAc                                       <  1:266460/93‑1 (MQ=255)

                                                     gCCTGCGGGTCGTCGTCGAAATCTTCCAGATCGCAAACCCACTGATGCATATCGGTGAATCGAACCGTTTTCGGATCAAGATCGGGATACGCa                                    >  2:1010889/1‑93 (MQ=255)

                                                      ccTGCGGGGCGTCGCCGAAATCTTCCAGATCGCAAACCCACTGATGCATATCGGTGAATCGAACCGTTTTCGGATCAAGATCGGGATACGCAt                                   <  1:3997601/93‑1 (MQ=255)

                                                       cTGCGGGTCGTCGTCGAAATCTTCCAGATCGCAAACCCACTGATGCATATCGGTGAATCGAACCGTTTTCGGATCAAGATCGGGATACGCATc                                  >  2:2397406/1‑93 (MQ=255)

                                                       cTGCGGGTCGTCGTCGAAATCTTCCAGATCGCAAACCCACTGATGCATATCGGTGAATCGAACCGGTTTCGGATCAAGATCGGGATACGCATc                                  >  2:2893243/1‑93 (MQ=255)

                                                        tGCGGGTCGTCGTCGAAATCTTCCAGATCGCAAACCCACTGATGCATATCGGTGAATCGAACCGTTTTCGGATCAAGATCGGGATACGCATCg                                 >  2:4026118/1‑93 (MQ=255)

                                                        tGCGGGTCGTCGTCGAAATCTTCCAGATCGCAAACCCACTGATGCATATCGGTGAATCGAACCGTTTTCGGATCAAGATCGGGATACGCATCg                                 >  1:848659/1‑93 (MQ=255)

                                                          cGGGTCGTCGTCGAAATCTTCCAGATCGCAAACCCACTGATGCATATCGGTGAATCGAACCGTTCTCGGATCAAGATCGGGATACGCATCGTa                               <  1:3240907/93‑1 (MQ=255)

                                                            ggTCGTCGTCGAAATCTTCCAGATCGCAAACCCACTGATGCATATCGGTGAATCGAACCGTTTTCGGATCAAGATCGGGATACGCATCGTACa                             <  2:1126106/93‑1 (MQ=255)

                                                            ggTCGTCGTCGAAATCTTCCAGATCGCAAACCCACTGATGCATATCGGTGAATCGAACCGTTTTCGGATCAAGATCGGGATACGCATCGTACa                             <  2:2442104/93‑1 (MQ=255)

                                                            ggTCGTCGTCGAAATCTTCCAGATCGCAAACCCACTGATGCATATCGGTGAATCGAACCGTTTTCGGATCAAGATCGGGATACGCATCGTACa                             >  2:3055148/1‑93 (MQ=255)

                                                            ggTCGTCGTCGAAATCTTCCAGATCGCAAACCCACTGATGCATATCGGTGAATCGAACCGTTTTCGGATCAAGATCGGGATACGCATCGTACa                             >  1:2973933/1‑93 (MQ=255)

                                                              gtggcgtcgAAATCTTCCAGATCGCAAACCCACTGATGCATATCGGTGAATCGAACCGTTTTCGGATCAAGATCGGGATACGCATCGTACAGt                           <  1:1334512/89‑1 (MQ=255)

                                                                gtcgtcgAAATCTTCCAGATCGCAAACCCACTGATGCATATCGGTGAATCGAACCGTTTTCGGATCAAGATCGGGATACGCATCGTACAGTGc                         >  1:2175822/1‑93 (MQ=255)

                                                                 tcgtcgAAATCTTCCAGATCGCAAACCCACTGATGCATATCGGTGAATCGAACCGTTTTCGGATCAAGATCGGGATACGCATCGTACAGTGCt                        >  1:4046397/1‑93 (MQ=255)

                                                                 tcgtcgAAATCTTCCAGATCGCAAACCCACTGATGCATATCGGTGAATCGAACCGTTTTCGGATCAAGATCGGGATACGCATCGTACAGTGCt                        >  1:268384/1‑93 (MQ=255)

                                                                 tcgtcgAAATCTTCCAGATCGCAAACCCACTGATGCATATCGCTGAATCGAACCGTTTTCGGATCAAGATCGGGATACGCATCGTACAGTGCt                        >  1:3296896/1‑93 (MQ=255)

                                                                   gtcgAAATCTTCCAGATCGCAAACCCACTGATGCATATCGGTGAATCGAACCGTTTTCGGATCAAGATCGGGATACGCATCGTACAGTGCTTc                      >  1:1010553/1‑93 (MQ=255)

                                                                   gtcgAAATCTTCCAGATCGCAAACCCACTGATGCATATCGGTGAATCGAACCGTTTTCGGATCAAGATCGGGATACGCATCGTACAGTGCTTc                      >  1:723395/1‑93 (MQ=255)

                                                                    tcgAAATCTTCCAGATCGCAAACCCACTGATGCATATCGGTGAATCGAACCGTTTTCGGATCAAGATCGGGATACGCATCGTACAGTGCTTCg                     >  1:2375092/1‑93 (MQ=255)

                                                                    tcgAAATCTTCCAGATCGCAAACCCACTGATGCATATCGGTGAATCGAACCGTTTTCGGATCAAGATCGGGATACGCATCGTACAGTGCTTCg                     >  1:651488/1‑93 (MQ=255)

                                                                      gAAATCTTCCAGATCGCAAACCCACTGATGCATATCGGTGAATCGAACCGTTTTCGGATCAAGATCGGGATACGCATCGTACAGTGCTTCGcc                   >  2:2580970/1‑93 (MQ=255)

                                                                      gAAATCTTCCAGATCGCAAACCCACTGATGCATATCGGTGAATCGAACCGTTTTCGGATCAAGATCGGGATACGCATCGTACAGTGCTTCGcc                   >  2:1273772/1‑93 (MQ=255)

                                                                       aaaTCTTCCAGATCGCAAACCCACTGATGCATATCGGTGAATCGAACCGTTTTCGGATCAAGATCGGGATACGCATCGTACAGTGCTTCGCCa                  >  2:1685569/1‑93 (MQ=255)

                                                                        aagcTTCCAGATCGCAAACCCACTGATGCATATCGGTGAATCGAACCGTTTTCGGATCAAGATCGGGATACGCATCGTACAGTGCTTCGCCaa                 <  2:2277829/90‑1 (MQ=255)

                                                                        aaTCTTCCAGATCGCAAACCCACTGATGCATATCGGTGAATCGAACCGTTTTCGGATCAAGATCGGGATACGCATCGTACAGTGCTTCGCCaa                 <  2:3138860/93‑1 (MQ=255)

                                                                        aaTCTTCCAGATCGCAAACCCACTGATGCATATCGGTGAATCGAACCGTTTTCGGATCAAGATCGGGATACGCATCGTACAGTGCTTCGCCaa                 <  2:1775581/93‑1 (MQ=255)

                                                                        aaTCTTCCAGATCGCAAACCCACTGATGCATATCGGTGAATCGAACCGTTTTCGGATCAAGATCGGGATACGCATCGTACAGTGCTTCGCCaa                 <  2:2109513/93‑1 (MQ=255)

                                                                          tCTTCCAGATCGCAAACCCACTGATGCATATCGGTGAATCGAACCGTTTTCGGATCAAGATCGGGATACGCATCGTACAGTGCTTCGCCAAtt               >  2:4227188/1‑93 (MQ=255)

                                                                            ttCCAGATCGCAAACCCACTGATGCATATCGGTGAATCGAACCGTTTTCGGATCAAGATCGGGATACGCATCGTACAGTGCTTCGCCAATTTc             >  1:3967605/1‑93 (MQ=255)

                                                                            ttCCAGATCGCAAACCCACTGATGCATATCGGTGAATCGAACCGTTTTCGGATCAAGATCGGGATACGCATCGTACAGTGCTTCGCCAATTTc             <  1:4662714/93‑1 (MQ=255)

                                                                            ttCCAGATCGCAAACCCACTGATGCATATCGGTGAATCGAACCGTTTTCGGATCAAGATCGGGATACGCATCGTACAGTGCTTCGCCAATTTc             <  2:2262793/93‑1 (MQ=255)

                                                                             tCCAGATCGCAAACCCACTGATGCATATCGGTGAATCGAACCGTTTTCGGATCAAGATCGGGATACGCATCGTACAGTGCTTCGCCAATTTcg            >  1:533983/1‑93 (MQ=255)

                                                                             tCCAGATCGCAAACCCACTGATGCATATCGGTGAATCGAACCGTTTTCGGATCAAGATCGGGATACGCATCGTACAGTGCTTCGCCAATTTcg            <  1:1158203/93‑1 (MQ=255)

                                                                              ccAGATCGCAAACCCACTGATGCATATCGGTGAATCGAACCGTTTTCGGATCAAGATCGGGATACGCATCGTACAGTGCTTCGCCAATTTcgc           <  2:3057533/93‑1 (MQ=255)

                                                                               cAGATCGCAAACCCACTGATGCATATCGGTGAATCGAACCGTTTTCGGATCAAGATCGGGATACGCATCGTACAGTGCTTCGCCAATTTcgcg          <  2:4275455/93‑1 (MQ=255)

                                                                                 gATCGCAAACCCACTGATGCATATCGGTGAATCGAACCGTTTTCGGATCAAGATCGGGATACGCATCGTACAGTGCTTCGCCAATTTCGCGGc        <  1:279810/93‑1 (MQ=255)

                                                                                  aTCGCAAACCCACTGATGCATATCGGTGAATCGAACCGTTTTCGGATCAAGATCGGGATACGCATCGTACAGTGCTTCGCCAATTTCGCGGCt       <  1:1578171/93‑1 (MQ=255)

                                                                                    cGCAAACCCACTGATGCATATCGGTGAATCGAACCGTTTTCGGATCAAGATCGGGATACGCATCGTACAGTGCTTCGCCAATTTCGCGGCTAt     >  2:1625895/1‑93 (MQ=255)

                                                                                     gCAAACCCACTGATGCATATCGGTGAATCGAACCGTTTTCGGATCAAGATCGGGATACGCATCGTACAGTGCTTCGCCAATTTCGCGGCTATc    <  1:4660113/93‑1 (MQ=255)

                                                                                     gCAAACCCACTGATGCATATCGGTGAATCGAACCGTTTTCGGATCAAGATCGGGATACGCATCGTACAGTGCTTCGCCAATTTCGCGGCTATc    <  1:1795701/93‑1 (MQ=255)

                                                                                      cAAACCCACTGATGCATATCGGTGAATCGAACCGTTTTCGGATCAAGATCGGGATACGCATCGTACAGTGCTTCGCCAATTTCGCGGCTATCg   >  2:829152/1‑93 (MQ=255)

                                                                                       aaaCCCACTGATGCATATCGGTGAATCGAACCGTTTTCGGATCAAGATCGGGATACGCATCGTACAGTGCTTCGCCAATTTCGCGGCTATCgg  >  2:2805976/1‑93 (MQ=255)

CGGCCTCGTCCAGCCAGACTAACAAAATCGCTTCGAGGATTTTCTCGTTGGATGCCTGCGGGTCGTCGTCGAAATCTTCCAGATCGCAAATCCACTGATGCATATCGGTGAATCGAACCGTTTTCGGATCAAGATCGGGATACGCATCGTACAGTGCTTCGCCAATTTCGCGGCTATCGG  >  NC_000913/2654561‑2654740

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: `‑`