BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,158,977	T→C	V111A (GTC→GCC)	yijC →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,158,977	0	T	C	98.8%	501.2 / ‑4.9	161	V111A (GTC→GCC)	yijC	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (1/1); new base C (62/96); total (64/97)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 170 total aligned reads displayed.

TGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGTTTCGT  >  NC_000913/4158887‑4159065

tGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGcct                                                                                        >  1:1636780/1‑93 (MQ=255)

 ggTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGcctc                                                                                       >  1:3912945/1‑93 (MQ=255)

  gTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTcc                                                                                      <  1:3460755/93‑1 (MQ=255)

   ttGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCca                                                                                     >  1:3284337/1‑93 (MQ=255)

    tgatgaGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCcac                                                                                    <  2:1592838/93‑1 (MQ=255)

     gatgaGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCcaca                                                                                   <  1:1101902/93‑1 (MQ=255)

       tgaGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACAtt                                                                                 <  1:1910756/93‑1 (MQ=255)

       tgaGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACAtt                                                                                 <  2:1982409/93‑1 (MQ=255)

       tgaGAGCGGTTTAATGCTACGCCAACGCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACAtt                                                                                 <  1:2317279/93‑1 (MQ=255)

        gagaGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACAttt                                                                                >  2:3337009/1‑93 (MQ=255)

         agagCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTa                                                                               <  1:3404916/93‑1 (MQ=255)

         agagCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTa                                                                               <  1:3161067/93‑1 (MQ=255)

          gagCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTAt                                                                              <  2:3022466/93‑1 (MQ=255)

          gagCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTAt                                                                              <  1:2085338/93‑1 (MQ=255)

           agCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATg                                                                             >  2:854618/1‑93 (MQ=255)

             cGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGa                                                                           <  2:248208/93‑1 (MQ=255)

             cGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGa                                                                           <  1:3963611/93‑1 (MQ=255)

              ggTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAg                                                                          >  1:1113819/1‑93 (MQ=255)

               gTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGt                                                                         >  2:2429438/1‑93 (MQ=255)

               gTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGt                                                                         >  1:301008/1‑93 (MQ=255)

               gTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCAACATTTATGGAGt                                                                         >  1:301009/1‑93 (MQ=255)

                tttAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGtt                                                                        >  1:2455611/1‑93 (MQ=255)

                tttAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGtt                                                                        <  2:4392774/93‑1 (MQ=255)

                  tAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCa                                                                      <  1:1815001/93‑1 (MQ=255)

                   aaTGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGGCTCCACATTTATGGAGTTCAt                                                                     >  2:2206471/1‑93 (MQ=255)

                   aaTGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCAt                                                                     <  2:956963/93‑1 (MQ=255)

                   aaTGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCATCCACATTTATGGAGTTCAt                                                                     <  1:3930808/93‑1 (MQ=255)

                    aTGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTGATc                                                                    <  2:1019543/93‑1 (MQ=255)

                    aTGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATc                                                                    <  2:2615490/93‑1 (MQ=255)

                     tGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCg                                                                   >  1:1378023/1‑93 (MQ=255)

                      gCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCgg                                                                  >  1:1256141/1‑93 (MQ=255)

                       cTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGt                                                                 <  2:1568098/93‑1 (MQ=255)

                        tACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGta                                                                >  2:1363815/1‑93 (MQ=255)

                         aCGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGtaa                                                               >  1:2108644/1‑93 (MQ=255)

                             cAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTaataat                                                           <  2:1162928/93‑1 (MQ=255)

                               aCTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATcc                                                         >  2:297593/1‑93 (MQ=255)

                                cTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCt                                                        >  1:4723789/1‑93 (MQ=255)

                                cTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCt                                                        <  1:2493152/93‑1 (MQ=255)

                                  cATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTaa                                                      <  1:3860548/93‑1 (MQ=255)

                                  cATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTaa                                                      <  2:1496755/93‑1 (MQ=255)

                                  cATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTaa                                                      <  2:4532467/93‑1 (MQ=255)

                                       gcCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGGGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCt                                                 >  1:3608710/1‑93 (MQ=255)

                                        ccAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCtt                                                >  1:126501/1‑93 (MQ=255)

                                         cAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCGCATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTc                                               >  2:3326501/1‑93 (MQ=255)

                                         cAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTc                                               >  2:2863273/1‑93 (MQ=255)

                                         cAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTc                                               >  1:917875/1‑93 (MQ=255)

                                          aGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTcc                                              <  2:1278544/93‑1 (MQ=255)

                                           ggCGCGTCAGCGTATCGCCAAAGGCGGGCGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAACCCTAACGCCTTCCg                                             >  1:1527178/1‑93 (MQ=255)

                                           ggCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCg                                             >  2:3441014/1‑93 (MQ=255)

                                            gcgcGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCgg                                            >  2:3753899/1‑93 (MQ=255)

                                            gcgcGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCgg                                            <  2:2461082/93‑1 (MQ=255)

                                             cgcgTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGt                                           <  2:3918150/93‑1 (MQ=255)

                                                gTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGttat                                        >  2:533083/1‑93 (MQ=255)

                                                gTCAGCGTATCGCCAAAGCCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGttat                                        >  2:723675/1‑93 (MQ=255)

                                                 tCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGttatt                                       <  2:4296232/93‑1 (MQ=255)

                                                 tCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGttatt                                       <  2:3916416/93‑1 (MQ=255)

                                                 tCAGCGTATCGCCAAAGGCGAGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGttatt                                       <  2:3143846/93‑1 (MQ=255)

                                                  cAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGttatta                                      <  2:802159/93‑1 (MQ=255)

                                                    gCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTattatt                                    >  1:261966/1‑93 (MQ=255)

                                                    gCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTattatt                                    >  1:4154582/1‑93 (MQ=255)

                                                    gCGTATCGCCAAAGGCGGGAGCGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTattatt                                    >  2:4322534/1‑93 (MQ=255)

                                                    gCGTATCCCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGATTTCATCGGTAATAATCCTAACGCCTTCCGGTTattatt                                    <  1:3385773/93‑1 (MQ=255)

                                                    gCGGATCGCCAAAGGGGGGAGGGTGATCCGGACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTattatt                                    <  2:3756905/93‑1 (MQ=255)

                                                      gTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGc                                  >  1:2563957/1‑93 (MQ=255)

                                                      gTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGc                                  <  2:4126181/93‑1 (MQ=255)

                                                         tCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTACGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCggg                               <  2:2548422/93‑1 (MQ=255)

                                                          cGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGa                              <  2:4703285/93‑1 (MQ=255)

                                                           gCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGaa                             <  1:2222937/93‑1 (MQ=255)

                                                            ccAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAAc                            <  2:3631915/93‑1 (MQ=255)

                                                             cAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACg                           >  2:2236095/1‑93 (MQ=255)

                                                             cAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACg                           <  1:1138969/93‑1 (MQ=255)

                                                             cAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACg                           >  2:3808184/1‑93 (MQ=255)

                                                              aaaGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGc                          <  1:631837/93‑1 (MQ=255)

                                                               aaGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCt                         >  1:1918031/1‑93 (MQ=255)

                                                               aaGGCGGGAGGGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGcc                         >  2:1251830/1‑92 (MQ=255)

                                                                 ggCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTcc                       <  1:1769932/93‑1 (MQ=255)

                                                                 ggCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTcc                       <  1:420537/93‑1 (MQ=255)

                                                                    gggAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCGTCCGGTTATTATTGCGGGAACGCTCCgca                    >  1:3225735/1‑91 (MQ=255)

                                                                     ggAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCa                   <  2:2122732/93‑1 (MQ=255)

                                                                     ggAGTGTGATCCGCACCTCGGCCTCCACATTTAAGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCa                   <  1:4586975/93‑1 (MQ=255)

                                                                      gAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCAc                  >  2:1099276/1‑93 (MQ=255)

                                                                        gtgtGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCt                >  2:4066731/1‑93 (MQ=255)

                                                                         tgtgATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTc               <  1:474500/93‑1 (MQ=255)

                                                                         tgtgATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTc               <  2:3386672/93‑1 (MQ=255)

                                                                           tgATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCg             >  1:2205990/1‑93 (MQ=255)

                                                                            gATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGc            >  1:2844026/1‑93 (MQ=255)

                                                                              tCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTg          >  1:3737368/1‑93 (MQ=255)

                                                                                 gCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGt       <  2:1842317/93‑1 (MQ=255)

                                                                                 gCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGt       >  1:148449/1‑93 (MQ=255)

                                                                                 gCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGt       <  1:4495302/93‑1 (MQ=255)

                                                                                  cACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGtt      <  2:1942205/93‑1 (MQ=255)

                                                                                  cACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGtt      <  1:2263903/93‑1 (MQ=255)

                                                                                   cccTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGttt     <  1:2699857/92‑1 (MQ=255)

                                                                                   aCCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGttt     >  2:4573923/1‑93 (MQ=255)

                                                                                   aCCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGttt     >  1:4301945/1‑93 (MQ=255)

                                                                                   aCCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGttt     <  1:1422217/93‑1 (MQ=255)

                                                                                   aCCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGttt     <  2:3150528/93‑1 (MQ=255)

                                                                                    ccTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGTTTc    <  2:42585/93‑1 (MQ=255)

                                                                                    ccTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGTTTc    <  2:4192519/93‑1 (MQ=255)

                                                                                      tCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGTTTCGt  <  1:4394718/93‑1 (MQ=255)

TGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGTTTCGT  >  NC_000913/4158887‑4159065

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