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	99.6%	820.5 / NA	259	V111A (GTC→GCC)	yijC	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/0); major base C (123/135); minor base G (1/0); total (124/135)
Fisher's exact test for biased strand distribution p-value = 4.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 270 total aligned reads displayed.

TGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGTTTCGTG  >  NC_000913/4158887‑4159066

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

  gTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTcc                                                                                       >  1:371499/1‑93 (MQ=255)

  gTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTcc                                                                                       >  2:421035/1‑93 (MQ=255)

   ttGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCca                                                                                      <  2:137645/93‑1 (MQ=255)

   ttGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCca                                                                                      <  1:319043/93‑1 (MQ=255)

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

      atgaGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACAt                                                                                   <  2:5845993/93‑1 (MQ=255)

       tgaGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACAtt                                                                                  >  1:5945567/1‑93 (MQ=255)

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

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

       tgaGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGGGTGATCCGCACCTCGGCCTCCATAtt                                                                                  >  2:2175043/1‑93 (MQ=255)

        gagaGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACAttt                                                                                 <  2:5362066/93‑1 (MQ=255)

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

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

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

         agagCAGTTGAATGCAACGCCAACTCATGCGCCAGGCGCGTCAGCTTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTa                                                                                <  2:564845/93‑1 (MQ=255)

            gCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATgg                                                                             <  1:1886903/93‑1 (MQ=255)

            gCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATgg                                                                             >  1:3828117/1‑93 (MQ=255)

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

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

             aggTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGa                                                                            <  1:3798914/92‑1 (MQ=255)

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

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

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

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

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

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

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

                  tAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCa                                                                       >  2:294270/1‑93 (MQ=255)

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

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

                    aTGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATc                                                                     >  1:2133568/1‑93 (MQ=255)

                     tGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCTGCACCTCGGCCTCCACATTTATGGAGTTCATCg                                                                    <  1:687565/93‑1 (MQ=255)

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

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

                         aCGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGtaa                                                                <  2:97527/93‑1 (MQ=255)

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

                         aCGCCAACTCAGGCGCCAGGCGCGTCAGCGTAGCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGtaa                                                                <  2:2854894/93‑1 (MQ=255)

                            ccAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGtaataa                                                             <  1:2319422/93‑1 (MQ=255)

                            ccAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGtaataa                                                             >  2:5184042/1‑93 (MQ=255)

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

                                cTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGGGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCt                                                         >  1:1534401/1‑93 (MQ=255)

                                 tCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTa                                                        <  1:2936661/93‑1 (MQ=255)

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

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

                                  cATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTaa                                                       >  1:5805023/1‑93 (MQ=255)

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

                                   aggcgcGAGGCGCGTCATCGTATCGCGAAAGGCGGGAGGGTGATCCGCACCGCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAAc                                                      <  2:2636550/91‑1 (MQ=255)

                                   aTGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAAc                                                      <  2:4973255/93‑1 (MQ=255)

                                   aTGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGCCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAAc                                                      <  1:3752429/93‑1 (MQ=255)

                                    tGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACg                                                     >  1:6031510/1‑93 (MQ=255)

                                       gcCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCt                                                  >  1:2185528/1‑93 (MQ=255)

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

                                        ccAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCtt                                                 <  2:446295/93‑1 (MQ=255)

                                         cAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTc                                                <  1:3020242/93‑1 (MQ=255)

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

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

                                           ggCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCAGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAAAGACTTCCg                                              >  2:921630/1‑93 (MQ=255)

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

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

                                              gcgTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGtt                                           >  2:2001597/1‑93 (MQ=255)

                                              gcgTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGtt                                           >  2:114810/1‑93 (MQ=255)

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

                                                 tCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGttatt                                        >  1:5487713/1‑93 (MQ=255)

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

                                                  cAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGttatta                                       >  1:4932275/1‑93 (MQ=255)

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

                                                   aGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTtattat                                      <  1:4004161/93‑1 (MQ=255)

                                                   aGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTtattat                                      >  1:3953425/1‑93 (MQ=255)

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

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

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

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

                                                       tATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCg                                  <  2:4910245/93‑1 (MQ=255)

                                                         tCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCggg                                >  2:5495897/1‑93 (MQ=255)

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

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

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

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

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

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

                                                                aGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTc                         >  1:4791307/1‑93 (MQ=255)

                                                                aGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTc                         <  1:3402139/93‑1 (MQ=255)

                                                                aGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTc                         <  2:25521/93‑1 (MQ=255)

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

                                                                  gCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCg                       <  1:2510387/93‑1 (MQ=255)

                                                                   cGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCgg                      <  1:1111942/93‑1 (MQ=255)

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

                                                                          gtgATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTcc               <  1:5080420/93‑1 (MQ=255)

                                                                          gtgATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTcc               >  1:3266967/1‑93 (MQ=255)

                                                                          gtgATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTcc               <  2:3222502/93‑1 (MQ=255)

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

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

                                                                               ccGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGc          >  2:2769429/1‑93 (MQ=255)

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

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

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

                                                                                     cTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGTTTCg    >  1:4500887/1‑93 (MQ=255)

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

                                                                                       cGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGTTTCGTg  >  2:660530/1‑93 (MQ=255)

TGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGTTTCGTG  >  NC_000913/4158887‑4159066

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